Next Article in Journal
Air Combat Intention Recognition with Incomplete Information Based on Decision Tree and GRU Network
Next Article in Special Issue
On Constructing Informationally Complete Covariant Positive Operator-Valued Measures
Previous Article in Journal
A Novel Underwater Acoustic Target Identification Method Based on Spectral Characteristic Extraction via Modified Adaptive Chirp Mode Decomposition
Previous Article in Special Issue
Generalized Toffoli Gate Decomposition Using Ququints: Towards Realizing Grover’s Algorithm with Qudits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Semi-Empirical Satellite-to-Ground Quantum Key Distribution Model for Realistic Receivers

by
Aleksandr V. Khmelev
1,2,3,*,
Egor I. Ivchenko
1,2,3,
Alexander V. Miller
3,4,
Alexey V. Duplinsky
3,5,
Vladimir L. Kurochkin
2,3,6,* and
Yury V. Kurochkin
7
1
Landau Phystech School of Physics and Research, Moscow Institute of Physics and Technology, Dolgoprudny, Moscow 141701, Russia
2
Russian Quantum Center, Skolkovo, Moscow 143025, Russia
3
QSpace Technologies, Moscow 121205, Russia
4
Moscow Centre for Quantum Technologies, Moscow 127204, Russia
5
Moscow Institute of Electronics and Mathematics, HSE University, Moscow 101000, Russia
6
NTI Center for Quantum Communications, NUST MISiS, Moscow 115419, Russia
7
Quantum Research Centre, Technology Innovation Institute, Masdar City, Abu Dhabi 9639, United Arab Emirates
*
Authors to whom correspondence should be addressed.
Entropy 2023, 25(4), 670; https://doi.org/10.3390/e25040670
Submission received: 21 February 2023 / Revised: 7 April 2023 / Accepted: 11 April 2023 / Published: 17 April 2023

Abstract

:
Satellite-based link analysis is valuable for efficient and secure quantum communication, despite seasonal limits and restrictions on transmission times. A semi-empirical quantum key distribution model for satellite-based systems was proposed that simplifies simulations of communication links. Unlike other theoretical models, our approach was based on the experimentally-determined atmospheric extinction coefficient typical for mid-latitude ground stations. The parameter was measured for both clear and foggy conditions, and it was validated using published experimental data from the Micius satellite. Using this model, we simulated secure QKD between the Micius satellite and ground stations with 300 mm and 600 mm aperture telescopes.

1. Introduction

Quantum communication technology provides a theoretically unhackable communication channel for encryption key exchange and long-term data security [1]. Based on quantum physics principles, the method, also known as a quantum key distribution (QKD), allows users to detect eavesdropping and estimate key leakage caused by hacker-induced errors in quantum states [2,3].
Typically, quantum communications rely on photon states being transferred across an optical fiber or through the atmosphere. However, the amount of light in the medium decreases exponentially with distance, and no replication or amplification is possible. As a result, the feasible range of QKD in optical fibers [4,5] and in terrestrial free space [6,7] is limited to a few hundred kilometers. Satellite-based QKD was proposed to expand the range of quantum communication because most of the signal transmission path between the satellite and the Earth passes in a vacuum. Several scientific organizations initiated satellite-based QKD research projects and demonstrated certain technological concepts [8,9].
Significant accomplishments were reported in Japan, where results on the conservation polarization characteristics of the laser beam from a satellite and a link parameters estimation for satellite-to-ground QKD were provided [10]. The first quantum-limited experiment between a low-Earth orbit (LEO) satellite and a ground station was conducted [11]. Quantum-limited tests were carried out on a geostationary satellite in 2017 [12]. In 2016, researchers from USTC (China) launched Micius, the first operating quantum science satellite, and successfully demonstrated satellite-to-ground QKD via polarization photon states [13,14].
As a promising concept, the first quantum science satellite paved the way for a future global quantum network, enabling secure worldwide data transfer [15,16,17]. Meanwhile, as a communication network capable of providing secret keys as a service, industrial quantum key distribution networks must meet the core requirements of flexible extension, cost effectiveness, component compatibility, and standardization [18,19,20,21]. As a result, a strict and comprehensive analysis of a quantum communication channel is needed to build an efficient quantum network [22,23], despite seasonal limits and restrictions on transmission times between satellites and ground stations.
Satellite-based uplink–downlink [24] and measurement-device-independent QKD models [25] were proposed recently. Additionally, a theoretical analysis of the quantum communication link between the Micius satellite and the Delingha receiving ground station was performed [22,26]. The existing theoretical models are valuable for a comprehensive analysis of free-space QKD. However, these models take into consideration many factors not known with sufficient accuracy, such as air turbulence, background noise, and pointing and tracking losses [8], which significantly limit their practical usage.
In this paper, we present a semi-empirical satellite-to-ground QKD model relying on phenomenological data by employing atmospheric extinction measurements. This allows for a realistic analysis of developed ground stations with 300 and 600 mm telescope apertures in terms of secure key rate and error rate. Considering our simulation of experimental results obtained with Micius [13] and a theoretical model [8,22], a more relevant and real range of sifted key rate for the satellite-based QKD was provided.
This paper is organized as follows: in Section 2, we describe the communication channel geometry for a series of satellite passages over Zvenigorod observatory and computed the time-dependent distance between the satellite and ground stations; in Section 3, we calculated overall link losses based on the measured characteristics of ground stations [27] and the photon source parameters of a satellite [13]; Section 4 examined the secure key generation and error rates of QKD systems using the standard BB84 decoy-state protocol [28,29]; finally, in Section 5, we discuss the deployment of our model to future space-to-ground quantum communication experiments.

2. Satellite-to-Ground Communication Geometry

