3. 量子多体动力学的基本方程

量子多体动力学模拟的第一步，不是选算法，而是确定“用哪个方程描述系统”。四个基本方程并非彼此替代，而是对应四种不同的物理情境：态矢是否已知、系统是否封闭、是否需要统计混合、以及我们关心的是态还是算符。选错方程，后续的数值方法、边界条件和误差分析都会失去意义。

一、四个方程的适用边界

薛定谔方程描述封闭纯态：iℏ ∂_t|ψ(t)⟩ = Ĥ|ψ(t)⟩。它要求系统与环境完全隔离，且初态是单一态矢。对 N 个自旋或 N 个格点，态矢维度随 N 指数增长，因此它只适合小规模精确模拟或作为基准。

冯·诺依曼方程把描述对象从态矢换成密度算符：iℏ ∂_tρ = [Ĥ, ρ]。它仍描述封闭系统，但允许初态是统计混合（如有限温度热态 ρ ∝ e^−βĤ），也便于计算任意可观测量 ⟨Â⟩ = Tr(ρÂ)。代价是 ρ 的矩阵规模为态矢维度的平方。

林德布拉德方程在冯·诺依曼方程上加入耗散项，描述开放系统：
∂_tρ = −(i/ℏ)[Ĥ, ρ] + Σ_k γ_k(L_kρL_k^† − ½{L_k^†L_k, ρ})。
其中 L_k 是跳跃算符，γ_k 是对应速率。它适用于系统与马尔可夫环境耦合、需要描述退相干与耗散的情形，例如超导量子比特的弛豫、冷原子中的粒子损失。

海森堡运动方程把时间演化放到算符上：dÂ/dt = (i/ℏ)[Ĥ, Â] + ∂_tÂ。态矢固定，算符随时间变化。它适合关心少数几个可观测量（如自旋关联、序参量）而不需要完整态的情形，也是推导关联函数和响应函数的自然起点。

二、从方程到可执行模拟：一个两自旋示例

考虑两个耦合自旋-1/2，哈密顿量取 Ĥ = J(σ_x⁽¹⁾σ_x⁽²⁾ + σ_y⁽¹⁾σ_y⁽²⁾) + B(σ_z⁽¹⁾ + σ_z⁽²⁾)，取 J = 1.0、B = 0.5（示例参数，单位任意）。初态设为 |↑↓⟩。若用薛定谔方程，态矢维度为 4，可直接对 4×4 矩阵做指数演化；若用冯·诺依曼方程，ρ 为 4×4 矩阵，演化后可直接读出 ⟨σ_z⁽¹⁾⟩ = Tr(ρ σ_z⁽¹⁾)。

若给第一个自旋加上振幅阻尼通道，跳跃算符 L = σ₋⁽¹⁾，速率 γ = 0.1，则必须用林德布拉德方程。此时 ρ 的对角元不再守恒，⟨σ_z⁽¹⁾⟩ 会从初值逐渐衰减到 −1 附近，而非在 ±1 之间振荡。这个差异是判断“是否需要林德布拉德”的直接依据：如果实验上观察到关联函数单调衰减而非持续振荡，就说明环境耗散不可忽略。

海森堡方程则适合只追踪 ⟨σ_z⁽¹⁾(t)σ_z⁽²⁾(t)⟩ 这类关联。对上述 Ĥ，可以写出 dσ_x⁽¹⁾/dt = −2B σ_y⁽¹⁾ − 2J σ_y⁽¹⁾σ_z⁽²⁾ 等一组耦合方程。若只关心少数算符，这组方程规模远小于 4×4 密度矩阵，但会引入更高阶关联，需要截断近似。

三、数值实现中的关键选择

实际编码时，四个方程对应不同的数据结构与积分器。薛定谔方程用态矢数组，配合四阶龙格-库塔或 Chebyshev 多项式展开；冯·诺依曼与林德布拉德用密度矩阵，林德布拉德还需把耗散项写成 Lindblad 超算符 ℒ，再对 ℒ 做指数或分段积分。海森堡方程则用算符的矩阵表示，或直接对易关系做符号推导后数值积分。

四、方程选择对照表

四个方程在数学上通过取迹、部分迹、对偶变换相互联系：对冯·诺依曼方程取 Tr(·) 即得海森堡方程；对林德布拉德方程令 γ_k = 0 即退化为冯·诺依曼方程；对纯态 ρ = |ψ⟩⟨ψ| 代入冯·诺依曼方程即回到薛定谔方程。理解这层关系，才能在模拟中根据物理情境灵活切换，而不是把某一个方程当作唯一正确的出发点。