Open AccessArticle

Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer

Huicun Yu

^1,2

Bangying Tang

³,

Haolin Ding

⁴

Yang Xue

⁵,

Jie Tang

¹,

Xingyu Wang

Bo Liu

^2,*

and

Lei Shi

^1,*

Information and Navigation College, Air Force Engineering University, Xi’an 710077, China

College of Advanced Interdisciplinary Studies, National University of Defense Technology, Changsha 410073, China

College of Computer and Science, National University of Defense Technology, Changsha 410073, China

⁴

College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China

⁵

Academy of Military Sciences, Beijing 100864, China

Authors to whom correspondence should be addressed.

Entropy 2023, 25(3), 472; https://doi.org/10.3390/e25030472

Submission received: 30 September 2022 / Revised: 27 February 2023 / Accepted: 6 March 2023 / Published: 8 March 2023

(This article belongs to the Collection Quantum Information)

Download

Browse Figures

Figure 1
Schematic diagram of airborne QKD downlink. "> Figure 2
Two-dimensional schematic diagram of the ray tracing method. "> Figure 3
The sectional drawing of refractive index field of (a) the airfoil boundary layer and (b) the supersonic boundary layer. The coordinate axis represents the size of the boundary layer, in millimeters. The incident points are respectively (x = 340 mm, y = 0) and (x = 9 mm, y = 0). And the color represents the refractive index distribution. "> Figure 4
(a) The deflection angle of transmitted photons. (b) The drifted offset of the beam, which reaches to the ground station. The speed of airborne is 0.7 Ma. "> Figure 5
(a) The deflection angle of transmitted photons. (b) The drifted offset of the beam, which via the supersonic boundary layer and reaches the ground station. "> Figure 6
(a) The photons distribution probability when the photons via the supersonic BL or not, which are denoted by a red triangle and blue point respectively. (b) The photons distribution probability when the photons via the supersonic BL. The azimuth angle is 10°. "> Figure 7
The receiving photons distribution probability ratio over the azimuth angle, when the photon via the supersonic BL. "> Figure 8
(a) The receiving photons distribution when the azimuth angle is 10° with and without the effect of supersonic BL, which are denoted by a red triangle and blue point respectively. (b) The receiving photons distribution when the azimuth angle is 20° with and without the effect of supersonic BL. The green dotted line indicates the aperture of the receiving telescope. "> Figure 9
The photons distribution probability with different azimuth when the photon via the boundary layer or not. (a–c) show the photons distribution probability of the Gaussian beam without BL. (d–f) show the photons distribution probability of the Gaussian beam that passes through the BL. (a,d) show the situation when the azimuth angle is 10°. (b,e) show the situation that the azimuth angle is 30°. (c,f) show the situation the azimuth angle is 50°. "> Figure 10
The secure key rate over the azimuth angle when the photon via subsonic BL, supersonic BL, or not. ">

Versions Notes

Abstract

Airborne quantum key distribution (QKD) that can synergize with terrestrial networks and quantum satellite nodes is expected to provide flexible and relay links for the large-scale integrated communication network. However, the photon transmission rate would be randomly reduced, owing to the random distributed boundary layer that surrounding to the surface of the aircraft when the flight speed larger than Mach 0.3. Here, we investigate the airborne QKD performance with the BL effects. Furthermore, we take experimental data of supersonic BL into the model and compare the airborne QKD performance under different conditions. Simulation results show that, owing to the complex small-scale turbulence structures in the supersonic boundary layer, the deflection angle and correspondingly drifted offset of the beam varied obviously and randomly, and the distribution probability of photons are redistributed. And the subsonic and supersonic boundary layer would decrease ~35.8% and ~62.5% of the secure key rate respectively. Our work provides a theoretical guidance towards a possible realization of high-speed airborne QKD.

Keywords:

boundary layer; quantum key distribution; airborne

1. Introduction

Quantum key distribution (QKD) is a promising technology to realize secure transmission, the fundamental principles of quantum mechanics guarantee the safety of QKD, which ensures the generation of information-theoretical-secure keys for distant users [1,2,3,4,5]. Furthermore, QKD can resist eavesdropping detection and tamper relying on physical principles, which has tremendous utilization potential in terms of finance, government, and military [6,7]. Up to now, remarkable progress has been achieved both in fiber links [8,9,10,11] and quantum satellites [12,13,14,15,16,17,18,19] QKD systems, and gradually transferred from the laboratory to realistic applications. However, with constant orbits, limited communication time window and night-only quantum satellites, airborne [20,21,22,23,24,25,26,27] platforms are ideal mobile nodes that can synergize with terrestrial links and quantum satellites to build the mobile, on-demand, and real-time coverage quantum network. For example, airborne quantum nodes could serve as temporary relays to solve the last-mile quantum key exchange for an inner-city or a field network benefiting from their rapid deployment capabilities.

