1. 课程导入

量子多体动力学研究的是大量相互作用的量子自由度在时间中的演化规律。这里的“多体”不是指三五个粒子，而是指自由度数目随体系尺寸指数增长、无法用单粒子图像近似描述的强关联体系。典型对象包括：冷原子光晶格中的玻色–哈伯德模型、量子自旋链中的横场伊辛模型、凝聚态中的费米–哈伯德模型、分子体系中的电子关联，以及量子计算平台上的量子比特阵列。它们共同的特征是：哈密顿量由局域相互作用项构成，但希尔伯特空间维度随粒子数指数增长，例如 N 个自旋-1/2 的维度为 2^N，N=40 时已达约 10^12，直接存储态矢量就已超出常规内存。

因此，本课程的核心问题可以概括为一句话：在指数大的希尔伯特空间中，如何以可控的误差和可接受的计算代价，获得体系随时间的演化信息，并从中提取可观测的物理量。围绕这一目标，需要同时处理三类困难。第一是表示困难：态矢量不能全量存储，必须用张量网络、稀疏基矢、采样表示或量子线路等压缩结构来近似。第二是演化困难：时间演化算符 e^{-iHt} 不能直接对指数维矩阵做指数运算，需要借助 Trotter 分解、Krylov 子空间、时间依赖变分原理或量子线路编译来近似实现。第三是测量困难：动力学关联函数、谱函数、纠缠熵等目标量往往涉及非局域算符或两点时间关联，需要设计专门的测量与重构流程。

上图给出了本课程的整体方法链路：左侧是物理对象，中间是表示与演化，右侧是测量、分析与验证。三者不是单向流水线，而是存在反馈——测量结果会反过来指导表示截断和演化步长的选择。例如，当纠缠熵增长过快导致 MPS 截断误差超标时，需要减小时间步长或改用更高阶 Trotter 分解，甚至切换到 PEPS 或量子线路表示。

核心挑战的具体表现

以横场伊辛模型 H = -J Σ σ^z_i σ^z_{i+1} - h Σ σ^x_i 为例，取 J=1、h=1、N=20，采用 MPS 表示并做时间演化。初始态设为所有自旋沿 x 方向极化，即 |+⟩^⊗N。在 t=0 到 t=5 的演化中，可以观察到三个典型现象。第一，纠缠熵从 0 开始近似线性增长，斜率约为 0.3～0.5，这意味着 MPS 的键维数 χ 需要随时间指数增长才能维持固定精度。第二，若固定 χ=64，截断误差在 t≈2 之后开始显著上升，从 10^{-10} 量级升至 10^{-4} 量级，此时 ⟨σ^z_i(t)⟩ 的曲线会出现可见偏差。第三，Trotter 步长 δt=0.05 与 δt=0.01 的结果在 t=5 时差异约为 10^{-3}，说明二阶 Trotter 的误差累积不可忽略。这些数字是教学示例中的合理量级，实际值取决于具体参数和实现。

这些现象对应三个核心挑战：纠缠墙（χ 随 t 指数增长）、时间步长与精度权衡（δt 越小越准但步数越多）、可观测量的测量代价（两点关联函数需要多次测量或重构）。量子模拟平台还额外面临门保真度限制和测量采样噪声，例如当前超导量子比特的单门保真度约 99.9%、双门约 99.5%，演化 100 个双门后总保真度可能降至 60% 以下，因此线路编译需要尽量压缩门数。

课程路线图

本课程共 30 章，按“基础—方法—实现—分析—前沿”五段推进。第 1～5 章建立量子多体动力学的基本语言：希尔伯特空间、算符代数、时间演化算符、关联函数与谱表示。第 6～12 章系统讲解张量网络方法，包括 MPS/PEPS 的规范形式、TEBD、TDVP、DMRG 与动力学 DMRG，并给出截断误差与键维数的控制策略。第 13～18 章转向精确对角化与 Krylov 方法，覆盖稀疏矩阵存储、Lanczos 迭代、Chebyshev 展开和有限温度动力学。第 19～24 章进入量子模拟实现，包括 Trotter 线路编译、变分量子模拟、测量方案与误差缓解。第 25～30 章聚焦分析、验证与前沿，涵盖纠缠动力学、热化与局域化、谱函数重构、方法交叉验证以及经典–量子混合方案。

需要特别强调的是，本课程不追求覆盖所有方法，而是以“能复现、能判读、能验证”为标准。每学完一种方法，应当能够回答三个问题：该方法适用于什么参数区间？关键误差来源是什么？如何判断结果是否可信？只有把这三个问题落实到具体计算中，模拟结果才具有物理意义，而不只是一条好看的曲线。