3. 量子开放系统的基本概念

前两章我们把量子态写成密度矩阵 ρ，并看到纯态退相干后会出现非对角元衰减。但真实的量子比特从来不是孤立的：它总与周围的电磁场、声子、控制线路甚至测量仪器耦合。把“我们关心的那部分”称为系统 S，把其余一切称为环境 E，整体 S+E 仍是封闭的，服从幺正演化；而系统自身的演化一般不再是幺正的，这正是开放系统的核心特征。

系统-环境耦合：从联合幺正到约化动力学

设初始时刻系统与环境无关联，联合态为 ρ_SE(0)=ρ_S(0)⊗ρ_E(0)。整体演化由幺正算符 U(t) 给出，系统约化密度矩阵为

ρ_S(t) = Tr_E[ U(t) (ρ_S(0) ⊗ ρ_E(0)) U†(t) ]

这个式子说明：环境被“求迹”掉，但它留下的影响全部编码在映射 Λ_t 中。耦合强度、环境谱密度、温度以及系统-环境相互作用的形式（如 σ_z⊗B、σ_x⊗B）共同决定 Λ_t 的具体形状。以量子比特退相位为例，取相互作用 H_int=σ_z⊗B，环境算符 B 的涨落会让 |0⟩ 与 |1⟩ 之间积累随机相位，非对角元按 e^−Γ(t) 衰减，而布居数保持不变——这就是纯退相位（pure dephasing）的物理图像。

马尔可夫近似：无记忆环境的代价与便利

若环境关联时间 τ_E 远小于系统弛豫时间 τ_S，环境“忘记”自己过去的状态，系统演化只依赖当前 ρ_S(t)，于是可写成微分方程形式（Lindblad 主方程）：

dρ/dt = -i[H, ρ] + Σ_k ( L_k ρ L_k† - ½{L_k†L_k, ρ} )

其中 L_k 是跳跃算符。对纯退相位，取 L=√(γ_φ/2) σ_z，可得非对角元 ρ₀₁(t)=ρ₀₁(0)e^−γ_φt，布居数不变。马尔可夫近似把复杂的记忆核替换为瞬时衰减率，代价是丢失短时非马尔可夫反弹；当耦合很强或环境有窄带结构时，这一近似会失效。

动力学映射与完全正定保迹（CPTP）

更一般地，任意时刻的演化可写成算符和表示（Kraus 形式）：

ρ_S(t) = Σ_i K_i(t) ρ_S(0) K_i†(t),   Σ_i K_i†(t) K_i(t) = I

保迹条件 ΣK_i^†K_i=I 保证概率守恒；完全正定要求对任意辅助系统维数扩展后映射仍保持正定，这排除了转置这类“正但非完全正”的映射。CPTP 映射是开放系统动力学的合法数学骨架，Lindblad 演化正是其连续时间生成元。

量子比特退相位：把四个概念串成一条物理链

现在用纯退相位把上述概念落到具体数值上。设量子比特初态为 |+⟩=(|0⟩+|1⟩)/√2，密度矩阵 ρ₀₁(0)=0.5。环境为马尔可夫型，退相位率 γ_φ=0.2 ns⁻¹。则 t 时刻相干性为 ρ₀₁(t)=0.5e^−0.2t（t 以 ns 计）。取几个时间点：t=0 时 0.500；t=5 ns 时约 0.184；t=10 ns 时约 0.068；t=20 ns 时约 0.009。可见相干性以时间常数 T₂*=1/γ_φ=5 ns 指数衰减，而 ρ₀₀=ρ₁₁=0.5 始终不变。若同时存在能量弛豫（T₁），则布居数也会向热平衡值演化，此时需用完整的 Lindblad 方程处理。

从 CPTP 角度看，纯退相位的 Kraus 算符可写为 K₀=√((1+e^−γ_φt)/2) I，K₁=√((1−e^−γ_φt)/2) σ_z，直接验证 ΣK_i^†K_i=I 成立，且映射完全正定。这说明退相位不是简单的“乘以一个因子”，而是一个合法的量子通道。

小结与判断依据

本节建立了开放系统的四层逻辑：系统-环境耦合给出相互作用的物理来源；联合幺正加求迹给出精确但通常难以求解的约化动力学；马尔可夫近似在 τ_E≪τ_S 时把动力学简化为 Lindblad 主方程；CPTP 映射则从数学上刻画了所有合法量子通道。对量子比特退相位，相干性按 e^−γ_φt 衰减、布居数不变，是检验这些概念最直观的示例。实际建模时，先估计 τ_E 与 τ_S 的比值，再决定用 Lindblad 还是非马尔可夫方法，最后用 CPTP 条件校验映射的合法性。