Several groups have reported their research in airborne quantum communications in recent years, the first demonstration was implemented in 2013, with the platform flying at the speed of 290 km/h and a height of 1.1 km [20].

In 2020, Nanjing University accomplished the first drone-based entanglement distribution [23] and achieved 200 m coverage and duration of 40 min, and in the next year, they expand the distance to 1000 m by means of one more drone [21] as an optical-relay. Compared to satellite quantum communication, airborne QKD features in high-speed maneuverability and suffers from the complicate atmosphere environment, including atmospheric turbulence, background noise and attitude disturbance, worse still, the boundary layer (BL) would also severely impact the quantum signal transmission. When the aircraft is flying at a high speed, usually larger than 0.3 Ma, a very thin layer of air will stick over the surface of the aircraft with high velocity, resulting in the BL, and it would introduce random disturbance to the transmitted photons, which would reduce the coupling efficiency and fidelity of quantum states [28], and decrease the performance of aircraft-based QKD. Currently, most implementations of airborne QKD only analyzed the noise from atmospheric turbulence and molecular scattering [25,26], but ignored the boundary layer effects. In 2021, our group proposed an airborne QKD performance evaluation scheme that takes the boundary layer effect into account. The analysis results show that the aerodynamic optical effects caused by the boundary layer should be attended to, which will greatly reduce the final key rate [29].

In this article, an airborne QKD performance evaluation scheme with supersonic BL effect is proposed. We first simulate the density and refractive index field of the BL. Then, the photon propagation model in the boundary layer is established by performing the ray tracing method that uses the Adams linear multistep method. In the end, we present the analysis results about the key rate of photons transmitted through the measured supersonic BL. Owing to the complex small-scale turbulence structures in the supersonic boundary layer, the deflection angle and correspondingly drifted offset of the beam varied obviously and randomly, and the distribution probability of photons is redistributed. The result shows that the subsonic and supersonic boundary layer would decrease ~35.8% and ~62.5% of the secure key rate respectively. With the increase of speed v, the key-rate curve is obviously jittering, and the QKD performance is continuously reduced. Our work provides a theoretical guidance towards a possible realization of high-speed airborne QKD.

2. Method

2.1. Background

In the air-to-ground communication scenario of airborne QKD, we assume

d

is the shortest projection distance to the horizon plane, between the aircraft (Alice) and the receiving ground station (Bob). Alice flies at a constant speed in a certain path obliquely above Bob, as shown in Figure 1.

h

is the relative height of Alice and Bob which is constant.

l

is the distance between Alice and Bob.

φ

is zenith angle. During the entire flight, the azimuth angle

α

of the aircraft is varying with a range of [−80°, 80°].

2.2. Photon Scattering under the Boundary Layer

Given the density field distribution of the boundary layer, the refractive index distribution of the flow field can be calculated by

n = 1 + K_{G D} \cdot ρ,

(1)

where ρ is the density, and n is the refractive index of BL. K_GD is the Gladstone-Dale constent with unit of m³/kg. K_GD is usually calculated by the empirical formula

K_{G D} = 2.23 \times 10^{- 4} (1 + \frac{7.52 \times 10^{- 3}}{λ^{2}}),

(2)

where λ is the wavelengthof photons with unit of micron. With the settings in Table 1, the G-D coefficient is 2.237 × 10⁻⁴ m³/kg.

The scattered photon path through the boundary layer can be calculated by ray tracing methods [30]. With the varying refractive index, the boundary layer can be divided into several sufficiently small cube units. And the internal refractive index of each unit

(i, j, k)

is uniform which is

n (i, j, k)

. Then, the photon will be refracted at the boundary of the two grids, as shown in Figure 2, and the new direction of communication can be calculated by Snell law

θ_{2} = \arcsin [\frac{n (x_{1}, y_{1}, z_{1}) \sin (θ_{1})}{n (x_{2}, y_{2}, z_{2})}],

(3)

where

n (x_{1}, y_{1}, z_{1})

is the refractive index of unit

(x_{1}, y_{1}, z_{1})

θ_{1}

is the incident angle,

θ_{2}

is the refraction angle. In the air-to-ground communication scenario, the incident angle

θ_{1}

is equal to zenith angle

