1.1

Problem 1.1

srednickiChapter 1

习题 1.1

来源: 第1章, PDF第29页


1.1 Show that the Dirac matrices must be even dimensional. Hint: show that the eigenvalues of β\beta are all ±1\pm 1, and that Tr β=0\text{Tr } \beta = 0. To show that Tr β=0\text{Tr } \beta = 0, consider, e.g., Tr α12β\text{Tr } \alpha_1^2 \beta. Similarly, show that Tr αi=0\text{Tr } \alpha_i = 0.

习题 1.1 - 解答


物理背景与分析

在构建相对论性量子力学波动方程(狄拉克方程)时,为了保证方程对时间和空间坐标是一阶导数,且满足相对论协变性与能量-动量关系 E2=p2c2+m2c4E^2 = \mathbf{p}^2 c^2 + m^2 c^4,狄拉克引入了四个矩阵 α1,α2,α3\alpha_1, \alpha_2, \alpha_3β\beta。这些矩阵必须满足以下克利福德代数(Clifford Algebra)的反对易关系: {αi,αj}=αiαj+αjαi=2δijI(i,j=1,2,3)\{\alpha_i, \alpha_j\} = \alpha_i \alpha_j + \alpha_j \alpha_i = 2\delta_{ij} I \quad (i, j = 1, 2, 3) {αi,β}=αiβ+βαi=0\{\alpha_i, \beta\} = \alpha_i \beta + \beta \alpha_i = 0 β2=I\beta^2 = I 其中 II 是单位矩阵。本题旨在通过矩阵的迹(Trace)和本征值性质,从上述代数关系中严格证明狄拉克矩阵的维数 NN 必须是偶数。

推导过程

第一步:证明 β\beta 的本征值只能是 ±1\pm 1

λ\lambda 为矩阵 β\beta 的任意一个本征值,对应的非零本征矢为 v|v\rangle,即: βv=λv\beta |v\rangle = \lambda |v\rangle 在等式两边再次作用 β\beta 矩阵: β2v=β(λv)=λ(βv)=λ2v\beta^2 |v\rangle = \beta (\lambda |v\rangle) = \lambda (\beta |v\rangle) = \lambda^2 |v\rangle 根据狄拉克矩阵的定义,β2=I\beta^2 = I,因此: Iv=λ2v    λ2=1I |v\rangle = \lambda^2 |v\rangle \implies \lambda^2 = 1 解得 β\beta 的本征值为: λ=±1\lambda = \pm 1

第二步:证明 Tr(β)=0\text{Tr}(\beta) = 0

利用矩阵迹的循环性质 Tr(ABC)=Tr(CAB)\text{Tr}(ABC) = \text{Tr}(CAB) 以及 α12=I\alpha_1^2 = I,我们可以将 β\beta 的迹改写为: Tr(β)=Tr(Iβ)=Tr(α12β)=Tr(α1α1β)\text{Tr}(\beta) = \text{Tr}(I \beta) = \text{Tr}(\alpha_1^2 \beta) = \text{Tr}(\alpha_1 \alpha_1 \beta) 应用迹的循环性质,将最左侧的 α1\alpha_1 移至最右侧: Tr(α1α1β)=Tr(α1βα1)\text{Tr}(\alpha_1 \alpha_1 \beta) = \text{Tr}(\alpha_1 \beta \alpha_1) 根据反对易关系 {α1,β}=0\{\alpha_1, \beta\} = 0,有 α1β=βα1\alpha_1 \beta = -\beta \alpha_1。将其代入上式: Tr(α1βα1)=Tr(βα1α1)=Tr(βα12)\text{Tr}(\alpha_1 \beta \alpha_1) = \text{Tr}(-\beta \alpha_1 \alpha_1) = -\text{Tr}(\beta \alpha_1^2) 再次利用 α12=I\alpha_1^2 = ITr(βα12)=Tr(βI)=Tr(β)-\text{Tr}(\beta \alpha_1^2) = -\text{Tr}(\beta I) = -\text{Tr}(\beta) 由此得到等式 Tr(β)=Tr(β)\text{Tr}(\beta) = -\text{Tr}(\beta),移项可得 2Tr(β)=02\text{Tr}(\beta) = 0,即: Tr(β)=0\text{Tr}(\beta) = 0

第三步:证明 Tr(αi)=0\text{Tr}(\alpha_i) = 0

同理,利用 β2=I\beta^2 = I 和迹的循环性质,可以证明任意 αi\alpha_i 的迹也为零: Tr(αi)=Tr(β2αi)=Tr(ββαi)=Tr(βαiβ)\text{Tr}(\alpha_i) = \text{Tr}(\beta^2 \alpha_i) = \text{Tr}(\beta \beta \alpha_i) = \text{Tr}(\beta \alpha_i \beta) 利用反对易关系 βαi=αiβ\beta \alpha_i = -\alpha_i \betaTr(βαiβ)=Tr(αiββ)=Tr(αiβ2)=Tr(αiI)=Tr(αi)\text{Tr}(\beta \alpha_i \beta) = \text{Tr}(-\alpha_i \beta \beta) = -\text{Tr}(\alpha_i \beta^2) = -\text{Tr}(\alpha_i I) = -\text{Tr}(\alpha_i) 因此 2Tr(αi)=02\text{Tr}(\alpha_i) = 0,即: Tr(αi)=0\text{Tr}(\alpha_i) = 0

第四步:证明矩阵维数必须为偶数

在线性代数中,矩阵的迹等于其所有本征值之和。 设狄拉克矩阵的维数为 NN。由第一步可知,β\beta 的本征值只能是 +1+11-1。 假设在 β\betaNN 个本征值中,有 n+n_++1+1nn_-1-1(其中 n+,nn_+, n_- 均为非负整数)。 矩阵的总维数 NN 为本征值的个数之和: N=n++nN = n_+ + n_- 矩阵 β\beta 的迹为本征值之和: Tr(β)=n+(+1)+n(1)=n+n\text{Tr}(\beta) = n_+(+1) + n_-(-1) = n_+ - n_- 由第二步已证 Tr(β)=0\text{Tr}(\beta) = 0,因此: n+n=0    n+=nn_+ - n_- = 0 \implies n_+ = n_- 将此结果代入维数表达式中: N=n++n+=2n+N = n_+ + n_+ = 2n_+ 因为 n+n_+ 是一个非负整数,所以 NN 必定是一个偶数。

N 必须为偶数\boxed{N \text{ 必须为偶数}}

1.2

Problem 1.2

srednickiChapter 1

习题 1.2

来源: 第1章, PDF第29页


1.2 With the hamiltonian of eq. (1.32), show that the state defined in eq. (1.33) obeys the abstract Schrödinger equation, eq. (1.1), if and only if the wave function obeys eq. (1.30). Your demonstration should apply both to the case of bosons, where the particle creation and annihilation operators obey the commutation relations of eq. (1.31), and to fermions, where the particle creation and annihilation operators obey the anticommutation relations of eq. (1.38).


Referenced Equations:

Equation (1.1):

itψ,t=Hψ,t,(1.1)i\hbar \frac{\partial}{\partial t}|\psi, t\rangle = H|\psi, t\rangle , \tag{1.1}

Equation (1.30):

itψ=[j=1n(22mj2+U(xj))+j=1nk=1j1V(xjxk)]ψ,(1.30)i\hbar \frac{\partial}{\partial t} \psi = \left[ \sum_{j=1}^n \left( -\frac{\hbar^2}{2m} \nabla_j^2 + U(\mathbf{x}_j) \right) + \sum_{j=1}^n \sum_{k=1}^{j-1} V(\mathbf{x}_j - \mathbf{x}_k) \right] \psi , \tag{1.30}

Equation (1.31):

[a(x),a(x)]=0,[a(x),a(x)]=0,[a(x),a(x)]=δ3(xx),(1.31)\begin{aligned} [a(\mathbf{x}), a(\mathbf{x}')] &= 0 , \\ [a^{\dagger}(\mathbf{x}), a^{\dagger}(\mathbf{x}')] &= 0 , \\ [a(\mathbf{x}), a^{\dagger}(\mathbf{x}')] &= \delta^3(\mathbf{x} - \mathbf{x}') , \end{aligned} \tag{1.31}

Equation (1.32):

H=d3x a(x)(22m2+U(x))a(x)+12d3x d3y V(xy)a(x)a(y)a(y)a(x).(1.32)\begin{aligned} H = &\int d^3x \ a^{\dagger}(\mathbf{x}) \left( -\frac{\hbar^2}{2m} \nabla^2 + U(\mathbf{x}) \right) a(\mathbf{x}) \\ &+ \frac{1}{2} \int d^3x \ d^3y \ V(\mathbf{x} - \mathbf{y}) a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) . \end{aligned} \tag{1.32}

Equation (1.33):

ψ,t=d3x1d3xn ψ(x1,,xn;t)a(x1)a(xn)0,(1.33)|\psi, t\rangle = \int d^3x_1 \dots d^3x_n \ \psi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) a^{\dagger}(\mathbf{x}_1) \dots a^{\dagger}(\mathbf{x}_n) |0\rangle , \tag{1.33}

Equation (1.38):

{a(x),a(x)}=0,{a(x),a(x)}=0,{a(x),a(x)}=δ3(xx),(1.38)\begin{aligned} \{a(\mathbf{x}), a(\mathbf{x}')\} &= 0 , \\ \{a^\dagger(\mathbf{x}), a^\dagger(\mathbf{x}')\} &= 0 , \\ \{a(\mathbf{x}), a^\dagger(\mathbf{x}')\} &= \delta^3(\mathbf{x} - \mathbf{x}') , \end{aligned} \tag{1.38}

习题 1.2 - 解答


为了统一处理玻色子和费米子,我们引入符号 η\eta:对于玻色子 η=1\eta = 1,对于费米子 η=1\eta = -1。根据题意中的对易关系 (1.31) 和反对易关系 (1.38),产生与湮灭算符满足统一的代数关系: a(x)a(x)ηa(x)a(x)=δ3(xx)a(\mathbf{x}) a^{\dagger}(\mathbf{x}') - \eta a^{\dagger}(\mathbf{x}') a(\mathbf{x}) = \delta^3(\mathbf{x} - \mathbf{x}') a(x)a(x)=ηa(x)a(x),a(x)a(x)=ηa(x)a(x)a(\mathbf{x}) a(\mathbf{x}') = \eta a(\mathbf{x}') a(\mathbf{x}), \quad a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{x}') = \eta a^{\dagger}(\mathbf{x}') a^{\dagger}(\mathbf{x})

定义 nn 粒子基矢为 x1,,xn=a(x1)a(xn)0|\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = a^{\dagger}(\mathbf{x}_1) \dots a^{\dagger}(\mathbf{x}_n) |0\rangle。 利用上述关系,湮灭算符 a(x)a(\mathbf{x}) 作用在基矢上时,需要依次与 a(xj)a^{\dagger}(\mathbf{x}_j) 交换直到作用于真空态 0|0\rangle 得到 00。每次交换跨过一个算符会产生因子 η\eta,因此有: a(x)x1,,xn=j=1nηj1δ3(xxj)x1,,xj^,,xna(\mathbf{x}) |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = \sum_{j=1}^n \eta^{j-1} \delta^3(\mathbf{x} - \mathbf{x}_j) |\mathbf{x}_1, \dots, \widehat{\mathbf{x}_j}, \dots, \mathbf{x}_n\rangle 其中 xj^\widehat{\mathbf{x}_j} 表示该位置的算符被移除。

我们将哈密顿量 (1.32) 分为单体部分 H1H_1 和两体部分 H2H_2。记单体微分算符为 h(x)=22mx2+U(x)h(\mathbf{x}) = -\frac{\hbar^2}{2m} \nabla_{\mathbf{x}}^2 + U(\mathbf{x})

1. 单体哈密顿量 H1H_1 的作用

H1=d3x a(x)h(x)a(x)H_1 = \int d^3x \ a^{\dagger}(\mathbf{x}) h(\mathbf{x}) a(\mathbf{x}) 作用于态 ψ,t|\psi, t\rangleH1ψ,t=d3x1d3xn ψ(x1,,xn;t)j=1nd3x a(x)[h(x)δ3(xxj)]ηj1x1,,xj^,,xnH_1 |\psi, t\rangle = \int d^3x_1 \dots d^3x_n \ \psi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) \sum_{j=1}^n \int d^3x \ a^{\dagger}(\mathbf{x}) \left[ h(\mathbf{x}) \delta^3(\mathbf{x} - \mathbf{x}_j) \right] \eta^{j-1} |\mathbf{x}_1, \dots, \widehat{\mathbf{x}_j}, \dots, \mathbf{x}_n\rangle 注意到 a(x)ηj1x1,,xj^,,xn=x1,,xj1,x,xj+1,,xna^{\dagger}(\mathbf{x}) \eta^{j-1} |\mathbf{x}_1, \dots, \widehat{\mathbf{x}_j}, \dots, \mathbf{x}_n\rangle = |\mathbf{x}_1, \dots, \mathbf{x}_{j-1}, \mathbf{x}, \mathbf{x}_{j+1}, \dots, \mathbf{x}_n\rangle。 利用 δ\delta 函数的性质,微分算符作用满足 h(x)δ3(xxj)=h(xj)δ3(xjx)h(\mathbf{x}) \delta^3(\mathbf{x} - \mathbf{x}_j) = h(\mathbf{x}_j) \delta^3(\mathbf{x}_j - \mathbf{x})。代入后,我们在第 jj 项的积分中交换积分变量名 xxj\mathbf{x} \leftrightarrow \mathbf{x}_jH1ψ,t=j=1nd3x1d3xn d3x ψ(x1,,x,,xn;t)[h(xj)δ3(xjx)]x1,,xj,,xnH_1 |\psi, t\rangle = \sum_{j=1}^n \int d^3x_1 \dots d^3x_n \ d^3x \ \psi(\mathbf{x}_1, \dots, \mathbf{x}, \dots, \mathbf{x}_n; t) \left[ h(\mathbf{x}_j) \delta^3(\mathbf{x}_j - \mathbf{x}) \right] |\mathbf{x}_1, \dots, \mathbf{x}_j, \dots, \mathbf{x}_n\rangle 由于 h(xj)h(\mathbf{x}_j) 仅包含对 xj\mathbf{x}_j 的导数,它可以从对 x\mathbf{x} 的积分中提出。完成对 x\mathbf{x}δ\delta 函数积分后,得到: H1ψ,t=d3x1d3xn [j=1nh(xj)ψ(x1,,xn;t)]x1,,xnH_1 |\psi, t\rangle = \int d^3x_1 \dots d^3x_n \ \left[ \sum_{j=1}^n h(\mathbf{x}_j) \psi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) \right] |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle

2. 两体哈密顿量 H2H_2 的作用

首先计算两个湮灭算符的作用: a(y)a(x)x1,,xn=j=1nηj1δ3(xxj)a(y)x1,,xj^,,xna(\mathbf{y}) a(\mathbf{x}) |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = \sum_{j=1}^n \eta^{j-1} \delta^3(\mathbf{x} - \mathbf{x}_j) a(\mathbf{y}) |\mathbf{x}_1, \dots, \widehat{\mathbf{x}_j}, \dots, \mathbf{x}_n\rangle =j=1nηj1δ3(xxj)(k=1j1ηk1δ3(yxk),xk^,,xj^,+k=j+1nηk2δ3(yxk),xj^,,xk^,)= \sum_{j=1}^n \eta^{j-1} \delta^3(\mathbf{x} - \mathbf{x}_j) \left( \sum_{k=1}^{j-1} \eta^{k-1} \delta^3(\mathbf{y} - \mathbf{x}_k) |\dots, \widehat{\mathbf{x}_k}, \dots, \widehat{\mathbf{x}_j}, \dots\rangle + \sum_{k=j+1}^n \eta^{k-2} \delta^3(\mathbf{y} - \mathbf{x}_k) |\dots, \widehat{\mathbf{x}_j}, \dots, \widehat{\mathbf{x}_k}, \dots\rangle \right)H2=12d3xd3y V(xy)a(x)a(y)a(y)a(x)H_2 = \frac{1}{2} \int d^3x d^3y \ V(\mathbf{x} - \mathbf{y}) a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) 作用于基矢,对 x,y\mathbf{x}, \mathbf{y} 的积分会将它们分别替换为 xj,xk\mathbf{x}_j, \mathbf{x}_kH2x1,,xn=12j=1n(k=1j1ηj+k2V(xjxk)a(xj)a(xk),xk^,,xj^,H_2 |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = \frac{1}{2} \sum_{j=1}^n \left( \sum_{k=1}^{j-1} \eta^{j+k-2} V(\mathbf{x}_j - \mathbf{x}_k) a^{\dagger}(\mathbf{x}_j) a^{\dagger}(\mathbf{x}_k) |\dots, \widehat{\mathbf{x}_k}, \dots, \widehat{\mathbf{x}_j}, \dots\rangle \right. +k=j+1nηj+k3V(xjxk)a(xj)a(xk),xj^,,xk^,)\left. + \sum_{k=j+1}^n \eta^{j+k-3} V(\mathbf{x}_j - \mathbf{x}_k) a^{\dagger}(\mathbf{x}_j) a^{\dagger}(\mathbf{x}_k) |\dots, \widehat{\mathbf{x}_j}, \dots, \widehat{\mathbf{x}_k}, \dots\rangle \right) 接下来将产生算符移回原位以恢复基矢 x1,,xn|\mathbf{x}_1, \dots, \mathbf{x}_n\rangle

  • 对于 k<jk < j 的项:a(xk)a^{\dagger}(\mathbf{x}_k) 移回第 kk 位需跨过 k1k-1 个算符(产生 ηk1\eta^{k-1}),随后 a(xj)a^{\dagger}(\mathbf{x}_j) 移回第 jj 位需跨过 j1j-1 个算符(产生 ηj1\eta^{j-1})。总符号为 ηj+k2\eta^{j+k-2}。与原系数相乘得 η2j+2k4=1\eta^{2j+2k-4} = 1
  • 对于 k>jk > j 的项:a(xk)a^{\dagger}(\mathbf{x}_k) 移回第 kk 位需跨过 k2k-2 个算符(因为 xj\mathbf{x}_j 缺失,产生 ηk2\eta^{k-2}),随后 a(xj)a^{\dagger}(\mathbf{x}_j) 移回第 jj 位需跨过 j1j-1 个算符(产生 ηj1\eta^{j-1})。总符号为 ηj+k3\eta^{j+k-3}。与原系数相乘得 η2j+2k6=1\eta^{2j+2k-6} = 1

由于 V(xjxk)=V(xkxj)V(\mathbf{x}_j - \mathbf{x}_k) = V(\mathbf{x}_k - \mathbf{x}_j),上述两项求和完全相同,合并后得到: H2x1,,xn=j=1nk=1j1V(xjxk)x1,,xnH_2 |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = \sum_{j=1}^n \sum_{k=1}^{j-1} V(\mathbf{x}_j - \mathbf{x}_k) |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle 因此,将 H2H_2 作用于 ψ,t|\psi, t\rangle 得到: H2ψ,t=d3x1d3xn [j=1nk=1j1V(xjxk)ψ(x1,,xn;t)]x1,,xnH_2 |\psi, t\rangle = \int d^3x_1 \dots d^3x_n \ \left[ \sum_{j=1}^n \sum_{k=1}^{j-1} V(\mathbf{x}_j - \mathbf{x}_k) \psi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) \right] |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle

3. 薛定谔方程的等价性

将抽象薛定谔方程 (1.1) 的左边展开: itψ,t=d3x1d3xn (itψ(x1,,xn;t))x1,,xni\hbar \frac{\partial}{\partial t} |\psi, t\rangle = \int d^3x_1 \dots d^3x_n \ \left( i\hbar \frac{\partial}{\partial t} \psi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) \right) |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle 结合前面求得的 Hψ,t=(H1+H2)ψ,tH |\psi, t\rangle = (H_1 + H_2) |\psi, t\rangle,方程 itψ,tHψ,t=0i\hbar \frac{\partial}{\partial t} |\psi, t\rangle - H |\psi, t\rangle = 0 可写为: d3x1d3xn Φ(x1,,xn;t)x1,,xn=0\int d^3x_1 \dots d^3x_n \ \Phi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) |\mathbf{x}_1, \dots, \mathbf{x}_n\rangle = 0 其中 Φ=iψt[j=1n(22mj2+U(xj))+j=1nk=1j1V(xjxk)]ψ\Phi = i\hbar \frac{\partial \psi}{\partial t} - \left[ \sum_{j=1}^n \left( -\frac{\hbar^2}{2m} \nabla_j^2 + U(\mathbf{x}_j) \right) + \sum_{j=1}^n \sum_{k=1}^{j-1} V(\mathbf{x}_j - \mathbf{x}_k) \right] \psi

由于波函数 ψ\psi 对于玻色子是全对称的,对于费米子是全反对称的,且方括号内的 nn 体哈密顿算符在粒子交换下是对称的,因此 Φ\Phi 具有与 ψ\psi 相同的交换对称性。基矢 x1,,xn|\mathbf{x}_1, \dots, \mathbf{x}_n\rangle 同样具有该对称性,积分非零当且仅当 Φ\Phi 本身非零。因此,上述积分恒为零的充要条件是 Φ(x1,,xn;t)0\Phi(\mathbf{x}_1, \dots, \mathbf{x}_n; t) \equiv 0

这直接给出了波函数必须满足的方程:

itψ=[j=1n(22mj2+U(xj))+j=1nk=1j1V(xjxk)]ψ\boxed{ i\hbar \frac{\partial}{\partial t} \psi = \left[ \sum_{j=1}^n \left( -\frac{\hbar^2}{2m} \nabla_j^2 + U(\mathbf{x}_j) \right) + \sum_{j=1}^n \sum_{k=1}^{j-1} V(\mathbf{x}_j - \mathbf{x}_k) \right] \psi }

此即为方程 (1.30),证明完毕。

1.3

Problem 1.3

srednickiChapter 1

习题 1.3

来源: 第1章, PDF第29页


1.3 Show explicitly that [N,H]=0[N, H] = 0, where HH is given by eq. (1.32) and NN by eq. (1.35).


Referenced Equations:

Equation (1.32):

H=d3x a(x)(22m2+U(x))a(x)+12d3x d3y V(xy)a(x)a(y)a(y)a(x).(1.32)\begin{aligned} H = &\int d^3x \ a^{\dagger}(\mathbf{x}) \left( -\frac{\hbar^2}{2m} \nabla^2 + U(\mathbf{x}) \right) a(\mathbf{x}) \\ &+ \frac{1}{2} \int d^3x \ d^3y \ V(\mathbf{x} - \mathbf{y}) a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) . \end{aligned} \tag{1.32}

Equation (1.35):

N=d3x a(x)a(x)(1.35)N = \int d^3x \ a^{\dagger}(\mathbf{x}) a(\mathbf{x}) \tag{1.35}

习题 1.3 - 解答


物理背景与分析

本题要求证明总粒子数算符 NN 与系统的哈密顿量 HH 对易,即 [N,H]=0[N, H] = 0。在二次量子化表象中,这对应于系统总粒子数守恒的物理事实。哈密顿量 HH 包含单体相互作用项(动能与外势)和两体相互作用项,每一项中产生算符 aa^\dagger 与湮灭算符 aa 的数目总是相等的。这反映了系统具有全局 U(1)U(1) 规范对称性,从而导致粒子数守恒。

为了保证证明的普适性,以下推导将同时适用于玻色子和费米子。我们引入统一定义的交换/反对易括号 [A,B]±=AB±BA[A, B]_\pm = AB \pm BA(取 ++ 对应费米子的反对易,取 - 对应玻色子的对易)。场算符满足以下基本关系:

[a(x),a(y)]±=δ(3)(xy)[a(\mathbf{x}), a^\dagger(\mathbf{y})]_\pm = \delta^{(3)}(\mathbf{x} - \mathbf{y})
[a(x),a(y)]±=0,[a(x),a(y)]±=0[a(\mathbf{x}), a(\mathbf{y})]_\pm = 0, \quad [a^\dagger(\mathbf{x}), a^\dagger(\mathbf{y})]_\pm = 0

1. 基本对易关系推导

首先计算总粒子数算符 N=d3z a(z)a(z)N = \int d^3z \ a^{\dagger}(\mathbf{z}) a(\mathbf{z}) 与单个场算符的对易子。利用算符恒等式 [AB,C]=A[B,C]±[A,C]±B[AB, C] = A[B, C]_\pm \mp [A, C]_\pm B

对于湮灭算符 a(x)a(\mathbf{x})

[a(z)a(z),a(x)]=a(z)[a(z),a(x)]±[a(z),a(x)]±a(z)=a(z)(0)(δ(3)(zx))a(z)=δ(3)(zx)a(z)\begin{aligned} [a^{\dagger}(\mathbf{z}) a(\mathbf{z}), a(\mathbf{x})] &= a^{\dagger}(\mathbf{z}) [a(\mathbf{z}), a(\mathbf{x})]_\pm \mp [a^{\dagger}(\mathbf{z}), a(\mathbf{x})]_\pm a(\mathbf{z}) \\ &= a^{\dagger}(\mathbf{z}) (0) \mp \left( \mp \delta^{(3)}(\mathbf{z} - \mathbf{x}) \right) a(\mathbf{z}) \\ &= -\delta^{(3)}(\mathbf{z} - \mathbf{x}) a(\mathbf{z}) \end{aligned}

z\mathbf{z} 积分得到:

[N,a(x)]=d3z(δ(3)(zx)a(z))=a(x)[N, a(\mathbf{x})] = \int d^3z \left( -\delta^{(3)}(\mathbf{z} - \mathbf{x}) a(\mathbf{z}) \right) = -a(\mathbf{x})

对于产生算符 a(x)a^\dagger(\mathbf{x})

[a(z)a(z),a(x)]=a(z)[a(z),a(x)]±[a(z),a(x)]±a(z)=a(z)δ(3)(zx)0=δ(3)(zx)a(z)\begin{aligned} [a^{\dagger}(\mathbf{z}) a(\mathbf{z}), a^{\dagger}(\mathbf{x})] &= a^{\dagger}(\mathbf{z}) [a(\mathbf{z}), a^{\dagger}(\mathbf{x})]_\pm \mp [a^{\dagger}(\mathbf{z}), a^{\dagger}(\mathbf{x})]_\pm a(\mathbf{z}) \\ &= a^{\dagger}(\mathbf{z}) \delta^{(3)}(\mathbf{z} - \mathbf{x}) \mp 0 \\ &= \delta^{(3)}(\mathbf{z} - \mathbf{x}) a^{\dagger}(\mathbf{z}) \end{aligned}

z\mathbf{z} 积分得到:

[N,a(x)]=d3z(δ(3)(zx)a(z))=a(x)[N, a^{\dagger}(\mathbf{x})] = \int d^3z \left( \delta^{(3)}(\mathbf{z} - \mathbf{x}) a^{\dagger}(\mathbf{z}) \right) = a^{\dagger}(\mathbf{x})

这两个关系式表明,产生算符使系统粒子数加1,湮灭算符使系统粒子数减1。

2. 与单体哈密顿量 H1H_1 的对易

将哈密顿量分为单体部分 H1H_1 和两体部分 H2H_2。单体部分为:

H1=d3x a(x)H0(x)a(x)H_1 = \int d^3x \ a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) a(\mathbf{x})

其中 H0(x)=22m2+U(x)\mathcal{H}_0(\mathbf{x}) = -\frac{\hbar^2}{2m} \nabla^2 + U(\mathbf{x}) 是作用在空间坐标上的微分算符,它与 NN 对易。利用对易子求导法则 [A,BC]=[A,B]C+B[A,C][A, BC] = [A, B]C + B[A, C]

[N,a(x)H0(x)a(x)]=[N,a(x)]H0(x)a(x)+a(x)H0(x)[N,a(x)]=(a(x))H0(x)a(x)+a(x)H0(x)(a(x))=a(x)H0(x)a(x)a(x)H0(x)a(x)=0\begin{aligned} [N, a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) a(\mathbf{x})] &= [N, a^{\dagger}(\mathbf{x})] \mathcal{H}_0(\mathbf{x}) a(\mathbf{x}) + a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) [N, a(\mathbf{x})] \\ &= \left( a^{\dagger}(\mathbf{x}) \right) \mathcal{H}_0(\mathbf{x}) a(\mathbf{x}) + a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) \left( -a(\mathbf{x}) \right) \\ &= a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) a(\mathbf{x}) - a^{\dagger}(\mathbf{x}) \mathcal{H}_0(\mathbf{x}) a(\mathbf{x}) \\ &= 0 \end{aligned}

因此,对空间积分后有:

[N,H1]=0[N, H_1] = 0

3. 与两体哈密顿量 H2H_2 的对易

两体相互作用部分为:

H2=12d3x d3y V(xy)a(x)a(y)a(y)a(x)H_2 = \frac{1}{2} \int d^3x \ d^3y \ V(\mathbf{x} - \mathbf{y}) a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x})

同样利用对易子的展开法则 [N,ABCD]=[N,A]BCD+A[N,B]CD+AB[N,C]D+ABC[N,D][N, ABCD] = [N,A]BCD + A[N,B]CD + AB[N,C]D + ABC[N,D],计算被积函数中算符乘积与 NN 的对易子:

[N,a(x)a(y)a(y)a(x)]=[N,a(x)]a(y)a(y)a(x)+a(x)[N,a(y)]a(y)a(x)+a(x)a(y)[N,a(y)]a(x)+a(x)a(y)a(y)[N,a(x)]=a(x)a(y)a(y)a(x)+a(x)a(y)a(y)a(x)a(x)a(y)a(y)a(x)a(x)a(y)a(y)a(x)=(1+111)a(x)a(y)a(y)a(x)=0\begin{aligned} &[N, a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x})] \\ &= [N, a^{\dagger}(\mathbf{x})] a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) + a^{\dagger}(\mathbf{x}) [N, a^{\dagger}(\mathbf{y})] a(\mathbf{y}) a(\mathbf{x}) \\ &\quad + a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) [N, a(\mathbf{y})] a(\mathbf{x}) + a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) [N, a(\mathbf{x})] \\ &= a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) + a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) \\ &\quad - a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) - a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) \\ &= (1 + 1 - 1 - 1) a^{\dagger}(\mathbf{x}) a^{\dagger}(\mathbf{y}) a(\mathbf{y}) a(\mathbf{x}) \\ &= 0 \end{aligned}

由于被积算符与 NN 的对易子严格为零,积分后必然有:

[N,H2]=0[N, H_2] = 0

4. 结论

综合单体与两体部分的计算结果,总哈密顿量 H=H1+H2H = H_1 + H_2 与总粒子数算符 NN 的对易关系为:

[N,H]=[N,H1]+[N,H2]=0+0=0[N, H] = [N, H_1] + [N, H_2] = 0 + 0 = 0

证明完毕。

[N,H]=0\boxed{[N, H] = 0}