2. 多体哈密顿量建模
多体动力学模拟的第一步，是把一个物理问题写成可计算的哈密顿量。这一步决定了后续用哪种算法、能算多大体系、能提取哪些物理量。本节围绕四种最常用的建模路径展开：二次量子化、格点模型、连续空间模型，以及自旋-玻色子耦合模型。它们并非互斥，而是同一物理问题的不同抽象层次。
一、二次量子化：从粒子标签到占据数
第一量子化把 N 个粒子写成波函数 ψ(r₁,…,r_N)，粒子数固定但坐标维度随 N 增长。二次量子化改用产生算符 a†ᵢ 和湮灭算符 aᵢ 描述单粒子态 i 的占据情况，多体态写成 Fock 态 |n₁n₂…⟩。其核心代数关系为：
[aᵢ, aⱼ†] = δᵢⱼ （玻色子）
{aᵢ, aⱼ†} = δᵢⱼ （费米子）费米子用反对易关系，自动保证交换两个粒子波函数反号，从而满足泡利不相容原理。构建哈密顿量的通用形式是：
H = Σᵢⱼ tᵢⱼ aᵢ† aⱼ + ½ Σᵢⱼₖₗ Vᵢⱼₖₗ aᵢ† aⱼ† aₖ aₗ其中 tᵢⱼ 是单粒子项（动能加外势），Vᵢⱼₖₗ 是两体相互作用矩阵元。实际建模时，先选定单粒子基（如平面波、谐振子基或格点瓦尼尔基），再计算这些矩阵元。若基组截断到 M 个态，Fock 空间维度对费米子是 C(M,N)，对玻色子则无上限，需要额外截断最大占据数。
二、格点模型：把空间离散成可计算的图
格点模型把连续空间替换为离散节点，粒子只能位于格点或格点间的键上。最典型的是紧束缚模型：
H = -t Σ⟨i,j⟩,σ (cᵢσ† cⱼσ + h.c.) + U Σᵢ nᵢ↑ nᵢ↓ - μ Σᵢ nᵢ这里 t 是最近邻跃迁振幅，U 是同一格点上的相互作用能，μ 是化学势。当 U/t 很大时，系统趋向莫特绝缘体；U/t 很小时，趋向金属。格点模型的关键参数包括：格点几何（链、方格子、三角格子、六角格子）、边界条件（开放或周期）、填充数 ν=N/L，以及是否包含次近邻跃迁 t′。
构建格点哈密顿量时，通常先生成格点坐标与邻接表，再按邻接关系填充稀疏矩阵。以 1D 链为例，若 L=12、N=6、U/t=4，可用精确对角化或密度矩阵重整化群（DMRG）求解。判断结果是否合理，可检查基态能量随 L 的收敛趋势、单粒子谱函数在费米面附近的权重，以及电荷能隙 Δc=E(N+1)+E(N−1)−2E(N) 是否随 U 增大而打开。
三、连续空间模型：保留真实坐标与长程相互作用
当体系具有连续平移对称性、长程库仑相互作用或需要精确描述散射长度时，连续空间模型更合适。其哈密顿量写作：
H = Σσ ∫ dr ψσ†(r) [−ℏ²∇²/2m + Vext(r)] ψσ(r)
+ ½ Σσσ′ ∫∫ dr dr′ ψσ†(r) ψσ′†(r′) V(r−r′) ψσ′(r′) ψσ(r)实际计算中，连续空间通常用平面波基或实空间网格离散化。平面波基的截断由动能上限 Ecut 控制，对应最大波矢 kcut=√(2mEcut)/ℏ。对冷原子体系，接触相互作用可用赝势 V(r)=gδ(r) 近似，g=4πℏ²a/m，a 为 s 波散射长度。对电子体系，则保留 1/(4πε₀|r−r′|) 的库仑形式。
连续空间建模的常见陷阱是紫外发散和有限尺寸效应。判断离散化是否足够，可检查总能量随 Ecut 或网格间距 Δx 的收敛曲线；若能量变化小于目标精度（例如 1 meV/粒子），则认为截断合理。此外，长程相互作用在周期边界下需要 Ewald 求和，否则会引入虚假的有限尺寸偏差。
四、自旋-玻色子耦合模型：开放量子系统的标准范式
自旋-玻色子模型描述一个量子两能级系统（自旋）与大量玻色模（声子、光子或环境模）的耦合，是研究退相干、耗散和量子相变的经典框架。其哈密顿量为：
H = Δ σx /2 + ε σz /2 + Σₖ ωₖ bₖ† bₖ + σz Σₖ gₖ (bₖ† + bₖ)其中 Δ 是隧穿振幅，ε 是偏置，ωₖ 是玻色模频率，gₖ 是耦合强度。环境的影响完全由谱密度 J(ω)=Σₖ gₖ² δ(ω−ωₖ) 刻画。常用的 Ohmic 谱密度为 J(ω)=2παω exp(−ω/ωc)，α 为耗散强度，ωc 为截断频率。当 α<1/2 时系统处于相干区，α>1/2 时进入局域化相，这是量子相变的典型判据。
构建该模型时，需要先确定谱密度的离散化方案。常见做法是把连续 J(ω) 采样为有限个玻色模，例如在 [0,ωc] 上取 M 个频点，每个频点的耦合强度由 gₖ²=J(ωₖ)Δω 给出。M 的选取需保证低温下热占据数 n(ωₖ)=1/(e^{ℏωₖ/kBT}−1) 被充分覆盖。若 M 太小，低频模缺失会导致退相干时间被高估；若 M 太大，希尔伯特空间维度指数增长，需要改用层级运动方程（HEOM）或矩阵乘积态方法。
五、四条路径的衔接与选择
实际建模时，四条路径常常组合使用。例如，研究冷原子在光晶格中的超流-莫特相变，可先用连续空间模型写出真实哈密顿量，再通过瓦尼尔函数局域化导出格点模型，最后用二次量子化形式在 Fock 空间中求解。研究超导量子比特的退相干，则从自旋-玻色子模型出发，把环境谱密度离散化为有限玻色模，再用 HEOM 或 MPS 计算约化密度矩阵的演化。
选择路径时，可参考以下判断：若体系具有强关联和短程相互作用，格点模型通常最经济；若需要精确描述长程库仑力或散射性质，连续空间模型更合适；若关注开放系统动力学和退相干，自旋-玻色子模型是标准起点；而二次量子化则是所有路径的共同语言，负责把物理图像翻译成算符代数。
总结来说，多体哈密顿量建模的核心不是写出一个公式，而是明确“基是什么、截断在哪里、参数如何对应物理量、结果如何判断收敛”。把这四点落实到具体模型上，后续的动力学模拟才有可靠起点。