φ

which can be calculated as

φ = \arctan (\frac{d}{h \cos α}) .

(4)

When the azimuths angle is

0 °

, the incident angle

θ_{1} = φ = \arctan (d / h)

. It is worth noting that when the azimuth angle is

0 °

, the incident angle is not

0 °

And the new incident angle can be calculated by

θ_{2}

when the photon propagated to the next grid. Then, recording the whole location

(x_{1}, y_{1}, z_{1}), (x_{2}, y_{2}, z_{2}) \dots (x_{n}, y_{n}, z_{n})

by reciprocating the Snell law, the scattered photon path X via the boundary layer can be calculated.

Then, the wavefront information could be calculated. Optical path length (OPL) is calculated by integrating the refractive index n along the propagation path X [30].

OPL = \int_{X} n d x .

(5)

And the optical path difference (OPD) is defined as follows,

OPD = OPL - 〈 OPL 〉 .

(6)

Afterward, assumed that the photon source conforms to the distribution of the Gaussian beam, the normalized intensity of the Gaussian beam can be expressed as

I (r, l) = \frac{2}{π \cdot ω_{L P}} \exp (\frac{- 2 r^{2}}{ω_{L P}^{2}}),

(7)

where r is the radial distance from the center axis of the beam, l is the axial distance from the beam’s focus, and

ω_{L P}

is the effective beam waist of the downlink signal at the ground station [31]

ω_{L P} = \sqrt{ω_{L}^{2} + {(σ_{T} \cdot l)}^{2}} .

(8)

Here, we assumed that the Transmitter pointing precision

σ_{T}

= 150 μrad that is given in [21].

ω_{L}

is the beam waist at the Ground station prior to pointing errors:

ω_{L} = l \frac{λ}{π \cdot ω_{0}} {[1 + 0.83 \cdot \sec (φ) {(\frac{D_{T}}{r_{0}})}^{5 / 3}]}^{3 / 5} .

(9)

where r₀ = 0.4 m is the Fried parameter in zenith, and

φ

is zenith angle as shown in Figure 1. ω₀ is the waist radius of the Gaussian beam

ω_{0} = 0.316 D_{T} .

(10)

The angled brackets denote the spatial average over the optical aperture.

2.3. Transmission Efficiency Analysis

Assuming that the ATP (Acquisition, Tracking, and Pointing) technique is perfect, the center axis of the beam can be aligned with the center of the receiver telescope. And then, the transmission efficiency from aircraft to the ground station can be calculated when the light is via the boundary layer or not.

Without the boundary layer effects, the detected events fulfil the Gaussian distribution. When the beam illuminates the receiving telescope, the transmission efficiency η₀ can be calculated by geometric optics as [31]

η_{0} = [1 - \exp [- 0.5 {(\frac{D_{R}}{ω_{L P}})}^{2}]] \cdot \exp [- β \cdot \sec (φ)],

(11)

where D_R is the diameter of the receiving telescope, and β is the extinction optical thickness between sea level and altitude.

After passing through the boundary layer, the intensity of the beam will be redistributed. Then, by weighting each beam with the Gaussian distribution in Equation (7), the photons distribution probability eventually reaches the receiver telescope will be calculated, when the light is via the boundary layer or not. Then, the photons distribution probability ratio

E_{b o u n d a r y}

can be calculated as

E_{b o u n d a r y} = \frac{ε_{b o u n d a r y}}{ε_{0}},

(12)

where

ε_{b o u n d a r y}

is the statistics of photons distribution probability within the range of receiving telescope per unit time with the effect of BL. While

ε_{0}

is the value without the effect of BL.

After passing through the boundary layer, the transmission efficiency

η_{b o u n d a r y}

can be calculated as

η_{b o u n d a r y} = η_{0} \cdot SR \cdot E_{b o u n d a r y},

(13)

where SR is the Strehl ratio [32]

SR \approx \exp [- {(\frac{2 π {OPD}_{rms}}{λ})}^{2}] .

(14)

where

{OPD}_{rms}

is the RMS of the

OPD

on the optical aperture.

2.4. Secure Key Rate Estimation

The decoy state method [33], as an important weapon to combat photon number splitting attack, is proposed to use a weakly coherent light source in the QKD protocol to replace the ideal single-photon source that cannot be achieved at present. Therefore, the Vacuum + weak BB84 protocol is selected in the airborne QKD system, and a formula for secure key rate is

R \geq q {Q_{1} [1 - H_{2} (e_{1})] - Q_{μ} f (E_{μ}) H_{2} (E_{μ})}

