噪声模型理论

量子噪声建模的理论基础，包括 CPTP 映射、Kraus 表示、量子误差通道以及 NISQ 设备的实际噪声考虑。本指南将数学形式化与 pyqpanda3 的噪声仿真能力联系起来。

---

噪声建模的必要性

真实的量子计算机在噪声环境中运行。与经典比特不同，量子比特是脆弱的量子系统，会通过与环境的相互作用而失去其量子特性。噪声建模对于以下方面至关重要：

1. 精确仿真：在运行实验之前预测真实硬件行为
2. 错误缓解：设计降低噪声影响的策略
3. 错误纠正：理解噪声以构建合适的量子纠错码
4. 基准测试：将硬件性能与理论噪声模型进行比较

在 pyqpanda3 中，NoiseModel 类和相关的误差通道为含噪声仿真提供了完整的框架。

---

量子操作作为映射

密度矩阵形式化

在讨论噪声之前，我们需要用于描述量子态的密度矩阵形式化：

纯态（Pure State）：可以写成 ket $|\psi ⟩$ 的态。其密度矩阵为：

$$\rho =|\psi ⟩⟨\psi |$$

混合态（Mixed State）：纯态的统计混合 $\{p_i,|{\psi }_i⟩\}$：

$$\rho =\sum _i{p}_i|{\psi }_i⟩⟨{\psi }_i|$$

其中 $p_i\ge 0$ 且 $\sum _i{p}_i=1$。

有效密度矩阵的性质：

• 厄米性（Hermitian）：$\rho ={\rho }^{†}$
• 半正定性：$\rho \ge 0$（所有特征值 $\ge 0$）
• 迹为 1：$\text{Tr}(\rho )=1$

纯态与混合态

当且仅当 $\text{Tr}({\rho }^2)=1$ 时态为纯态，$\text{Tr}({\rho }^2)<1$ 时为混合态。纯度（Purity） $\gamma =\text{Tr}({\rho }^2)$ 量化态的"纯度"：

• $\gamma =1$：纯态
• $1/d<\gamma <1$：混合态（$d$ 是希尔伯特空间维数）
• $\gamma =1/d$：最大混合态 $\rho =I/d$

pyqpanda3 的 DensityMatrix 类提供了 purity() 方法用于此计算。

---

完全正保迹（CPTP）映射

定义

量子通道（噪声模型）在数学上由完全正保迹（Completely Positive Trace-Preserving, CPTP）映射 $\mathcal{E}$ 描述，它将输入密度矩阵变换为输出密度矩阵：

$$\mathcal{E}:{\rho }_{\text{in}}\mapsto {\rho }_{\text{out}}$$

两个要求是：

完全正性（CP）：对于任意维度的任意辅助系统 $d_A$，扩展映射 $(\mathcal{E}\otimes I_{d_A})$ 将正算符映射为正算符：

$$\rho \ge 0⟹(\mathcal{E}\otimes I)(\rho )\ge 0\mathrm{\forall }\text{ 扩展}$$

完全正性比单纯正性更强 — 一个正但不完全正的映射在应用于纠缠态的子系统时可能产生非物理结果。

保迹性（TP）：密度矩阵的迹被保持：

$$\text{Tr}(\mathcal{E}(\rho ))=\text{Tr}(\rho )=1$$

这确保总概率守恒。

为什么是 CPTP？

CPTP 映射是量子噪声的正确数学框架，因为：

1. 物理性：CPTP 映射恰好是自然界中可以发生的变换集合（根据 Stinespring 扩张定理，任何 CPTP 映射都可以通过更大系统上的幺正演化实现）
2. 可组合性：两个 CPTP 映射的组合也是 CPTP 的 — 噪声模型可以链接
3. 凸性：CPTP 映射的概率混合也是 CPTP 的 — 不同噪声源正确组合

Stinespring 扩张

系统 ${\mathcal{H}}_S$ 上的每个 CPTP 映射 $\mathcal{E}$ 都可以实现为更大系统上的幺正演化：

$$\mathcal{E}(\rho )={\text{Tr}}_E[U_{SE}(\rho \otimes |e_0⟩⟨e_0|)U_{SE}^{†}]$$

