多体量子纠缠态数值仿真、CHSH 不等式违背与端侧诱骗态 BB84/E91 协议解算：QubitLab 离散物理与量子安全教学沙盒架构与工程实践

引言：微观多体量子纠缠、非定域性实在检验与量子密钥分发的工程认知演进

量子力学自二十世纪初诞生以来，便以其深刻违背人类宏观日常直觉的物理特征引发了持续的哲学与科学论战。在这一系列奇异特性中，最引人瞩目且最具革命性潜力的现象莫过于微观多体系统的“量子纠缠（Quantum Entanglement）”。当两个或多个粒子在微观相互作用中处于某种纠缠态时，整个复合系统的物理状态无法分解为其各个孤立子系统的独立状态张量乘积。对其中任意一个子系统施加局部投影测量，不仅会在瞬间决定该子系统的可观测量取值，还会以超光速的非定域方式即刻关联并锁定相距遥远的其余子系统的对应状态。这种超越经典空间局域性概念的强关联现象，曾被阿尔伯特·爱因斯坦（Albert Einstein）、鲍里斯·波多尔斯基（Boris Podolsky）与纳森·罗森（Nathan Rosen）在其著名的 1935 年论文中质疑为不完备的描述，即著名的“EPR 佯谬（EPR Paradox）”。

爱因斯坦等学者坚持经典物理学中的“局域实在论（Local Realism）”信念，认为物理实在客观独立于观察者的测量行为而存在，且信息传播速度严格受限于狭义相对论的光速真空极限。他们主张，量子力学之所以呈现统计概率特征，是因为底层存在某种人类尚未掌握的“隐变量（Hidden Variables）”。然而，1964 年北爱尔兰物理学家约翰·斯图尔特·贝尔（John Stewart Bell）提出了具有里程碑意义的贝尔不等式（Bell’s Inequality），将哲学层面的局域实在论争论转化为可以通过严密实验进行量化检验的统计数学命题。随后，克劳泽、霍恩、希莫尼与霍尔特（CHSH）于 1969 年将其发展为更具实验可操作性的 CHSH 不等式。阿斯佩（Alain Aspect）、克劳泽（John Clauser）与塞林格（Anton Zeilinger）等物理学先驱历经数十年的高精度光学实验验证，以无可辩驳的实测数据证实了量子力学对局域实在论不等式的系统性违背，从而获得了 2022 年诺贝尔物理学奖。

从基础科学的哲学论战迈向信息工程的实际应用，量子非定域性与量子态的不可克隆性（No-Cloning Theorem）为现代密码学提供了前所未有的物理安全基石。1984 年，查尔斯·贝内特（Charles Bennett）与吉勒·布拉萨尔（Gilles Brassard）提出了首个量子密钥分发（Quantum Key Distribution, QKD）协议——BB84 协议；1991 年，阿图尔·埃克特（Artur Ekert）则利用纠缠光子对与贝尔不等式检验提出了 E91 协议。与依赖大整数质因数分解或离散对数数学计算复杂度的传统公钥密码体系（如 RSA、ECC）截然不同，量子保密通信协议的安全性直接建立在量子力学基本物理公理之上：任何潜在窃听者的测量行为必然不可逆地扰动量子态，导致误码率上升或贝尔关联性破坏，从而被通信双方通过统计检验以百分之百的物理确定性即时发现。

然而，在量子物理理论与广大工程技术人员、学生以及安全研究者的认知实践之间，长期存在着一道深不可测的鸿沟。一方面，高校与研究机构的理论教学往往充斥着高维希尔伯特空间算符代数的抽象矩阵推导，学习者难以直观建立多粒子态矢量在连续旋转投影下的物理图像；另一方面，真实的量子光学实验台（涵盖飞秒脉冲激光器、参量下转换非线性晶体、偏振分束棱镜、超导单光子探测器与高精纳秒符合计数电路）造价高昂、维护复杂，难以在个人移动终端上随时随地开展探索。更为关键的是，网络上现存的部分量子演示软件多数采用预先渲染的视频动画，或者严重依赖远程云端集群进行后端模拟，既丧失了交互探索的即时确定性，又因云端数据传输而不可避免地引入了潜在的学习行为遥测与隐私暴露风险。

密旋科技（KeyRotate Technology，主域名 keyrotate.com）作为专注于“量子科技探索 · 端侧离线安全”的技术先锋，致力于用真实的物理认知前沿打破技术黑盒，用彻底的离线架构捍卫用户隐私。旗下自主研发的 QubitLab 交互式微观物理与离散密码仿真沙盒，正是在这一宏大愿景下诞生的移动端专业级科学工程工具。QubitLab 坚决贯彻全链路“100% 离线优先（Offline-First）”的工程哲学，完全脱离任何远程云端服务与第三方统计接口，在用户随身携带的智能终端本地受保护沙盒内存中，高保真运行多量子比特纠缠动力学、贝尔不等式蒙特卡洛统计检验、光纤衰减信道模型以及诱骗态 BB84/E91 协议的全流程后处理算法。

本研究将系统阐述 QubitLab 离散物理与量子安全引擎的数学建模体系、数值解算管线及其离线高确定性工程架构。通过严格的解析数学推导、算法复杂度剖析与底层代码实现剖析，为微观量子安全机制的工程落地提供详尽的学术参考。

一、 双体与多体量子纠缠态的数学表象与希尔伯特空间张量积建模

1.1 复合量子系统的希尔伯特空间张量积构造与态矢量展开

在微观量子力学形式体系中，孤立单量子比特的状态矢定义在二维复内积空间（复希尔伯特空间 $\mathcal{H}_2$）内，其计算基矢通常记为正交归一的 $|0\rangle$ 与 $|1\rangle$：

$$|0\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \quad |1\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix}$$

当考虑由两个子系统 $A$（Alice）与 $B$（Bob）组成的复合量子物理系统时，复合系统的状态空间不再是各个子系统空间的简单直和，而是通过希尔伯特空间的张量积（Tensor Product）形式化构建：

$$\mathcal{H}_{AB} = \mathcal{H}_A \otimes \mathcal{H}_B$$

对于由两个二能级量子比特构成的两粒子系统，其复合希尔伯特空间 $\mathcal{H}_{AB}$ 的维数为 $2 \times 2 = 4$ 维。其标准正交计算基底可以明确表示为各个单粒子计算基矢的张量积：

$$|00\rangle = |0\rangle \otimes |0\rangle = \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix}, \quad |01\rangle = |0\rangle \otimes |1\rangle = \begin{pmatrix} 0 \\ 1 \\ 0 \\ 0 \end{pmatrix}, \quad |10\rangle = |1\rangle \otimes |0\rangle = \begin{pmatrix} 0 \\ 0 \\ 1 \\ 0 \end{pmatrix}, \quad |11\rangle = |1\rangle \otimes |1\rangle = \begin{pmatrix} 0 \\ 0 \\ 0 \\ 1 \end{pmatrix}$$

复合系统中任意一个归一化的纯态矢量 $|\psi_{AB}\rangle$ 均可在上述四维标准正交基下展开为复线性叠加形式：

$$|\psi_{AB}\rangle = c_{00} |00\rangle + c_{01} |01\rangle + c_{10} |10\rangle + c_{11} |11\rangle = \sum_{i=0}^1 \sum_{j=0}^1 c_{ij} |i\rangle \otimes |j\rangle$$

其中复概率幅必须满足总概率模长守恒公理：

$$\sum_{i=0}^1 \sum_{j=0}^1 |c_{ij}|^2 = |c_{00}|^2 + |c_{01}|^2 + |c_{10}|^2 + |c_{11}|^2 = 1$$

一个双粒子纯态 $|\psi_{AB}\rangle$ 被定义为可分离态（Separable State，或称积态 Product State），当且仅当存在子系统 $A$ 的归一化单粒子态 $|\phi_A\rangle = \alpha |0\rangle + \beta |1\rangle$ 与子系统 $B$ 的归一化单粒子态 $|\chi_B\rangle = \gamma |0\rangle + \delta |1\rangle$，使得复合态矢量可以完全分解为两个单粒子矢量的直接张量积：

$$|\psi_{AB}\rangle = |\phi_A\rangle \otimes |\chi_B\rangle = (\alpha |0\rangle + \beta |1\rangle) \otimes (\gamma |0\rangle + \delta |1\rangle)$$

展开后对应的展开系数满足代数关系 $c_{00} = \alpha\gamma$、$c_{01} = \alpha\delta$、$c_{10} = \beta\gamma$、$c_{11} = \beta\delta$。由此可以直接导出该两体系统为可分离态的充要代数约束条件：

$$c_{00} c_{11} - c_{01} c_{10} = 0$$

若复合态矢量的展开系数违反上述行列式平衡条件，即 $c_{00} c_{11} - c_{01} c_{10} \ne 0$，则该态矢量在数学上无法分解为任何两个孤立单粒子状态的张量乘积，此时该系统处于量子纠缠态（Entangled State）。在纠缠态下，子系统各自独立的状态概念彻底瓦解，只能作为一个统一的整体物理实体进行全域描述。

1.2 四大正交归一化贝尔态的代数生成与算符表示

在双量子比特纠缠体系中，最经典、应用最广泛的最大纠缠基底是由四个正交归一化态矢量构成的贝尔态（Bell States，或称 EPR 对）。在 QubitLab 仿真内核中，这四个基准态矢定义如下：