(15)

where Q₁ is the gain of single-photon states, e₁ is the error rate of single-photon states, f(x) is the bidirectional error correction efficiency as a function of error rate, H₂(x) is the binary Shannon information function and μ is the intensity of the signal state. Qμ and Eμ respectively represent the gain of signal states and the overall quantum bit error rate.

In free-space quantum communication, it is necessary to consider the reduction of efficiency caused by the diffusion of the light spot at the receiving terminal. With the influence of the boundary layer, the total transmission efficiency is

η = η_{b o u n d a r y} η_{s} η_{d} .

(16)

where η_s is receiving optical module efficiency, and η_d is detector efficiency. Substituting Equations (16) into the formula about Vacuum + weak BB84 protocol in literature [33], and the secure key rate R can be calculated.

3. Evaluation Result and Discussion

Here, the typical airfoil that named “NACA0015” is chosen for the performance analysis of our specified air-to-ground QKD system. The specific parameter settings of the aircraft, source, ground station and protocols are shown in Table 1. The boundary layer will be generated around the airfoil and its density field distribution can be simulated by the computational fluid dynamics software (Ansys Fluent, Canonsburg, PA, USA) as shown in Figure 3a, with the parameters shown in Table 1. Usually, the boundary layer thickness based on the velocity boundary layer concept that the region in which flow adjusts from zero velocity at the wall to a maximum in the main stream of the flow is termed the boundary layer. According to the simulation results, the airfoil velocity boundary layer thickness is about 10 mm. Here, in order to ensure that the path lengths in the boundary layer are approximately equal and the results are more accurate under different incident angles, we expend the concept of the boundary layer as the field that caused by the motion of aircraft. So, the generalized boundary layer thickness is 400 mm.

Literature [34] used the nano-particle-based planar laser scattering technique to measure the density distribution of the supersonic (Ma = 3.0) turbulent boundary layers, as shown in Figure 3b. The velocity boundary layer thickness is about 10 mm [33], and the generalized boundary layer thickness is about 13 mm. Although this kind of boundary layer is not the actual aircraft boundary layer, the turbulence structure is similar to the actual situation. Due to the limited amount of data, we only analyze the results at a certain moment, without considering an average over enough runs. It is also meaningful to research the influence of this kind of boundary layer on QKD performance.

When the photon trajectories of different incident angles pass through the boundary layer and reached the ground station, the evaluated deflection angle and the drifted offset of the beam are shown in Figure 2. The deflection angle can be calculated as

| θ_{1} - θ_{n} |

which is absolute value of the angle difference between the incident angle

θ_{1}

and the last refracted angle

θ_{n}

when the photon propagated via the boundary layer. In both cases, the deflection angle and correspondingly drifted offset of the beam were varied obviously in Figure 4 and Figure 5. Because of the complex small-scale turbulence structures in the supersonic boundary layer, the spatial distribution of the density field is anisotropic and random, and the results vary randomly as shown in Figure 5. When the azimuth angle is near zero, the incident photon will pass through the region with a large refractive index gradient, as shown in Figure 3b, which could influence the deflection angle obviously. In a practical airborne quantum communication system, owing to the ATP technique, the center axis of the beam can always be aligned with the center of the receiver telescope. Therefore, the deflection angle and the drifted offset were not analyzed emphatically in this paper, which could been ignored.

Afterward, the photons distribution probability that eventually reaches the receiver telescope at azimuth angle is 10°,was calculated as shown in Figure 6, when the photon via the supersonic BL or not. It shows that the photons distribution probability is redistributed when the scattered photon passes through the supersonic BL. But the photons distribution probability is fundamentally invariant when the photon via the BL when flight speed is 0.7 Ma, which is not displayed.

Then, the receiving photons distribution probability ratio

E_{b o u n d a r y}

in Equation (12) over the azimuth angle can be calculated, as shown in Figure 7, when the photon via the supersonic BL. The curve random jitter is more obvious, and

E_{b o u n d a r y} > 1

at some azimuth angle which shows that photons bunching has happened. Moreover, when the refractive index gradient is less than zero, the refraction angle is larger than incident angle. In this case, the photon spot dispersion caused by the small incident angle is more serious due to the limited number of grids. As shown in Figure 4, The small incidence angle will cause the beam to separate to the two sides, intensify the diffusion of light spots, and reduce the receiving photons distribution probability ratio

E_{b o u n d a r y}

. The case of the small incident angle with the negative refractive index gradient is indicated by the black box in Figure 3b. So, it is possible that significant dispersion phenomenon occurs around the zero azimuthal angle.

