1. 课程导入与密度矩阵基础
在量子力学入门阶段，我们习惯用态矢量 \(|\psi\rangle\) 描述系统状态。态矢量简洁、直观，能完美刻画一个“确定知道”的量子态——我们称之为纯态。但真实量子系统几乎从不孤立：它与环境耦合、经历退相干、存在制备误差，或者我们只掌握系统的部分信息。此时，态矢量不再够用，必须引入更一般的描述工具——密度矩阵（density matrix），记作 \(\rho\)。
本课程的目标，是建立一套从密度矩阵出发、能够模拟开放量子系统的完整方法。第 1 章先解决最基础的问题：密度矩阵是什么、它如何统一纯态与混合态、如何用布洛赫球直观理解它，以及整门课的知识地图长什么样。
为什么态矢量不够用
设想一个场景：某台设备制备量子比特时，有 70% 的概率输出 \(|0\rangle\)，30% 的概率输出 \(|1\rangle\)。我们并不知道这一次拿到的是哪一个，只知道这个概率分布。若强行写成一个态矢量，比如 \(\alpha|0\rangle+\beta|1\rangle\)，那描述的是 \(|0\rangle\) 与 \(|1\rangle\) 的相干叠加，会出现干涉项；而上述情形是经典概率混合，两者物理上完全不同。混合态的测量统计无法用单个态矢量复现，这正是需要密度矩阵的原因。
密度矩阵的定义正是为处理这种“系综”而生。若系统以概率 \(p_i\) 处于纯态 \(|\psi_i\rangle\)，则
ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|, 其中 p_i ≥ 0, Σ_i p_i = 1纯态是它的特例：只有一个 \(p_i=1\)，此时 \(\rho=|\psi\rangle\langle\psi|\)。注意 \(\rho\) 是算符而非矢量，它把“量子叠加”与“经典不确定性”统一在同一个数学对象里。
密度矩阵的三条基本性质
任何合法的密度矩阵必须满足以下条件，这也是模拟中检查数值结果是否合理的第一道关卡：
- 厄米性：\(\rho^\dagger=\rho\)，保证测量期望值为实数。
- 单位迹：\(\mathrm{Tr}(\rho)=1\)，对应总概率归一。
- 半正定：\(\rho\ge 0\)，即所有本征值非负，保证任意测量的概率非负。
由此可定义纯度 \(\gamma=\mathrm{Tr}(\rho^2)\)。纯态满足 \(\gamma=1\)；混合态满足 \(\gamma<1\)。对单量子比特，\(\gamma=1\) 对应布洛赫矢量长度 1，\(\gamma<1\) 对应长度小于 1。纯度是判断“这个态有多混合”的最直接指标。
布洛赫球：单量子比特的几何图像
对单量子比特，密度矩阵可写成
ρ = (I + r·σ) / 2, r = (⟨σx⟩, ⟨σy⟩, ⟨σz⟩)其中 \(\sigma=(\sigma_x,\sigma_y,\sigma_z)\) 是泡利矩阵，\(r\) 称为布洛赫矢量。它的几何意义非常清晰：
- \(|r|=1\)：球面上的点，纯态；
- \(|r|<1\)：球内部的点，混合态；
- \(r=0\)：球心，完全混合态 \(I/2\)，纯度最低。
布洛赫球的三个坐标轴分别对应 \(\langle\sigma_x\rangle,\langle\sigma_y\rangle,\langle\sigma_z\rangle\)。例如 \(|0\rangle\) 位于 \(+z\) 极点，\(|1\rangle\) 位于 \(-z\) 极点，\(|+\rangle=(|0\rangle+|1\rangle)/\sqrt2\) 位于 \(+x\) 极点。退相干过程在几何上表现为布洛赫矢量向球心收缩——这正是开放系统模拟中最常观察到的现象之一。
开放系统模拟的整体学习地图
有了密度矩阵这个统一语言，后续课程将沿着一条清晰的路线展开。可以把它理解为四个层次：
- 描述层：密度矩阵、约化密度矩阵、纯度与纠缠熵，回答“系统现在是什么状态”。
- 演化层：封闭系统用幺正演化 \(\rho(t)=U\rho(0)U^\dagger\)；开放系统用主方程，如 Lindblad 形式 \(\dot\rho=-i[H,\rho]+\sum_k(L_k\rho L_k^\dagger-\tfrac12\{L_k^\dagger L_k,\rho\})\)。
- 模拟层：把连续演化离散化，常用方法包括量子轨迹（quantum jump）、蒙特卡洛波函数、矩阵乘积算符（MPO）与时间演化块消减（TEBD）等。
- 验证层：检查迹守恒、厄米性、半正定性，比较解析解或小系统精确解，观察纯度与纠缠熵的物理趋势。
本课程会反复用到几个关键参数：耦合强度 \(\gamma\)、环境温度对应的热占据数 \(n_{th}\)、时间步长 \(dt\)、截断维数 \(\chi\)。它们分别控制退相干速率、热噪声强度、数值精度和计算代价。理解这些参数如何影响布洛赫矢量长度和主方程解，是后续动手模拟的核心。