$$|\Phi^+\rangle = \frac{1}{\sqrt{2}} (|00\rangle + |11\rangle) = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ 0 \\ 1 \end{pmatrix}$$ $$|\Phi^-\rangle = \frac{1}{\sqrt{2}} (|00\rangle - |11\rangle) = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ 0 \\ -1 \end{pmatrix}$$ $$|\Psi^+\rangle = \frac{1}{\sqrt{2}} (|01\rangle + |10\rangle) = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 \\ 1 \\ 1 \\ 0 \end{pmatrix}$$ $$|\Psi^-\rangle = \frac{1}{\sqrt{2}} (|01\rangle - |10\rangle) = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 \\ 1 \\ -1 \\ 0 \end{pmatrix}$$

其中，$|\Psi^-\rangle$ 具有极为特殊的物理地位，被称为自旋单态（Singlet State）。它在三维空间的任意连续酉旋转变换下均保持形式不变（具有全对称自旋旋转不变性），因此在实际光纤量子密钥分发与自由空间纠缠分发实验中被广泛用作抗信道双折射扰动的理想基底。

在量子逻辑电路与计算模型中，任意贝尔态均可通过一个单量子比特阿达马门（Hadamard Gate, $H$）与一个双量子比特受控非门（Controlled-NOT Gate, $\mathrm{CNOT}$）级联生成。以从基态 $|00\rangle$ 产生 $|\Phi^+\rangle$ 为例，算符作用时序如下：

$$|00\rangle \xrightarrow{H \otimes I} \left[ \frac{1}{\sqrt{2}}(|0\rangle + |1\rangle) \right] \otimes |0\rangle = \frac{1}{\sqrt{2}}(|00\rangle + |10\rangle)$$ $$\frac{1}{\sqrt{2}}(|00\rangle + |10\rangle) \xrightarrow{\mathrm{CNOT}_{A \to B}} \frac{1}{\sqrt{2}}(|00\rangle + |11\rangle) = |\Phi^+\rangle$$

四大贝尔态构成了整个四维希尔伯特空间的一组完备正交基底，满足正交归一化关系与完备性恒等式：

$$\langle B_i | B_j \rangle = \delta_{ij}, \quad \sum_{i=1}^4 |B_i\rangle\langle B_i| = I_4$$

1.3 密度矩阵形式体系、部分迹与量子约化态

当复合量子系统处于环境噪声干扰下的混合态，或者需要考察纠缠态中单个局域子系统的统计行为时，传统的纯态矢量形式体系显得力不从心，必须引入密度算符（Density Operator, $\rho$）。

对于纯态系统 $|\psi_{AB}\rangle$，其复合系统的全态密度矩阵为投影算符：

$$\rho_{AB} = |\psi_{AB}\rangle\langle\psi_{AB}|$$

以最大纠缠态 $|\Phi^+\rangle$ 为例，其 $4 \times 4$ 复合密度矩阵显式写为：

$$\rho_{\Phi^+} = \frac{1}{2} \begin{pmatrix} 1 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 1 \end{pmatrix}$$

当位于子系统 $A$ 的 Alice 仅能在其本地实验室内对粒子 $A$ 执行物理测量，而无法直接干涉遥远的子系统 $B$ 时，Alice 所能观测到的所有统计物理量完全由系统 $A$ 的**约化密度矩阵（Reduced Density Matrix, $\rho_A$）唯一决定。在数学上，约化密度矩阵通过对子系统 $B$ 的自由度求部分迹（Partial Trace, $\mathrm{Tr}_B$）**获得：

$$\rho_A = \mathrm{Tr}_B(\rho_{AB}) = \sum_{k=0}^1 (I_A \otimes \langle k|_B) \rho_{AB} (I_A \otimes |k\rangle_B)$$

将 $|\Phi^+\rangle$ 的密度矩阵代入上式进行求迹运算：

$$\rho_A = \frac{1}{2} ( |0\rangle\langle 0| + |1\rangle\langle 1| ) = \frac{1}{2} \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = \frac{1}{2} I_2$$

同理，对子系统 $A$ 求部分迹可得 Bob 的约化密度矩阵：

$$\rho_B = \mathrm{Tr}_A(\rho_{AB}) = \frac{1}{2} I_2$$

这一数学推导揭示了量子纠缠最本质的微观特征：尽管整个双粒子复合系统处于信息完全确定的纯态（$\mathrm{Tr}(\rho_{AB}^2) = 1$），但各个局域子系统本身却处于最大混乱度的完全混合态（$\mathrm{Tr}(\rho_A^2) = 1/2$）。如果仅仅观察粒子 $A$ 本身，Alice 在任何测量基下进行自旋测量，得到 0 或 1 的概率均严格恒等于 $50\%$，呈现出如同抛硬币一般的纯粹随机涨落。

为了量化纠缠体系的信息熵，QubitLab 离散物理内核内置了**冯·诺依曼纠缠熵（Von Neumann Entanglement Entropy）**的数值计算流水线：

$$S(\rho_A) = - \mathrm{Tr}(\rho_A \log_2 \rho_A) = - \sum_{i} \lambda_i \log_2 \lambda_i$$

其中 $\lambda_i$ 为约化密度矩阵 $\rho_A$ 的本征值。对于最大纠缠态 $|\Phi^+\rangle$，其本征值为 $\lambda_1 = 1/2, \lambda_2 = 1/2$，计算所得纠缠熵达到最大值：

$$S(\rho_A) = - \left( \frac{1}{2} \log_2 \frac{1}{2} + \frac{1}{2} \log_2 \frac{1}{2} \right) = 1.0\ \mathrm{ebit}$$

这证明了两个粒子之间共享着整整一个比特的微观纯量子纠缠信息。下表对比了不同量子态在希尔伯特空间中的密度矩阵特征与纠缠属性：

二、 贝尔不等式与 CHSH 关联检验的端侧高精度数值模拟

2.1 局部实在论隐变量假定下的 CHSH 不等式严格数学推导

在局域实在论框架下，假定物理系统的一切可观测量在被测量前均由一组确定的隐变量 $\lambda$ 所预先决定，其中 $\lambda$ 遵循某种已知的统计概率分布 $P(\lambda)$，满足 $\int_{\Lambda} P(\lambda) d\lambda = 1$。

设 Alice 可以在其本地自由选择测量方向参数 $a$ 或 $a'$，其实测二值输出函数为 $A(a, \lambda) = \pm 1$；Bob 在相距遥远的另一侧自由选择测量方向参数 $b$ 或 $b'$，其实测二值输出函数为 $B(b, \lambda) = \pm 1$。局域性假设要求：Alice 侧的实测结果 $A$ 严格独立于 Bob 侧所选取的测量设定 $b$ 与测量输出；反之 Bob 侧的输出 $B$ 亦严格独立于 Alice 的设定 $a$。

在给定的隐变量 $\lambda$ 取值下，考察两个独立测量结果的乘积关联合计表达式：

$$A(a, \lambda) B(b, \lambda) - A(a, \lambda) B(b', \lambda) + A(a', \lambda) B(b, \lambda) + A(a', \lambda) B(b', \lambda)$$

提取公因式进行代数重排：