When Figure 6 is converted into a top view, as shown in Figure 8, the distribution of photons can be well displayed with and without the effect of supersonic BL. By comparing the two cases where the azimuth is 10° and 20°, the probability of photon appearing in the range of the receiving telescope after passing through the boundary layer is different. When the azimuth is 10°, the photons passe through the boundary layer and are partially refracted out of the scope of the telescope. When the azimuth is 20°, some photons that are not in the scope of the telescope are refracted into the telescope after passing through the boundary layer. By contrast, the case of

E_{b o u n d a r y} > 1

can be explained.

And the photons distribution probability with several different azimuth angles can be calculated, when the photons via the supersonic BL or not. As shown in Figure 9, different from the subsonic boundary layer, the supersonic boundary layer can diffuse the probability of photons distribution. It can cause the photon of receiving terminal to be so irregular that the transmission efficiency and secure key rate decrease.

By taking the total transmission efficiency into the formula of the key rate, the QKD performance with the BL effects is evaluated and the result is shown in Figure 10. When the photon via the BL, the key rate curve drops obviously. Compared with the subsonic surface layer, the effect of the supersonic surface layer is more obvious. The curve varies even more dramatically, and the effect of supersonic BL at individual points is very obvious. Therefore, the estimated average secure key rate is around 943 bit/s and 551 bit/s respectively when the photon via the subsonic and supersonic surface layer. If there’s no boundary layer surrounding the aircraft, the estimated average secure key rate would be around 1468 bit/s. It shows that the actual supersonic BL can also have an impact on the QKD performance, and the effect is more obvious. In addition, due to the higher speed, the effective communication time between the vehicle and the ground station will be shorter, and the total key rates will be smaller. Further, although the impact of the deflection angle and the drifted offset could be ignored, some potential factors will also affect the performance of QKD which from ATP system. Therefore, it is feasible to implement QKD by supersonic aircraft, but the results would be unsatisfactory.

4. Conclusions

Airborne quantum key distribution (QKD) that can synergize with terrestrial networks and quantum satellite nodes is expected to provide flexible and relay links for the large-scale integrated communication network. However, the photon transmission rate would be randomly reduced, owing to the random distributed boundary layer that surrounding to the surface of the aircraft when the flight speed larger than Mach 0.3, which would change the local refractive index and energy flux density drastically. Different from satellite based implementations, airborne platforms need to consider the influence of boundary layer due to the atmosphere. In this article, an airborne QKD performance evaluation scheme with supersonic BL effect is proposed. Through modeling and analysis, owing to the complex small-scale turbulence structures in the supersonic boundary layer, the deflection angle and correspondingly drifted offset of the beam varied obviously and randomly, and the distribution probability of photons are redistributed. The result shows that the subsonic and supersonic boundary layer would decrease ~35.8% and ~62.5% of the secure key rate respectively. Our work provides a theoretical guidance towards a possible realization of high-speed airborne QKD.

Author Contributions

