2. 量子开放系统基础
封闭量子系统用态矢量 \(|\psi\rangle\) 描述，演化由薛定谔方程给出，是幺正的、可逆的。但任何真实量子比特都不可避免地与周围环境耦合：超导量子比特会通过电容耦合到控制线路和衬底声子，离子阱会受激光相位噪声和背景气体碰撞影响，光子会因光纤损耗而泄漏。这类系统称为开放量子系统，其核心特征是：我们只关心系统本身，环境自由度被"积分掉"，系统演化不再是幺正的，而是出现退相干、耗散和弛豫。
密度矩阵：从纯态到混合态的统一语言
当系统状态不完全已知，或与另一系统纠缠时，态矢量不再够用。密度矩阵定义为 \(\rho=\sum_i p_i|\psi_i\rangle\langle\psi_i|\)，其中 \(p_i\) 是系统处于 \(|\psi_i\rangle\) 的经典概率。它把"量子叠加"与"经典不确定性"统一在一个算符里。几个必须记住的性质：\(\rho\) 是厄米的、半正定的、迹为 1；纯态满足 \(\mathrm{Tr}(\rho^2)=1\)，混合态则 \(\mathrm{Tr}(\rho^2)<1\)。对单量子比特，\(\rho=\frac{1}{2}(I+\vec r\cdot\vec\sigma)\)，其中 \(\vec r=(r_x,r_y,r_z)\) 是布洛赫矢量，\(|\vec r|=1\) 对应纯态，\(|\vec r|<1\) 对应混合态，\(|\vec r|=0\) 是完全混态。
密度矩阵最重要的用途是取偏迹。若系统与环境整体处于纯态 \(|\Psi\rangle_{SE}\)，系统约化密度矩阵为 \(\rho_S=\mathrm{Tr}_E(|\Psi\rangle\langle\Psi|)\)。这一步正是"环境被忽略"的数学表达，也是退相干的来源：系统与环境的纠缠使 \(\rho_S\) 从纯态变为混合态，相干项衰减。
主方程：把环境效应压缩成演化方程
直接求解"系统+环境"的联合薛定谔方程通常不可行，因为环境自由度太多。主方程的思路是：在若干近似下，写出只含系统自由度的演化方程。最常用的近似有三条：
- 玻恩近似：系统与环境耦合弱，环境状态几乎不受系统影响，联合态近似为 \(\rho_{SE}(t)\approx\rho_S(t)\otimes\rho_E\)。
- 马尔可夫近似：环境关联时间远短于系统弛豫时间，环境"忘记"过去，演化只依赖当前 \(\rho_S(t)\)。
- 旋转波近似：忽略快速振荡项，保留能量守恒的跃迁。
在这些近似下得到Lindblad 主方程（也称 GKSL 方程）：
dρ/dt = -i/ħ [H, ρ] + Σ_k γ_k ( L_k ρ L_k† - 1/2 {L_k† L_k, ρ} )第一项是幺正演化，由系统哈密顿量 \(H\) 决定；第二项是耗散项，由一组跳跃算符 \(L_k\) 和速率 \(\gamma_k\) 描述。\(L_k\) 的物理含义是"环境通过哪种方式作用于系统"：例如 \(L=\sigma_-=\begin{pmatrix}0&1\\0&0\end{pmatrix}\) 描述从 \(|1\rangle\) 到 \(|0\rangle\) 的自发辐射或弛豫，速率 \(\gamma_1=1/T_1\)；\(L=\sigma_z\) 描述纯退相位，速率 \(\gamma_\phi=1/T_\phi\)。对超导量子比特，典型参数如 \(T_1\approx 50\,\mu s\)、\(T_2\approx 30\,\mu s\)，则 \(\gamma_1=2\times10^4\,s^{-1}\)，\(\gamma_\phi=1/T_2-1/(2T_1)\approx 2.3\times10^4\,s^{-1}\)。
Lindblad 形式的关键价值在于：它保证 \(\rho(t)\) 始终是合法密度矩阵（迹为 1、半正定），且演化是完全正定的。这使它成为量子光学、量子信息和量子热力学的标准工具。
量子轨迹：把主方程"展开"成随机历史
主方程给出的是系综平均 \(\rho(t)=\mathbb{E}[|\psi(t)\rangle\langle\psi(t)|]\)。但单次实验观测到的往往是一条具体轨迹：量子比特在某随机时刻突然从 \(|1\rangle\) 跳到 \(|0\rangle\)，发出一个光子；或者相位在连续扩散。量子轨迹方法正是把主方程"解构"为大量随机纯态轨迹的系综。
以自发辐射为例，在时间步 \(dt\) 内：
- 以概率 \(dp=\gamma_1 dt\langle\psi|\sigma_+\sigma_-|\psi\rangle\) 发生跳跃，态变为 \(|\psi\rangle\to\sigma_-|\psi\rangle/\|\sigma_-|\psi\rangle\|\)；
- 以概率 \(1-dp\) 不发生跳跃，态按非厄米有效哈密顿量 \(H_{\mathrm{eff}}=H-\frac{i}{2}\sum_k\gamma_k L_k^\dagger L_k\) 连续演化，再归一化。
对 \(N\) 条独立轨迹取平均，\(\rho(t)\approx\frac{1}{N}\sum_{n=1}^N|\psi_n(t)\rangle\langle\psi_n(t)|\)，当 \(N\to\infty\) 时收敛到主方程解。这就是量子轨迹的物理起源：环境对系统的每一次"测量"都迫使系统选择一条随机历史，而主方程只是这些历史的统计平均。这也解释了为什么量子轨迹天然适合模拟单次实验记录，例如超导量子比特的连续零差测量、离子阱的量子跳跃荧光信号。
三者的关系与实操要点
把上面三块串起来：密度矩阵是描述开放系统的状态语言；Lindblad 主方程是它的演化方程；量子轨迹是主方程的随机展开。在数值模拟中，如果只关心系综平均量（如布洛赫矢量衰减、保真度），直接积分主方程即可；如果要模拟单次测量记录、条件反馈或量子跳跃统计，则用轨迹方法更自然。
一个常见的实操陷阱是：Lindblad 方程要求跳跃算符和速率满足完全正定条件，不能随意把 \(T_2\) 写成 \(2T_1\) 之外的组合。对纯退相位，正确写法是 \(\gamma_\phi=1/T_2-1/(2T_1)\)，而不是 \(1/T_2\)。若忽略这一点，模拟出的相干衰减会偏快。另一个要点是时间步长：轨迹模拟中 \(dt\) 必须远小于 \(1/\gamma_1\) 和 \(1/\gamma_\phi\)，否则跳跃概率 \(dp\) 可能超过 1，导致非物理结果。通常取 \(dt\le 0.01/\max(\gamma_k)\) 作为经验准则。