3. 开放量子系统框架

封闭量子系统的演化由薛定谔方程描述,是幺正的、可逆的。但任何真实量子器件——超导量子比特、离子阱、光子线路——都不可能完全孤立,它始终与周围的电磁场、声子、控制线路、测量仪器耦合。退相干的物理起源,正藏在这种“系统 + 环境”的耦合结构之中。本节建立开放量子系统的基本框架:把总系统写成系统与环境的直积,写出总哈密顿量,用幺正演化描述整体,再通过对环境求偏迹得到系统的约化密度矩阵,从而看清相干性是如何“泄漏”到环境中的。

一、总哈密顿量的三段式结构

设系统 S 的希尔伯特空间为 ℋS,环境 E 的为 ℋE,总空间为 ℋS ⊗ ℋE。总哈密顿量通常写成三部分之和:

H_total = H_S ⊗ I_E  +  I_S ⊗ H_E  +  H_int

H_S   : 系统自身哈密顿量(如 qubit 的 (ħω/2)σ_z)
H_E   : 环境哈密顿量(如热库、传输线模的集合)
H_int : 系统-环境耦合项,一般可展开为
        H_int = Σ_k  S_k ⊗ B_k
        S_k 为系统算符,B_k 为环境算符(如位移、相位噪声算符)

这个结构的关键在于 H_int 的“算符张量积”形式:系统算符 Sk 与 Bath 算符 Bk 相乘,意味着环境对系统的作用不是简单的标量扰动,而是通过特定系统算符通道施加影响。例如纯退相位耦合常取 H_int = σ_z ⊗ B,此时系统能量本征态不被翻转,但相位会随环境算符 B 的涨落而随机累积;而能量弛豫(T₁)则来自 σ_x ⊗ B 或 σ_- ⊗ B 型耦合,允许能级跃迁并交换能量。因此,耦合算符的代数形式直接决定了退相干是“丢相位”还是“丢能量”,这是后续建模的出发点。

二、整体幺正演化与约化密度矩阵

总系统是封闭的,其密度矩阵 ρSE(t) 满足冯·诺依曼方程,形式解为幺正演化:

ρ_SE(t) = U(t) ρ_SE(0) U†(t),   U(t) = exp(-i H_total t / ħ)

系统约化密度矩阵:
ρ_S(t) = Tr_E[ ρ_SE(t) ] = Tr_E[ U(t) (ρ_S(0) ⊗ ρ_E(0)) U†(t) ]

这里通常假设初始时刻系统与环境无关联,即 ρSE(0) = ρS(0) ⊗ ρE(0),环境处于某个热平衡态 ρE(0)。约化密度矩阵 ρS(t) 才是我们在实验上能通过层析测量到的对象。注意:ρS(t) 的演化一般不再是幺正的,它可能变成混合态,纯度 Tr(ρS²) 下降,这正是退相干的数学表现。物理上,信息没有消失,而是通过纠缠扩散到了系统与环境之间的关联中;对 E 求偏迹,等于我们主动“看不见”这部分信息。

判断退相干是否发生,最直接的指标是 ρS 的非对角元模长。以单 qubit 为例,若 ρS = [[ρ₀₀, ρ₀₁],[ρ₁₀, ρ₁₁]],纯退相位使 |ρ₀₁| 按 exp(-t/T₂*) 衰减,而布居数 ρ₀₀、ρ₁₁ 保持不变;能量弛豫则同时改变布居数与非对角元,非对角元按 exp(-t/T₂) 衰减,且 T₂ ≤ 2T₁。

三、从总框架到主方程:一个可操作的推导路径

直接对 U(t) 求偏迹通常无法解析计算,工程上采用两条等效路径。第一条是算符和表示(Kraus 表示):把环境基矢 {|en⟩} 插入,得到

ρ_S(t) = Σ_{m,n} K_{mn}(t) ρ_S(0) K_{mn}†(t)
K_{mn} = ⟨e_m| U(t) |e_n⟩   (作用在系统空间上的算符)

Kraus 算符满足 Σ K†K ≤ I,等号成立时演化为保迹的完全正映射(量子信道)。第二条路径是在马尔可夫近似(环境关联时间远短于系统演化时间)下,对时间做二阶微扰并取粗粒化,得到 Lindblad 主方程:

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

L_k : 跳跃算符(如 σ_- 对应能量弛豫,σ_z 对应纯退相位)
γ_k : 对应速率,与耦合谱密度在系统频率处的取值相关

以纯退相位为例,取 L = σz、速率 γφ,可解得非对角元 ρ₀₁(t) = ρ₀₁(0) exp(-γφt),即 T₂* = 1/γφ。若同时存在能量弛豫速率 γ₁,则总退相干速率满足 1/T₂ = 1/(2T₁) + 1/T₂*,这一关系是区分两类噪声、反推耦合参数的重要依据。

四、框架的物理含义与建模要点

下图把本节的核心结构串起来:左侧是系统与环境的张量积结构及总哈密顿量三段式,中间是整体幺正演化,右侧是求偏迹后得到的约化密度矩阵及其非幺正表现,底部标出两条等效计算路径与可观测判据。

开放量子系统:从总哈密顿量到约化密度矩阵 总空间 ℋ_S ⊗ ℋ_E 系统 S H_S 环境 E H_E ⊗ H_int = Σ_k S_k ⊗ B_k σ_z⊗B → 纯退相位 σ_-⊗B → 能量弛豫 U(t) 整体幺正演化 ρ_SE(t) = U ρ_SE(0) U† U = exp(−i H_total t / ħ) 初始:ρ_S(0) ⊗ ρ_E(0) 演化中:S 与 E 产生纠缠 信息流向环境关联 Tr_E 约化密度矩阵 ρ_S(t) ρ_S(t) = Tr_E[ ρ_SE(t) ] 一般非幺正、可能变混合态 退相干判据 |ρ₀₁| 衰减:纯退相位 T₂* 布居数变化:能量弛豫 T₁ 1/T₂ = 1/(2T₁) + 1/T₂* 路径一:Kraus 算符和表示 ρ_S(t) = Σ K_{mn} ρ_S(0) K_{mn}† K_{mn} = ⟨e_m| U(t) |e_n⟩ 完全正、保迹映射(量子信道) 适用于任意强度耦合与短时演化 需已知环境基或可数值求解 U(t) 输出:ρ_S(t) 的完整时间序列 路径二:Lindblad 主方程 dρ/dt = −i[H_S,ρ]/ħ + Σ γ_k 𝒟[L_k]ρ 𝒟[L]ρ = LρL† − ½{L†L, ρ} 马尔可夫近似:环境关联时间 ≪ 系统时间 L = σ_-(T₁)、L = σ_z(T₂*) 速率 γ_k 由谱密度在系统频率处取值决定 输出:可直接拟合 T₁、T₂ 的指数衰减

建模时有三点需要特别注意。第一,初始关联假设:ρSE(0) = ρS(0) ⊗ ρE(0) 是推导主方程的前提,若系统制备后立即与强关联环境耦合(如低温下 qubit 与滤波器的近场耦合),该假设失效,需引入关联修正。第二,马尔可夫近似的适用边界:当环境具有长记忆(如 1/f 噪声、高 Q 谐振腔)时,Lindblad 方程会低估相干性的部分恢复(非马尔可夫回流),此时应改用记忆核主方程或层级运动方程(HEOM)。第三,参数的可辨识性:T₁、T₂、T₂* 并非独立,实验上通过拉比振荡、Ramsey 干涉、自旋回波分别提取,再与耦合谱密度模型对照,才能反推出 Sk、Bk 的具体形式与强度。

把这三层结构记牢——总哈密顿量决定耦合通道,幺正演化保证信息守恒,偏迹操作制造了“表观不可逆”——就抓住了退相干物理起源的主线:相干性并未被破坏,而是从系统自由度转移到了系统与环境的纠缠关联中,约化描述让我们只看到衰减的痕迹。