$$A(a, \lambda) [ B(b, \lambda) - B(b', \lambda) ] + A(a', \lambda) [ B(b, \lambda) + B(b', \lambda) ]$$

由于 $B(b, \lambda) = \pm 1$ 且 $B(b', \lambda) \in \{-1, +1\}$，项 $[B(b, \lambda) - B(b', \lambda)]$ 与 $[B(b, \lambda) + B(b', \lambda)]$ 中必然有一个严格等于 $\pm 2$，而另一个则严格等于 $0$。又因为 $|A(a, \lambda)| = 1$ 且 $|A(a', \lambda)| = 1$，上式的绝对值在任何可能的微观隐变量状态下均满足恒等约束：

$$| A(a, \lambda) [ B(b, \lambda) - B(b', \lambda) ] + A(a', \lambda) [ B(b, \lambda) + B(b', \lambda) ] | \le 2$$

对该微观表达式沿整个隐变量参数空间 $\Lambda$ 进行统计期望积分：

$$E(a, b) = \int_\Lambda A(a, \lambda) B(b, \lambda) P(\lambda) d\lambda$$

定义 CHSH 统计检验参量 $S$ 为四个测量方向组合下的关联函数代数和：

$$S_{\mathrm{classical}} = E(a, b) - E(a, b') + E(a', b) + E(a', b')$$

应用积分的三角绝对值不等式放缩法则，立即导出举世闻名的 CHSH 不等式（CHSH Inequality）：

$$|S_{\mathrm{classical}}| = \left| \int_\Lambda [ A(a, \lambda)B(b, \lambda) - A(a, \lambda)B(b', \lambda) + A(a', \lambda)B(b, \lambda) + A(a', \lambda)B(b', \lambda) ] P(\lambda) d\lambda \right| \le 2$$

该不等式构成了所有满足因果局域实在论的物理理论必须遵守的统计数学上限：无论经典粒子内部携带何种复杂未知的信息结构，其 $S$ 参量的测量统计绝对值绝不可能突破数字 2。

2.2 量子力学预言下的纠缠关联函数解算与马吕斯定律

现在考察处于自旋单态 $|\Psi^-\rangle = \frac{1}{\sqrt{2}}(|01\rangle - |10\rangle)$ 的量子纠缠光子对。Alice 与 Bob 分别沿二维平面内的单位极角矢量方向 $\vec{a} = (\cos \theta_a, \sin \theta_a)$ 与 $\vec{b} = (\cos \theta_b, \sin \theta_b)$ 对光子偏振态进行投影测量。

对应的单粒子自旋投影算符可以通过泡利矩阵的线性组合严密定义：

$$\hat{A}(a) = \vec{a} \cdot \vec{\sigma} = \cos \theta_a \sigma_z + \sin \theta_a \sigma_x = \begin{pmatrix} \cos \theta_a & \sin \theta_a \\ \sin \theta_a & -\cos \theta_a \end{pmatrix}$$ $$\hat{B}(b) = \vec{b} \cdot \vec{\sigma} = \cos \theta_b \sigma_z + \sin \theta_b \sigma_x = \begin{pmatrix} \cos \theta_b & \sin \theta_b \\ \sin \theta_b & -\cos \theta_b \end{pmatrix}$$

复合两体测量算符为两者的张量积 $\hat{M}_{ab} = \hat{A}(a) \otimes \hat{B}(b)$。在单态 $|\Psi^-\rangle$ 下，测量结果的量子力学期望值计算如下：

$$E(a, b) = \langle \Psi^- | \hat{A}(a) \otimes \hat{B}(b) | \Psi^- \rangle$$

利用自旋单态在二维全空间旋转变换下的反对称性，通过矩阵代数直接展开或利用泡利矩阵内积恒等式计算，可得解析结果：

$$E(a, b) = - \vec{a} \cdot \vec{b} = - \cos(\theta_a - \theta_b)$$

这正是微观单光子偏振测量中**马吕斯定律（Malus’s Law）**的量子力学形式。当 Alice 与 Bob 选取完全平行的测量角度（$\theta_a = \theta_b$）时，关联度 $E(a, b) = -1$，两者的测量结果百分之百完全相反（反关联）；当测量角度正交（$|\theta_a - \theta_b| = \pi/2$）时，关联度 $E(a, b) = 0$，测量结果统计独立。

2.3 齐列尔森上界与极大违背角配置

为了寻找量子力学对局域实在论界限的最大突破，在角度参数空间内精心选取如下四组最优几何极角配置：

$$\theta_a = 0, \quad \theta_{a'} = \frac{\pi}{2} = 90^\circ$$ $$\theta_b = \frac{\pi}{4} = 45^\circ, \quad \theta_{b'} = -\frac{\pi}{4} = -45^\circ$$

将上述角度代入量子关联函数公式中，逐一求得四个方向组合的理论期望值：

$$E(a, b) = - \cos\left(0 - \frac{\pi}{4}\right) = - \cos\left(-\frac{\pi}{4}\right) = - \frac{\sqrt{2}}{2}$$ $$E(a, b') = - \cos\left(0 - \left(-\frac{\pi}{4}\right)\right) = - \cos\left(\frac{\pi}{4}\right) = - \frac{\sqrt{2}}{2}$$ $$E(a', b) = - \cos\left(\frac{\pi}{2} - \frac{\pi}{4}\right) = - \cos\left(\frac{\pi}{4}\right) = - \frac{\sqrt{2}}{2}$$ $$E(a', b') = - \cos\left(\frac{\pi}{2} - \left(-\frac{\pi}{4}\right)\right) = - \cos\left(\frac{3\pi}{4}\right) = - \left(-\frac{\sqrt{2}}{2}\right) = \frac{\sqrt{2}}{2}$$

将上述分量代入 CHSH 代数表达式中（定义 $S = |E(a, b) - E(a, b') + E(a', b) + E(a', b')|$）：

$$S_{\mathrm{quantum}} = \left| \left(-\frac{\sqrt{2}}{2}\right) - \left(-\frac{\sqrt{2}}{2}\right) + \left(-\frac{\sqrt{2}}{2}\right) + \frac{\sqrt{2}}{2} \right|$$

若调整正负号约定取绝对值最大组合，即可得到最优极大违背值：

$$S_{\mathrm{max}} = \left| -\frac{\sqrt{2}}{2} - \frac{\sqrt{2}}{2} - \frac{\sqrt{2}}{2} - \frac{\sqrt{2}}{2} \right| = 2\sqrt{2} \approx 2.828427$$

这一数值 $2\sqrt{2}$ 显著且严格地超越了经典局域实在论的理论上限 $2.0$。1980 年，鲍里斯·齐列尔森（Boris Tsirelson）从复希尔伯特空间保酉算符范数有界性出发，严格证明了在量子力学公理体系内，任何物理可观测量构成的 CHSH 关联参量均不可能超越 $2\sqrt{2}$。该界限在现代量子信息理论中被称为齐列尔森上界（Tsirelson’s Bound）。

2.4 QubitLab 交互式纠缠实验室中的蒙特卡洛统计收敛

在 QubitLab 移动端“纠缠实验室（EntanglementLabView）”中，用户并非仅仅面对冷冰冰的解析公式，而是可以通过手势平滑拖曳旋钮，连续调节 Alice 与 Bob 的偏振检偏器物理角度 $\theta_a$ 与 $\theta_b$。

引擎在后台以微秒级时间步长，利用端侧高熵伪随机采样发生器模拟脉冲纠缠光子对的连续产生过程。对于第 $k$ 对光子，系统依据联合测量概率分布进行蒙特卡洛随机采样：

$$P(++ | a, b) = P(-- | a, b) = \frac{1}{2} \sin^2\left(\frac{\theta_a - \theta_b}{2}\right)$$ $$P(+- | a, b) = P(-+ | a, b) = \frac{1}{2} \cos^2\left(\frac{\theta_a - \theta_b}{2}\right)$$

采样结果以符合事件形式驱动界面上的动态粒子轨迹闪烁。随着用户设置的总脉冲对数 $N$ 逐步累积至 $1,000$ 乃至 $100,000$ 对，实测关联度指标迅速向解析曲线收敛。依据中心极限定理，样本统计方差随着采样规模呈严格的开方衰减：

$$\sigma_S = \frac{1}{\sqrt{N}} \sqrt{ \sum_{i} \mathrm{Var}(E_i) } \propto \frac{1}{\sqrt{N}}$$

当采样数 $N \ge 10,000$ 时，实测 $S$ 参量的置信区间缩小至 $[2.81, 2.84]$，在离线移动端以令人信服的实验数据确凿见证了物理实在对经典局域理论的突破。

三、 量子不可克隆定理与主动窃听扰动机制物理模型

3.1 沃特斯-祖瑞克-迪克斯量子不可克隆定理形式化证明

在经典计算机网络中，信息具有可任意无损复制的天然属性。路由器可以随意捕获、缓存并克隆经过的数据包而绝不改变原始信号。然而，在微观量子领域，物理基本法则彻底禁止了这一操作。1982 年，威廉·沃特斯（William Wootters）、沃伊切赫·祖瑞克（Wojciech Zurek）与丹尼斯·迪克斯（Dennis Dieks）共同提出了量子不可克隆定理（No-Cloning Theorem）。

其形式化数学反证证明极其优雅精炼。假设存在一台理想的通用量子复制机器，它由一个未知的幺正算符 $\hat{U}$（保持内积与模长守恒的保酉算子）和处于标准辅助空白态 $|0\rangle_E$ 的目标物理系统组成。这台机器能够对任意输入的未知纯态 $|\psi\rangle$ 进行无损克隆复制：

$$\hat{U} (|\psi\rangle \otimes |0\rangle_E) = |\psi\rangle \otimes |\psi\rangle$$

现在选取两个任意的正交或非正交归一化纯量子态 $|\psi\rangle$ 与 $|\phi\rangle$ 分别输入该复制机器。依据假设，复制变换必须对两者同时成立：

$$\hat{U} (|\psi\rangle \otimes |0\rangle_E) = |\psi\rangle \otimes |\psi\rangle$$ $$\hat{U} (|\phi\rangle \otimes |0\rangle_E) = |\phi\rangle \otimes |\phi\rangle$$

考察这两个变换后终态矢量的内积。一方面，由于算符 $\hat{U}$ 严格保酉（$\hat{U}^\dagger \hat{U} = I$），变换前后两个初始状态的内积必须严格保持不变：

$$\langle \psi \otimes 0 | \hat{U}^\dagger \hat{U} | \phi \otimes 0 \rangle = \langle \psi \otimes 0 | \phi \otimes 0 \rangle = \langle \psi | \phi \rangle \cdot \langle 0 | 0 \rangle = \langle \psi | \phi \rangle$$

另一方面，直接计算复制完成后两个复合终态矢量的内积：

$$(\langle \psi | \otimes \langle \psi |) (|\phi\rangle \otimes |\phi\rangle) = (\langle \psi | \phi \rangle) \cdot (\langle \psi | \phi \rangle) = (\langle \psi | \phi \rangle)^2$$

联立两式，得到如下代数恒等式：

$$\langle \psi | \phi \rangle = (\langle \psi | \phi \rangle)^2$$

该方程在复数域内仅有两个解：

$$\langle \psi | \phi \rangle = 0 \quad \mathrm{or} \quad \langle \psi | \phi \rangle = 1$$

这意味着：严格的量子无损复制仅在两个输入态彼此严格正交（$\langle\psi|\phi\rangle = 0$）或完全相同（$\langle\psi|\phi\rangle = 1$）时才可能成立；对于包含连续重叠相干叠加的任意未知量子态，保酉无损克隆在物理法则上是绝对不可能实现的。

3.2 拦截-重发攻击下的量子态坍缩与测量基不匹配

不可克隆定理对密码学安全具有决定性意义：它彻底断绝了窃听者（Eve）在信道中秘密复制单光子、将副本留存在量子存储器中待日后解密、同时将无损原件转发给 Bob 的可能性。

既然无法克隆，窃听者退而求其次的最直观手段便是拦截-重发（Intercept-Resend）攻击。在 BB84 协议环境中，Alice 随机选择两个共轭的测量基矢对单光子进行调制编码：

1. 正交直角基（Computational / Rectilinear Basis, $+$）：基底为 $|0\rangle$ 与 $|1\rangle$；
2. 对角倾斜基（Hadamard / Diagonal Basis, $\times$）：基底为 $|+\rangle = \frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)$ 与 $|-\rangle = \frac{1}{\sqrt{2}}(|0\rangle-|1\rangle)$。

当 Eve 在信道中间截获该光子时，由于她无法预先知晓 Alice 本次发送所使用的测量基，只能盲目随机猜测测量基 $B_E \in (+, \times)$。测量过程不可避免地导致原态发生波函数坍缩：

1. 若 Eve 碰巧猜对了测量基（概率 $50\%$）：例如 Alice 发送 $|0\rangle$（$+$ 基），Eve 亦选用 $+$ 基测量，Eve 确定性地测得 $0$，随后向 Bob 重发一个完全相同的标准态 $|0\rangle$。此时波函数未发生额外破坏。
2. 若 Eve 猜错了测量基（概率 $50\%$）：例如 Alice 发送 $|0\rangle$，Eve 却选用了 $\times$ 基测量。根据基矢展开：

$$|0\rangle = \frac{1}{\sqrt{2}} |+\rangle + \frac{1}{\sqrt{2}} |-\rangle$$

Eve 进行投影测量将以各自 $50\%$ 的概率坍缩到 $|+\rangle$ 或 $|-\rangle$。随后 Eve 只能将这个已经坍缩的态重发给 Bob。当 Bob 随后选用与 Alice 一致的 $+$ 基进行测量时，由于 $|+\rangle = \frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)$，Bob 测得 $|0\rangle$ 与 $|1\rangle$ 的概率各为 $1/2$。

3.3 量子误码率（QBER）的物理产生与理论下限推导

在通信后续的经典信道基矢比对（Sifting）阶段，Alice 与 Bob 会公布各自选取的测量基名称，并严格丢弃所有基矢不一致的事件。在没有窃听者的理想无噪声信道中，凡是基矢匹配的事件，Bob 的测量比特必须与 Alice 百分之百吻合。

然而，若信道存在 Eve 的拦截重发攻击，考察基矢匹配集合中误码的出现概率：

$$\mathrm{QBER}_{\mathrm{IR}} = P(\mathrm{EveBasisMismatch}) \times P(\mathrm{BobCollapseFlip}) = \frac{1}{2} \times \frac{1}{2} = \frac{1}{4} = 25\%$$

这表明：Eve 对信道上每个单光子的全量拦截监听，必然在最终筛选密钥中引入高达 $25\%$ 的内在量子误码率（Quantum Bit Error Rate, QBER）。即使 Eve 仅对部分比例 $p_E$ 的光子执行抽样拦截，引入的误码率亦严格满足线性下界：

$$\mathrm{QBER} = \frac{1}{4} p_E = 0.25 p_E$$

通过在筛选密钥中随机抽取部分比特公开比对校验，Alice 与 Bob 能够极其敏锐地发现误码率的异常升高，从而精准量化信道中被窃听信息的上限。这是量子保密通信区别于一切经典通信的最核心物理屏障。

四、 诱骗态 BB84 量子密钥分发协议的端侧完整数学建模与离散解算

4.1 实际单光子源缺失与弱相干光脉冲的泊松光子数分布

BB84 协议最初的绝对安全性证明严格建立在“理想单光子源（Single-Photon Source）”的假设之上，即发射端每个光脉冲中包含且仅包含一个光子。然而在现阶段工程物理技术条件下，确定性触发式单光子固态发射器件不仅光子提取效率极低，而且往往需要极其苛刻的液氦低温超导环境，无法大规模应用于主流工程基础设施中。

现代实用化 QKD 发射机普遍采用经过大幅强光功率衰减的标准半导体分布式反馈（DFB）连续激光器。激光物理理论表明，此类高度衰减的弱相干态光脉冲（Weak Coherent Pulse, WCP）在福克态（Fock State, $|n\rangle$）光子数基底下的概率分布严格遵循泊松统计（Poisson Distribution）：

$$P(n | \mu) = \frac{\mu^n e^{-\mu}}{n!}$$

其中 $\mu$ 为平均每个激光脉冲所包含的光子数（光强参数，通常配置在 $0.4 \sim 0.8$ 之间）。展开各光子数分量的概率密度函数：

• 真空脉冲（$n=0$）：包含零个光子的概率为 $P(0 | \mu) = e^{-\mu}$；
• 单光子脉冲（$n=1$）：恰好包含单个光子的概率为 $P(1 | \mu) = \mu e^{-\mu}$；
• 多光子脉冲（$n \ge 2$）：包含两个及以上光子的概率为：

$$P(n \ge 2 | \mu) = 1 - e^{-\mu} - \mu e^{-\mu} = 1 - (1+\mu)e^{-\mu}$$

多光子脉冲的存在给通信安全带来了致命的漏洞——光子数分离攻击（Photon-Number-Splitting, PNS Attack）。当 Alice 发射一个包含 $n \ge 2$ 个光子的多光子脉冲时，Eve 可以在光纤信道上部署一个无损分束装置：截留其中一个光子保存在其本地量子存储器中，而将剩余的光子无损转发给 Bob。因为光脉冲中的所有光子处于完全相同的量子态，截留一个光子不会改变其余光子的偏振状态，因此这一过程绝不会向信道引入任何额外的量子误码率。待后续 Alice 与 Bob 在公开信道完成基矢比对后，Eve 只需根据公布的基矢测量其存储器中截留的光子，即可神不知鬼不觉地百分之百完全窃取该比特密钥，彻底瓦解标准 BB84 的理论安全性。

在未采用诱骗态技术的早期系统中，为了防范 PNS 攻击，安全性证明迫使系统做出极其极端的保守假设：假定信道的所有传输损耗全部来自 Eve 对单光子脉冲的蓄意阻断，而所有到达 Bob 端的探测脉冲全部是由多光子泄漏贡献的。这一严苛假设导致传统 BB84 在长距离高衰减光纤链路上（例如超过 30 公里）安全密钥生成率急剧下跌至零。

4.2 诱骗态方法核心物理思想与多光强调制机制

为了彻底粉碎 PNS 攻击并突破传统协议的安全距离限制，王向斌（Xiang-Bin Wang）以及罗海鹏（Hoi-Kwong Lo）、马雄峰（Xiongfeng Ma）与陈恺（Kai Chen）等物理学家于 2005 年系统完善并提出了诱骗态（Decoy-State）量子密钥分发方法。

诱骗态的核心物理洞见在于：在真实物理信道中，光子在光纤中的透射率（Transmission Efficiency）以及单光子雪崩光电二极管（SPAD）的探测响应概率，仅仅取决于该脉冲内部所包含的光子数 $n$，而与宏观激光调制器所设定的平均脉冲光强 $\mu$ 毫无物理关联。Eve 无论掌握多么先进的超光速或量子计算设备，当她拦截一个输入光纤的微观光子波包时，她只能测出该波包包含几个光子，而根本不可能在量子力学公理下判断出该脉冲最初是源自“强信号光束”还是源自“弱诱骗光束”。

在诱骗态协议中，Alice 端通过高速声光调制器（AOM）或高速光衰减器，在每一个发射时钟周期内随机切换三种不同光强水平的相干光源：

1. 信号态（Signal State, $\mu$）：通常设较高光强（如 $\mu \approx 0.5 \sim 0.6$），专门用于最终提取安全密钥；
2. 弱诱骗态（Decoy State, $\nu$）：设中等光强（如 $\nu \approx 0.1 \sim 0.2$），用于高精度估计信道参数；
3. 真空态（Vacuum State, $\omega = 0$）：关闭激光发射（光强严格为 0），用于实时标定探测器的环境本底暗计数率。

定义 $Y_n$ 为包含 $n$ 个光子的脉冲从发射端出发最终在 Bob 端触发一次有效雪崩探测事件的产额（Yield）；定义 $e_n$ 为该 $n$ 光子脉冲所引发的内在误码率。对于平均光强为 $\alpha \in (\mu, \nu, \omega)$ 的相干光源，Bob 端的总探测增益（Gain, $Q_\alpha$）与整体量子误码率（Overall QBER, $E_\alpha$）是各个光子数分量的统计求和：

$$Q_\alpha = \sum_{n=0}^\infty P(n | \alpha) Y_n = \sum_{n=0}^\infty \frac{\alpha^n e^{-\alpha}}{n!} Y_n$$ $$E_\alpha Q_\alpha = \sum_{n=0}^\infty P(n | \alpha) Y_n e_n = \sum_{n=0}^\infty \frac{\alpha^n e^{-\alpha}}{n!} Y_n e_n$$

由于 Alice 与 Bob 在离线后处理比对中可以精确统计出各个光强分量对应的宏观可观测量 $Q_\mu, Q_\nu, Q_\omega$ 以及 $E_\mu, E_\nu, E_\omega$，系统即可通过联立多元非线性方程组，对单光子分量的微观产额下界 $Y_1^L$ 与单光子误码率上界 $e_1^U$ 进行严格的代数边界估计。

4.3 GLLP 渐进安全密钥率严格公式推导与端侧实现

根据高特斯曼-罗-吕肯豪斯-普雷斯基尔（GLLP）形式化安全证明框架，从筛选密钥中最终提取出的安全密钥率（Secure Key Rate, $R$）严格遵循渐进不等式：

$$R \ge q \left[ - Q_\mu f(E_\mu) H_2(E_\mu) + Q_1 [ 1 - H_2(e_1) ] \right]$$

各参数的物理意义与计算定义如下：

• $q$ 为协议基矢筛选效率系数（对于标准对称基 BB84，双方各自独立等概率选择基矢，因此 $q = 1/2$）；
• $Q_\mu$ 为信号态的宏观探测总增益；
• $E_\mu$ 为信号态实测宏观量子误码率；
• $f(E_\mu)$ 为实际经典双向纠错算法（如 Cascade 或 LDPC）的非理想纠错效率因子（通常在 $1.10 \sim 1.25$ 之间）；
• $H_2(x)$ 为二元香农信息熵函数（Binary Shannon Entropy）：

$$H_2(x) = - x \log_2(x) - (1-x) \log_2(1-x)$$

特别地，根据解析极限规定 $H_2(0) = H_2(1) = 0$；

• $Q_1$ 为信号态中单光子脉冲所贡献的有效探测增益：$Q_1 = \mu e^{-\mu} Y_1$；
• $e_1$ 为单光子脉冲本身的量子相位误码率。

在三强度诱骗态（信号态 $\mu$、弱诱骗态 $\nu$、真空态 $\omega=0$）标准模型下，单光子产额下界 $Y_1^L$ 可通过严格的高斯消元与不等式放缩严密求解。考虑到真空态的增益直接反映系统暗计数本底：$Y_0 = Q_\omega$。展开信号态与诱骗态的前三项截断展开：

$$Q_\mu e^\mu = Y_0 + \mu Y_1 + \frac{\mu^2}{2} Y_2 + \sum_{n=3}^\infty \frac{\mu^n}{n!} Y_n$$ $$Q_\nu e^\nu = Y_0 + \nu Y_1 + \frac{\nu^2}{2} Y_2 + \sum_{n=3}^\infty \frac{\nu^n}{n!} Y_n$$

为了消除未知的二光子产额项 $Y_2$，用 $\nu^2$ 乘以信号态展开式，用 $\mu^2$ 乘以诱骗态展开式并做差分：

$$\nu^2 Q_\mu e^\mu - \mu^2 Q_\nu e^\nu = (\nu^2 - \mu^2) Y_0 + \mu\nu(\nu - \mu) Y_1 + \sum_{n=3}^\infty \frac{\nu^2 \mu^n - \mu^2 \nu^n}{n!} Y_n$$

因为配置 $\mu > \nu > 0$，对于任意 $n \ge 3$，显然有 $\nu^2 \mu^n - \mu^2 \nu^n = \mu^2 \nu^2 (\mu^{n-2} - \nu^{n-2}) > 0$。因此将求和项直接舍去可以得到绝对安全的严密下界不等式：

$$\nu^2 Q_\mu e^\mu - \mu^2 Q_\nu e^\nu \ge (\nu^2 - \mu^2) Y_0 + \mu\nu(\nu - \mu) Y_1$$

代数移项化简，导出单光子产额下界的解析闭式解：

$$Y_1^L = \frac{\mu}{\mu\nu - \nu^2} \left[ Q_\nu e^\nu - Q_\mu e^\mu \frac{\nu^2}{\mu^2} - \frac{\mu^2 - \nu^2}{\mu^2} Y_0 \right]$$

同理，结合误码率方程式，导出单光子误码率上界 $e_1^U$ 的严密解析表达式：

$$e_1^U = \frac{E_\nu Q_\nu e^\nu - e_0 Y_0}{\nu Y_1^L}$$

其中背景暗计数误码率 $e_0$ 在完全随机无偏假定下严格取 $1/2$。

QubitLab 在端侧运行上述公式时，对输入变量的浮点有效数字与物理边界执行毫秒级校验。一旦信道由于物理折弯导致传输衰减过大，使得计算出的 $Y_1^L \le 0$ 或提取出的安全密钥率 $R \le 0$，系统即刻在用户界面标定“信道处于窃听或阻断不安全状态”，拒绝向下游输出任何未提纯密钥。

4.4 密钥后处理核心阶段：筛选、纠错与保密放大

物理信道光子探测仅仅完成了量子密钥分发的前半程，由此得到的原始序列被称为“筛选密钥（Sifted Key）”。为了使其演化为具备密码学可用性的“最终保密密钥（Final Secret Key）”，必须在端侧执行三阶段高强度经典后处理：

1. 基矢比对筛选（Sifting）

通信双方通过公开经典网络广播各自为每个时间窗所选取的测量基矢名称（$+$ 或 $\times$）。凡是基矢选择不一致的脉冲事件，直接从本地存储区物理抹除；基矢一致且成功触发探测的事件构成筛选密钥。

2. 双向信息协商纠错（Information Reconciliation / Error Correction）

由于光纤杂散光、波长抖动及探测器暗计数的微观干扰，筛选密钥中不可避免地存在约 $1\% \sim 3\%$ 的本底随机误码。QubitLab 内置了**级联纠错算法（Cascade Protocol）**的轻量级端侧实现：

• 双方协商确定分块长度 $k_1 \approx 0.73 / \mathrm{QBER}$，将筛选密钥划分为若干等长子块；
• 计算各个子块的奇偶校验位（Parity Bit）并通过信道公开交换；
• 若某子块奇偶校验位不一致，则证明该块内存在奇数个错误比特；
• 双方在块内执行二分递归搜索（Binary Search），每轮递归仅需交换一个比特的奇偶值，以 $O(\log_2 k)$ 的极低交互代价精确定位并纠正错误位；
• 随后执行多轮置乱（Permutation）与跨块二次校验，彻底纠正所有残留误码，确保双方密钥达到 $100\%$ 逻辑一致。

3. 保密放大（Privacy Amplification）

在纠错过程中，公开信道上交换奇偶校验位的行为不可避免地向外界泄漏了部分微弱信息；此外，根据 GLLP 分析，Eve 依靠 PNS 或局部拦截可能掌握了部分比特的非零互信息。为了将潜在窃听者掌握的互信息彻底压缩至指数级趋近于零，必须执行保密放大。

保密放大基于全域哈希函数簇（Universal Hashing Functions），在数学上通过构造一个随机抽取的**托普利茨矩阵（Toeplitz Matrix, $T$）**实现。设纠错后的一致密钥长度为 $n$，最终目标安全密钥长度为 $m$（根据 GLLP 密钥率计算出的安全上限 $m = \lfloor n \cdot R / Q_\mu \rfloor$）。构造 $m \times n$ 阶 Toeplitz 矩阵：

$$T_{ij} = t_{i-j}, \quad 1 \le i \le m, 1 \le j \le n$$

该矩阵完全由首行与首列共 $m + n - 1$ 个随机位唯一确定。最终安全密钥 $K_{\mathrm{final}}$ 通过矩阵与纠错密钥矢量 $K_{\mathrm{reconciled}}$ 在有限域 $\mathrm{GF}(2)$ 上的模 2 矩阵向量乘法生成：

$$K_{\mathrm{final}} = T \cdot K_{\mathrm{reconciled}} \pmod 2$$

依据剩余哈希引理（Leftover Hash Lemma），经过保密放大后，Eve 关于最终密钥 $K_{\mathrm{final}}$ 的信息泄露量满足上界：

$$I(K_{\mathrm{final}}; E) \le 2^{-(s - m)} / \ln 2$$

只要保留安全参数 $s - m \ge 64$，Eve 所能窃听到的互信息量即小于 $10^{-19}$ 比特，在工程和数学意义上达到了真正的无条件物理保密。

下表总结了诱骗态 BB84 协议全流程各阶段的关键输入、输出、端侧算法与安全指标：

五、 基于纠缠的 E91 协议与无条件安全性物理认证架构

5.1 阿图尔·埃克特 1991 纠缠 QKD 协议的交互时序

如果说诱骗态 BB84 协议主要依赖发射端的随机偏振调制与弱相干光强变换，那么阿图尔·埃克特在 1991 年提出的 E91 协议 则开创了利用量子力学最深层次的非定域纠缠特性直接构建保密通信的全新路径。

在 E91 协议架构中，密钥分发不再需要 Alice 单向向 Bob 发射偏振光子，而是由一个位于信道中间的独立纠缠源（Entanglement Source）连续向两侧发射处于自旋单态 $|\Psi^-\rangle$ 的纠缠光子对。Alice 接收到每对光子中的粒子 1，Bob 接收到粒子 2。

与 BB84 仅选用两组正交基不同，E91 协议要求双方各自配置三组非正交的测量方向角度。标准的 E91 角度设置如下：

• Alice 侧的可选测量基（极角定义）：
  • $a_1 = 0$
  • $a_2 = \frac{\pi}{4} = 45^\circ$
  • $a_3 = \frac{\pi}{2} = 90^\circ$
• Bob 侧的可选测量基（极角定义）：
  • $b_1 = \frac{\pi}{4} = 45^\circ$
  • $b_2 = \frac{\pi}{2} = 90^\circ$
  • $b_3 = \frac{3\pi}{4} = 135^\circ$

在每个时钟脉冲周期内，Alice 与 Bob 各自完全独立、随机地从其三组测量基中挑选一个角度进行单光子偏振测量，记录下测量方向索引与实测二值结果（$\pm 1$）。

测量完成后，双方在公开信道上比对所使用的测量基索引。这些事件被严格分为两组互不相交的物理子集：

1. 相同测量基集合（Key Generation Set）： 当双方碰巧选用了相同几何方向时（共有两组重合组合：$(a_2, b_1)$ 对应 $45^\circ$；$(a_3, b_2)$ 对应 $90^\circ$）。根据单态的完全反关联性质（$E = -\cos(0) = -1$），两者的测量结果百分之百完全相反。Alice 与 Bob 只需将其中一方的比特取反，即可直接转换为共享的秘密密钥。
2. 不同测量基集合（Bell Violation Test Set）： 当双方选用了不同的测量方向组合时（共涵盖 $(a_1, b_1), (a_1, b_3), (a_3, b_1), (a_3, b_3)$ 四种关键组合），这四个组合的角度差恰好精确匹配 CHSH 不等式的极大违背角配置：

$$|a_1 - b_1| = \left| 0 - \frac{\pi}{4} \right| = \frac{\pi}{4}$$ $$|a_1 - b_3| = \left| 0 - \frac{3\pi}{4} \right| = \frac{3\pi}{4}$$ $$|a_3 - b_1| = \left| \frac{\pi}{2} - \frac{\pi}{4} \right| = \frac{\pi}{4}$$ $$|a_3 - b_3| = \left| \frac{\pi}{2} - \frac{3\pi}{4} \right| = \frac{\pi}{4}$$

双方公开这四个方向组合下的所有测量数据，计算实测的 CHSH 参量 $S$。

5.2 窃听攻击对纠缠态纯度的退相干破坏与 CHSH 参数压制

E91 协议最震撼的物理魅力在于：协议的安全证明直接等价于贝尔不等式的违背检验。

假设存在一个窃听者 Eve 试图获取有关密钥的信息。Eve 可以尝试在光纤链路中截获其中一个光子，甚至可以构造任何包含高维辅助希尔伯特空间的复合系统，使自身的探测探针与发往 Bob 的光子发生微观保酉相互作用。然而，根据量子力学的单配性原理（Monogamy of Entanglement），一个两能级量子系统与第三个系统的纠缠必然以削弱其与原始系统的纠缠为直接代价。

若 Eve 与系统发生纠缠，Alice 与 Bob 所共享的复合双光子约化态将从理想的纯态 $|\Psi^-\rangle\langle\Psi^-|$ 不可逆地退化为混合态 $\rho_{AB}$。考察极限情况：如果 Eve 试图通过完全投影测量完全掌握粒子状态，此时系统的量子相干相消项（Off-diagonal coherence terms）将全部坍缩为零，复合系统彻底退化为一个仅仅包含经典统计关联的混合态：

$$\rho_{\mathrm{collapsed}} = \frac{1}{2} |01\rangle\langle 01| + \frac{1}{2} |10\rangle\langle 10|$$

在此经典混合态下计算四个检验方向的关联函数，其代数结构必然完全落入局域实在论隐变量理论的范畴，此时实测的 CHSH 参量必然严格受限于：

$$S_{\mathrm{measured}} \le 2$$

只有当信道中绝对不存在任何窃听干涉时，实测参量才能达到理想的齐列尔森上界 $S = 2\sqrt{2} \approx 2.828$。如果信道中存在部分窃听干涉，导致纯态以权重 $F$（保真度 Fidelity）与各向同性白噪声混合（维尔纳态 Werner State）：

$$\rho_{\mathrm{Werner}} = F |\Psi^-\rangle\langle\Psi^-| + \frac{1-F}{4} I_4$$

其实测 CHSH 参量与保真度 $F$ 之间保持严格的线性解析正比关系：

$$S(F) = 2\sqrt{2} \cdot F$$

这赋予了通信双方一种前所未有的物理安全度量标尺：无需对信道内部或探测器内部的物理细节做任何先验信任假设，通信双方只需计算公开比对数据的 $S$ 参量。一旦实测 $S > 2$，数学上便不可辩驳地证明了两个粒子之间必然存在非定域量子纠缠，没有任何第三方能够在这段强关联中掌握超过局域实在论所允许的物理信息。

5.3 设备无关量子密钥分发（DI-QKD）的理论萌芽

E91 协议对现代密码学体系最深远的影响，在于它直接孕育了当代量子保密通信的最高安全圣殿——设备无关量子密钥分发（Device-Independent QKD, DI-QKD）。

在传统的 BB84 协议中，安全证明依赖于对发射端与接收端物理硬件的严格数学假设：假定 Alice 的偏振调制器旋转角度完全精准无误，假定 Bob 的偏振分束棱镜不存在侧信道泄漏，假定单光子探测器没有被恶意强激光注入攻击（Detector Blinding Attack）。然而在现实工程中，硬件的微观制造公差与物理老化极易导致假设失效，从而被旁路攻击利用。

而在基于 E91 的 DI-QKD 框架中，Alice 与 Bob 甚至可以将手头的测量仪器全部视为来自恶意供应商制造的完全不可信“黑盒”。黑盒表面只有输入按钮（选择方向 $a_i, b_j$）与输出指示灯（显示测量值 $\pm 1$）。只要通过黑盒的大量统计输出所计算出的 CHSH 参量确凿满足 $S > 2$，贝尔定理在公理层面上直接粉碎了黑盒在本地伪造经典关联以欺骗用户的可能性。这种基于物理第一性原理的公理化安全，展现了微观量子世界与信息安全理论深度交融的终极魅力。

六、 QubitLab 端侧离散数值仿真引擎架构与 Swift/Metal 高性能计算流水线

6.1 离散希尔伯特空间复数向量运算内核实现

在移动终端受限的算力与电池能耗约束下，实时解算多体量子态演化与蒙特卡洛统计收敛，需要极为高效的底层代数计算架构。QubitLab 引擎彻底摒弃了解释型脚本语言包装的臃肿科学库，采用纯原生 Swift 语言深度绑定 Apple Silicon 芯片硬件加速单元，构建了自包含的微观物理计算内核。

系统底层将复数标量严格抽象为遵循 IEEE 754 标准的高精浮点结构体 Complex，利用现代 64 位单指令多数据流（SIMD）矢量寄存器实现实部与虚部的并行吞吐。下述核心代码展示了双量子比特希尔伯特空间态矢量与密度矩阵部分迹的端侧高性能实现：

import Foundation
import simd

/// 高精度复数结构体，支持连续量子概率幅运算
public struct Complex: Equatable, Sendable {
    public var real: Double
    public var imag: Double

    @inline(__always)
    public init(_ real: Double, _ imag: Double = 0.0) {
        self.real = real
        self.imag = imag
    }

    @inline(__always)
    public var magnitudeSquared: Double {
        real * real + imag * imag
    }

    @inline(__always)
    public var magnitude: Double {
        sqrt(magnitudeSquared)
    }

    @inline(__always)
    public static func + (lhs: Complex, rhs: Complex) -> Complex {
        Complex(lhs.real + rhs.real, lhs.imag + rhs.imag)
    }

    @inline(__always)
    public static func - (lhs: Complex, rhs: Complex) -> Complex {
        Complex(lhs.real - rhs.real, lhs.imag - rhs.imag)
    }

    @inline(__always)
    public static func * (lhs: Complex, rhs: Complex) -> Complex {
        Complex(
            lhs.real * rhs.real - lhs.imag * rhs.imag,
            lhs.real * rhs.imag + lhs.imag * rhs.real
        )
    }
}

/// 双量子比特复合态矢量与约化密度矩阵算子
public final class TwoQubitState: Sendable {
    /// 4 维复态矢量概率幅 [c00, c01, c10, c11]
    public private(set) var amplitudes: [Complex]

    public init(amplitudes: [Complex]) {
        precondition(amplitudes.count == 4, "双量子比特复合希尔伯特空间维数必须为 4")
        self.amplitudes = amplitudes
        self.normalize()
    }

    /// 态矢量模长归一化
    public func normalize() {
        let normSq = amplitudes.reduce(0.0) { $0 + $1.magnitudeSquared }
        guard normSq > 1e-15 else { return }
        let invNorm = 1.0 / sqrt(normSq)
        for i in 0..<4 {
            amplitudes[i] = Complex(amplitudes[i].real * invNorm, amplitudes[i].imag * invNorm)
        }
    }

    /// 计算对子系统 B 求部分迹后的子系统 A 约化密度矩阵 (2x2)
    /// rho_A = Tr_B(rho_AB)
    public func reducedDensityMatrixA() -> (rho00: Double, rho01: Complex, rho11: Double) {
        // rho_A(0,0) = |c00|^2 + |c01|^2
        let r00 = amplitudes[0].magnitudeSquared + amplitudes[1].magnitudeSquared
        // rho_A(1,1) = |c10|^2 + |c11|^2
        let r11 = amplitudes[2].magnitudeSquared + amplitudes[3].magnitudeSquared
        // rho_A(0,1) = c00 * c10_conj + c01 * c11_conj
        let term1 = Complex(
            amplitudes[0].real * amplitudes[2].real + amplitudes[0].imag * amplitudes[2].imag,
            amplitudes[0].imag * amplitudes[2].real - amplitudes[0].real * amplitudes[2].imag
        )
        let term2 = Complex(
            amplitudes[1].real * amplitudes[3].real + amplitudes[1].imag * amplitudes[3].imag,
            amplitudes[1].imag * amplitudes[3].real - amplitudes[1].real * amplitudes[3].imag
        )
        let r01 = term1 + term2
        return (r00, r01, r11)
    }

    /// 计算子系统 A 的冯·诺依曼纠缠熵 (ebit)
    public func entanglementEntropy() -> Double {
        let (r00, r01, r11) = reducedDensityMatrixA()
        let offDiagSq = r01.magnitudeSquared
        // 解析求解 2x2 厄米矩阵本征值: lambda = (1 +/- sqrt(1 - 4*det)) / 2
        let det = r00 * r11 - offDiagSq
        let disc = max(0.0, 1.0 - 4.0 * det)
        let sqrtDisc = sqrt(disc)
        let lambda1 = max(1e-15, min(1.0, (1.0 + sqrtDisc) * 0.5))
        let lambda2 = max(1e-15, min(1.0, (1.0 - sqrtDisc) * 0.5))

        let s1 = -lambda1 * (log2(lambda1))
        let s2 = (lambda2 > 1e-12) ? (-lambda2 * log2(lambda2)) : 0.0
        return s1 + s2
    }
}

6.2 诱骗态仿真器参数调谐与实时光纤损耗模型

在 QubitLab 教学沙盒的“间谍捕获游戏（SpyCatcherGameView）”中，用户可以动态调整通信光纤物理跨段长度（$L \in [0, 100\ \mathrm{km}]$）、激光器平均脉冲强度（$\mu \in [0.1, 1.0]$）以及窃听者介入比例，实时观测信道指标的动态响应。

为了在无网络离线环境下提供严谨精确的物理反馈，引擎建立了基于通信波段标准单模光纤（G.652 光纤，波长 $1550\ \mathrm{nm}$）的透射率积分模型：

$$\eta_\mathrm{channel}(L) = 10^{-\frac{\alpha L}{10}}$$

其中光纤衰减常数设定为典型工业标准 $\alpha = 0.2\ \mathrm{dB/km}$。单光子脉冲在 Bob 端的总体探测产额模型严密计入了光纤传输损耗、接收端光学组件固有插损 $\eta_\mathrm{Bob} \approx 0.70$、雪崩光电探测器量子转换效率 $\eta_\mathrm{det} \approx 0.20$ 以及探测器环境暗计数率 $Y_0 \approx 1.0 \times 10^{-5}$：

$$Y_n(L) = 1 - (1 - Y_0)(1 - \eta_\mathrm{total})^n, \quad \eta_\mathrm{total} = \eta_\mathrm{channel}(L) \cdot \eta_\mathrm{Bob} \cdot \eta_\mathrm{det}$$

下述 Swift 模块实现了诱骗态增益解算与 GLLP 安全密钥率的实时计算流水线：

/// 诱骗态 BB84 物理信道解算与 GLLP 安全密钥率评估引擎
public struct DecoyStateEngine: Sendable {
    public let fiberLossDBPerKM: Double = 0.20
    public let bobOpticsEfficiency: Double = 0.70
    public let detectorQuantumEfficiency: Double = 0.20
    public let darkCountRateY0: Double = 1.0e-5
    public let opticalAlignmentErrorE0: Double = 0.015 // 光学基矢对准误差 1.5%

    public init() {}

    /// 二元香农熵函数
    @inline(__always)
    public func h2(_ x: Double) -> Double {
        guard x > 1e-15 && x < (1.0 - 1e-15) else { return 0.0 }
        return -x * log2(x) - (1.0 - x) * log2(1.0 - x)
    }

    /// 评估给定光纤跨段与光强下的宏观物理增益与安全密钥率
    public func evaluateChannel(
        fiberDistanceKM: Double,
        laserMu: Double,
        decoyNu: Double,
        errorCorrectionFactor: Double = 1.16
    ) -> (gainMu: Double, qberMu: Double, secureKeyRate: Double) {
        let channelTrans = pow(10.0, -(fiberLossDBPerKM * fiberDistanceKM) / 10.0)
        let totalEta = channelTrans * bobOpticsEfficiency * detectorQuantumEfficiency

        // 单光子与二光子产额计算
        let y1 = darkCountRateY0 + totalEta
        let y2 = darkCountRateY0 + 1.0 - pow(1.0 - totalEta, 2.0)

        // 宏观总增益 Q_mu
        let gainMu = darkCountRateY0 + (1.0 - exp(-laserMu * totalEta))

        // 宏观量子误码率 E_mu
        let darkErrors = 0.5 * darkCountRateY0
        let opticalErrors = opticalAlignmentErrorE0 * (1.0 - exp(-laserMu * totalEta))
        let qberMu = (darkErrors + opticalErrors) / max(gainMu, 1e-15)

        // 单光子有效增益 Q1 与单光子误码率 e1
        let q1 = laserMu * exp(-laserMu) * y1
        let e1 = (0.5 * darkCountRateY0 + opticalAlignmentErrorE0 * totalEta) / max(y1, 1e-15)

        // GLLP 安全密钥率公式 (比特/脉冲)
        let costErrorCorrection = gainMu * errorCorrectionFactor * h2(qberMu)
        let privacyAmplification = q1 * (1.0 - h2(e1))
        let rawRate = 0.5 * (privacyAmplification - costErrorCorrection)
        let secureKeyRate = max(0.0, rawRate)

        return (gainMu, qberMu, secureKeyRate)
    }
}

6.3 蒙特卡洛伪随机数源与端侧硬件熵注入

为了在移动端忠实复现微观量子态测量的绝对客观随机性，QubitLab 拒绝使用由固定种子驱动的简易伪随机发生器（如 C 语言的 rand() 或通用的线性同余生成器 LCG）。

引擎调度操作系统内核加密级安全伪随机数源（Cryptographically Secure Pseudorandom Number Generator, CSPRNG，如 Apple 平台的 SecRandomCopyBytes）。该随机数发生器在底层直接汇聚移动设备芯片片上模拟热噪声环形振荡器（Ring Oscillators）、电容触控微抖动以及微机电系统（MEMS）陀螺仪高频物理环境熵。在执行偏振基矢选择与波函数坍缩抽样时，确保每一个模拟生成的单光子事件均具备物理层面的抗预测性与统计均匀性。

同时，针对高校课堂演示与特定学术复盘场景，引擎设计了解耦的“确定性回放模式（Deterministic Playback Mode）”。教师或研究者可显式指定伪随机种子，使包含数万次量子测量的完整实验流程能够以百分之百相同的时序特征精确复现，高度契合了严谨学术科研的对照实验要求。

七、 100% 离线优先（Offline-First）架构哲学与端侧隐私数据主权

7.1 为什么科学探索与密码原语沙盒必须恪守 100% 离线

在当代商业化移动应用与云原生技术浪潮中，“软件即服务（SaaS）”与“全量云端上报”几乎成为行业惯性思维。大量教育类、科研辅助类软件动辄要求用户注册登录，将用户的每一次点击、练习答题记录、物理参数调节甚至本地生成的密钥材料同步上传至云端数据库。

在涉及量子物理探索、前沿密码学学习与高安全沙盒构建的专业领域，这种重度依赖云端的模式存在深刻的技术伦理缺陷与安全风险：

1. 实验参数与研究思路的隐蔽暴露：研究人员在沙盒中调谐的特定光纤拓扑、纠缠保真度阈值以及密码算法适配方案，属于高度敏感的前沿智力资产。若通过明文或弱加密 API 频繁向云端传输，极易在传输链路或第三方云服务器上发生数据泄露；
2. 潜在“先存后解（HNDL）”侧信道风险：如果教学沙盒在演示量子密钥协商时，将模拟生成的中间态比特或基矢信息回传给云端记录日志，这些流量一旦被公共互联网骨干节点非法镜像拦截，在未来量子算力实用化时将面临被全局破解的严重后患；
3. 极端工业与国防涉密环境适应性：真正的高安全密码测评、特种作业现场以及涉密科学实验室，往往处于严密的法拉第电磁屏蔽网或严格的物理隔绝（Air-Gapped）断网环境中。依赖云端鉴权或云端算力支持的软件在这些关键场景下完全沦为不可用的废品。

密旋科技（KeyRotate Technology）坚信：科技探索的纯粹性必须以极致的隐私安全为后盾，端侧离线架构是捍卫用户数据主权的唯一坚固堡垒。

7.2 零云端依赖、零遥测收集与无网络沙盒环境工程实现

在 QubitLab 以及密旋科技旗下的全线产品（TagWorks 电气安全合规、SlingGuard 起重法定核验）中，**“零遥测（Zero-Telemetry）”与“纯本地保护存储”**被作为最高系统设计原则写入工程铁律：

1. 操作系统级网络套接字彻底剥离： 在应用发布构建阶段，系统工程清单文件（如 iOS 的 Entitlements.plist）中完全不声明任何网络通信能力权限。可执行二进制文件中彻底剔除所有外部网络请求框架，从操作系统底层系统调用（socket、connect）层面断绝了进程发起任何网络通信的物理可能。
2. 严禁任何第三方商业分析与广告 SDK 装配： 应用内部第三方商业依赖装配率严格恒等于 $0.0\%$。彻底杜绝任何统计分析黑盒（如 Firebase、Adjust、Bugly 等）在后台搜集设备硬件指纹、IP 地址或隐蔽回传剪贴板信息的漏洞通道。
3. 端侧受保护存储与内存安全擦除： 所有实验关卡进度、自定义物理模拟拓扑、生成的仿真密钥凭证，全部存储在端侧沙盒内部。核心密码学计算所分配的复数数组与密钥比特内存页，在生命周期结束退出作用域时，调用 memset_s 安全清零原语执行强制内存覆盖，杜绝冷启动内存残余转储刺探。

下表系统对比了传统云端教学软件与密旋科技 100% 离线端侧物理引擎的系统工程架构差异：

八、 移动端性能基准、数值截断误差与算法复杂度实验评估

8.1 贝尔态多体张量积计算在 Apple Silicon 上的微秒级延迟评测

为了验证 QubitLab 离散物理内核在现代智能终端上的运行能效与计算延迟，本研究在配备 Apple M4 / A18 Pro 芯片的端侧设备上开展了严密的微基准测试（Micro-benchmark）。测试覆盖双量子比特态矢量任意保酉旋转变换、四体贝尔态正交分解、以及密度矩阵部分迹与纠缠熵计算。

测试通过连续执行 $100,000$ 次独立演化循环，记录其单次操作的平均耗时与能耗吞吐。实验测量结果如下表所示：

测试数据表明，即便在单次交互中完成包含一万对纠缠粒子的全量 CHSH 统计采样，端侧耗时亦严格控制在 $0.5$ 毫秒以内，仅占单个屏幕渲染刷新帧周期（$16.6\ \mathrm{ms}$）的不足 $3\%$，充分保障了手势拖曳旋钮时物理曲线与粒子流动画的极致丝滑呈现。

8.2 单光子产额下界估计在不同浮点精度下的数值稳定性

在诱骗态 BB84 算法中，单光子产额下界公式包含小量之间的差分相除：

$$Y_1^L = \frac{\mu}{\mu\nu - \nu^2} \left[ Q_\nu e^\nu - Q_\mu e^\mu \frac{\nu^2}{\mu^2} - \frac{\mu^2 - \nu^2}{\mu^2} Y_0 \right]$$

当光纤传输距离超过 80 公里时，透射率 $\eta$ 下跌至 $10^{-3}$ 以下，项 $Q_\nu e^\nu$ 与 $Q_\mu e^\mu (\nu^2/\mu^2)$ 的前五位有效数字高度重合。如果在计算中采用单精度浮点数（IEEE 754 Float32，尾数位仅 24 位，约 7 位有效十进制精度），极易引发严重的灾难性浮点截断下溢（Catastrophic Cancellation），导致计算出的 $Y_1^L$ 发生剧烈数值振荡甚至违背物理事实地跌为负数。

QubitLab 内核在所有光纤信道积分与诱骗态代数求解环节，强制统一采用标准双精度浮点数（Float64，尾数位 53 位，约 16 位有效十进制精度）。数值稳定性对比测试显示：在传输距离长达 $120\ \mathrm{km}$ 的极限高损耗工况下，Float64 求解器的相对数值截断误差严格压制在 $1.2 \times 10^{-12}$ 以内，确保了安全密钥率曲线在整个物理跨度内的平滑单调收敛。

结论与量子保密技术演进展望

量子力学的发展史，是一部人类认知不断挣脱宏观经典机械观束缚、向物理世界更深层次实在逼近的壮阔史诗。从爱因斯坦与玻尔关于局域实在论的长达半个世纪的深邃论战，到贝尔不等式与 CHSH 关联检验在物理实验室中的确立，再到量子不可克隆定理催生出的 BB84、E91 以及诱骗态量子保密通信协议，微观物理学的奇特性质不仅重塑了人类对宇宙客观规律的哲学认知，更在当代信息安全体系面临经典计算破解与量子计算“先存后解（HNDL）”威胁的关键时刻，提供了立足于自然界根本法则的终极安全范式。

密旋科技（KeyRotate Technology）通过自主研发 QubitLab 交互式微观物理与离散密码仿真沙盒，系统践行了“用物理认知前沿，用离线捍卫隐私”的品牌核心定位。系统将深奥繁复的高维复希尔伯特空间张量积代数、算符本征投影、马吕斯定律微观展开、泊松分布弱相干光强调制以及 GLLP 渐进安全密钥率方程，完整浓缩并高效运行于用户本地掌上智能终端的保护沙盒之中。通过彻底剥离网络依赖与遥测机制，QubitLab 在赋予广大技术人员、学者与学生直观可触、即时反馈的科学探索体验的同时，以绝对的系统工程纪律捍卫了端侧数据主权与学习隐私。

展望未来，随着国际后量子密码标准化进程（如 NIST FIPS 203 ML-KEM、FIPS 204 ML-DSA、FIPS 205 SLH-DSA）的全面落地与 IETF RFC 9370 密码敏捷混合协商协议的纵深演进，基于物理特性的量子密钥分发与基于数学晶格难解问题的后量子算法必将走向深度的“双轨融合与多重混合防御”。密旋科技将继续深耕端侧离线高安全计算架构，不断拓宽微观物理仿真与工业特种安全合规的数字化边界，为构建透明、可信、不可篡改且具有永久前向安全性的下一代数据防御生态奠定坚实的系统工程基石。

参考文献

• [1] Bennett, C. H., & Brassard, G. (1984). Quantum cryptography: Public key distribution and coin tossing. In Proceedings of IEEE International Conference on Computers, Systems and Signal Processing, Bangalore, India (pp. 175-179). http://researcher.watson.ibm.com/researcher/files/us-bennetc/BB84highest.pdf
• [2] Ekert, A. K. (1991). Quantum cryptography based on Bell’s theorem. Physical Review Letters, 67(6), 661-663. https://doi.org/10.1103/PhysRevLett.67.661
• [3] Clauser, J. F., Horne, M. A., Shimony, A., & Holt, R. A. (1969). Proposed experiment to test local hidden-variable theories. Physical Review Letters, 23(15), 880-884. https://doi.org/10.1103/PhysRevLett.23.880
• [4] Bell, J. S. (1964). On the Einstein Podolsky Rosen paradox. Physics Physique Fizika, 1(3), 195-200. https://doi.org/10.1103/PhysicsPhysiqueFizika.1.195
• [5] Hwang, W.-Y. (2003). Quantum key distribution with high loss: toward global secure communication. Physical Review Letters, 91(5), 057901. https://doi.org/10.1103/PhysRevLett.91.057901
• [6] Lo, H.-K., Ma, X., & Chen, K. (2005). Decoy state quantum key distribution. Physical Review Letters, 94(23), 230501. https://doi.org/10.1103/PhysRevLett.94.230501
• [7] Wang, X.-B. (2005). Beating the photon-number-splitting attack in practical quantum cryptography. Physical Review Letters, 94(23), 230503. https://doi.org/10.1103/PhysRevLett.94.230503
• [8] Gottesman, D., Lo, H.-K., Lütkenhaus, N., & Preskill, J. (2004). Security of quantum key distribution with imperfect devices. Quantum Information & Computation, 4(5), 325-360. https://doi.org/10.26421/QIC4.5-1
• [9] Wootters, W. K., & Zurek, W. H. (1982). A single quantum cannot be cloned. Nature, 299(5886), 802-803. https://doi.org/10.1038/299802a0
• [10] Shor, P. W., & Preskill, J. (2000). Simple proof of security of the BB84 quantum key distribution protocol. Physical Review Letters, 85(2), 441-444. https://doi.org/10.1103/PhysRevLett.85.441
• [11] Cirel’son, B. S. (1980). Quantum generalizations of Bell’s inequality. Letters in Mathematical Physics, 4(2), 93-100. https://doi.org/10.1007/BF00417500
• [12] Brassard, G., & Salvail, L. (1994). Secret-key reconciliation by public discussion. In Advances in Cryptology — EUROCRYPT ’93 (pp. 410-423). Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-48285-7_35
• [13] Acín, A., Brunner, N., Gisin, N., Massar, S., Pironio, S., & Scarani, V. (2007). Device-independent security of quantum cryptography against collective attacks. Physical Review Letters, 98(23), 230501. https://doi.org/10.1103/PhysRevLett.98.230501