To determine the channel efficiency in a satellite-based QKD model, one needs to calculate the time-dependent distance d ( t ) between the transmitter and the receivers, as well as the satellite elevation angle θ E l above the horizon. We considered several typical satellite passages and calculated one passing through the zenith in a circular orbit.
Three non-zenith actual passages of the Micius satellite were simulated for 31 October 2021, 9 March 2022, and 23 March 2022, with the maximum elevation angles of 32.5°, 57.5°, and 42°, respectively. To compute the time dependency of the satellite’s elevation angle for the (55° 42 N, 36° 45 E) location, we used the simplified general perturbations model SGP4 [30], given the satellite’s two-line elements (TLE) orbital data [31,32].
We also considered a simple model to calculate the trace-time function of the elevation angle for the limit case corresponding to the satellite passing through the zenith point over the ground station.Here, we assumed that the Micius satellite moves in a perfectly circular orbit (the eccentricity is 0) at an altitude of 500 km and an inclination of 97.33° [13].
The angular velocity ω of a low-Earth orbit satellite may be approximated by a constant in the Earth-centered fixed frame, according to Ref. [33]:
ω ω S ω E cos i ,
where ω S g / R is the angular velocity of the satellite in the Earth-centred inertial coordinate system with acceleration due to gravity of g and an orbit radius of R, the orbit inclination is i, and ω E is the angular velocity of the Earth’s rotation. The numerical values of the orbital parameters and angular velocity data of the Micius satellite for our model are listed in Table 1.
We approximated the magnitude of the tangential velocity by its minimum value for a given orbit in the Earth-centered fixed frame. The absolute fluctuation of the tangential velocity magnitude for the LEO satellite is quite negligible for any latitude, according to Ref. [34]. Therefore, angular velocity is independent of the ground station latitude in this model.
The elevation angle for the zenith satellite passage was calculated using a set of equations derived from simple geometric considerations (see Figure 1 and Appendix A.1), as follows:
θ E l = arccos sin 2 α 1 + R E R 2 2 R E R cos α ,
where the central angle can be expressed as a linear function of time as a α = ω t .
Next, we determined the distance between a satellite and a ground station using the elevation angle from TLE data over time for real satellite passage as well as a computed satellite pass above the zenith. Appendix A.2 has a detailed derivation of the distance, which is given by:
d 2 + d · 2 R E sin θ E l + R E 2 R 2 = 0 .
Figure 2 shows the dependence of the satellite elevation angle and the distance between the satellite and the ground station on time obtained from Equations (1)–(3) for the typical satellite passages, by using TLE data and the calculated zenith satellite passage. The calculations were carried out when the satellite was in the ground station’s field of view, i.e., when the elevation angle was above 20 degrees above the horizon.
The time duration of quantum communication increases with a peak elevation angle in Figure 2, whereas the minimal distance between the satellite and receiver decreases. The satellite passage through the zenith offers the longest communication time and and the shortest communication channel length. For comparison, the communication time of four satellite passages with the maximum elevation angles of 32.5°, 42°, 57.5°, and 90° are equal to 221 s, 260 s, 283 s, and 285 s, respectively.

3. Link Efficiency for 300-mm and 600-mm Aperture Ground Stations

Channel losses between a satellite and a ground station dynamically change during quantum communication, which significantly distinguishes this method from fiber optic or terrestrial QKD. In general, diffraction losses vary with the change of the distance between a satellite and a ground station. Additionally, varying the elevation angle results in various effective thicknesses of the atmosphere during light transmission.
The diffraction loss can be calculated using the laser source divergence γ , communication channel length d, and effective receiving area with telescope diameter D T and obstruction ε of the telescope secondary mirror, as follows ε D T 2 / ( γ d ) 2 . However, the atmospheric extinction should be described in more detail.
In astronomy, the magnitude of a star or another visible celestial object, such as a satellite, indicates its brightness class. The magnitudes m 1 and m 2 of two objects are related to the corresponding flux densities F 1 and F 2 by equation [35]:
m 1 m 2 = 2.5 log 10 ( F 1 / F 2 ) .
Given that starlight travels a greater distance through the atmosphere near the horizon than it does at the zenith, it follows that a star near the horizon will be less luminous. Therefore, we can relate the magnitude of the observed object outside the atmosphere m 0 (apparent magnitude) to the magnitude of the object at the Earth’s surface m, as follows:
m m 0 = ϰ · f ( θ E l ) ,
where ϰ is the atmospheric extinction coefficient and f ( θ E l ) is the air mass function depending on the elevation angle of the object. The amount of air at the zenith corresponds to one air mass.
Young and Irvine [36] proposed the relationship between the air mass function and elevation angles as:
f ( θ E l ) = csc θ E l · ( 1 0.0012 cot 2 θ E l ) .
The atmospheric extinction coefficient was determined from the measurements of stars’ fluxes with known apparent magnitudes outside the atmosphere in the spectral range 845 nm–855 nm [37] at several elevation angles by means of a 600-mm ground station. We assumed that, for the same ground station, the star fluxes were equally proportional to the photon count rate. Appendix B contains additional information and data.
Figure 3 depicts the dependence of atmospheric extinction ϰ f ( θ E l ) on air mass up to a constant, according to Equations (A13) and (6). The measurements of star fluxes were performed on 9 March 2022, with clear weather conditions (humidity < 20 % ) and on 24 June 2021, with hazy, foggy weather conditions (humidity < 80 % near the Earth’s surface). As a result, the slope of atmospheric extinction equals the atmospheric extinction coefficient and ranges from 0.23 ± 0.08 for clear nights to 0.41 ± 0.09 for foggy nights.
Hence, link efficiency, or the overall transmission of photons at 850 nm between the Micius satellite and a ground station, including detection efficiency, is given by:
η ( t ) = ε D T 2 ( γ d ) 2 · 10 0.4 ϰ csc θ E l · ( 1 0.0012 cot 2 θ E l ) η o p t η d e t ,
where η o p t is an optical efficiency of the respective ground station and η d e t is a quantum efficiency of a single photon detector.
Table 2 summarizes the baseline parameters of the setups required for further modeling. The 600 mm ground station is a fixed telescope equipped with a quantum states decoder, as detailed in Ref. [27]. The 300 mm ground station is a portable receiving telescope equipped with a quantum communications optical system [38], which has a better optical efficiency than a 600 mm telescope. The detailed description of the experimental setups are provided in Appendix C.
Figure 4 presents calculated link loss η l i n k = 10 log 10 η ( t ) for our satellite passages in the satellite QKD experiment under clear weather conditions. For the given system parameters, the peak link losses at zenith for the ground stations with the telescope apertures of 300 mm and 600 mm are equal to 32 dB and 29 dB, respectively. The characteristic difference of channel losses between elevation angles of 20° and 90° are 9 dB for both receivers.