The model design and simulations were done by H.Y. and B.T. The experimental data of supersonic boundary layer was obtained by H.D. The background was analysed by Y.X., J.T. and X.W. are responsible for the review. The effort was conceived and supervised by B.L. and co-supervised by L.S., H.Y. and B.L. wrote the draft, and all the authors reviewed the script. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China under Grant No. 61971436, No. 61972410 and No. 61803382, the Natural Science Basic Research Plan in Shaanxi Province of China (No. 2018020JQ6020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Arbekov, I.M.; Molotkov, S.N. Secret keys agreement in communication networks with quantum key distribution and trusted nodes. Laser Phys. Lett. 2020, 17, 055202. [Google Scholar] [CrossRef]
Erven, C.; Heim, B.; Meyer-Scott, E.; Bourgoin, J.P.; Laflamme, R.; Weihs, G.; Jennewein, T. Studying free-space transmission statistics and improving free-space quantum key distribution in the turbulent atmosphere. New J. Phys. 2012, 14, 123018. [Google Scholar] [CrossRef]
Kaur, E.; Wilde, M.M.; Winter, A. Fundamental limits on key rates in device-independent quantum key distribution. New J. Phys. 2020, 22, 023039. [Google Scholar] [CrossRef]
Lütkenhaus, N.; Shields, A.J. Focus on Quantum Cryptography: Theory and Practice. New J. Phys. 2009, 11, 045005. [Google Scholar] [CrossRef]
Stucki, D.; Legré, M.; Buntschu, F.; Clausen, B.; Felber, N.; Gisin, N.; Henzen, L.; Junod, P.; Litzistorf, G.; Monbaron, P.; et al. Long-term performance of the SwissQuantum quantum key distribution network in a field environment. New J. Phys. 2011, 13, 123001. [Google Scholar] [CrossRef]
Chen, H.; Wang, J.; Tang, B.; Li, Z.; Liu, B.; Sun, S. Field demonstration of time-bin reference-frame-independent quantum key distribution via an intracity free-space link. Opt. Lett. 2020, 45, 3022–3025. [Google Scholar] [CrossRef]
Xue, Y.; Shi, L.; Chen, W.; Yin, Z.; Fan-Yuan, G.-j.; Fu, H.; Lu, Q.; Wei, J. Improving the performance of reference-frame-independent quantum key distribution through a turbulent atmosphere. Phys. Rev. A 2020, 102, 062602. [Google Scholar] [CrossRef]
Chen, W.; Han, Z.; Zhang, T.; Wen, H.; Yin, Z.; Xu, F.; Wu, Q.; Liu, Y.; Zhang, Y.; Mo, X.; et al. Field Experiment on a “Star Type” Metropolitan Quantum Key Distribution Network. IEEE Photonics Technol. Lett. 2009, 21, 575–577. [Google Scholar] [CrossRef]
Wang, S.; Chen, W.; Yin, Z.-Q.; He, D.-Y.; Hui, C.; Hao, P.-L.; Fan-Yuan, G.-J.; Wang, C.; Zhang, L.-J.; Kuang, J.; et al. Practical gigahertz quantum key distribution robust against channel disturbance. Opt. Lett. 2018, 43, 2030–2033. [Google Scholar] [CrossRef]
Wang, X.-B.; Yu, Z.-W.; Hu, X.-L. Twin-field quantum key distribution with large misalignment error. Phys. Rev. A 2018, 98, 062323. [Google Scholar] [CrossRef] [Green Version]
Kozlowski, W.; Wehner, S. Towards Large-Scale Quantum Networks. In Proceedings of the Sixth Annual ACM International Conference on Nanoscale Computing and Communication, Dublin, Ireland, 25–27 September 2019; p. 3. [Google Scholar]
Hughes, R.J.; Nordholt, J.E. Quantum space race heats up. Nat. Photonics 2017, 11, 456–458. [Google Scholar] [CrossRef]
Hughes, R.J.; Nordholt, J.E.; Derkacs, D.; Peterson, C.G. Practical free-space quantum key distribution over 10 km in daylight and at night. New J. Phys. 2002, 4, 43. [Google Scholar] [CrossRef] [Green Version]
Bedington, R.; Arrazola, J.M.; Ling, A. Progress in satellite quantum key distribution. npj Quantum Inf. 2017, 3, 30. [Google Scholar] [CrossRef] [Green Version]
Liao, S.-K.; Cai, W.-Q.; Handsteiner, J.; Liu, B.; Yin, J.; Zhang, L.; Rauch, D.; Fink, M.; Ren, J.-G.; Liu, W.-Y.; et al. Satellite-Relayed Intercontinental Quantum Network. Phys. Rev. Lett. 2018, 120, 030501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Schmitt-Manderbach, T.; Weier, H.; Fürst, M.; Ursin, R.; Tiefenbacher, F.; Scheidl, T.; Perdigues, J.; Sodnik, Z.; Kurtsiefer, C.; Rarity, J.G.; et al. Experimental Demonstration of Free-Space Decoy-State Quantum Key Distribution over 144 km. Phys. Rev. Lett. 2007, 98, 010504. [Google Scholar] [CrossRef] [PubMed]
Vallone, G.; Bacco, D.; Dequal, D.; Gaiarin, S.; Luceri, V.; Bianco, G.; Villoresi, P. Experimental Satellite Quantum Communications. Phys. Rev. Lett. 2015, 115, 040502. [Google Scholar] [CrossRef] [Green Version]
Yin, H.-L.; Chen, T.-Y.; Yu, Z.-W.; Liu, H.; You, L.-X.; Zhou, Y.-H.; Chen, S.-J.; Mao, Y.; Huang, M.-Q.; Zhang, W.-J.; et al. Measurement-Device-Independent Quantum Key Distribution Over a 404 km Optical Fiber. Phys. Rev. Lett. 2016, 117, 190501. [Google Scholar] [CrossRef] [Green Version]
Calderaro, L.; Agnesi, C.; Dequal, D.; Vedovato, F.; Schiavon, M.; Santamato, A.; Luceri, V.; Bianco, G.; Vallone, G.; Villoresi, P. Towards quantum communication from global navigation satellite system. Quantum Sci. Technol. 2018, 4, 015012. [Google Scholar] [CrossRef] [Green Version]
Nauerth, S.; Moll, F.; Rau, M.; Fuchs, C.; Horwath, J.; Frick, S.; Weinfurter, H. Air-to-ground quantum communication. Nat. Photonics 2013, 7, 382–386. [Google Scholar] [CrossRef]
Liu, H.-Y.; Tian, X.-H.; Gu, C.; Fan, P.; Ni, X.; Yang, R.; Zhang, J.-N.; Hu, M.; Guo, J.; Cao, X.; et al. Optical-Relayed Entanglement Distribution Using Drones as Mobile Nodes. Phys. Rev. Lett. 2021, 126, 020503. [Google Scholar] [CrossRef]
Pugh, C.J.; Kaiser, S.; Bourgoin, J.-P.; Jin, J.; Sultana, N.; Agne, S.; Anisimova, E.; Makarov, V.; Choi, E.; Higgins, B.L.; et al. Airborne demonstration of a quantum key distribution receiver payload. Quantum Sci. Technol. 2017, 2, 024009. [Google Scholar] [CrossRef] [Green Version]
Liu, H.; Tian, X.-H.; Gu, C.; Fan, P.; Ni, X.; Yang, R.; Zhang, J.-N.; Hu, M.; Guo, J.; Cao, X.; et al. Drone-based entanglement distribution towards mobile quantum networks. Natl. Sci. Rev. 2020, 7, 921–928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
DeCesare, A.; Snyder, R.; Carvalho, D.; Miller, W.; Alsing, P.; Ahn, D. Toward mobile free-space optical QKD: Characterization of a polarization-based receiver. In Proceedings of the SPIE Defense + Commercial Sensing, Virtual, 27 April–9 May 2020; Volume 11391. [Google Scholar]
Zhang, M.; Zhang, L.; Wu, J.; Yang, S.; Wan, X.; He, Z.; Jia, J.; Citrin, D.S.; Wang, J. Detection and compensation of basis deviation in satellite-to-ground quantum communications. Opt. Express 2014, 22, 9871–9886. [Google Scholar] [CrossRef]
Moll, F.; Botter, T.; Marquardt, C.; Pusey, D.; Shrestha, A.; Reeves, A.; Jaksch, K.; Gunthner, K.; Bayraktar, O.; Mueller-Hirschkorn, C.; et al. Stratospheric QKD: Feasibility analysis and free-space optics system concept. In Proceedings of the SPIE Security + Defence, Strasbourg, France, 9–12 September 2019; Volume 11167. [Google Scholar]
Xue, Y.; Chen, W.; Wang, S.; Yin, Z.; Shi, L.; Han, Z. Airborne quantum key distribution: A review. Chin Opt Lett. 2021, 19, 122702. [Google Scholar] [CrossRef]
Gordeyev, S.; Rennie, M.R.; Cain, A.B.; Hayden, T. Aero-optical measurements of high-Mach supersonic boundary layers. In Proceedings of the 46th AIAA Plasmadynamics and Lasers Conference, Dallas, TX, USA, 22–26 June 2015; p. 3246. [Google Scholar]
Yu, H.-C.; Tang, B.-Y.; Chen, H.; Xue, Y.; Tang, J.; Yu, W.-R.; Liu, B.; Shi, L. Airborne quantum key distribution with boundary layer effects. EPJ Quantum Technol. 2021, 8, 26. [Google Scholar] [CrossRef]
Guo, G.; Liu, H.; Zhang, B. Aero-optical effects of an optical seeker with a supersonic jet for hypersonic vehicles in near space. Appl. Opt. 2016, 55, 4741–4751. [Google Scholar] [CrossRef] [PubMed]
Neumann, S.P.; Joshi, S.K.; Fink, M.; Scheidl, T.; Blach, R.; Scharlemann, C.; Abouagaga, S.; Bambery, D.; Kerstel, E.; Barthelemy, M.; et al. Q3Sat: Quantum communications uplink to a 3U CubeSat—Feasibility & design. EPJ Quantum Technol. 2018, 5, 4. [Google Scholar] [CrossRef] [Green Version]
Ross, T.S. Limitations and applicability of the Maréchal approximation. Appl. Opt. 2009, 48, 1812–1818. [Google Scholar] [CrossRef]
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] [Green Version]
Ding, H.; Yi, S.; Zhu, Y.; He, L. Experimental investigation on aero-optics of supersonic turbulent boundary layers. Appl. Opt. 2017, 56, 7604–7610. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Schematic diagram of airborne QKD downlink.

Figure 2. Two-dimensional schematic diagram of the ray tracing method.

Figure 3. The sectional drawing of refractive index field of (a) the airfoil boundary layer and (b) the supersonic boundary layer. The coordinate axis represents the size of the boundary layer, in millimeters. The incident points are respectively (x = 340 mm, y = 0) and (x = 9 mm, y = 0). And the color represents the refractive index distribution.

Figure 4. (a) The deflection angle of transmitted photons. (b) The drifted offset of the beam, which reaches to the ground station. The speed of airborne is 0.7 Ma.

Figure 5. (a) The deflection angle of transmitted photons. (b) The drifted offset of the beam, which via the supersonic boundary layer and reaches the ground station.

Figure 6. (a) The photons distribution probability when the photons via the supersonic BL or not, which are denoted by a red triangle and blue point respectively. (b) The photons distribution probability when the photons via the supersonic BL. The azimuth angle is 10°.

Figure 7. The receiving photons distribution probability ratio over the azimuth angle, when the photon via the supersonic BL.

Figure 8. (a) The receiving photons distribution when the azimuth angle is 10° with and without the effect of supersonic BL, which are denoted by a red triangle and blue point respectively. (b) The receiving photons distribution when the azimuth angle is 20° with and without the effect of supersonic BL. The green dotted line indicates the aperture of the receiving telescope.

Figure 9. The photons distribution probability with different azimuth when the photon via the boundary layer or not. (a–c) show the photons distribution probability of the Gaussian beam without BL. (d–f) show the photons distribution probability of the Gaussian beam that passes through the BL. (a,d) show the situation when the azimuth angle is 10°. (b,e) show the situation that the azimuth angle is 30°. (c,f) show the situation the azimuth angle is 50°.

Figure 10. The secure key rate over the azimuth angle when the photon via subsonic BL, supersonic BL, or not.

Table 1. List of parameters and values for which we assigned fixed values.

	Symbol	Parameter	Value
Aircraft	v	Flight speed	0.7 Ma
	h	Relative flying height (h = h_T − h_R)	10 km
	ρh	Air density	0.41271 kg/m³
	d	The shortest projection distance to the ground of the aircraft and the ground station	10 km
Source	h_T	Altitude of the aircraft	11 km
	D_T	Transmitter telescope diameter	0.1 m
	δ_T	Transmitter pointing precision	2.4 μrad
	λ	Transmitter wavelength	1550 mm
	ω₀	Waist radius	0.0316 m
Ground station	h_R	Altitude of the ground station	1 km
	D_R	Receiver telescope diameter	0.3 m
	e_d	Detection error rate	1%
	p_d	Dark count	2 × 10⁻⁶
	η_d	Detector efficiency	15%
	η_s	Receiving optical module efficiency	60%
Protocols	μ	Expected photon number of signal states	0.1
	ν	Expected photon number of decoy states	0.05
	N	System repetition rate	100 MHz
	Ps	Sent Probability of signal states	50%
	Pd	Sent Probability of decoy states	25%

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.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Yu, H.; Tang, B.; Ding, H.; Xue, Y.; Tang, J.; Wang, X.; Liu, B.; Shi, L. Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer. Entropy 2023, 25, 472. https://doi.org/10.3390/e25030472

AMA Style

Yu H, Tang B, Ding H, Xue Y, Tang J, Wang X, Liu B, Shi L. Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer. Entropy. 2023; 25(3):472. https://doi.org/10.3390/e25030472

Chicago/Turabian Style

Yu, Huicun, Bangying Tang, Haolin Ding, Yang Xue, Jie Tang, Xingyu Wang, Bo Liu, and Lei Shi. 2023. "Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer" Entropy 25, no. 3: 472. https://doi.org/10.3390/e25030472

APA Style

Yu, H., Tang, B., Ding, H., Xue, Y., Tang, J., Wang, X., Liu, B., & Shi, L. (2023). Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer. Entropy, 25(3), 472. https://doi.org/10.3390/e25030472

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

Article Menu

Airborne Quantum Key Distribution Performance Analysis under Supersonic Boundary Layer

Abstract

1. Introduction

2. Method

2.1. Background

2.2. Photon Scattering under the Boundary Layer

2.3. Transmission Efficiency Analysis

2.4. Secure Key Rate Estimation

3. Evaluation Result and Discussion

4. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI