2. 二次量子化与格点表示

在量子格点模型中，我们关心的不是单个粒子的波函数，而是大量全同粒子在离散格点上的集体量子态。二次量子化提供了描述这类系统的标准语言：把“粒子”变成算符作用在态上的激发，把格点上的跃迁、相互作用写成产生湮灭算符的组合。本节围绕四个核心问题展开：产生湮灭算符的定义与代数、福克空间的构造、哈密顿量在格点上的表达，以及玻色子与费米子算符代数的区别。

一、产生湮灭算符与占据数表示

设格点系统有 N 个单粒子态（例如一维链的 N 个格点，或每个格点上的轨道）。对每个单粒子态 i，定义一对算符：湮灭算符 a_i 和产生算符 a_i^†。它们作用在福克态上，改变第 i 个态的粒子数。

对玻色子，算符满足对易关系：

[a_i, a_j†] = δ_ij
[a_i, a_j] = [a_i†, a_j†] = 0

对费米子，算符满足反对易关系：

{a_i, a_j†} = δ_ij
{a_i, a_j} = {a_i†, a_j†} = 0

这里的 δ_ij 是克罗内克符号。对易与反对易的差别，直接导致两种统计：玻色子可以多个占据同一态，费米子则受泡利不相容原理限制，每个态占据数只能是 0 或 1。

算符作用在福克态上的规则可以统一写成：

a_i† |n_0, n_1, ..., n_i, ...> 
  = √(n_i + 1) |..., n_i + 1, ...>   (玻色子)
  = (-1)^(Σ_{j<i} n_j) |..., n_i + 1, ...>  (费米子，n_i = 0)

a_i |n_0, n_1, ..., n_i, ...>
  = √(n_i) |..., n_i - 1, ...>   (玻色子)
  = (-1)^(Σ_{j<i} n_j) |..., n_i - 1, ...>  (费米子，n_i = 1)

费米子表达式中的符号因子来自反对易关系，它保证交换两个费米子算符时态矢出现负号。这个符号在格点模型中非常关键，尤其在处理跃迁项和超导配对项时，必须与格点排序约定保持一致。

二、福克空间的结构

福克空间由所有可能的占据数组态张成。对玻色子系统，每个格点的占据数 n_i 可以取 0, 1, 2, …；对费米子系统，n_i 只能取 0 或 1。因此，一个 N 格点的费米子系统，其福克空间维数是 2^N；玻色子系统的维数则是无穷维的，实际计算中需要截断到某个最大占据数。

在格点模型中，我们通常选择一组基矢，例如按格点顺序排列的占据数基：

|n_0 n_1 ... n_{N-1}>

所有物理算符都可以表示为这组基上的矩阵。对于费米子系统，由于维数随 N 指数增长，精确对角化通常只能处理十几个格点；更大系统需要借助密度矩阵重整化群、量子蒙特卡洛或张量网络等方法。

三、哈密顿量在格点上的表达

格点模型中最常见的哈密顿量由两部分组成：单粒子跃迁项和相互作用项。以紧束缚模型为例，一维链上的跃迁项写作：

H_hop = -t Σ_{<i,j>} (a_i† a_j + a_j† a_i)

其中 t 是跃迁振幅，<i,j> 表示最近邻格点对。这一项描述粒子从一个格点跳到相邻格点。若格点之间有相位或规范场，跃迁项会带上相位因子 e^iθ_ij。

相互作用项以 Hubbard 模型为例：

H_int = U Σ_i n_{i↑} n_{i↓}

其中 n_iσ = a_iσ^†a_iσ 是格点 i 上自旋 σ 的粒子数算符，U 是同一格点上两个自旋相反电子的库仑排斥能。完整的 Hubbard 哈密顿量为：

H = -t Σ_{<i,j>,σ} (a_{iσ}† a_{jσ} + h.c.) + U Σ_i n_{i↑} n_{i↓} - μ Σ_{i,σ} n_{iσ}

其中 μ 是化学势，用于控制总粒子数。在格点表示下，哈密顿量被写成算符的二次型加四次型组合，所有物理量都可以通过算符代数计算。

四、不同统计的算符代数与格点排序

玻色子与费米子的算符代数差异，在格点模型中会直接影响数值实现。对玻色子，算符可以按任意顺序排列，对易关系简单；对费米子，必须固定一个格点排序，并在跨越多个格点的算符乘积中引入符号因子。

例如，在费米子系统中，算符乘积 a_i^†a_j 当 i ≠ j 时，若 i < j，则不需要额外符号；若 i > j，则需要交换算符顺序并引入负号。这一规则在编写程序时通常通过 Jordan-Wigner 变换或直接构造费米子符号表来实现。

下面这张图概括了从单粒子态到格点哈密顿量的构造关系，以及玻色子与费米子算符代数的核心区别。

五、从算符到可计算形式

要在计算机上处理格点模型，需要把算符代数转化为矩阵或位串操作。对费米子系统，常用 Jordan-Wigner 变换把费米子算符映射为自旋算符：

a_i = (Π_{j<i} σ_j^z) σ_i^-
a_i† = (Π_{j<i} σ_j^z) σ_i^+

这样，费米子符号被编码进 σ^z 串中，格点排序自然体现在乘积顺序里。对玻色子系统，则可以直接用截断的占据数基构造矩阵，每个格点的最大占据数取 n_max，总维数为 (n_max+1)^N。

实际计算中，判断算符代数是否正确，可以检查几个基本关系：对易/反对易关系是否满足、粒子数算符是否守恒、哈密顿量是否厄米。例如，对 Hubbard 模型，总粒子数算符 N = Σ_iσ n_iσ 与 H 对易，因此总粒子数是好量子数，可以在固定粒子数扇区中求解，从而大幅降低矩阵维数。

总结来说，二次量子化把格点模型写成产生湮灭算符的多项式，福克空间提供了所有可能占据数组态的线性空间，哈密顿量则通过跃迁项和相互作用项在格点上展开。玻色子与费米子的区别集中体现在对易与反对易代数上，而费米子的符号规则是数值实现中最容易出错、也最需要严格检查的部分。