2. 量子模拟的数学与物理基础
量子模拟器输出的从来不是“物理量本身”，而是一组按量子力学规则组织起来的数学对象。要把这些对象翻译成可验证的物理结论，必须先厘清四个层次：态矢量描述纯态、密度矩阵描述一般态、可观测量期望值把态与测量结果连接起来、哈密顿量演化决定态如何随时间变化。本节围绕这四层展开，并给出模拟结果到真实物理量的映射路径。
一、态矢量：纯态的完备描述
一个 n 量子比特系统的纯态可写成计算基的线性叠加：
|ψ⟩ = Σ_i c_i |i⟩, i = 0,1,...,2^n − 1
归一化条件：Σ_i |c_i|² = 1
测量概率：P(i) = |c_i|²
相位：c_i = |c_i| e^{iφ_i}，全局相位不可观测，相对相位可观测态矢量是复希尔伯特空间中的单位向量，其维度随比特数指数增长：10 比特对应 1024 个复振幅，30 比特对应约 10⁹ 个。这正是经典计算机难以直接存储全态矢量的原因，也是量子模拟必须借助张量网络、稀疏表示或采样方法的原因。在模拟中，态矢量通常以复数数组保存，实部与虚部分开存储；读取时需注意归一化误差，一般要求 |Σ|c_i|² − 1| < 10⁻¹⁰。
二、密度矩阵：混合态与退相干的载体
当系统与环境耦合、或我们只掌握部分信息时，纯态描述不再充分，需用密度矩阵 ρ：
ρ = Σ_k p_k |ψ_k⟩⟨ψ_k|, Σ_k p_k = 1, p_k ≥ 0
纯态判据：Tr(ρ²) = 1；混合态：Tr(ρ²) < 1
约化密度矩阵：ρ_A = Tr_B(ρ_AB)
期望值：⟨O⟩ = Tr(ρ O)密度矩阵是厄米、半正定、迹为 1 的算符。它的对角元给出各计算基的布居概率，非对角元（相干项）反映量子叠加是否保留。退相干过程正是非对角元衰减的过程。在模拟中，若采用含噪声模型，密度矩阵维度为 4ⁿ，比态矢量更昂贵，因此常用蒙特卡洛轨迹或量子信道表示来近似。
三、可观测量期望值：从态到测量结果
可观测量由厄米算符 O 表示，其期望值 ⟨O⟩ = ⟨ψ|O|ψ⟩（纯态）或 Tr(ρO)（混合态）。模拟中最常测量的是泡利算符组合，例如：
Z 期望值：⟨Z⟩ = P(0) − P(1)
X 期望值：⟨X⟩ = 2·Re(⟨0|ρ|1⟩)
关联：⟨Z_i Z_j⟩ = P(00)+P(11) − P(01) − P(10)实际模拟中，期望值通过多次采样估计：若测量 N 次得到 m 次“0”，则 ⟨Z⟩ ≈ (2m − N)/N，统计误差约为 1/√N。因此 1000 次采样对应约 3% 的统计不确定度。判断模拟结果是否可信，首先要区分统计误差与模型误差：前者随采样数下降，后者来自哈密顿量近似、截断或噪声建模。
四、哈密顿量演化：时间如何改变态
封闭量子系统的演化由薛定谔方程决定：
iℏ ∂|ψ⟩/∂t = H|ψ⟩
形式解：|ψ(t)⟩ = e^{−iHt/ℏ} |ψ(0)⟩
若 H = Σ_k E_k |k⟩⟨k|，则 |ψ(t)⟩ = Σ_k c_k e^{−iE_k t/ℏ} |k⟩模拟中通常取 ℏ = 1，时间以能量倒数计。演化算符 U(t) = e^{−iHt} 是酉算符，保持态矢量归一化。对于多比特系统，直接指数化不可行，需用 Trotter 分解：
e^{−i(H_A + H_B)t} ≈ (e^{−iH_A t/n} e^{−iH_B t/n})^n + O(t²/n)Trotter 步数 n 越大，误差越小，但线路深度增加。这是模拟精度与资源之间的核心权衡。若 H 含非对易项，Trotter 误差随 t²/n 下降；若项对易，则分解精确。
五、模拟结果到真实物理量的映射
模拟输出的是量子态或测量统计，而真实物理量往往是能量、磁化强度、关联函数或谱函数。映射关系如下：
- 能量：E = ⟨ψ|H|ψ⟩，对基态可用变分量子本征求解器（VQE）估计，误差来源包括 ansatz 表达能力和测量噪声。
- 磁化强度：M = (1/N) Σ_i ⟨Z_i⟩，直接由泡利 Z 期望值平均得到。
- 关联函数：C(r) = ⟨Z_i Z_{i+r}⟩ − ⟨Z_i⟩⟨Z_{i+r}⟩，用于判断相变与关联长度。
- 谱函数：A(ω) 需通过实时演化后做傅里叶变换，模拟中常以时间步长 Δt 采样，频率分辨率受总演化时间限制。
映射过程中必须保持参照约定一致：例如自旋向上对应 |0⟩ 还是 |1⟩、哈密顿量符号约定、时间演化方向。若模拟采用 |0⟩ 为自旋向下，而实验数据以自旋向上为正，则磁化强度符号会整体翻转。这类约定错误在验证阶段极易被误判为物理效应。
六、验证时的关键检查点
拿到模拟结果后，建议按以下顺序检查，而不是直接与实验值比较：
- 归一化与厄米性：态矢量模长是否为 1，密度矩阵是否满足 Tr(ρ)=1、ρ=ρ†、本征值非负。
- 极限行为：当相互作用强度 J→0 时，磁化强度是否回到自由自旋值；当演化时间 t→0 时，态是否回到初态。
- 对称性：若哈密顿量具有平移对称性，关联函数是否只依赖距离 r 而非绝对位置。
- 收敛性：增加 Trotter 步数 n 或采样数 N，结果是否收敛到稳定值；若未收敛，误差棒是否覆盖变化范围。
- 符号一致性：能量是否为负（对吸引相互作用）、磁化方向是否与约定一致。