高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于GMM的幅度相位联合编码CVQKD安全性分析

赵常兰 王天一

引用本文:
Citation:

基于GMM的幅度相位联合编码CVQKD安全性分析

    通讯作者: 王天一, ty.wang@gzu.edu.cn
  • 基金项目:

    贵州省科技计划资助项目 黔科合基础-ZK一般304

    贵州大学科技计划资助项目 贵大培育56号

  • 中图分类号: O431.2

GMM-based amplitude-phase joint coding CVQKD security analysis

    Corresponding author: WANG Tianyi, ty.wang@gzu.edu.cn
  • CLC number: O431.2

  • 摘要: 为了提高离散调制连续变量量子密钥分发协议性能,采用幅度相移键控(APSK)联合调制格式方法,在接收端采用高斯混合模型分类算法识别量子态来提升系统的性能。将密钥传输系统分为状态学习和状态预测两个阶段,在状态学习阶段基于高斯混合模型的分类器对已知类别的量子态进行训练,学习不同类别量子态的幅度相位分布情况;在状态预测阶段则采用最小欧氏距离计算出待测量子态属于每个类别的后验概率,从而判定待测量子态的类别,并通过参数估计、反向协调和保密增强生成最终密钥。结果表明,在反向协调和集体攻击下128-APSK离散调制连续变量量子密钥分发协议能够有效生成安全密钥,当安全码率为10-6 bit/symbol时,传输距离可接近60 km。该研究为进一步提高离散调制连续变量量子密钥分发协议的系统性能提供了参考。
  • 图 1  M-APSK编码CVQKD协议在相空间的分布示意图

    Figure 1.  Schematic distribution of M-APSK encoded CVQKD protocol in phase space

    图 2  16-APSK调制相干态经过量子信道传输在相空间的分布

    Figure 2.  Distribution of 16-APSK modulated coherent states in phase space through quantum channel transmission

    图 3  基于高斯混合模型的离散调制CVQKD协议

    Figure 3.  Discrete modulation CVQKD protocol based on GMM

    图 4  16-APSK调制相干态的特征提取

    Figure 4.  Feature extraction of 16-APSK modulated coherent states

    图 5  不同M-APSK调制格式下识别精度与信噪比的变化曲线

    Figure 5.  Variation cur ves of recognition accuracy and SNR under different M-APSK modulation formats

    图 6  考虑有限码长效应M-APSK调制下CVQKD协议的密钥速率

    Figure 6.  Key rate of CVQKD protocol with M-APSK modulation considering finite code length effect

    图 7  不同M-APSK调制格式下系统性能与过量噪声的变化曲线

    Figure 7.  Variation curves of system performance and excess noise under different M-APSK modulation formats

    图 8  不同码长对系统性能的变化曲线

    Figure 8.  Variation curve of system performance for different code lengths

    图 9  基于GMM的M-APSK调制格式CVQKD协议密钥速率

    Figure 9.  GMM-based M-APSK modulation format CVQKD protocol key rate

    表 1  APSK调制类型的星座点

    Table 1.  Constellation points for APSK modulation types

    constellation points of 16-APSK, 32-APSK, 64-APSK, and 128-APSK
    $C_k \in R_k \exp \left[\mathrm{j}\left(\frac{2 {\rm{ \mathsf{ π} }}}{n_k} i_k+\theta_k\right)\right], $
    $\left\{{array}{l} k \in[1, 2], n_k \in[4, 12], i_k \in[0, \cdots, 11], \text { if } M=16 \\ k \in[1, 2, 3], n_k \in[4, 12, 16], i_k \in[0, \cdots, 15], \text { if } M=32 \\ k \in[1, 2, 3, 4], n_k \in[4, 12, 16, 32], i_k \in[0, \cdots, 31], \text { if } M=64 \\ k \in[1, 2, 3, 4, 5], n_k \in[4, 12, 16, 32, 64], i_k \in[0, \cdots, 63], \text { if } M=128 {array}\right.$
    下载: 导出CSV
  • [1]

    GUO Y, DING J Z, MAO Y, et al. Quantum catalysis-based discrete modulation continuous variable quantum key distribution with eight states[J]. Physics Letters, 2020, A384(12): 126340.
    [2]

    MIKHAIL E, PINAKI M. Continuous-variable quantum key distribution with discretized modulations in the strong noise regime[J]. Phy-sical Review, 2020, A101(6): 062313.
    [3]

    ZHAO H, LI H, XU Y, et al. Simple continuous-variable quantum key distribution scheme using a Sagnac-based Gaussian modulator[J]. Optics Letters, 2022, 47(12): 2939-2942. doi: 10.1364/OL.458443
    [4]

    PAN Y, WANG H, SHAO Y, et al. Experimental demonstration of high-rate discrete-modulated continuous-variable quantum key distribution system[J]. Optics Letters, 2022, 47(13): 3307-3310. doi: 10.1364/OL.456978
    [5]

    WANG H, PAN Y, SHAO Y, et al. Performance analysis for OFDM-based multi-carrier continuous-variable quantum key distribution with an arbitrary modulation protocol[J]. Optics Express, 2023, 31(4): 5577-5592. doi: 10.1364/OE.482136
    [6]

    ALMEIDA M, PEREIRA D, MUGA N J, et al. Secret key rate of multiring M-APSK continuous variable quantum key distribution. Optics Express, 2021, 29(23): 38669-38682. doi: 10.1364/OE.439992
    [7]

    WANG Q, QUAN Z, SONG T Y, et al. M-APSK constellation optimization in the presence of phase reference error[J]. IEEE Wireless Communications Letters, 2020, 9(12): 2154-2158. doi: 10.1109/LWC.2020.3015890
    [8]

    LIAO Q, XIAO G, ZHONG H, et al. Multi-label learning for improving discretely-modulated continuous-variable quantum key distribution[J]. New Journal of Physics, 2020, 22(8): 083086. doi: 10.1088/1367-2630/abab3c
    [9]

    LIU W Q, HUANG P, PENG J Y, et al. Integrating machine learning to achieve an automatic parameter prediction for practical conti-nuous-variable quantum key distribution[J]. Physical Review, 2018, A97(2): 02231.
    [10]

    LIU Z P, ZHOU M G, LIU W B, et al. Automated machine learning for secure key rate in discrete-modulated continuous-variable quantum key distribution[J]. Optics Express, 2022, 30(9): 15024-15036. doi: 10.1364/OE.455762
    [11]

    WANG X, WANG H, ZHOU C, et al. Continuous-variable quantum key distribution with low-complexity information reconciliation[J]. Optics Express, 2022, 30(17): 30455-30465. doi: 10.1364/OE.461665
    [12]

    LI J W, GUO Y, WANG X D, et al. Discrete-modulated continuous-variable quantum key distribution with a machine-learning-based detector[J]. Optical Engineering, 2018, 57(6): 066109.
    [13] 吴晓东, 黄端, 黄鹏, 等. 基于实际探测器补偿的离散调制连续变量测量设备无关量子密钥分发方案[J]. 物理学报, 2022, 71(24): 56-65.

    WU X D, HUANG D, HUANG P, et al. Device independent quantum key distribution scheme for discrete modulated continuous variable measurements based on real detector compensation[J]. Journal of Physics, 2022, 71(24): 56-65(in Chinese). 
    [14]

    XIA Y H, WANG Q W, FENG H T, et al. Finite volume effects on the QCD chiral phase transition in the finite size de-pendent Nambu-Jona-Lasinio model[J]. Chinese Physics, 2019, C43(3): 43-48.
    [15]

    XU Sh J, LI Y, WANG Y J, et al. Security analysis of a passive continuous-variable quantum key distribution by considering finite-size effect[J]. Entropy, 2021, 23(12): 1698. doi: 10.3390/e23121698
    [16] 吴晓东, 黄端. 基于非高斯态区分探测的往返式离散调制连续变量量子密钥分发方案[J]. 物理学报, 2023, 72(5): 187-198.

    WU X D, HUANG D. Round-trip discrete modulated continuous variable quantum key distribution scheme based on non-Gaussian state distinguished detection[J]. Journal of Physics, 2023, 72(5): 187-198(in Chinese). 
    [17]

    SAADEH R, BURQAN A, ElAJOU A. Reliable solutions to fractional Lane-Emden equations via Laplace transform and residual error function[J]. Alexandria Engineering Journal, 2022, 61(12): 10551-10562. doi: 10.1016/j.aej.2022.04.004
    [18] 陶玉婷. 机器学习与深度学习[M]. 北京: 电子工业出版社, 2022: 146-167.

    TAO Y T. Machine learning and deep learning[M]. Beijing: Electronic Industry Press, 2022: 146-167(in Chinese).
    [19] 西班扬·达斯, 乌米特·卡卡马克. 自动机器学习入门与实践: 使用Python[M]. 武汉: 华中科技大学出版社, 2019: 50-66.

    DAS S, CAKMAK U. Introduction to and practice of automated machine learning: Using Python[M]. Wuhan: Huazhong University of Science and Technology Press, 2019: 50-66(in Chinese).
    [20]

    PEREIRA D, ALMEIDA M, FACAO M, et al. Probabilistic shaped 128-APSK CV-QKD transmission system over optical fibres[J]. Optics Letters, 2022, 47(15): 3948-3951.
    [21]

    JAIN N, CHIN H M, MANI H, et al. Practical continuous-variable quantum key distribution with composable security[J]. Nature Communications, 2022, 13: 4740.
  • [1] 徐强沈思谢修敏吴鹏周强邓光伟王浟宋海智 . 可用于激光雷达的量子光学技术. 激光技术, 2021, 45(1): 44-47. doi: 10.7510/jgjs.issn.1001-3806.2021.01.008
    [2] 董颖娣彭进业 . 量子密钥分配系统中非线性偏振耦合的研究. 激光技术, 2015, 39(6): 756-760. doi: 10.7510/jgjs.issn.1001-3806.2015.06.005
    [3] 徐强谢修敏张伟袁菲胡卫英邓杰赵新华宋海智 . 半导体量子点量子光源研究进展. 激光技术, 2020, 44(5): 575-586. doi: 10.7510/jgjs.issn.1001-3806.2020.05.009
    [4] 金锐博杨子祥郝向英李百宏 . 量子色散消除的研究进展. 激光技术, 2022, 46(1): 69-78. doi: 10.7510/jgjs.issn.1001-3806.2022.01.006
    [5] 成然黄帅徐强张伟邓光伟周强王浟宋海智 . 铌酸锂量子器件研究进展. 激光技术, 2022, 46(6): 722-728. doi: 10.7510/jgjs.issn.1001-3806.2022.06.002
    [6] 李城敬波廖金宇陈钰洁宋日尧张天乐宋海智周强 . 通信波段稀土离子掺杂固态量子存储进展. 激光技术, 2022, 46(1): 45-57. doi: 10.7510/jgjs.issn.1001-3806.2022.01.003
    [7] 张宝华刘鹤 . 基于区域定位与轮廓分割的红外目标检测. 激光技术, 2015, 39(6): 840-844. doi: 10.7510/jgjs.issn.1001-3806.2015.06.023
    [8] 王中结李聪 . 未知原子态经纠缠腔场信道的隐形传态. 激光技术, 2013, 37(6): 812-815. doi: 10.7510/jgjs.issn.1001-3806.2013.06.024
    [9] 王叶荟唐丽倪重文是度芳 . 一维光子晶体非线性微腔的缺陷模和双稳态. 激光技术, 2006, 30(5): 462-464,468.
    [10] 王中结方旭 . 在去相环境中原子态的确定性隐形传态. 激光技术, 2015, 39(2): 270-273. doi: 10.7510/jgjs.issn.1001-3806.2015.02.026
    [11] 李佳佳丁志超汪之国肖光宗胡绍民 . 铷原子磁力仪最佳抽运光强的研究. 激光技术, 2016, 40(5): 691-694. doi: 10.7510/jgjs.issn.1001-3806.2016.05.015
    [12] 王中结方旭 . 非线性囚禁离子振动模式互关联函数的研究. 激光技术, 2015, 39(1): 109-113. doi: 10.7510/jgjs.issn.1001-3806.2015.01.022
    [13] 胡小龙胡南邹锴孟赟许亮冯一帆 . SNSPD二十年:回顾与展望. 激光技术, 2022, 46(1): 1-37. doi: 10.7510/jgjs.issn.1001-3806.2022.01.001
    [14] 刘剑黄典贺青王轶文韦联福 . 基于光子数可分辨探测器的单脉冲光子数检测. 激光技术, 2022, 46(1): 58-63. doi: 10.7510/jgjs.issn.1001-3806.2022.01.004
    [15] 王中结李聪 . 在抗退相干子空间中原子态的隐形传态. 激光技术, 2013, 37(5): 682-685. doi: 10.7510/jgjs.issn.1001-3806.2013.05.026
    [16] 申静李俊奇 . 基于光纤差分相移键控的色散补偿方案的研究. 激光技术, 2019, 43(5): 641-645. doi: 10.7510/jgjs.issn.1001-3806.2019.05.010
    [17] 张开银王树春赵丽娟黄晖张光寅许京军 . 量子物理的基础及其光学实验. 激光技术, 2001, 25(3): 232-237.
    [18] 周贤韬江英华 . 基于n粒子GHZ态和单光子混合的量子安全直接通信. 激光技术, 2022, 46(1): 79-84. doi: 10.7510/jgjs.issn.1001-3806.2022.01.007
    [19] 饶金同王家驷梁军 . 强度依赖Jaynes-Cummings模型中的量子态保真度. 激光技术, 2012, 36(5): 708-712. doi: 10.3969/j.issn.1001-3806.2012.05.034
    [20] 王东方王玉德王景武 . 基于改进的混合高斯模型的运动目标检测方法. 激光技术, 2014, 38(6): 776-779. doi: 10.7510/jgjs.issn.1001-3806.2014.06.011
  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  1851
  • HTML全文浏览量:  821
  • PDF下载量:  57
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-06-05
  • 录用日期:  2023-08-17
  • 刊出日期:  2024-05-25

基于GMM的幅度相位联合编码CVQKD安全性分析

    通讯作者: 王天一, ty.wang@gzu.edu.cn
  • 贵州大学 大数据与信息工程学院, 贵阳 550025, 中国
基金项目:  贵州省科技计划资助项目 黔科合基础-ZK一般304贵州大学科技计划资助项目 贵大培育56号

摘要: 为了提高离散调制连续变量量子密钥分发协议性能,采用幅度相移键控(APSK)联合调制格式方法,在接收端采用高斯混合模型分类算法识别量子态来提升系统的性能。将密钥传输系统分为状态学习和状态预测两个阶段,在状态学习阶段基于高斯混合模型的分类器对已知类别的量子态进行训练,学习不同类别量子态的幅度相位分布情况;在状态预测阶段则采用最小欧氏距离计算出待测量子态属于每个类别的后验概率,从而判定待测量子态的类别,并通过参数估计、反向协调和保密增强生成最终密钥。结果表明,在反向协调和集体攻击下128-APSK离散调制连续变量量子密钥分发协议能够有效生成安全密钥,当安全码率为10-6 bit/symbol时,传输距离可接近60 km。该研究为进一步提高离散调制连续变量量子密钥分发协议的系统性能提供了参考。

English Abstract

    • 量子密钥分发(quantum key distribution, QKD)能够实现两个合法用户Alice和Bob之间的无条件安全密钥分发,其安全性由量子力学的基本原理保证[1]。首个QKD协议是BB84协议,通过随机选择偏振基使得两个非正交的偏振态携带同样的信息,因此其安全性来源于非正交态不可区分[2]

      虽然BB84协议的安全性已得到严格证明,但由于单光子源难以制备且单光子探测技术实现难度大、成本较高,研究者逐渐将研究方向从单光子转向了连续变量。连续变量量子密钥分发(continuous variable quantum key distribution,CVQKD)具有更高的密钥速率和与标准通信组件的兼容性等独特的优势,自诞生以来迅速受到学术界关注[3]。在高斯调制CVQKD中,信息被编码在量化电磁场的正交上:Alice在相空间上对相干态进行调制,而Bob执行零差或外差检测以建立与Alice之间的相关性[4],这些相关性可以通过后处理过程进一步被转变成安全密钥。

      由于理论上的高斯调制在物理上无法实现,基于离散调制格式的CVQKD协议近年来吸引了越来越多的关注,其优点包括更容易的实验实现和更高的纠错效率,这些特性为量子安全网络中的大规模部署带来了潜在优势[5]。相移键控(phase-shift keying,PSK)调制是常用的CVQKD离散调制方案,但其密钥速率与高斯调制相比仍存在差距[6]。相比之下,幅度相移键控(amplitude phase-shift keying,APSK)联合调制可以获得更高的密钥速率,更加接近真实的高斯调制的结果。

      为增强CVQKD的实用性,机器学习的方法被引入CVQKD领域,用于设计新型CVQKD方案、预测参数和检测量子攻击[7]。最近一些关于使用基于机器学习的方法改进离散调制CVQKD的研究已经证明了其可行性,如LIAO等人[8]提出的多标签分类算法作为嵌入式分类器应用于离散调制CVQKD系统;LIU等人[9]通过集成支持向量回归模型,以优化离散调制CVQKD系统的性能和实际的安全性。

      本文中提出了基于高斯混合模型(Gaussian mixture model,GMM)的幅度相位联合编码CVQKD方案。在该方案中,将量子通信系统分为两个阶段,即状态学习阶段和状态预测阶段,状态学习阶段基于高斯混合模型的分类器对已知类别的量子态进行训练,学习不同类别量子态的幅度相位分布情况;状态预测阶段则采用最小欧氏距离计算出待测量子态属于每个类别的后验概率,从而判定待测量子态的类别,并通过参数估计、反向协调和保密增强生成最终密钥。由已知相干态构成的训练集进行学习,再对未知量子态提取特征进行分类进而做出判决,设计一个量子态识别分类器,更精确地预测相干态调制提高了离散调制CVQKD的系统性能。

    • M阶APSK(M-APSK)调制方式的星座图通常由多个同心圆组成,每个同心圆上均匀分布着多个量子态,这些量子态构成的信号集可以表示为:

      $ C_k=R_k \exp \left[\mathrm{j}\left(\frac{2 {\rm{ \mathsf{ π} }}}{n_k} i_k+\theta_k\right)\right] $

      (1)

      式中: Rk为第k个同心圆的半径; 2πik/nk+θk表示相空间中量子态的相位; nk为第k个同心圆上的量子态数, n1+n2+…+nk=M; θk为第k个同心圆上的量子态的初始相位; ik(ik=0, 1, …, nk-1)为第k个同心圆上的一个量子态。

      本文中考虑了{16, 32, 64, 128}阶APSK调制的星座点,如表 1所示; 绘制了16-APSK和32-APSK编码的CVQKD在相空间的示意图,如图 1所示。

      表 1  APSK调制类型的星座点

      Table 1.  Constellation points for APSK modulation types

      constellation points of 16-APSK, 32-APSK, 64-APSK, and 128-APSK
      $C_k \in R_k \exp \left[\mathrm{j}\left(\frac{2 {\rm{ \mathsf{ π} }}}{n_k} i_k+\theta_k\right)\right], $
      $\left\{{array}{l} k \in[1, 2], n_k \in[4, 12], i_k \in[0, \cdots, 11], \text { if } M=16 \\ k \in[1, 2, 3], n_k \in[4, 12, 16], i_k \in[0, \cdots, 15], \text { if } M=32 \\ k \in[1, 2, 3, 4], n_k \in[4, 12, 16, 32], i_k \in[0, \cdots, 31], \text { if } M=64 \\ k \in[1, 2, 3, 4, 5], n_k \in[4, 12, 16, 32, 64], i_k \in[0, \cdots, 63], \text { if } M=128 {array}\right.$

      图  1  M-APSK编码CVQKD协议在相空间的分布示意图

      Figure 1.  Schematic distribution of M-APSK encoded CVQKD protocol in phase space

    • 在CVQKD系统中,发送方Alice和接收方Bob的设备是可信的, 并且不被潜在的窃听者Eve所操控[10],它们之间获得的密钥速率K为:

      $ K=\frac{m}{W}\left(\beta I_{\mathrm{BA}}-\chi_{\mathrm{BE}}-\Delta m\right) $

      (2)

      式中: m为用于信息协调的数据量; W为量子信道中传输的总量子态数,此处考虑参数估计的数据,将m/W设置为1/2;β为反向协调的协调效率; IBA为Bob和Alice之间的互信息[11]; χBE为Bob和Eve互信息的Holevo界; Δm为与隐私放大安全相关的参数。假定Bob使用外差探测同时测量Alice发送的两个正交分量,则互信息IBA可由下式计算:

      $ I_{\mathrm{BA}}=\operatorname{lb}\left(1+\frac{S}{N}\right)=\operatorname{lb}\left(1+\frac{2 T \eta\langle n\rangle}{2+T \eta \xi+2 \xi_{\mathrm{th}}}\right) $

      (3)

      式中: S/N为信噪比; 〈n〉是每个符号的平均光子数; T为信道传输率; ξ为过量噪声; ξth为检测系统的热噪声; η为光电二极管的检测效率。

      在集体攻击下,Bob和Eve互信息的Holevo界[12],即χBE可由下式给出:

      $ \chi_{\mathrm{BE}}=\sum\limits_{a=1}^2 G\left(\frac{\mu_a-1}{2}\right)-\sum\limits_{a=3}^5 G\left(\frac{\mu_a-1}{2}\right) $

      (4)

      其中,

      $ G(x)=(x+1) \operatorname{lb}(x+1)-x \operatorname{lb} x $

      (5)

      式中: a=1, 2时, μa为描述Alice和Bob共享状态的协方差矩阵γAB的辛特征值; a=3, 4, 5时, μa为Bob的投影测量的协方差矩阵的辛特征值[13],即:

      $ \left\{\begin{array}{l} \boldsymbol{\gamma}_{\mathrm{AB}}=\left[\begin{array}{cc} V \boldsymbol{I} & \sqrt{T} Z \boldsymbol{\sigma}_z \\ \sqrt{T} Z \boldsymbol{\sigma}_z & T\left(V+\chi_{\text {line }}\right) \boldsymbol{I} \end{array}\right] \\ \boldsymbol{I}=\left[\begin{array}{ll} 1 & 0 \\ 0 & 1 \end{array}\right] \\ \boldsymbol{\sigma}_z=\left[\begin{array}{cc} 1 & 0 \\ 0 & -1 \end{array}\right] \end{array}\right. $

      (6)

      $ \mu_a=\left\{\begin{array}{l} \sqrt{\frac{1}{2}\left(\varDelta \pm \sqrt{\varDelta^2-4 D}\right)}, (a=1, 2) \\ \sqrt{\frac{1}{2}\left(A \pm \sqrt{A^2-4 B}\right)}, (a=3, 4) \\ 1, (a=5) \end{array}\right. $

      (7)

      $ \left\{\begin{array}{l} \varDelta=V^2+T^2\left(V+\chi_{\text {line }}\right)^2-2 T Z^2 \\ D=T^2\left(V^2+V X_{\text {line }}-Z^2\right)^2 \end{array}\right. $

      (8)

      在接收端Bob执行外差探测,即:

      $ \begin{gathered} A= \\ \frac{\chi_{\text {hom }}{ }^2 \varDelta+D+1+2 \chi_{\text {hom }}\left[V \sqrt{D}+T\left(V+\chi_{\text {line }}\right)+2 T Z\right]}{T^2\left(V+\chi_{\text {tot }}\right)^2} \end{gathered} $

      (9)

      $ B=\left[\frac{V+\sqrt{D} \chi_{\text {het }}}{T\left(V+\chi_{\text {tot }}\right)}\right]^2 $

      (10)

      式中: χline=1/T-1+ξ为通道输入有关的量子信道过量噪声; χhet=[1+(1-η)+2υEN]/η为Bob外差检测输入相关的检测附加噪声; χhom=[(1-η)+υEN]/η为Bob零差检测输入相关的检测附加噪声; 信道输入端的等效总噪声则可以表示为χtot=χline+χhet/T; υEN为实际探测器电子器件产生的热噪声; V=VA+1, 且VA是发送方Alice的调制方差; 此外Z为Alice和Bob的相关性常数。

      Δm是与隐私放大安全相关的参数(由于有限码长效应)[14],由于Alice和Bob之间只共享有限数量的W个量子态,因此在安全性分析中,必须考虑有限码长效应导致的影响。有限码长效应降低了Alice和Bob对信道传输率和过量噪声的估计精度[15]。假设原始密钥以比特编码,则有限码长效应可以通过参数Δm表示,即:

      $ \Delta m=7 \sqrt{\frac{\operatorname{lb}\left(\frac{2}{\varepsilon}\right)}{m}}+\frac{2}{m} \mathrm{lb}\left(\frac{1}{\varepsilon_{\mathrm{PA}}}\right) $

      (11)

      式中: ε是平滑参数; εPA是隐私放大过程失败的概率。

      假定量子信道为线性信道,参数估计步骤需要通过分别考虑Alice和Bob的数据xy之间的线性模型来估计信道传输率和过量噪声[16]。在使用外差探测时,Alice和Bob的数据遵循标准线性模型y=tx+z,其中参数$t=\sqrt{T \eta / 2}$,z服从方差为σ2=ηTξ/2+1+ξth的正态分布,探测系统的热噪声ξth可在校准阶段测量。因为估计的好坏取决于参数估计步骤中考虑的状态数W-m,因此必须在有限尺寸效应中考虑,计算概率至少为1-εPE的密钥率的下限。密钥速率的下界由t下界tminσ2的上界σmax2计算,但概率为εPE/2的情况除外,由下式给出:

      $ t_{\min } \approx \sqrt{\frac{T \eta}{2}}-z_{\varepsilon_{\mathrm{PE}} / 2} \sqrt{\frac{\frac{\eta T}{2} \xi+1+\xi_{\mathrm{th}}}{2(W-m)(2\langle n\rangle)}} $

      (12)

      $ \begin{gathered} \sigma_{\max }^2 \approx \frac{\eta T}{2} \xi+1+\xi_{\mathrm{th}}+ \\ z_{\varepsilon_{\mathrm{PE}} / 2} \frac{\left(\frac{\eta T}{2} \xi+1+\xi_{\mathrm{th}}\right) \sqrt{2}}{\sqrt{2(W-m)}} \end{gathered} $

      (13)

      式中: zεPE/2是标准正态分布的100(1-εPE/2)百分位数,置信区间为(1-εPE)×100%,并且使得$1-\operatorname{erf}\left(z_{\varepsilon_{\mathrm{PE}} / 2} /\sqrt{2}\right)=\varepsilon_{\mathrm{PE}} / 2$,其中$\operatorname{erf}(x)=\frac{2}{{\rm{ \mathsf{ π} }}} \int_0^x \exp \left(-t^2\right) \mathrm{d} t$是误差函数[17]; 2(Wm)中的2表示的是对于外差检测,两个正交分量都被测量,并在参数估计步骤中独立处理。最后,可以计算密钥速率的下限,除了概率为εPE/2之外,考虑信道传输率的最小值TminTmin=2tmin2/η,以及过量噪声的最大值$\xi_{\max }=\left[2 /\left(\eta T_{\min }\right)\right]\left(\sigma_{\max }^2-1-\xi_{\mathrm{th}}\right)$,下文中所有考虑有限码长效应的密钥速率都是通过Tminξmax计算的。

    • GMM是一种参数概率密度函数,表示为高斯分量密度的加权和,使用期望最大化(expectation-maximum,EM)算法估计高斯参数[18]。采用GMM的基本假设是所有样本都来自l个高斯分布的混合,每一个调制符号都是一个单独的高斯分布。使用EM算法,可以基于高斯潜在分量对样本进行聚类,该算法能够确定每个量子态的类别。

      高斯混合模型,即:

      $ p(x)=\sum\limits_{j=1}^l w_j \times g\left(x \mid \mu_j, \varSigma_j\right) $

      (14)

      式中: wj>0是第j个高斯分量的权重,且$\sum\limits_{j=1}^l w_j=1$; $g\left(x \mid \mu_j, \varSigma_j\right)$是均值为μ和方差为Σ的高斯分量密度函数,即:

      $ \begin{gathered} g\left(x \mid \mu_j, \varSigma_j\right)=\frac{1}{\sqrt{(2 {\rm{ \mathsf{ π} }})^n\left|\varSigma_j\right|}} \times \\ \exp \left[-\frac{1}{2}\left(x-\mu_j\right)^{\mathrm{T}} \varSigma_j^{-1}\left(x-\mu_j\right)\right] \end{gathered} $

      (15)

      该算法的目的是找到参数$\lambda=\left\{w_1, \cdots, w_{l-1}; \mu_1, \cdots, \mu_l; \varSigma_1, \cdots, \varSigma_l\right\}$,使其对数似然最大化。

      EM算法是一种迭代优化策略,是概率模型依赖隐含变量寻找参数最大似然估计或者最大后验概率的算法[19]。它的计算方法中每一次迭代都分两步: 第1步为期望步,利用隐含变量计算其最大似然估计值; 第2步为极大步,通过一系列的参数使得似然不断地增加,即此次的似然高于上一次迭代产生的似然,通过上一步求得的最大似然值来计算参数,此过程反复迭代,直到收敛完成训练。

      (a) E-step计算每个样本对高斯分量的隶属度的后验概率:

      $ t_l\left(x_i\right)=\frac{w_j g\left(x_i \mid \mu_j, \varSigma_j\right)}{\sum\limits_{j=1}^l w_j \times g\left(x_i \mid \mu_j, \varSigma_j\right)} $

      (16)

      式中: xi是第i个数据点。

      (b) M-step应用最大似然估计来更新均值,协方差矩阵和混合比例,使用步骤(a)中获得的后验概率,即:

      $ L\left(\lambda \mid x_1, \cdots, x_m\right)=\prod\limits_{i=1}^m \sum\limits_{j=1}^l w_j \times g\left(x_i \mid \mu_j, \varSigma_j\right) $

      (17)

      式中:(λ|x1, …, xm)是某一个数据点。

    • 在多环M-APSK编码的CVQKD协议中,Alice制备量子态,通过不完美信道传输,因信道传输过程中存在损耗及其电子探测器噪声的影响,导致传输态会有相位漂移和能量衰减,以一定在概率分布在相空间,如图 2所示。16-APSK调制相干态经过100 km的量子信道传输后在相空间的表示,调制方差为0.35,接收端Bob执行外差探测[20],传输态与初始的调制量子态的不再相同,为提高系统的检测效率,本文作者提出在系统接收端添加基于最小欧氏距离的高斯混合模型算法的量子态分类器。

      图  2  16-APSK调制相干态经过量子信道传输在相空间的分布

      Figure 2.  Distribution of 16-APSK modulated coherent states in phase space through quantum channel transmission

      为了构造一个性能良好的分类器,将本文中提出的方案分为两个部分: 状态学习和状态预测,如图 3所示。状态学习阶段主要用于训练和评估分类器,由发送方Alice制备已知类别量子态,并通过不可信的量子传输信道发送给Bob,同时通过辅助信道将已知类别量子态的相关信息发送到接收端,对接收到的量子态Bob执行外差探测和特征提取,为提高接收端Bob对量子态的检测效率,为每一个相干态构造一组距离特征:

      $ \delta_{\min }=\min\limits_{x, x^{\prime} \in X, x \neq x^{\prime}}\left|x-x^{\prime}\right|^2 $

      (18)

      图  3  基于高斯混合模型的离散调制CVQKD协议

      Figure 3.  Discrete modulation CVQKD protocol based on GMM

      根据式(1),M-APSK调制相干态x, x′∈X可由极坐标幅度(半径)与相位的形式(Rk, 2πik/nk+θk)转为直角坐标的形式(xk, yk),x′为参考态,X为调制相干态的信号态集,ik为第k个同心圈圆周上的第i(i=0, 1, …, nk-1)个量子态,即:

      $ \left\{\begin{array}{l} x_{i_k}=R_k \cos \left(2 {\rm{ \mathsf{ π} }} \frac{i_k}{n_k}+\theta_k\right) \\ y_{i_k}=R_k \sin \left(2 {\rm{ \mathsf{ π} }} \frac{i_k}{n_k}+\theta_k\right) \end{array}\right. $

      (19)

      由式(18)可得:

      $ \delta_{\min }{ }^2=\min\limits_{x, x^{\prime} \in X, x \neq x^{\prime}}\left[\left(x_{i_k}-x_{i_k}{ }^{\prime}\right)^2+\left(y_{i_k}-y_{i_k}{ }^{\prime}\right)^2\right] $

      (20)

      16-APSK调制相干态的特征提取如图 4所示。橙色表示未知量子态,蓝色表示虚拟的参考态,即理想条件下的调制状态初始值。Bob接收传输的相干态并提取一组距离特征后,即通过判断可得到传输态与参考态之间的最小欧氏距离δmin,因此可识别出该信号态所属类别1。

      图  4  16-APSK调制相干态的特征提取

      Figure 4.  Feature extraction of 16-APSK modulated coherent states

      在提取出鲁棒性特征后用作状态学习分类器的输入数据,$Y=\left\{y_1, y_2, \cdots, y_l\right\}$是包含l个类别的标签空间,训练集为$U=\left\{\left(x_i, y_j\right) \mid 1 \leqslant i \leqslant m\right\}$,其中xiXd维属性向量(xi1, xi2, …, xid)T, yjYxi所属的一个类别(1≤jl)。状态学习阶段的目的是找到一个分类器函数h(·),给定一个阈值函数t: XR, 使得$h(x)=\{y \mid f(x, y)>t(x), y \in Y\}$。

      若$|x\rangle$表示未分类相干态,Cj计算是$|x\rangle$属于第jyj(1≤jl):

      $ C_j=\sum\limits_{\left(\left|x^*\right\rangle, Y^*\right) \in N(|x\rangle)}\left(y_j \in Y^*\right) $

      (21)

      式中: $\left(\left|x^*\right\rangle, Y^*\right)$表示训练集中已知类别的相干态。假设Hj表示相干态$|x\rangle$的类别为yj,则P(Hj|Cj)表示后验概率在Cj属于类别yj的条件下Hj为真,P(Hj|Cj)为后验概率在Cj属于类别yj的条件下Hj为假,让f($|x\rangle$, yj)=P(Hj|Cj)/P(Hj|Cj),则基于GMM的幅度相位联合编码CVQKD协议分类器为:

      $ \begin{aligned} h(|x\rangle)=y_j & \left\lvert\, \frac{P\left(H_j \mid C_j\right)}{P\left(\bar{H}_j \mid C_j\right)}>t(|x\rangle)\right., \\ & (1 \leqslant j \leqslant l) \end{aligned} $

      (22)

      P(Hj|Cj)后验概率大于t($|x\rangle$)·P(Hj|Cj)时, 未知类别的相干态$|x\rangle$属于yj类。

      根据贝叶斯定理可得到:

      $ f\left(|x\rangle, y_j\right)=\frac{P\left(H_j \mid C_j\right)}{P\left(\bar{H}_j \mid C_j\right)}=\frac{P\left(H_j\right) \cdot P\left(C_j \mid H_j\right)}{P\left(\bar{H}_j\right) \cdot P\left(C_j \mid \bar{H}_j\right)} $

      (23)

      式中: P(Hj)和P(Hj)分别表示Hj为真或假的先验概率; P(CjHj)和P(Cjj)分别代表事件Hj为真或假的条件下Cj属于类别yj的条件概率。

      图 5中表示了不同M-APSK调制格式下信噪比(signal-to-noise ratio, SNR)对系统分类精度的影响。分类器的效率在信噪比大于5 dB时增加,这种优越的性能可以归因于分类器的内部结构的鲁棒性和所提出的特征的优点; 在SNR低于10 dB时,调制状态的量子态数目越少精确度越高; 在SNR为9 dB时,4-APSK、16-APSK和32-APSK调制状态下系统的性能达到最佳; 在SNR为11 dB时,64-APSK调制状态下系统的性能达到最佳。

      图  5  不同M-APSK调制格式下识别精度与信噪比的变化曲线

      Figure 5.  Variation cur ves of recognition accuracy and SNR under different M-APSK modulation formats

      状态预测阶段,利用已构成的量子态识别分类器对未知量子态进行识别分类,经过多轮预测,Bob向Alice公开部分密钥信息对传输信道进行参数估计,反向协调和保密增强生成最终密钥。

    • 基于本文作者提出的高斯混合模型CVQKD协议,多环M-APSK连续变量量子密钥速率可由下式计算:

      $ K_{\text {GMM }}=\frac{m}{W}\left[\varLambda(T) \beta I_{\mathrm{BA}}-\chi_{\mathrm{E}, \text { GMM }}-\Delta m\right] $

      (24)

      式中: Λ(T)为分类器的效率,其取值将随着传输距离的变化而变化; χE, GMM表示在集体攻击下窃听者Eve通过与量子态相互作用获得的有用信息的Holevo量:

      $ \chi_{\mathrm{E}, \mathrm{GMM}}=S\left(\rho_{\mathrm{E}}\right)-\sum\limits_{y_j}^m P\left(y_j\right) S\left(\rho_{\mathrm{E}, y_j}\right) $

      (25)

      $ \rho_{\mathrm{E}}=\sum\limits_{y_j}^m P\left(y_j\right) \rho_{\mathrm{E}, y_j} $

      (26)

      式中: S(ρE)=-Tr(ρEln ρE)为冯诺依曼熵[21]; P(yj)为Bob测量得到的原始密钥yj的概率; ρE, yj是Bob测量原始密钥yj时窃听者Eve拦截的状态。

      在传统的离散调制CVQKD协议中,发送方包含m个有限且完整的编码事件A=Ai(i=1, 2, …, m),Alice随机选择发送给Bob,因离散调制相干态与其二进制的表示是固定的,Bob测量的原始密钥应服从均匀分布,在四态协议与八态协议中即P(yj)=1/4、P(yj)=1/8,因此窃听者在截获量子态后可恢复其对应编码;在本文作者提出的方案中,对于编码的对应关系有且仅有Alice知道,通过状态学习阶段Bob学习相对应的编码规则,因此对于窃听者Eve在截获到量子态时,即P(yj)=1/m→0(m→∞),很难获取有用的密钥信息。

    • 图 6中呈现了作为传输距离的函数密钥速率的变化曲线,考虑了有限码长效应。参数设置为ε=εPA=εPE=10-10; 协调效率β=0.96;过量噪声ξ=0.005SNU,其中, SNU(shot noise unit)为散粒噪声单元; 检测效率η=0.6;调制方差VA=0.35,且M-APSK星座调制态的连续环从内到外依次包含4、12、16、32、64、128和256个量子态。从图中可以观察到,有限码长效应明显限制了离散调制CVQKD系统在传输距离方面的性能。此外,超过4个环64-APSK的星座大小仅提供在密钥速率方面的低增益;通过增加环的数量可达到更远的传输距离,例如从1个环4-APSK增加到2个环16-APSK(内环4个点,外环12个点)增加了9.78 km,从2个环16-APSK增加到3环32-APSK(内环4个点,中环12个点,外环16个点),可以观察到增加了13.4 km。因此,通过增加环的数量,可以获得离散调制CVQKD系统的更高性能和效率。

      图  6  考虑有限码长效应M-APSK调制下CVQKD协议的密钥速率

      Figure 6.  Key rate of CVQKD protocol with M-APSK modulation considering finite code length effect

      图 7中显示了不同M-APSK调制格式下过量噪声对系统量子态识别精度的影响。从图中可以观察到,随着过量噪声的增加,系统量子态分类器的识别精确度呈下降趋势,特别是在过量噪声从0~0.1区间识别精确度急剧下降;过量噪声ξ=0.03SNU时系统的精确度低于0.6,说明本文作者提出的算法在低噪声场景下具备较好的性能。

      图  7  不同M-APSK调制格式下系统性能与过量噪声的变化曲线

      Figure 7.  Variation curves of system performance and excess noise under different M-APSK modulation formats

      图 8中显示出了不同码长效应对系统分类精度的影响。从图中可以观察到,使用较长的编码提取分类特征,提高了正确分类率。例如,当使用码长L=5000个符号时,与码长L=1000个符号相比,在大于5 dB的SNR下,分类准确度的改善几乎超过90%。结果表明,该方法是可行的, 并且分类准确率是信道的SNR和信号中的码长的函数,这种关系是由于基于这些参数所提取的特征值的波动,这可能表明误分类随着较低的SNR和较小的块长度而增加。

      图  8  不同码长对系统性能的变化曲线

      Figure 8.  Variation curve of system performance for different code lengths

      图 9是在考虑有限码长效应的情况下,在反向协调和集体攻击下基于高斯混合模型分类的离散调制CVQKD协议和传统外差探测下的离散调制CVQKD协议的安全密钥率与传输距离的仿真图。设置协调效率β=0.96,过量噪声ξ=0.005SNU,检测效率η=0.6。仿真结果表明,基于GMM算法的M-APSK调制协议的密钥速率得到提升,当密钥速率为10-6 bit/symbol时,针对基于GMM算法的32-APSK调制CVQKD协议的传输距离达到48 km,128-APSK调制CVQKD协议的传输距离接近60 km。

      图  9  基于GMM的M-APSK调制格式CVQKD协议密钥速率

      Figure 9.  GMM-based M-APSK modulation format CVQKD protocol key rate

    • 基于高斯混合模型的M-APSK离散调制CVQKD协议,可分为两个过程即状态学习和状态预测两个步骤,状态学习用于训练和测试分类器,而状态预测用于生成最终密钥。仿真结果表明,基于高斯混合模型的量子态识别分类器的识别精确度随着过量噪声的增加而下降;当过量噪声ξ=0.005SNU时,基于高斯混合模型的128-APSK离散调制CVQKD协议其传输距离和密钥率的性能随着调制方差的增加有所提升, 安全传输距离可接近60 km。

参考文献 (21)

目录

    /

    返回文章
    返回