其中 $U_{SE}$ 是联合系统-环境希尔伯特空间 ${\mathcal{H}}_S\otimes {\mathcal{H}}_E$ 上的幺正操作，$|e_0⟩$ 是固定的环境态，${\text{Tr}}_E$ 是对环境的偏迹。

这意味着所有物理噪声都可以理解为与未观测环境的幺正演化 — 随机性来自对（忽略）环境自由度的迹运算。

graph LR
    subgraph "Stinespring 扩张"
        S["系统 ρ"] --> US["U_SE"]
        E["环境 |e₀⟩"] --> US
        US --> SE["系统 ⊗ 环境"]
        SE --> TrE["Tr_E (求偏迹)"]
        TrE --> Out["ρ_out = ℰ(ρ)"]
    end

    style S fill:#e3f2fd,stroke:#1565c0
    style E fill:#fff3e0,stroke:#e65100
    style Out fill:#c8e6c9,stroke:#2e7d32

---

Kraus 表示

算符和表示

CPTP 映射最直观的表示是 Kraus（算符和）表示：

$$\mathcal{E}(\rho )=\sum _i{K}_i\rho K_i^{†}$$

其中 Kraus 算符 $\{K_i\}$ 满足完备关系：

$$\sum _i{K}_i^{†}{K}_i=I$$

这个条件确保保迹性：

$$\text{Tr}(\mathcal{E}(\rho ))=\text{Tr}\left(\sum _i{K}_i\rho K_i^{†}\right)=\text{Tr}\left(\rho \sum _i{K}_i^{†}{K}_i\right)=\text{Tr}(\rho )=1$$

Kraus 表示的性质

1. 非唯一性：同一通道可以有不同的 Kraus 表示，它们通过幺正混合关联：$K_i^{\prime }=\sum _j{u}_{ij}{K}_j$，其中 $U=(u_{ij})$ 是幺正矩阵。

2. 最小数量：所需 Kraus 算符的最小数量是通道的 Choi 秩。对于量子比特通道，该值在 1（幺正）到 4（完全去极化）之间。

3. 单个 Kraus 算符：如果只有一个 Kraus 算符，则通道是幺正的：$\mathcal{E}(\rho )=U\rho U^{†}$。

常见通道的 Kraus 算符

通道	Kraus 算符数量	Kraus 算符
幺正（$U$）	1	$K_0=U$
比特翻转	2	$K_0=\sqrt{1-p}I$, $K_1=\sqrt{p}X$
相位翻转	2	$K_0=\sqrt{1-p}I$, $K_1=\sqrt{p}Z$
去极化	4	$K_0=\sqrt{1-3p/4}I$, $K_{1,2,3}=\sqrt{p/4}\{X,Y,Z\}$
振幅阻尼	2	$K_0=\left(\begin{array}{cc}1& 0\\ 0& \sqrt{1-\gamma }\end{array}\right)$, $K_1=\left(\begin{array}{cc}0& \sqrt{\gamma }\\ 0& 0\end{array}\right)$
相位阻尼	2	$K_0=\left(\begin{array}{cc}1& 0\\ 0& \sqrt{1-\lambda }\end{array}\right)$, $K_1=\left(\begin{array}{cc}0& 0\\ 0& \sqrt{\lambda }\end{array}\right)$

pyqpanda3 的 quantum_info 模块中的 Kraus 类以此形式表示通道。通道还可以表示为 Chi、Choi、SuperOp 或 PTM，它们之间可以相互转换。

---

量子误差通道

去极化通道

去极化通道是最常用的噪声模型。以概率 $p$，态被替换为最大混合态 $I/2$：

$${\mathcal{E}}_{\text{depol}}(\rho )=(1-p)\rho +p\frac{I}{2}$$

在 Pauli 基中展开的等价形式为：

$${\mathcal{E}}_{\text{depol}}(\rho )=(1-\frac{3p}{4})\rho +\frac{p}{4}(X\rho X+Y\rho Y+Z\rho Z)$$

Kraus 算符（4 个算符）：

$$K_0=\sqrt{1-\frac{3p}{4}}I,K_1=\sqrt{\frac{p}{4}}X,K_2=\sqrt{\frac{p}{4}}Y,K_3=\sqrt{\frac{p}{4}}Z$$

性质：

• 对称性：平等对待三个轴
• 均匀缩小 Bloch 矢量：$\overrightarrow{r}\to (1-p)\overrightarrow{r}$
• 常用作最坏情况噪声模型
• $p=1$ 给出完全随机态 $I/2$

振幅阻尼通道

振幅阻尼通道建模能量耗散 — 激发态 $|1⟩$ 衰变到基态 $|0⟩$：

$${\mathcal{E}}_{\text{AD}}(\rho )=K_0\rho K_0^{†}+K_1\rho K_1^{†}$$

Kraus 算符：

$$K_0=\left(\begin{array}{cc}1& 0\\ 0& \sqrt{1-\gamma }\end{array}\right),K_1=\left(\begin{array}{cc}0& \sqrt{\gamma }\\ 0& 0\end{array}\right)$$

其中 $\gamma \in [0,1]$ 是阻尼概率（与 $T_1$ 弛豫时间相关）。

对基态的作用：

$$|0⟩⟨0|\to |0⟩⟨0|\text{（基态不变）}$$$$|1⟩⟨1|\to (1-\gamma )|1⟩⟨1|+\gamma |0⟩⟨0|\text{（衰变到基态）}$$

性质：

• 非幺正的：$\mathcal{E}(I)\ne I$（通道不保持恒等）
• 驱动态趋向 $|0⟩$
• 建模自发辐射、$T_1$ 衰变
• $\gamma =0$：恒等（无阻尼）；$\gamma =1$：完全弛豫到 $|0⟩$

与 $T_1$ 时间的关系：对于持续时间为 $t_g$ 的门：

$$\gamma =1-e^{-t_g/T_1}$$

相位阻尼（退相干）通道

相位阻尼通道建模相位相干性损失，不涉及能量损失：

$${\mathcal{E}}_{\text{PD}}(\rho )=K_0\rho K_0^{†}+K_1\rho K_1^{†}$$

Kraus 算符：

$$K_0=\left(\begin{array}{cc}1& 0\\ 0& \sqrt{1-\lambda }\end{array}\right),K_1=\left(\begin{array}{cc}0& 0\\ 0& \sqrt{\lambda }\end{array}\right)$$

其中 $\lambda \in [0,1]$ 是退相位概率。

使用 Pauli Z 的等价形式：

$${\mathcal{E}}_{\text{PD}}(\rho )=(1-\frac{p}{2})\rho +\frac{p}{2}Z\rho Z$$

其中 $p=1-e^{-t_g/T_{\varphi }}$。

性质：

• 幺正的：$\mathcal{E}(I)=I$
• 保持计算基：对角元素不变
• 抑制非对角元素（相干性）：${\rho }_{01}\to \sqrt{1-\lambda }{\rho }_{01}$
• 建模 $T_2$ 退相干
• $\lambda =0$：恒等；$\lambda =1$：完全退相位（经典混合）

与 $T_2$ 时间的关系：

$$\frac{1}{T_2}=\frac{1}{2T_1}+\frac{1}{T_{\varphi }}$$

其中 $T_{\varphi }$ 是纯退相位时间。

比特翻转通道

比特翻转通道随机施加 $X$ 门：

$${\mathcal{E}}_{\text{BF}}(\rho )=(1-p)\rho +pX\rho X$$

Kraus 算符：$K_0=\sqrt{1-p}I$, $K_1=\sqrt{p}X$

Bloch 球作用：以概率 $p$ 翻转 $z$ 分量：$r_z\to (1-2p)r_z$, $r_y\to (1-2p)r_y$, $r_x$ 不变。

关联：等价于 Hadamard 基中的相位翻转：$H\cdot {\mathcal{E}}_{\text{BF}}\cdot H={\mathcal{E}}_{\text{PF}}$。

相位翻转通道

相位翻转通道随机施加 $Z$ 门：

$${\mathcal{E}}_{\text{PF}}(\rho )=(1-p)\rho +pZ\rho Z$$

Kraus 算符：$K_0=\sqrt{1-p}I$, $K_1=\sqrt{p}Z$

Bloch 球作用：$r_z$ 不变, $r_x\to (1-2p)r_x$, $r_y\to (1-2p)r_y$。

比特-相位翻转通道

比特-相位翻转通道随机施加 $Y$ 门（同时进行比特翻转和相位翻转）：

$${\mathcal{E}}_{\text{BPF}}(\rho )=(1-p)\rho +pY\rho Y$$

Kraus 算符：$K_0=\sqrt{1-p}I$, $K_1=\sqrt{p}Y$

注意：$Y=iXZ$，因此这是比特翻转和相位翻转的组合（带有一个附加的全局相位）。

Pauli 噪声通道

通用 Pauli 噪声通道以各自的概率随机施加 Pauli 算符：

$${\mathcal{E}}_{\text{Pauli}}(\rho )=p_I\rho +p_{X}X\rho X+p_{Y}Y\rho Y+p_{Z}Z\rho Z$$

其中 $p_I+p_X+p_Y+p_Z=1$。

这是最通用的单量子比特 Pauli 通道，包括比特翻转、相位翻转和去极化作为特例。

热弛豫通道

热弛豫通道建模 $T_1$（振幅阻尼）和 $T_2$（退相位）过程在门持续时间 $t$ 内的综合效应：

$${\mathcal{E}}_{\text{thermal}}(\rho )=E_0\rho E_0^{†}+E_1\rho E_1^{†}+E_2\rho E_2^{†}$$

其中 Kraus 算符取决于弛豫时间 $T_1$、$T_2$ 和门时间 $t$：

$$p_{\text{reset}}=1-e^{-t/T_1}$$$$p_{\text{dephasing}}=1-e^{-t/T_2}$$

此通道对于超导量子比特硬件的真实仿真至关重要。

总结：误差通道比较

通道	参数	幺正？	物理来源
去极化	$p\in [0,1]$	是	通用噪声
振幅阻尼	$\gamma \in [0,1]$	否	$T_1$ 弛豫
相位阻尼	$\lambda \in [0,1]$	是	$T_2$ 退相位
比特翻转	$p\in [0,0.5]$	是	$X$ 错误
相位翻转	$p\in [0,0.5]$	是	$Z$ 错误
比特-相位翻转	$p\in [0,0.5]$	是	$Y$ 错误
Pauli 噪声	$p_I,p_X,p_Y,p_Z$	是	通用 Pauli 错误
热弛豫	$T_1,T_2,t$	否	综合弛豫

graph TD
    All["所有误差通道"] --> Unital["幺正通道"]
    All --> NonUnital["非幺正通道"]

    Unital --> PauliCh["Pauli 通道"]
    PauliCh --> BF["比特翻转"]
    PauliCh --> PF["相位翻转"]
    PauliCh --> BPF["比特-相位翻转"]
    PauliCh --> Depol["去极化"]
    PauliCh --> General["通用 Pauli"]

    Unital --> PD["相位阻尼"]

    NonUnital --> AD["振幅阻尼"]
    NonUnital --> TR["热弛豫"]
    AD --> TR

    style Unital fill:#e3f2fd,stroke:#1565c0
    style NonUnital fill:#fff3e0,stroke:#e65100

---

量子通道表示

单个量子通道可以用多种数学等价的方式表示。pyqpanda3 支持五种表示，均可相互转换：

Kraus 表示

$$\mathcal{E}(\rho )=\sum _i{K}_i\rho K_i^{†}$$

最直观的形式；直接展示物理操作。

超算符（Liouville）表示

通道表示为作用在矢量化密度矩阵上的 $d^2\times d^2$ 矩阵 $\mathcal{S}$：

$$\text{vec}(\mathcal{E}(\rho ))=\mathcal{S}\cdot \text{vec}(\rho )$$

其中 $\text{vec}(\cdot )$ 将矩阵列堆叠。在 pyqpanda3 中，这是 SuperOp 类。

Chi 矩阵表示

通道在 Pauli 基中展开：

$$\mathcal{E}(P_j)=\sum _i{\chi }_{ij}{P}_i$$

$\chi$ 矩阵（不要与 Choi 矩阵混淆）是基于 Pauli 基的表示。在 pyqpanda3 中，这是 Chi 类。

Choi 矩阵表示

Choi 矩阵定义为：

$$J(\mathcal{E})=(\mathcal{E}\otimes I)(|\mathrm{\Omega }⟩⟨\mathrm{\Omega }|)$$

其中 $|\mathrm{\Omega }⟩=\sum _i|i⟩\otimes |i⟩$ 是最大纠缠态。当且仅当 $J(\mathcal{E})\ge 0$（半正定）且 ${\text{Tr}}_1(J)=I$（偏迹条件）时，通道为 CPTP 的。在 pyqpanda3 中，这是 Choi 类。

Pauli 传输矩阵（PTM）

PTM 表示通道对广义 Pauli 基的作用：

$$R_{ij}=\frac{1}{d}\text{Tr}(P_i\mathcal{E}(P_j))$$

其中 $\{P_i\}$ 是归一化的 Pauli 算符。对于保持厄米性的通道，PTM 是实矩阵。在 pyqpanda3 中，这是 PTM 类。

表示之间的转换

在 pyqpanda3 中，所有五种表示都可以通过从一个构造另一个来相互转换：

graph TD
    Kraus["Kraus"]
    SuperOp["SuperOp"]
    Chi["Chi"]
    Choi["Choi"]
    PTM["PTM"]

    Kraus -->|constructor| SuperOp
    Kraus -->|constructor| Chi
    Kraus -->|constructor| Choi
    Kraus -->|constructor| PTM

    SuperOp -->|constructor| Kraus
    SuperOp -->|constructor| Chi
    SuperOp -->|constructor| Choi
    SuperOp -->|constructor| PTM

    Chi -->|constructor| Kraus
    Chi -->|constructor| SuperOp
    Chi -->|constructor| Choi
    Chi -->|constructor| PTM

    Choi -->|constructor| Kraus
    Choi -->|constructor| SuperOp
    Choi -->|constructor| Chi
    Choi -->|constructor| PTM

    PTM -->|constructor| Kraus
    PTM -->|constructor| SuperOp
    PTM -->|constructor| Chi
    PTM -->|constructor| Choi

    style Kraus fill:#e3f2fd,stroke:#1565c0
    style SuperOp fill:#c8e6d4,stroke:#2e7d32
    style Chi fill:#fff3e0,stroke:#e65100
    style Choi fill:#f3e5f5,stroke:#9c27b0
    style PTM fill:#ffebee,stroke:#c62828

---

NISQ 设备中的噪声

NISQ 噪声的特征

当前的含噪声中等规模量子（NISQ）设备表现出多种类型的噪声：

1. 门误差：不完美的门实现导致幺正误差和过度/不足旋转：

• 单量子比特门误差率：$\sim {10}^{-4}$ 到 ${10}^{-3}$
• 双量子比特门误差率：$\sim {10}^{-3}$ 到 ${10}^{-2}$
• 通常建模为去极化或相干旋转误差

2. 测量误差：读出（SPAM — 态制备和测量）误差：

• 分配误差率：$\sim 1\mathrm{\%}$ 到 $5\mathrm{\%}$
• $P(\text{读出 }0|\text{态 }1)$ 和 $P(\text{读出 }1|\text{态 }0)$ 可能不对称
• 建模为经典混淆矩阵

3. 退相干：量子比特随时间失去量子信息：

• $T_1$（弛豫时间）：能量衰变，通常 $20\mu s$ 到 $200\mu s$
• $T_2$（退相位时间）：相位相干性损失，通常 $10\mu s$ 到 $150\mu s$
• 门时间：$\sim 20ns$ 到 $200ns$

4. 串扰：操作期间量子比特之间不需要的耦合：

• 频谱串扰：同时门操作相互干扰
• 常开耦合：相邻量子比特之间的残余相互作用

5. 泄漏：量子比特态逃逸出计算子空间：

• 在 $|2⟩$ 或更高能级上的布居数
• 对超导 transmon 量子比特尤其成问题

使用 pyqpanda3 建模噪声

pyqpanda3 的 NoiseModel 类支持为真实仿真配置噪声：

from pyqpanda3.core import CPUQVM, NoiseModel, QProg, QGate
from pyqpanda3.core import H, CNOT, measure, depolarizing_error, GateType

# 创建噪声模型
noise = NoiseModel()

# 为所有单量子比特门添加去极化噪声
error_1q = depolarizing_error(0.001)
noise.add_all_qubit_quantum_error(error_1q, GateType.H)

# 为所有双量子比特门添加去极化噪声
error_2q = depolarizing_error(0.01)
noise.add_all_qubit_quantum_error(error_2q, GateType.CNOT)

# 运行含噪声仿真
qvm = CPUQVM()
prog = QProg()
prog << H(0) << CNOT(0, 1) << measure([0, 1], [0, 1])
qvm.run(prog, shots=1000, model=noise)
result = qvm.result()

噪声感知线路设计

设计抗噪声线路的原则：

1. 最小化双量子比特门：它们的误差率最高
2. 最小化线路深度：误差随每个时间步累积
3. 使用噪声感知转译：将线路映射到具有最佳连接性和相干性的量子比特
4. 利用对称性：设计噪声可以抵消的线路（如动态解耦）

错误缓解技术

错误缓解在不完全纠错的情况下减少噪声的影响：

零噪声外推（ZNE）：

• 在多个噪声水平运行线路（通过折叠门来放大噪声）
• 将结果外推到零噪声极限
• 通常使用线性或指数拟合

概率错误消除（PEC）：

• 表征噪声通道 $\mathcal{E}$
• 概率性地施加逆通道 ${\mathcal{E}}^{-1}$
• 以增加采样为代价产生无偏估计

测量误差缓解：

• 表征读出误差矩阵 $M$
• 将 $M^{-1}$ 应用于观测概率分布
• 需要为每个测量基进行校准运行

旋转读出误差消除（TREX）：

• 在测量前随机翻转量子比特（旋转）
• 将相关的读出误差转换为更简单的对角形式
• 更容易求逆

---

数学总结

关键关系

$$\boxed{{\displaystyle \text{CPTP 映射}}}⟷\boxed{{\displaystyle \text{Kraus 算符}}}⟷\boxed{{\displaystyle \text{Choi 矩阵}}}⟷\boxed{{\displaystyle \text{SuperOp}}}⟷\boxed{{\displaystyle \text{PTM}}}$$

Bloch 球图像

对于单量子比特通道，CPTP 映射对 Bloch 矢量 $\overrightarrow{r}$ 的作用为仿射变换：

$${\overrightarrow{r}}^{\prime }=M\overrightarrow{r}+\overrightarrow{t}$$

其中 $M$ 是 $3\times 3$ 实矩阵，$\overrightarrow{t}$ 是平移矢量。

通道类型	$M$	$\overrightarrow{t}$
去极化（$p$）	$(1-p)I$	$\overrightarrow{0}$
比特翻转（$p$）	$\text{diag}(1,1-2p,1-2p)$	$\overrightarrow{0}$
振幅阻尼（$\gamma$）	$\left(\begin{array}{ccc}\sqrt{1-\gamma }& 0& 0\\ 0& \sqrt{1-\gamma }& 0\\ 0& 0& 1-\gamma \end{array}\right)$	$(0,0,\gamma )$
相位阻尼（$\lambda$）	$\text{diag}(\sqrt{1-\lambda },\sqrt{1-\lambda },1)$	$\overrightarrow{0}$

幺正通道（$\overrightarrow{t}=\overrightarrow{0}$）保持最大混合态，仅缩小/旋转 Bloch 球。非幺正通道（$\overrightarrow{t}\ne \overrightarrow{0}$）同时平移 Bloch 球，驱动它趋向一个优先态。

---

另请参阅

• NoiseModel API — 噪声配置的 API 参考
• 量子通道表示 — Kraus、Chi、Choi、SuperOp、PTM 类
• DensityMatrix API — 密度矩阵操作
• 量子信息理论 — 态表示、保真度和距离度量
• 仿真工作流图 — 可视化噪声仿真流水线