4. Satellite-to-Ground QKD Analysis

In this section, we simulated a satellite-to-ground QKD experiment using the BB84 decoy-state protocol [28,39]. We determined the sifted key rate R s i f t and quantum bit error rate (QBER) of the designed receivers. Then, we estimated the secret key length by applying data post-processing for signal state QKD.

4.1. Satellite-QKD Model and Parameters Verification

We used the QKD model from Ref. [28], with the difference that the link transmission was time-dependent. Therefore, transmission of a k–photon state through the communication channel described in Section 3 is given by:
η k = 1 ( 1 η ( t ) ) k
for k = 0 , 1 , 2 , .
The photon number in each weak coherent pulse of the satellite transmitter has a Poisson distribution [13]. Pulses with a frequency f and intensities of signal μ , decoy ν , and vacuum λ =0 states are emitted at random with probabilities p s , p d , and p v , respectively.
Background noise Y 0 , which includes detector dark counts as well as additional background contributions, was considered to be 250 and 500 counts per second for ground stations with 300 mm and 600 mm apertures, respectively. Note that the values were assumed to be constant during communication and were obtained while tracking the switched-off satellite by the ground stations.
Additionally, intrinsic error e d e t caused by receiver measurement fidelity remained independent of distance and was set at 0.9 % and 0.5 % for the 300 mm and 600 mm ground stations, respectively. The values used to simulate our QKD setups are listed in Table 3.
The sifted key rate R s i f t and QBER were calculated from parameters of the optical link for the weak and vacuum decoy-state BB84 model provided by Xiongfeng Ma et al. [28]. The expressions of values for signal states are given by:
R s i f t = 1 2 f p s [ Y 0 + 1 e μ η ( t ) ] ,
Q B E R = 1 2 f p s [ e 0 Y 0 + e d e t ( 1 e μ η ( t ) ) ] / R s i f t .
Figure 5 presents the calculated R s i f t and QBER of signal states during the satellite passages with peak elevation angles of 32.5° and 90°. The sifted key rate has a strong dependence on the maximum elevation angle and changes non-uniformly over the communication time.
We list the sifted key rate at a maximum value are 3.3 kbit s 1 and 11.4 kbit s 1 for the 300 mm ground station, and 7.2 kbit s 1 and 25.3 kbit s 1 for the 600 mm one. The minimal sifted key generation corresponds to the elevation angle of 20° and equals 1.4 kbit s 1 and 3.1 kbit s 1 for the ground stations with the 300 mm and 600 mm apertures, respectively. The QBER of signal states E μ is also different for the ground stations [28] and is limited to 3 % for the 300 mm aperture ground station, and to 2.5% for the 600 mm aperture.
To verify our semi-empirical model with a phenomenological parameter ϰ , we compared the simulated results for the Xinglong ground station with published experimental data collected with the Micius satellite [13] and modelled findings by Vergoossen et al. [22]. Figure 6 illustrates the peak sifted key rate on the shortest distance between the Micius satellite and the Xinglong ground station for specific satellite passages between September 2016 and May 2017.
The Xinglong ground station is equipped with a 1000-mm-aperture telescope, where the obstruction of the telescope’s secondary mirror is approximated to be 80%. In order to compare the models, we calculated the atmospheric extinction coefficient ϰ V e r g o o s s e n = 0.80 , which is equivalent to a −3.2 dB Vergoossen atmospheric absorption loss at zenith. The other inputs for the key rate calculations are provided in Vergoossen et al. (Table 3 [22]).
Thus, the peak sifted key rate in our model with the atmospheric extinction parameter ϰ V e r g o o s s e n agrees well with the model prediction provided by Vergoossen et al. [22]. Meanwhile, the experimentally obtained atmospheric extinction coefficient for mid-latitude ground stations under favourable weather conditions provides a better evaluation of the Micius data with the shortest distances of 530.33 km, 591.56 km, 711.94 km, and 929.67 km. Moreover, the atmospheric extinction coefficient for foggy, hazy weather conditions, together with less than 3 dB loss due to pointing error [13], is suitable for the assessment of the lower bound of the peak sifted key rate.
Hence, we have verified our model with the real experimental data. As a result, our semi-empirical model with phenomenological parameters is robust and closely matches the satellite-to-ground QKD experiment.

4.2. Secret Keys for 300 mm and 600 mm Ground Stations

We stress that only single photon states are completely secure in practical QKD with an imperfect single photon source [40]. The communication time of Δ t corresponds to elevation angles greater than 20° and sifted key length l s i f t is given:
l s i f t = Δ t R s i f t d t .
Therefore, we estimated the proportion of secret bits from obtained sifted keys using the classic BB84 decoy-state protocol, according to Refs. [28,39]. The lower limit of the final secret key length l s e c is as follows:
l s e c N p s q Q 1 1 H 2 e 1 1.44 Q μ H 2 ( E μ ) ,
where N is the overall number of pulses emitted by satellite, q = 1 / 2 is the BB84 protocol due to sifting procedure, the subscript μ denotes the intensity of signal states, Q μ is the gain of signal states, E μ is the mean overall QBER across communication time, Q 1 is the gain of single photon states, e 1 is the error rate of single photon states, the H 2 ( x ) is the binary Shannon entropy function, and 1.44 is the efficiency of error correction [41].
Given the finite size of the sifted key length, the average values contain statistical fluctuations that can be estimated via Chernoff inequality [42] and the observed values. The failure probability for the Chernoff bound is equal to 10 9 .
Table 4 shows the calculated sifted key length and the estimated secret key length, using Equation (12) and Chernoff inequality, for four satellite passages with the maximum elevation angles of 32.5°, 42°, 57.5°, and 90°.
For typical satellite passages with various peak elevation angles, the secret key lengths differ by 3.2 to 4.6 times for the ground stations with 300 mm and 600 mm apertures, whereas telescope areas differ by four times. For some satellite passages, the 300 mm ground station is more cost-effective and efficient for space-based quantum communications.

5. Conclusions

A semi-empirical QKD model for satellite-based systems that simplifies the practical simulation of communication links has been proposed. In contrast to other pure theoretical models, our semi-empirical approach was based on atmospheric extinction coefficients common for ground stations at mid-latitude. These coefficients were determined experimentally by means of real ground-based receivers developed by the authors.
To verify our method, the Xinglong ground station’s peak sifted key rate was simulated and compared to the published experimental data from the Micius satellite for specific satellite passages. The lower and upper boundaries of the peak sifted key rate were appropriately estimated by our model, employing best-choice phenomenological parameters. Then, using this model for ground stations with 300 and 600 mm aperture telescopes, we simulated the secure key generation from the Micius satellite.
To conclude, we proposed and verified a simple practical method for estimating the performance of ground-based stations, which are one of the key elements of future QKD networks. We believe that the developed model will be useful for quick, right-on-the-spot estimation of satellite QKD links and the analysis of future satellite-to-ground QKD systems.

Author Contributions

Conceptualization, V.L.K. and A.V.K.; methodology, A.V.K., E.I.I. and A.V.M.; investigation, A.V.K. and A.V.D.; writing—original draft preparation and editing, A.V.K. and A.V.M.; supervision, V.L.K. and Y.V.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Education and Science of the Russian Federation in the framework of the Program of Strategic Academic Leadership “Priority 2030” (Strategic Project “Quantum Internet”).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank A. Tayduganov and O. Fat’yanov for helpful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Derivation of Equation (2)

In this section, we derive the elevation angle for the satellite passage to describe the quantum communication link geometry, as illustrated in Figure 1. The following equations hold for a triangle with vertices at the satellite (S), the ground station (G), and the centre of the earth (O), where all the quantities are defined in the Section 2:
d 2 = R E 2 + R 2 2 R E R cos α sin α d = sin ( 90 + θ E l ) R
By removing variable of the distance d between the ground station and the satellite from the equation system, we get:
R 2 sin 2 α = R E 2 + R 2 2 R E R cos α cos 2 θ E l
Then, we obtain the required Equation (2), which is given as:
cos θ E l = sin 2 α 1 + R E R 2 2 R E R cos α

Appendix A.2. Derivation of Equation (3)

In this section, we show how to obtain the relation between the elevation angle and the distance d, that we used in the communication link geometry study of Section 2.
Considering T G S (Figure 1) we can write the following equation:
s i n ( α / 2 + θ E l ) h = sin ( 90 + α / 2 ) d
The result of simplifying this relation is:
h d sin θ E l cos α 2 = cos θ E l · sin α 2
Inserting the second Equation (A1) in the Equation (A5) we obtain the following formula:
h d sin θ E l cos α 2 = R d sin α · sin α 2
That turns into:
h d sin θ E l = R · ( 1 cos α )
Using the relation R = R E + h according to Figure 1, it can be expressed as follows:
R E + d sin θ E l = R · cos α
Hence, the first Equation (A1) and the Equation (A8) can be reduced to the following expression:
d 2 + d · 2 R E sin θ E l + R E 2 R 2 = 0
Then, by solving the quadratic equation, we obtain the result:
d = R E sin θ E l + R E 2 sin 2 θ E l + h 2 + 2 R E h

Appendix B. Atmospheric Extinction Coefficient

In this section, we provide the experimental findings of the atmospheric extinction coefficient obtained from the measurements of star fluxes in the I wavelength band. To measure star fluxes, we point our 600-mm ground station at the chosen star and measure the photon count rate throughout a 20-s period (Table A1).
Using Equation (5), the difference in magnitude (m) for two stars at the Earth’s surface in the same weather conditions (i.e., with the same atmospheric extinction coefficient ϰ ) can be written as:
m i m 0 i m j + m 0 j = ϰ ( f ( θ E l i ) f ( θ E l j ) ) ,
where superscripts (i and j) denote the different stars.
According to Equation (4), it can be expressed as follows:
2.5 log 10 F i m 0 i = ϰ f ( θ E l i ) 2.5 log 10 F j m 0 j ϰ f ( θ E l j )
Then, for the stars observed under identical weather conditions, Equation (A12) turns into:
2.5 log 10 F i m 0 i = ϰ f ( θ E l i ) + c o n s t ,
Therefore, using Equations (A13) and (6), we found that the atmospheric extinction coefficient ϰ ranges from 0.23 on a clear night to 0.41 on a foggy, hazy night. Note that the spectral range for both ground stations is the same.
Table A1. The experimental data of photon count rate collected on clear and foggy nights.
Table A1. The experimental data of photon count rate collected on clear and foggy nights.
DateStar NameAltitude ( θ El )Apparent Magnitude ( m 0 ) I-Band [43]Count Rate, Kcps
9 March 2022
(clear night)
η Ursae Majoris49.6°2.161100 ± 20
η Ursae Majoris56.6°2.161090 ± 20
η Ursae Majoris62.7°2.161110 ± 20
ζ Draconis40.0°2.35344 ± 7
ζ Draconis44.1°2.35336 ± 7
ζ Draconis47.6°2.35344 ± 7
24 June 2021
(foggy night)
o Ursae Majoris30.8°2.31800 ± 15
η Boötis33.3°1.951350 ± 30
ϵ Cygni51.1°1.192850 ± 60
γ Boötis53.3°2.8738 ± 15
α Ursae Minoris55.0°1.223230 ± 60
α Cygni58.9°1.043940 ± 80
η Herculis70.6°2.351070 ± 20

Appendix C. 300 mm and 600 mm Aperture Ground Stations

The experimental setups with 300 mm (see Figure A1) [38] and 600 mm apertures (see Figure A2) [35] are suggested for satellite-to-ground QKD. The ground stations are designed to function with optical radiation at a wavelength of 850 nm (quantum channel), 532 nm for the synchronisation procedure (SCM) and guiding system (GT), and 671 nm for the pointing by the laser (3W beacon LD).
The system works as follows: once reaching the receiving telescope, the quantum signal enters the processing and analysis unit, followed by the polarisation analyzer, which functions as a quantum state decoder. After that, the signal is transmitted via an optical fibre and detected by single-photon detectors.
Figure A1. Portable ground station with a 300 mm aperture for satellite quantum key distribution (right) and its scheme (left) [38]. CPL—collimator; POL—polarizer; PBS—polarizing beam splitter cube; HWP—half wave plate; BS—50:50 beam splitter cube; IF1—interference filter (central wavelength 850 nm; half-width 10 nm); HWP-m—motorized half-wave plate; M1—mirror; L1—lens; RT—receiving telescope with an aperture of 300 mm; GT—guiding telescope; SCM—synchronisation module; CAM—camera; IF2–—interference filter (central wavelength 532 nm; half-width 10 nm); 5 ch SPDs—5-channel single photon detector module; TDC—time-to-digital converter; PC—personal computer; PA—polarization analyzer (decoder); APS—module for analysis and processing of optical signals.
Figure A1. Portable ground station with a 300 mm aperture for satellite quantum key distribution (right) and its scheme (left) [38]. CPL—collimator; POL—polarizer; PBS—polarizing beam splitter cube; HWP—half wave plate; BS—50:50 beam splitter cube; IF1—interference filter (central wavelength 850 nm; half-width 10 nm); HWP-m—motorized half-wave plate; M1—mirror; L1—lens; RT—receiving telescope with an aperture of 300 mm; GT—guiding telescope; SCM—synchronisation module; CAM—camera; IF2–—interference filter (central wavelength 532 nm; half-width 10 nm); 5 ch SPDs—5-channel single photon detector module; TDC—time-to-digital converter; PC—personal computer; PA—polarization analyzer (decoder); APS—module for analysis and processing of optical signals.
Entropy 25 00670 g0a1
Figure A2. Fixed ground station with a 600 mm aperture for satellite quantum key distribution (right) and its scheme (left) [27]. CPL—collimator; POL—polarizer; PBS—polarizing beam splitter cube; HWP—half wave plate; BS—50:50 beam splitter cube; IF1—interference filter (central wavelength 850 nm; half-width 10 nm); HWP-m—motorized half-wave plate; BE—beam expander; DM is a dichroic mirror; L1—lens; FSM—fast steering mirror; RCT—Ritchey-Chrétien telescope with an aperture of 600 mm; GT—guiding telescope; SCM—synchronisation module; CAM—camera; IF2—interference filter (central wavelength 532 nm; half-width 10 nm); 5 ch SPDs—5-channel single photon detector module; TDC—time-to-digital converter; PC—personal computer; PA—polarization analyzer (decoder); APS—module for analysis and processing of optical signals.
Figure A2. Fixed ground station with a 600 mm aperture for satellite quantum key distribution (right) and its scheme (left) [27]. CPL—collimator; POL—polarizer; PBS—polarizing beam splitter cube; HWP—half wave plate; BS—50:50 beam splitter cube; IF1—interference filter (central wavelength 850 nm; half-width 10 nm); HWP-m—motorized half-wave plate; BE—beam expander; DM is a dichroic mirror; L1—lens; FSM—fast steering mirror; RCT—Ritchey-Chrétien telescope with an aperture of 600 mm; GT—guiding telescope; SCM—synchronisation module; CAM—camera; IF2—interference filter (central wavelength 532 nm; half-width 10 nm); 5 ch SPDs—5-channel single photon detector module; TDC—time-to-digital converter; PC—personal computer; PA—polarization analyzer (decoder); APS—module for analysis and processing of optical signals.
Entropy 25 00670 g0a2
The ground station with the 300 mm aperture has a smaller design than the fixed station (Figure A2), because of the compact optical parts of the optical signal analysis and processing module and quantum state decoder. The optical elements for the portable ground station were picked up with a higher transmittance in comparison to the fixed one. In particular, a hard-coated bandpass filter (>90% transmission; Figure A1—IF1) was used instead of the laser line filter (60% transmission; Figure A2—IF1). As a result, the ground station’s optical efficiency was improved.

References

  1. Scarani, V.; Bechmann-Pasquinucci, H.; Cerf, N.J.; Dušek, M.; Lütkenhaus, N.; Peev, M. The security of practical quantum key distribution. Rev. Mod. Phys. 2009, 81, 1301–1350. [Google Scholar] [CrossRef]
  2. Bennett, C.H.; Brassard, G. Quantum cryptography: Public key distribution and coin tossing. Theor. Comput. Sci. 2014, 560, 7–11. [Google Scholar] [CrossRef]
  3. Bloom, Y.; Fields, I.; Maslennikov, A.; Rozenman, G.G. Quantum Cryptography—A Simplified Undergraduate Experiment and Simulation. Physics 2022, 4, 104–123. [Google Scholar] [CrossRef]
  4. Chen, J.P.; Zhang, C.; Liu, Y.; Jiang, C.; Zhang, W.; Hu, X.L.; Guan, J.Y.; Yu, Z.W.; Xu, H.; Lin, J.; et al. Sending-or-Not-Sending with Independent Lasers: Secure Twin-Field Quantum Key Distribution over 509 km. Phys. Rev. Lett. 2020, 124, 070501. [Google Scholar] [CrossRef] [PubMed]
  5. Pittaluga, M.; Minder, M.; Lucamarini, M.; Sanzaro, M.; Woodward, R.; Li, M.J.; Yuan, Z.; Shields, A. 600-km repeater-like quantum communications with dual-band stabilization. Nat. Photonics 2021, 15, 530–535. [Google Scholar] [CrossRef]
  6. Ursin, R.; Tiefenbacher, F.; Schmitt-Manderbach, T.; Weier, H.; Scheidl, T.; Lindenthal, M.; Blauensteiner, B.; Jennewein, T.; Perdigues, J.; Trojek, P.; et al. Entanglement-based quantum communication over 144 km. Nat. Phys. 2007, 3, 481–486. [Google Scholar] [CrossRef]
  7. Bisztray, T.; Bacsardi, L. The Evolution of Free-Space Quantum Key Distribution. Infocommun. J. 2018, 10, 22–30. [Google Scholar] [CrossRef]
  8. Bourgoin, J.; Meyer-Scott, E.; Higgins, B.L.; Helou, B.; Erven, C.; Huebel, H.; Kumar, B.; Hudson, D.; D’Souza, I.; Girard, R.; et al. A comprehensive design and performance analysis of low Earth orbit satellite quantum communication. New J. Phys. 2013, 15, 023006. [Google Scholar] [CrossRef]
  9. Bedington, R.; Mantilla, J.; Ling, A. Progress in satellite quantum key distribution. Npj Quantum Inf. 2017, 3, 30. [Google Scholar] [CrossRef]
  10. Toyoshima, M.; Takenaka, H.; Shoji, Y.; Takayama, Y.; Koyama, Y.; Kunimori, H. Polarization measurements through space-to-ground atmospheric propagation paths by using a highly polarized laser source in space. Opt. Express 2009, 17, 22333–22340. [Google Scholar] [CrossRef]
  11. Takenaka, H.; Carrasco-Casado, A.; Fujiwara, M.; Kitamura, M.; Sasaki, M.; Toyoshima, M. Satellite-to-ground quantum-limited communication using a 50-kg-class microsatellite. Nat. Photonics 2017, 11, 502–508. [Google Scholar] [CrossRef]
  12. Günthner, K.; Khan, I.; Elser, D.; Stiller, B.; Bayraktar, Ö.; Müller, C.; Saucke, K.; Tröndle, D.; Heine, R.; Seel, S.; et al. Quantum-limited measurements of optical signals from a geostationary satellite. Optica 2017, 4, 611–616. [Google Scholar] [CrossRef]
  13. Shengkai, L.; Cai, W.Q.; Liu, W.; Zhang, L.; Li, Y.; Wang, J.; Yin, J.; Shen, Q.; Cao, Y.; Li, Z.P.; et al. Satellite-to-ground quantum key distribution. Nature 2017, 549, 43–47. [Google Scholar] [CrossRef]
  14. Lu, C.Y.; Cao, Y.; Peng, C.Z.; Pan, J.W. Micius quantum experiments in space. Rev. Mod. Phys. 2022, 94, 035001. [Google Scholar] [CrossRef]
  15. Chen, Y.A.; Zhang, Q.; Chen, T.; Cai, W.Q.; Shengkai, L.; Zhang, J.; Chen, K.; Yin, J.; Wang, J.; Chen, Z.; et al. An integrated space-to-ground quantum communication network over 4,600 kilometres. Nature 2021, 589, 214–219. [Google Scholar] [CrossRef]
  16. Kurochkin, V.L.; Khmelev, A.V.; Mayboroda, V.F.; Bakhshaliev, R.M.; Duplinsky, A.V.; Kurochkin, Y.V. Elements of satellite quantum network. In Proceedings of the International Conference on Micro- and Nano-Electronics 2021, Zvenigorod, Russia, 4–9 October 2021; Lukichev, V.F., Rudenko, K.V., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, WA, USA, 2022; Volume 12157, p. 121571T. [Google Scholar] [CrossRef]
  17. Khatri, S.; Brady, A.J.; Desporte, R.A.; Bart, M.P.; Dowling, J.P. Spooky action at a global distance: Analysis of space-based entanglement distribution for the quantum internet. Npj Quantum Inf. 2021, 7, 4. [Google Scholar] [CrossRef]
  18. Shirko, O.; Askar, S. A Novel Security Survival Model for Quantum Key Distribution Networks Enabled by Software-Defined Networking. IEEE Access 2023, 11, 21641–21654. [Google Scholar] [CrossRef]
  19. Liu, R.; Rozenman, G.G.; Kundu, N.K.; Chandra, D.; De, D. Towards the industrialisation of quantum key distribution in communication networks: A short survey. IET Quantum Commun. 2022, 3, 151–163. [Google Scholar] [CrossRef]
  20. Cao, Y.; Zhao, Y.; Wang, Q.; Zhang, J.; Ng, S.X.; Hanzo, L. The Evolution of Quantum Key Distribution Networks: On the Road to the Qinternet. IEEE Commun. Surv. Tutor. 2022, 24, 839–894. [Google Scholar] [CrossRef]
  21. De Forges de Parny, L.; Alibart, O.; Debaud, J.; Gressani, S.; Lagarrigue, A.; Martin, A.; Metrat, A.; Schiavon, M.; Troisi, T.; Diamanti, E.; et al. Satellite-based quantum information networks: Use cases, architecture, and roadmap. Commun. Phys. 2023, 6, 12. [Google Scholar] [CrossRef]
  22. Vergoossen, T.; Loarte, S.; Bedington, R.; Kuiper, H.; Ling, A. Modelling of satellite constellations for trusted node QKD networks. Acta Astronaut. 2020, 173, 164–171. [Google Scholar] [CrossRef]
  23. Ecker, S.; Liu, B.; Handsteiner, J.; Fink, M.; Rauch, D.; Steinlechner, F.; Scheidl, T.; Zeilinger, A.; Ursin, R. Strategies for achieving high key rates in satellite-based QKD. Npj Quantum Inf. 2021, 7, 5. [Google Scholar] [CrossRef]
  24. Pirandola, S. Satellite quantum communications: Fundamental bounds and practical security. Phys. Rev. Res. 2021, 3, 023130. [Google Scholar] [CrossRef]
  25. Wang, X.; Dong, C.; Zhao, S.; Liu, Y.; Liu, X.; Zhu, H. Feasibility of space-based measurement-device-independent quantum key distribution. New J. Phys. 2021, 23, 045001. [Google Scholar] [CrossRef]
  26. Sidhu, J.; Brougham, T.; Mcarthur, D.; Pousa, R.; Oi, D. Finite key effects in satellite quantum key distribution. Npj Quantum Inf. 2022, 8, 18. [Google Scholar] [CrossRef]
  27. Khmelev, A.; Duplinsky, A.; Mayboroda, V.; Bakhshaliev, R.; Balanov, M.; Kurochkin, V.; Kurochkin, Y. Recording of a Single-Photon Signal from Low-Flying Satellites for Satellite Quantum Key Distribution. Tech. Phys. Lett. 2021, 47, 858–861. [Google Scholar] [CrossRef]
  28. Ma, X.; Qi, B.; Zhao, Y.; Lo, H.K. Practical decoy state for quantum key distribution. Phys. Rev. A 2005, 72, 012326. [Google Scholar] [CrossRef]
  29. Hwang, W.Y. Quantum Key Distribution with High Loss: Toward Global Secure Communication. Phys. Rev. Lett. 2003, 91, 057901. [Google Scholar] [CrossRef]
  30. Vallado, D.; Crawford, P.; Hujsak, R.; Kelso, T. Revisiting Spacetrack Report# 3: Rev. In Proceedings of the AIAA/AAS Astrodynamics Specialist Conference and Exhibi, Keystone, CO, USA, 21–24 August 2006. [Google Scholar] [CrossRef]
  31. Heavens-Above. Available online: https://www.heavens-above.com/ (accessed on 10 March 2023).
  32. Space-Track.org. Available online: https://www.space-track.org/ (accessed on 10 March 2023).
  33. Seyedi, Y.; Rahimi, F. A Trace-Time Framework for Prediction of Elevation Angle Over Land Mobile LEO Satellites Networks. Wirel. Pers. Commun. 2012, 62, 793–804. [Google Scholar] [CrossRef]
  34. Ali, I.; Al-Dhahir, N.; Hershey, J.E. Doppler characterization for LEO satellites. IEEE Trans. Commun. 1998, 46, 309–313. [Google Scholar] [CrossRef]
  35. Karttunen, H.; Kröger, P.; Oja, H.; Poutanen, M.; Donner, K.J. (Eds.) Photometric Concepts and Magnitudes. In Fundamental Astronomy; Springer New York: New York, NY, USA, 1987; pp. 87–102. [Google Scholar] [CrossRef]
  36. Young, A.T.; Irvine, W.M. Multicolor photoelectric photometry of the brighter planets. I. Program and Procedure. Astron. J. 1967, 72, 945. [Google Scholar] [CrossRef]
  37. Khmelev, A.V.; Duplinsky, A.V.; Kurochkin, V.L.; Kurochkin, Y.V. Stellar calibration of the single-photon receiver for satellite-to-ground quantum key distribution. J. Phys. Conf. Ser. 2021, 2086, 012137. [Google Scholar] [CrossRef]
  38. Duplinsky, A.; Khmelev, A.; Merzlinkin, V.; Pismeniuk, L.; Bakhshaliev, R.; Tikhonova, K.; Nesterov, I.; Kurochkin, V.; Kurochkin, Y. Cubsat based scalable quantum key distribution satellite network. In Proceedings of the The XIV International Conference on Computer-Aided Technologies in Applied Mathematics, Listvyanka, Russia, 19–24 September 2022; p. 64. Available online: https://www.elibrary.ru/item.asp?id=49449483&pff=1 (accessed on 10 March 2023).
  39. Lo, H.K.; Ma, X.; Chen, K. Decoy State Quantum Key Distribution. Phys. Rev. Lett. 2005, 94, 230504. [Google Scholar] [CrossRef] [PubMed]
  40. Shor, P.W.; Preskill, J. Simple Proof of Security of the BB84 Quantum Key Distribution Protocol. Phys. Rev. Lett. 2000, 85, 441–444. [Google Scholar] [CrossRef]
  41. Kiktenko, E.O.; Trushechkin, A.S.; Lim, C.C.W.; Kurochkin, Y.V.; Fedorov, A.K. Symmetric Blind Information Reconciliation for Quantum Key Distribution. Phys. Rev. Appl. 2017, 8, 044017. [Google Scholar] [CrossRef]
  42. Zhang, Z.; Zhao, Q.; Razavi, M.; Ma, X. Improved key-rate bounds for practical decoy-state quantum-key-distribution systems. Phys. Rev. A 2017, 95, 012333. [Google Scholar] [CrossRef]
  43. SIMBAD Astronomical Database. Available online: https://simbad.unistra.fr/simbad/ (accessed on 10 March 2023).
