1. 变分量子算法在酉扩张中的核心原理
量子计算中的酉扩张技术是实现非酉量子操作的关键方法。简单来说，酉扩张就像是为一个不完美的量子操作"搭建脚手架"——通过引入额外的量子比特（称为辅助比特），我们可以将这个不完美的操作嵌入到一个更大的、完美的酉操作中。这种技术的重要性在于，量子计算机本质上只能执行酉操作（可逆的量子门操作），而许多实际应用（如开放量子系统模拟）需要处理非酉操作。
传统Sz.-Nagy扩张算法虽然数学上完备，但在实际量子硬件实现时需要大量量子门，这在当前含噪声中等规模量子（NISQ）设备上几乎不可行。变分量子算法(VQA)的创新之处在于，它采用参数化量子电路来近似实现酉扩张，通过经典-量子混合优化找到最优参数，大幅降低了量子资源需求。
具体数学表述为：给定一个非酉算子M∈C^(d×d)，我们寻找一个作用于n+1个量子比特的酉算子U_SE(θ)∈U(2^(n+1))，使得其左上块近似αM：
˜M(θ) := (⟨0|_a ⊗ I) U_SE(θ) (|0⟩_a ⊗ I) ≈ αM
其中α是归一化常数，|0⟩_a表示辅助比特的初始状态。这个近似过程通过优化参数θ最小化代价函数实现。
2. 三种量子态测量方案对比
2.1 基于Choi态保真度的方案
Choi态是量子信道的一种表示方法，可以理解为将量子信道"编码"到一个量子态中。对于我们的非酉算子M，其对应的Choi态定义为：
Φ_M = (I ⊗ M)(|Φ^+⟩⟨Φ^+|)
其中|Φ^+⟩ = 1/√d Σ_j |j⟩_R|j⟩_S是最大纠缠态。变分量子电路产生的近似Choi态Φ_θ通过以下步骤制备：
- 准备两个n量子比特寄存器R2S2和辅助寄存器E
- 在R2S2上制备|Φ^+⟩态
- 对S2E应用参数化酉算子U_SE(θ)
- 后选择辅助寄存器E处于|0⟩态
代价函数采用保真度形式： L_Choi = 1 - F(Φ_θ, Φ_M)
保真度测量采用直接保真度估计(DFE)技术，通过采样Pauli字符串P∈P_2n进行估计：
F = Σ_(P∈P_2n) w_P μ_P
其中w_P = Tr(Φ_M P)^2 / Σ_Q Tr(Φ_M Q)^2是重要性权重，μ_P = Tr(Φ_θ P)/Tr(Φ_M P)是测量结果。
提示：在实际实验中，可以采用重要性采样策略——先用少量测量估计Tr(Φ_M P)的分布，然后对重要Pauli项进行精细测量，这能显著减少测量次数。
2.2 基于SWAP测试的方案
SWAP测试提供了一种完全相干的保真度测量方法。电路实现包含：
- 准备控制比特|+⟩态
- 在控制下执行寄存器对(R1↔R2)和(S1↔S2)的受控SWAP操作
- 测量控制比特的σ_x期望值
测量结果与保真度的关系为： 1/2 + Tr(Φ_θ Φ_M)/2
为降低电路深度，可采用破坏性SWAP测试变体：
- 对两个Choi态对应的量子比特对执行Bell基测量
- 通过经典后处理计算保真度
或者采用双拷贝随机测量方案：
- 对两个拷贝施加相同的随机局域Clifford门
- 在计算基下测量
- 使用标准重叠估计器计算Tr(ρσ)
2.3 基于状态回归的方案
当Choi态难以制备时，可采用回归型代价函数：
L_reg(θ) = 1/K Σ_(k=1)^K ||˜M(θ)|ψ_k⟩ - αM|ψ_k⟩||_2^2
其中{|ψ_k⟩}是一组探测态（如计算基态、乘积态等）。展开后包含三类矩阵元：
- ⟨ψ_j|˜M^† ˜M|ψ_j⟩：可通过量子电路直接测量
- ⟨ψ_j|αM^† αM|ψ_j⟩：经典已知或一次性估计
- ⟨ψ_j|˜M^† αM|ψ_j⟩：通过图S1(c)所示电路测量
3. 量子电路实现细节
3.1 二能级系统模拟
对于二能级开放系统，我们设计四层变分量子电路，每层包含：
- 参数化单量子比特旋转门R(θ_i)
- CNOT纠缠门
优化目标是最小化Frobenius范数： min_θ ||U^(2)(θ) - αM||_F
实验结果显示：
- 马尔可夫情况（γ_0t=6π）：最终误差2.6×10^-3
- 非马尔可夫情况（α=0.9）：最终误差3.3×10^-6
相比传统Sz.-Nagy扩张，门数量减少87.5%（单量子比特门）和90.5%（双量子比特门）。
3.2 多能级系统模拟
对于三能级和四能级系统，电路架构采用：
- 参数化R(θ_i)旋转门
- CZ纠缠门
稳态模拟结果（t→∞，α=0.5）：
- 三能级系统：误差5.0×10^-4
- 四能级系统：误差9.0×10^-5
资源对比（相比文献[S4,S5]）：
- 三能级系统：单量子比特门减少96.8%，双量子比特门减少97.7%
- 四能级系统：单量子比特门减少95.6%，双量子比特门减少96.8%
4. 超导量子处理器实验
实验在"百花"超导量子计算平台进行，主要参数：
- 量子比特数：>100
- 典型相干时间：T1≈50-85μs，T2≈40-50μs
- 门错误率：
- 单量子比特门：~5×10^-4
- 双量子比特门：~8×10^-3
- 读取错误率：~2×10^-2
关键优化策略：
- 利用硬件原生门集设计参数化电路
- 根据量子比特连通性优化门序列
- 采用动态解耦技术延长相干时间
5. Kraus算符方法扩展
酉扩张技术同样适用于Kraus算符描述的开放量子系统。动力学演化表示为：
ρ(t) = Σ_k M_k ρ M_k^†
每个Kraus算符M_k可通过独立量子电路实现。以二能级系统为例：
M_0 = [1 0; 0 √e^(-γt)] M_1 = [0 √(1-e^(-γt)); 0 0]
实验步骤：
- 将初始态ρ分解为纯态混合（如ρ=(|1⟩⟨1|+|+⟩⟨+|)/2）
- 对每个纯态和Kraus算符构建独立电路
- 通过测量结果重构ρ(t)
6. 实用技巧与注意事项
参数初始化策略：
- 对于浅层电路，采用Xavier初始化
- 对于深层电路，使用identity-block初始化保持初始酉性
梯度估计：
- 利用参数位移法则计算解析梯度
- 对于含噪声情况，采用随机梯度下降变体
测量优化：
- 对重要Pauli项进行分组测量
- 利用Clifford随机化降低测量次数
误差缓解：
- 采用零噪声外推技术
- 使用测量后选择减少辅助比特误差影响
实际限制与应对：
- 当前硬件限制：≈50-100个参数为实际可优化上限
- 解决方案：分层优化、参数共享、子空间方法
我在实际实验中发现，对于非马尔可夫系统，在优化后期引入小幅度的随机扰动有助于跳出局部最优。此外，将辅助比特的测量结果作为经典反馈信号调整后续电路参数，能提升约15%的收敛速度。