Figure 1. Schematic of the Earth-satellite geometry. The satellite passing over the ground station through the zenith (i.e., θ E l m a x = 90 °) has an orbit radius of R = R E + h , where R E is the Earth’s radius (6364 km) for the Zvenigorod location and h is the altitude of the satellite orbit. The distance between the satellite and the ground station is determined solely by an elevation angle θ E l , that is, the angle between the horizontal line of sight and the direction of the satellite from the site of the ground station. The central angle between a direction on the ground station and a direction on the satellite from the Earth’s center is denoted by α .
Figure 1. Schematic of the Earth-satellite geometry. The satellite passing over the ground station through the zenith (i.e., θ E l m a x = 90 °) has an orbit radius of R = R E + h , where R E is the Earth’s radius (6364 km) for the Zvenigorod location and h is the altitude of the satellite orbit. The distance between the satellite and the ground station is determined solely by an elevation angle θ E l , that is, the angle between the horizontal line of sight and the direction of the satellite from the site of the ground station. The central angle between a direction on the ground station and a direction on the satellite from the Earth’s center is denoted by α .
Entropy 25 00670 g001
Figure 2. Calculated characteristics of the satellite passages. Dashed curves: elevation angle for the actual passages using the satellite TLE data and for the limit zenith passage with a maximum elevation angle of 90°. Solid curves: communication channel length for the actual passages and for the satellite passage through the zenith. The satellite passages are symmetrically aligned with respect to the highest elevation angle in each one that corresponds to zero time.
Figure 2. Calculated characteristics of the satellite passages. Dashed curves: elevation angle for the actual passages using the satellite TLE data and for the limit zenith passage with a maximum elevation angle of 90°. Solid curves: communication channel length for the actual passages and for the satellite passage through the zenith. The satellite passages are symmetrically aligned with respect to the highest elevation angle in each one that corresponds to zero time.
Entropy 25 00670 g002
Figure 3. The experimental data and the trend lines of the relationship between atmospheric extinction increase and air mass as estimated by stars’ fluxes observations on clear and foggy nights. The raw data are listed in Table A1 and were acquired with the 600-mm ground station.
Figure 3. The experimental data and the trend lines of the relationship between atmospheric extinction increase and air mass as estimated by stars’ fluxes observations on clear and foggy nights. The raw data are listed in Table A1 and were acquired with the 600-mm ground station.
Entropy 25 00670 g003
Figure 4. Link losses calculated for real passes of the Micius satellite over Zvenigorod observatory, with max elevation angles of 32.5° (31 October 2021; black curves), 57.5° (9 March 2022; red curves) and 42° (23 March 2022; blue curves). The yellow curves represent theoretical zenith passage with a maximum elevation angle of 90°. Solid lines exhibit the 600 mm telescope data; dashed lines—300 mm. The computations are carried out for elevation angles greater than 20°. We assumed good weather conditions, that is, ϰ = 0.23 .
Figure 4. Link losses calculated for real passes of the Micius satellite over Zvenigorod observatory, with max elevation angles of 32.5° (31 October 2021; black curves), 57.5° (9 March 2022; red curves) and 42° (23 March 2022; blue curves). The yellow curves represent theoretical zenith passage with a maximum elevation angle of 90°. Solid lines exhibit the 600 mm telescope data; dashed lines—300 mm. The computations are carried out for elevation angles greater than 20°. We assumed good weather conditions, that is, ϰ = 0.23 .
Entropy 25 00670 g004
Figure 5. Modelled sifted key rate R s i f t and QBER of signal states during the satellite passage above a ground station on 31 October 2021 and during the zenith satellite passage. The key parameters of the model are listed in Table 3. The computations were carried out for clear weather conditions and for elevation angles greater than 20°.
Figure 5. Modelled sifted key rate R s i f t and QBER of signal states during the satellite passage above a ground station on 31 October 2021 and during the zenith satellite passage. The key parameters of the model are listed in Table 3. The computations were carried out for clear weather conditions and for elevation angles greater than 20°.
Entropy 25 00670 g005
Figure 6. Comparison of peak sifted key rates simulated by a semi-empirical model with Vergoossen’s model data as a function of the shortest distance between satellite and ground stations. The experimental data collected by the Micius satellite for specific satellite passages are presented [13].
Figure 6. Comparison of peak sifted key rates simulated by a semi-empirical model with Vergoossen’s model data as a function of the shortest distance between satellite and ground stations. The experimental data collected by the Micius satellite for specific satellite passages are presented [13].
Entropy 25 00670 g006
Table 1. The orbital parameters of the satellite.
Table 1. The orbital parameters of the satellite.
R (km)i (Degrees) ω E (Rad/s) ω S (Rad/s)
686497.3 7.3 × 10 5 111.4 × 10 5
Table 2. Experimental setup parameters for modeling channel between satellite and ground stations [27,37,38].
Table 2. Experimental setup parameters for modeling channel between satellite and ground stations [27,37,38].
D T  (mm) ε η opt η det γ  (Rad) ϰ
3000.810.440.55 10 5 0.23–0.41
6000.730.270.55 10 5 0.23–0.41
Table 3. Key parameters of the QKD experimental setups.
Table 3. Key parameters of the QKD experimental setups.
μ p s f (Hz) Y 0 300 (/pulse) e det 300
0.80.5 10 8 2.5 × 10 6 0.9 %
ν p d p v Y 0 600  (/pulse) e det 600
0.10.250.25 5 × 10 6 0.5 %
Note: The superscript of Y 0 and e d e t denotes the aperture of a ground station.
Table 4. Sifted key length and secret key length were simulated for the satellite-to-ground QKD with 300 mm/600 mm ground stations for four typical satellite passages. The communication time of Δ t corresponds to elevation angles greater than 20°.
Table 4. Sifted key length and secret key length were simulated for the satellite-to-ground QKD with 300 mm/600 mm ground stations for four typical satellite passages. The communication time of Δ t corresponds to elevation angles greater than 20°.
PassageSifted Key Length, l sift (kbit)Secret Key Length, l sec (kbit)
2021-10-31
( θ E l m a x = 32.5 °, Δ t = 221 s)
547/120740/182
2022-03-23
( θ E l m a x = 42 °, Δ t = 260 s)
859/1897101/357
2022-03-09
( θ E l m a x = 57.5 °, Δ t = 283 s)
1263/2790187/600
Zenith passage
( θ E l m a x = 90 °, Δ t = 285 s)
1586/3505260/803
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.

Share and Cite

MDPI and ACS Style

Khmelev, A.V.; Ivchenko, E.I.; Miller, A.V.; Duplinsky, A.V.; Kurochkin, V.L.; Kurochkin, Y.V. Semi-Empirical Satellite-to-Ground Quantum Key Distribution Model for Realistic Receivers. Entropy 2023, 25, 670. https://doi.org/10.3390/e25040670

AMA Style

Khmelev AV, Ivchenko EI, Miller AV, Duplinsky AV, Kurochkin VL, Kurochkin YV. Semi-Empirical Satellite-to-Ground Quantum Key Distribution Model for Realistic Receivers. Entropy. 2023; 25(4):670. https://doi.org/10.3390/e25040670

Chicago/Turabian Style

Khmelev, Aleksandr V., Egor I. Ivchenko, Alexander V. Miller, Alexey V. Duplinsky, Vladimir L. Kurochkin, and Yury V. Kurochkin. 2023. "Semi-Empirical Satellite-to-Ground Quantum Key Distribution Model for Realistic Receivers" Entropy 25, no. 4: 670. https://doi.org/10.3390/e25040670

APA Style

Khmelev, A. V., Ivchenko, E. I., Miller, A. V., Duplinsky, A. V., Kurochkin, V. L., & Kurochkin, Y. V. (2023). Semi-Empirical Satellite-to-Ground Quantum Key Distribution Model for Realistic Receivers. Entropy, 25(4), 670. https://doi.org/10.3390/e25040670

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop