图1

# 相容与不可相容观测量

# 相容观测量

\quad 算符可以当作是一个观测的仪器,它测量出来的量就是它的本征值,这一点再量子力学中就以及学过。若算符A^\hat{A}B^\hat{B} 满足对易关系[A^,B^]=0[\hat{A},\hat{B}]=0,则这两个算符称为是相容的 (compatible),则称它俩的观测值 (本征值) 是相容可观测量。

定理:假若A^\hat{A}B^\hat{B} 是相容的,而且A^\hat{A} 的本征值非简并,则矩阵B^\hat{B}(在算符A^\hat{A} 表象下所构成的矩阵)是全对角的。

证明:若右矢i\ket{i}A^\hat{A} 的本征态矢

A^i=aii(0.1.1)\hat{A} \ket{i} = a_i \ket{i} \tag{0.1.1}

由于A^\hat{A}B^\hat{B} 相容,满足对易关系,因此有

A^(B^i)=B^A^i=ai(B^i)(0.1.2)\hat{A} (\hat{B} \ket{i}) = \hat{B} \hat{A} \ket{i} = a_i (\hat{B} \ket{i}) \tag{0.1.2}

因此态矢i\ket{i}B^i\hat{B}\ket{i} 都是算符A^\hat{A} 对应本征值为aia_i 的本征态矢。由由于A^\hat{A} 是非简并的,因此i\ket{i}B^i\hat{B}\ket{i} 就是同一个态矢,因此它们只能相差一个常数倍:

B^i=bii(0.1.3)\hat{B}\ket{i}=b_i\ket{i} \tag{0.1.3}

因此,i\ket{i} 即是算符A^\hat{A} 对应本征值为aia_i 的本征态,也是算符B^\hat{B} 对应本征值为bib_i 的本征态。这就说明了很重要的一点:当两算符是对易时,它们拥有共同本征矢量完全集。 由此也不难证明矩阵B^\hat{B}(在算符A^\hat{A} 表象下所构成的矩阵)是全对角的:

jjBi^=biδi,j(0.1.4)\sum_j \braket{j|\hat{B|i}} = b_i \delta_{i,j} \tag{0.1.4}

证毕。

上述的定理同样可以推广到A^\hat{A} 可以是简并的情况,详情参考喀兴林的《高等量子力学》。而上面画红线的这句话,也可以推广到多个算符,即若干个算符如果都相互对易的话,那么它们有一组共同的本征矢量完全集。

\quad 以及,如果两个算符是相容 ,那么他俩的测量是互不干涉的,是可以同时测量并测准的。

图2

# 不可相容观测量

\quadA^\hat{A}B^\hat{B} 是不可相容的,也就是说它们不是对易的,那么它俩的测量是会互相干涉的。我们考虑两种测量过程,假若A^\hat{A}B^\hat{B} 是不对易的,且B^\hat{B}C^\hat{C} 是不对易的:

  1. 测量仪器A^\hat{A} 测量得到al\ket{a_l},那么继续做B^\hat{B} 的测量,再接着做C^\hat{C} 的测量,最终测量结果是的概率是cn\ket{c_n} 的概率为

mcnbm2bmal2=mcnbmbmalalbmbmcn(0.1.5)\sum_m |\braket{c_n|b_m}|^2 |\braket{b_m|a_l}|^2 = \sum_m \braket{c_n|b_m} \braket{b_m|a_l} \braket{a_l|b_m} \braket{b_m|c_n} \tag{0.1.5}

图3

  1. 测量仪器A^\hat{A} 测量得到al\ket{a_l},之后直接做C^\hat{C} 的测量,最终测量结果是的概率是cn\ket{c_n} 的概率为

cnal2=mcnbmbmal2=mmcnbmbmalalbmbmcn(0.1.6)\begin{aligned} |\braket{c_n|a_l}|^2 &= \mid \sum_m \braket{c_n|b_m} \braket{b_m|a_l} \mid^2 \\&= \sum_m \sum_{m'} \braket{c_n|b_m} \braket{b_m|a_l} \braket{a_l|b_{m'}} \braket{b_{m'}|c_n} \end{aligned} \tag{0.1.6}

图4

对于上述这两种测量过程,我们最终得到相同测量结构的概率是不同的

mcnbmbmalalbmbmcnmmcnbmbmalalbmbmcn\sum_m \braket{c_n|b_m} \braket{b_m|a_l} \braket{a_l|b_m} \braket{b_m|c_n} \ne \sum_m \sum_{m'} \braket{c_n|b_m} \braket{b_m|a_l} \braket{a_l|b_{m'}} \braket{b_{m'}|c_n}

这种情况就是因为算符B^\hat{B} 的与其他两个算符不可相容导致的。那么在什么样的条件下,这两种过程最终表达式是一样的呢?在不存在简并时,只要[A^,B^]=0[\hat{A},\hat{B}]=0[B^,C^]=0[\hat{B},\hat{C}]=0 两种结果就会相等啦!(那么存在简并的情况呢?(⊙_⊙)?)

# 纯系综、混合系综与密度算符

# 纯系综与混合系综

\quad 在斯特恩 - 盖拉赫实验中,考虑通过SGSG 过滤器的一束原子系综。束流中的每个原子其自旋指向相同的方向,即由该过滤器的非均匀磁场确定的方向 (zz 轴方向)。我们取过滤出来自旋方向是zz 轴正向的原子束流系综,这种系综,我们可以用态矢

α=sz,+(0.2.1)\ket{\alpha} = \ket{s_z,+} \tag{0.2.1}

来描述。因为在这个系综里面随机选一个原子出来,其自旋方向沿zz 轴正向是确定的。当然,上式我们也可以写成在xx 方向自旋表象的形式:

α=12sx,++12sx,(0.2.2)\ket{\alpha} = \frac{1}{\sqrt{2}} \ket{s_x,+} + \frac{1}{\sqrt{2}} \ket{s_x,-} \tag{0.2.2}

因为我们可以证明sz,+=12sx,++12sx,\ket{s_z,+}=\frac{1}{\sqrt{2}} \ket{s_x,+} + \frac{1}{\sqrt{2}} \ket{s_x,-}。这种能用态矢描述的系综称为 纯系综
\quad 而如果我们做这样的考虑,让通过SGSG 过滤器得到的sz,+\ket{s_z,+} 原子束和sz,\ket{s_z,-} 原子束混合。那么混合的原子束系综可以简单地用如下的叠加态矢表示吗?

mix=c+sz,++csz,(0.2.3)\ket{mix} = c_+ \ket{s_z,+} + c_- \ket{s_z,-} \tag{0.2.3}

答案是:不可以的!因为你在这个系综里面取一个原子,它的系综方向要么向上,要么向下,不可能存在向上、向下叠加态的原子。因为系综内的各个原子的自旋方向早已因通过SGSG 过滤器决定了,只是你在这个系综里面究竟取到那个原子有概率分布。像这种系综称为 混合系综,它可以通过引入权重来描述[1]

mix={sz,+:w1sz,:w2(0.2.4)\ket{mix} = \begin{cases} \ket{s_z,+} : w_1 \\ \ket{s_z,-} : w_2 \end{cases} \tag{0.2.4}

举个例子,假若有两套SGSG 设备。一套设备加热原子的热炉子只有 “30%30\% 的产生原子速率”,其原子束通过zz 轴方向的过滤器,取其sz,+\ket{s_z,+} 的原子束;另一套设备的炉子有 “70%70\% 的产生原子速率”,同样通过zz 轴方向的过滤器,但取其sz,\ket{s_z,-} 的原子束。将俩原子束混合,那么这个混合系综可表示为:

mix={sz,+:30%sz,:70%(0.2.5)\ket{mix} = \begin{cases} \ket{s_z,+} : 30 \% \\ \ket{s_z,-} : 70\% \end{cases} \tag{0.2.5}

或许此时就明白了,混合系综其实就是纯系综的混合。

再举个例子,假若有两套SGSG。一套 “加热率为30%30\%”,通过zz 轴方向的过滤器,取其sz,+\ket{s_z,+} 的原子束;另一套 “加热率为70%70\%”,通过xx 轴方向的过滤器,取其sx,+\ket{s_x,+} 的原子束。将俩原子束混合,那么这个混合系综可表示为:

mix={sz,+:30%sx,+:70%={sz,+:30%12sz,++12sz,:70%(0.2.6)\ket{mix} = \begin{cases} \ket{s_z,+} : 30\% \\ \ket{s_x,+} : 70\% \end{cases} = \begin{cases} \ket{s_z,+} : 30\% \\ \frac{1}{\sqrt{2}} \ket{s_z,+} + \frac{1}{\sqrt{2}} \ket{s_z,-} : 70\% \end{cases} \tag{0.2.6}

其中上式利用到了sx,+=12sz,++12sz,\ket{s_x,+}=\frac{1}{\sqrt{2}} \ket{s_z,+} + \frac{1}{\sqrt{2}} \ket{s_z,-}。通过这个例子,我们可以明白,混合态里面的纯态甚至可以不正交!

\quad 对于任何一个混合系综,比如以下这个多个纯系综混合的混合系综

mix={α(1):w1α(2):w2α(n):wn(0.2.7)\ket{mix} = \begin{cases} \ket{\alpha^{(1)}} : w_1 \\ \ket{\alpha^{(2)}} : w_2 \\ \cdots \\ \ket{\alpha^{(n)}} : w_n \end{cases} \tag{0.2.7}

其权重参数,或称为布居数wiw_i 应该要满足归一化条件:

i=1nwi=1(0.2.8)\sum_{i=1}^n w_i = 1 \tag{0.2.8}

# 系综的平均与密度算符

\quad 假定在一个混合系综上测量某可观测量A^\hat{A}。我们可以问:当测量的次数很多时,A^\hat{A} 的平均测量值是多少,答案有A^\hat{A} 的系综平均值来确定,它被定义为:

A^=iwiα(i)A^α(i)=iwiα(i)(jajaj)A^(kakak)α(i)=ijkwiajA^akα(i)ajakα(i)=ijkwiakδj,kα(i)ajakα(i)=ijwiajajα(i)2(0.2.9)\begin{aligned} \braket{\hat{A}} &= \sum_i w_i \braket{\alpha^{(i)}|\hat{A}|\alpha^{(i)}} \\&= \sum_i w_i \bra{\alpha^{(i)}} \left( \sum_j \ket{a_j}\bra{a_j} \right) \hat{A} \left( \sum_k \ket{a_k}\bra{a_k} \right) \ket{\alpha^{(i)}} \\&= \sum_i \sum_j \sum_k w_i \braket{a_j|\hat{A}|a_k} \braket{\alpha^{(i)}|a_j} \braket{a_k|\alpha^{(i)}} \\&= \sum_i \sum_j \sum_k w_i \ a_k \ \delta_{j,k} \braket{\alpha^{(i)}|a_j} \braket{a_k|\alpha^{(i)}} \\&= \sum_i \sum_j w_i \ a_j |\braket{a_j|\alpha^{(i)}}|^2 \end{aligned} \tag{0.2.9}

其中,aj\ket{a_j} 代表这A^\hat{A} 的对应本征值为aja_j 的本征态。我们仔细观测上式,不难注意到计算平均值是怎么引入了两次概率的概念:第一次是在A^\hat{A} 的一个本征态aj\ket{a_j} 中找到态α(i)\ket{\alpha^{(i)}} 的量子力学概率ajα(i)2|\braket{a_j|\alpha^{(i)}}|^2,而第二次是在系综中找到一个由α(i)\ket{\alpha^{(i)}} 表征的量子力学的概率因子wiw_i

\quad 现在,我们考虑不要A^\hat{A} 的本征态a\ket{a} 作为基,而是考虑利用任意一个更一般的基b\ket{b} 来改写系综的平均值(0.2.9)(0.2.9) 式:

A^=iwiα(i)A^α(i)=iwiα(i)(jbjbj)A^(kbkbk)α(i)=jk(iwibkα(i)α(i)bj)bjA^bk(0.2.10)\begin{aligned} \braket{\hat{A}} &= \sum_i w_i \braket{\alpha^{(i)}|\hat{A}|\alpha^{(i)}} \\&= \sum_i w_i \bra{\alpha^{(i)}} \left( \sum_j \ket{b_j}\bra{b_j} \right) \hat{A} \left( \sum_k \ket{b_k}\bra{b_k} \right) \ket{\alpha^{(i)}} \\&= \sum_j \sum_k \left( \sum_i w_i \braket{b_k|\alpha^{(i)}} \braket{\alpha^{(i)}|b_j} \right) \braket{b_j|\hat{A}|b_k} \end{aligned} \tag{0.2.10}

此时,我们定义一个密度算符ρ^\hat{\rho},如下:

ρ^iwiα(i)α(i)(0.2.11)\hat{\rho} \equiv \sum_i w_i \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} \tag{0.2.11}

这个密度算符包含了我们可能得到的,关于所论及系综的所有有物理意义的信息,将它带回(0.2.10)(0.2.10) 式,得:

A^=jkbkρ^bjbjA^bk=kbkρ^(jbjbj)A^bk=kbkρ^A^bk=tr(ρ^A^)(0.2.12)\begin{aligned} \braket{\hat{A}} &= \sum_j \sum_k \braket{b_k|\hat{\rho}|b_j} \braket{b_j|\hat{A}|b_k} \\&= \sum_k \bra{b_k} \hat{\rho} \left( \sum_j \ket{b_j} \bra{b_j} \right) \hat{A} \ket{b_k} \\&= \sum_k \braket{b_k|\hat{\rho}\hat{A}|b_k} \\&= tr(\hat{\rho}\hat{A}) \end{aligned} \tag{0.2.12}

上式就说明了:系综中,A^\hat{A} 的平均值为算符ρ^A^\hat{\rho}\hat{A} 的迹!因为矩阵的迹不依赖于表象,引入tr(ρ^A^)tr(\hat{\rho}\hat{A}) 可以利用任何方便的基来计算。

# 密度算符的性质

\quad 密度算符有几个性质值得记住。

  1. 密度算符是厄密的。 证明:

ρ^=(iwiα(i)α(i))=iwiα(i)α(i)=iwiα(i)α(i)=ρ^(0.2.13)\hat{\rho}^\dagger = \left(\sum_i w_i \ket{\alpha^{(i)}} \bra{\alpha^{(i)}}\right)^\dagger =\sum_i w_i \bra{\alpha^{(i)}}^\dagger \ket{\alpha^{(i)}}^\dagger = \sum_i w_i \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} =\hat{\rho} \tag{0.2.13}

证毕。

  1. 密度算符的迹为 1。证明:

tr(ρ^)=ijwibjα(i)α(i)bj=iwiα(i)(jbjbj)α(i)=iwiα(i)α(i)=iwi=1(0.2.14)\begin{aligned} tr(\hat{\rho}) &= \sum_i \sum_j w_i \braket{b_j|\alpha^{(i)}} \braket{\alpha^{(i)}|b_j} \\&= \sum_i w_i \bra{\alpha^{(i)}} \left( \sum_j \ket{b_j} \bra{b_j} \right) \ket{\alpha^{(i)}} \\&= \sum_i w_i \braket{\alpha^{(i)}|\alpha^{(i)}} = \sum_i w_i =1 \end{aligned} \tag{0.2.14}

证毕。这条性质也成为密度算符的归一化性质。

  1. 纯系综的密度算符满足tr(ρ^2)=1tr(\hat{\rho}^2)=1。证明:

ρ^=α(i)α(i)ρ^2=α(i)α(i)α(i)α(i)=ρ^\begin{aligned} \hat{\rho} &= \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} \\ \hat{\rho}^2 &= \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} = \hat{\rho} \end{aligned}

因此,有:

tr(ρ^2)=tr(ρ^)=1(0.2.15)tr(\hat{\rho}^2) = tr(\hat{\rho}) = 1 \tag{0.2.15}

证毕。

  1. 混合系综的密度算符满足tr(ρ^2)<1tr(\hat{\rho}^2)<1,且是个正数。证明:

ρ^=iwiα(i)α(i)\hat{\rho} = \sum_i w_i \ket{\alpha^{(i)}} \bra{\alpha^{(i)}}

ρ^2=(iwiα(i)α(i))(jwjα(j)α(j))=ijwiwjα(i)α(i)α(j)α(j)=iwi2α(i)α(i)(0.2.16)\begin{aligned} \hat{\rho}^2 &= (\sum_i w_i \ket{\alpha^{(i)}} \bra{\alpha^{(i)}}) (\sum_j w_j \ket{\alpha^{(j)}} \bra{\alpha^{(j)}}) \\&= \sum_i \sum_j w_i w_j \ket{\alpha^{(i)}} \braket{\alpha^{(i)}|\alpha^{(j)}} \bra{\alpha^{(j)}} \\&= \sum_i w_i^2 \ket{\alpha^{(i)}} \bra{\alpha^{(i)}} \end{aligned} \tag{0.2.16}

因此:

tr(ρ^2)=iwi2<iwi=1(0.2.17)tr(\hat{\rho}^2) = \sum_i w_i^2 < \sum_i w_i=1 \tag{0.2.17}

证毕。

\quad 下面给几个很有益的例子。

11:一束sz,+\ket{s_z,+} 的原子束纯系综

ρ^=sz,+sz,+=(10)(10)=(1000)(0.2.18)\hat{\rho} = \ket{s_z,+} \bra{s_z,+} = \left( \begin{matrix} 1 \\ 0 \end{matrix} \right) \left( \begin{matrix} 1 & 0 \end{matrix} \right) = \left( \begin{matrix} 1 & 0 \\ 0 & 0 \end{matrix} \right) \tag{0.2.18}

ρ^2=(1000)(1000)=(1000)(0.2.19)\hat{\rho}^2 = \left( \begin{matrix} 1 & 0 \\ 0 & 0 \end{matrix} \right) \left( \begin{matrix} 1 & 0 \\ 0 & 0 \end{matrix} \right) = \left( \begin{matrix} 1 & 0 \\ 0 & 0 \end{matrix} \right) \tag{0.2.19}

很显然,tr(ρ^2)=1tr(\hat{\rho}^2)=1,而它的确是一个纯系综。

22:一束sx,+\ket{s_x,+} 的原子束纯系综

ρ^=sx,+sx,+=12(11)12(11)=(1/21/21/21/2)(0.2.20)\hat{\rho} = \ket{s_x,+} \bra{s_x,+} = \frac{1}{\sqrt{2}} \left( \begin{matrix} 1 \\ 1 \end{matrix} \right) \frac{1}{\sqrt{2}} \left( \begin{matrix} 1 & 1 \end{matrix} \right) = \left( \begin{matrix} 1/2 & 1/2 \\ 1/2 & 1/2 \end{matrix} \right) \tag{0.2.20}

ρ^2=(1/21/21/21/2)(1/21/21/21/2)=(1/21/21/21/2)(0.2.21)\hat{\rho}^2 = \left( \begin{matrix} 1/2 & 1/2 \\ 1/2 & 1/2 \end{matrix} \right) \left( \begin{matrix} 1/2 & 1/2 \\ 1/2 & 1/2 \end{matrix} \right) = \left( \begin{matrix} 1/2 & 1/2 \\ 1/2 & 1/2 \end{matrix} \right) \tag{0.2.21}

很明显,tr(ρ^2)=1tr(\hat{\rho}^2)=1,而它的确是一个纯系综。

33:一个zz 轴自旋向上的系综和一个zz 轴自旋向下的系综以相等权重混合

ρ^=12sz,+sz,++12sz,sz,=(1/2001/2)(0.2.22)\hat{\rho} = \frac{1}{2} \ket{s_z,+} \bra{s_z,+} + \frac{1}{2} \ket{s_z,-} \bra{s_z,-} = \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) \tag{0.2.22}

ρ^2=(1/2001/2)(1/2001/2)=(1/4001/4)(0.2.23)\hat{\rho}^2 = \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) = \left( \begin{matrix} 1/4 & 0 \\ 0 & 1/4 \end{matrix} \right) \tag{0.2.23}

很明显,tr(ρ^2)=1/2<1tr(\hat{\rho}^2)=1/2<1,而它的确是一个混合系综。而且,若我们讨论这个系综各个方向自旋的平均值:

s^x=tr(ρ^s^x)=tr[2(1/2001/2)(0110)]=tr[2(01/21/20)]=0(0.2.24)\braket{\hat{s}_x} = tr(\hat{\rho} \hat{s}_x) = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) \left( \begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix} \right) \right] = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 0 & 1/2 \\ 1/2 & 0 \end{matrix} \right) \right] = 0 \tag{0.2.24}

s^y=tr(ρ^s^x)=tr[2(1/2001/2)(0ii0)]=tr[2(0i/2i/20)]=0(0.2.25)\braket{\hat{s}_y} = tr(\hat{\rho} \hat{s}_x) = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) \left( \begin{matrix} 0 & -i \\ i & 0 \end{matrix} \right) \right] = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 0 & -i/2 \\ i/2 & 0 \end{matrix} \right) \right] = 0 \tag{0.2.25}

s^z=tr(ρ^s^z)=tr[2(1/2001/2)(1001)]=tr[2(1/2001/2)]=0(0.2.26)\braket{\hat{s}_z} = tr(\hat{\rho} \hat{s}_z) = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 1/2 & 0 \\ 0 & 1/2 \end{matrix} \right) \left( \begin{matrix} 1 & 0 \\ 0 & -1 \end{matrix} \right) \right] = tr \left[ \frac{\hbar}{2} \left( \begin{matrix} 1/2 & 0 \\ 0 & -1/2 \end{matrix} \right) \right] = 0 \tag{0.2.26}

的确,我们凭直觉也不难想到这个系综的自旋平均值为零。

以上的例子默认采用了泡利表象。

# 约化密度算符

\quad 现在我们要考虑这样的情况:有一个大系统,而希望求平均值的物理量只与系统的一部分有关。例如,在两个粒子121,2 构成的系统中,希望求粒子11 的某一物理量F^(1)\hat{F}(1) 的平均值。这时上述所有内容当然仍然适用,不过可以做一些简化。
\quad 就以上述双粒子系统为例,设粒子11 有一组基底{φi}\{\ket{\varphi_i}\},粒子22 有一组基底{χm}\{\ket{\chi_m}\},则在121,2 粒子空间的直积空间中,系统的态矢的一般形式为:

ψ=imcimφiχm(0.2.27)\ket{\psi} = \sum_i \sum_m c_{im} \ket{\varphi_i} \otimes \ket{\chi_m} \tag{0.2.27}

为了使ψ\ket{\psi} 归一化,系数cimc_{im} 应满足:

imcim=1(0.2.28)\sum_i \sum_m c_{im} = 1 \tag{0.2.28}

那么,系统处于纯态 ψ\ket{\psi} 时的密度算符为:

ρ^=ψψ=iimmφiχmcimcimψiχm(0.2.29)\hat{\rho} = \ket{\psi}\bra{\psi} = \sum_{i'i} \sum_{m'm} \ket{\varphi_{i'}} \otimes \ket{\chi_{m'}} c_{i'm'}^* c_{im} \bra{\psi_i} \otimes \bra{\chi_m} \tag{0.2.29}

而密度矩阵的矩阵元则为:

ρ^im,im=cimcim(0.2.30)\hat{\rho}_{i'm',im} = c_{i'm'}^* c_{im} \tag{0.2.30}

\quad 现在求粒子11 的某物理量F^(1)\hat{F}(1) 的平均值,其实可以转化为求直积空间中,物理量F^(1)I^(2)\hat{F}(1)\otimes\hat{I}(2) 的平均值,利用(0.2.12)(0.2.12) 式,可得:

F^(1)=F^(1)I^(2)=tr[(F^(1)I^(2))ρ^]=imφiχm(F^(1)I^(2))ρ^φiχm=iimmφiχm(F^(1)I^(2))φiχmψiχmρ^φiχm=iimmφiF^(1)φiδm,mψiχmρ^φiχm=iiφiF^(1)φimφiχmρ^χmφi=iiφiF^(1)φiφitr2(ρ^)φi=iφiF^(1)tr2(ρ^)φi=tr(F^(1)tr2(ρ^))(0.2.31)\begin{aligned} \braket{\hat{F}(1)} &= \braket{\hat{F}(1)\otimes\hat{I}(2)} = tr \left[\left(\hat{F}(1)\otimes\hat{I}(2)\right) \hat{\rho}\right] \\&= \sum_i\sum_m \bra{\varphi_i} \otimes \bra{\chi_m} \left(\hat{F}(1)\otimes\hat{I}(2)\right) \hat{\rho} \ket{\varphi_i} \otimes \ket{\chi_m} \\&= \sum_{i'i} \sum_{m'm} \bra{\varphi_i} \otimes \bra{\chi_m} \left(\hat{F}(1)\otimes\hat{I}(2)\right) \ket{\varphi_{i'}} \otimes \ket{\chi_{m'}} \bra{\psi_{i'}} \otimes \bra{\chi_{m'}} \hat{\rho} \ket{\varphi_i} \otimes \ket{\chi_m} \\&= \sum_{i'i} \sum_{m'm} \braket{\varphi_i|\hat{F}(1)|\varphi_{i'}} \delta_{m,m'} \bra{\psi_{i'}} \otimes \bra{\chi_{m'}} \hat{\rho} \ket{\varphi_i} \otimes \ket{\chi_m} \\&= \sum_{i'i} \braket{\varphi_i|\hat{F}(1)|\varphi_{i'}} \sum_m \bra{\varphi_{i'}} \braket{\chi_m|\hat{\rho}|\chi_m} \ket{\varphi_i} \\&= \sum_{i'i} \braket{\varphi_i|\hat{F}(1)|\varphi_{i'}} \braket{\varphi_{i'}|tr_2(\hat{\rho})|\varphi_i} \\&= \sum_i \braket{\varphi_i|\hat{F}(1)tr_2(\hat{\rho})|\varphi_{i'}} \\&= tr \left(\hat{F}(1)tr_2(\hat{\rho})\right) \end{aligned} \tag{0.2.31}

上式中,tr2(ρ^)tr_2(\hat{\rho}) 的意思是只对粒子22 取迹,取迹完之后,它就从双粒子直积空间中的算符变成了粒子11 空间中的算符,记作ρ^(1)\hat{\rho}(1)

ρ^(1)=tr2(ρ^)=mχmρ^χm(0.2.32)\hat{\rho}(1) = tr_2(\hat{\rho}) = \sum_m \braket{\chi_m|\hat{\rho}|\chi_m} \tag{0.2.32}

这种ρ^(1)\hat{\rho}(1) 算符就称为描写粒子11约化密度算符。这样一来,F^(1)\hat{F}(1) 的平均值(0.2.31)(0.2.31) 式就可以改写为:

F^(1)=tr[F^(1)ρ^(1)](0.2.33)\braket{\hat{F}(1)} = tr \left[\hat{F}(1)\hat{\rho}(1)\right] \tag{0.2.33}

\quad 同理,粒子22 的约化密度算符ρ^(2)\hat{\rho}(2) 为:

ρ^(2)=tr1(ρ^)=iφiρ^φi(0.2.34)\hat{\rho}(2) = tr_1(\hat{\rho}) = \sum_i \braket{\varphi_i|\hat{\rho}|\varphi_i} \tag{0.2.34}

粒子22 的物理量F^(2)\hat{F}(2) 平均值为:

F^(2)=tr[F^(2)ρ^(2)](0.2.35)\braket{\hat{F}(2)} = tr \left[\hat{F}(2)\hat{\rho}(2)\right] \tag{0.2.35}

# 连续谱

\quad 大多数时候所考虑的可观测量都是被假定具有分立的本征值谱。然而,在量子力学中有一些可观测量具有连续的本征值。例如,位置算符r^\hat{r} 的本质值就是连续的,它能取(一定范围内的)任何实数。
\quad 对于一些分立本征值谱对应的本征态,可以把它们作为基,张成一个表象空间,这可能是有限维的(例如单粒子自旋算符s^z\hat{s}_z 的本征值就只有两个),也有可能是无穷维的(例如简谐振子的本质值就有分立的无穷多个)。而对于连续本征谱对应的本征态。它所张成的空间明显是无穷维的。很多时候有限维空间的结论是可以直接推广到无穷维空间中去的,但也有不可以的时候,这是值得注意的。

\quad 下面稍微来讲以下连续谱的一些写法。它的本征方程和分立谱的区别不大,都是

ξ^ξ=ξξ(0.3.1)\hat{\xi} \ket{\xi'} = \xi' \ket{\xi'} \tag{0.3.1}

对于某一态矢α\ket{\alpha},若用一组分立谱的基底叠加而成,用求和的形式;但若用一组连续谱的基底叠加,则要把求和改为积分:

α=iaiαaiα=ξαξdξ(0.3.2)\ket{\alpha} = \sum_i \braket{a_i|\alpha} \ket{a_i} \ \longrightarrow \ \ket{\alpha} = \int \braket{\xi'|\alpha} \ket{\xi'} d\xi' \tag{0.3.2}

其内积也要从求和改为积分:

iaiα2=1ξα2dξ(0.3.3)\sum_i |\braket{a_i|\alpha}|^2 = 1 \ \longrightarrow\ \int |\braket{\xi'|\alpha}|^2 d\xi' \tag{0.3.3}

总之,抓住求和改积分就行啦!ᖘ ❛‿˂̵✧

βα=iβaiaiαβα=βξaiξαdξ(0.3.4)\braket{\beta|\alpha} = \sum_i \braket{\beta|a_i} \braket{a_i|\alpha} \ \longrightarrow\ \braket{\beta|\alpha} = \int \braket{\beta|\xi'} \braket{a_i\xi'|\alpha} d\xi' \tag{0.3.4}

ajA^ai=aiδi,jξξ^ξ=ξδ(ξξ)(0.3.5)\braket{a_j|\hat{A}|a_i} = a_i \delta_{i,j} \ \longrightarrow\ \braket{\xi'|\hat{\xi}|\xi''} = \xi'' \delta(\xi'-\xi'') \tag{0.3.5}

# 平移

# 无穷小平移的具体形式

\quad 现在,我们引入非常重要的平移或空间位移的概念。我们定义一个无穷小平移算符T^(dr)\hat{\mathcal{T}}(d\vec{r})[2],它的作用就是

T^(dr)r=r+dr(0.4.1)\hat{\mathcal{T}}(d\vec{r}) \ket{\vec{r}} = \ket{\vec{r}+d\vec{r}} \tag{0.4.1}

上式中,理论上结果应该会有一个相因子,但我们约定,这个相因子为 1。
\quad 用位置本征态把任意一个态α\ket{\alpha} 展开,可以考查无穷小平移算符对于α\ket{\alpha} 的影响,如下:

T^(dr)α=T^(dr)rαrdr=rαr+drdr=rdrαrdr(0.4.2)\begin{aligned} \hat{\mathcal{T}}(d\vec{r}) \ket{\alpha} &= \hat{\mathcal{T}}(d\vec{r}) \int \braket{\vec{r}|\alpha} \ket{\vec{r}} d\vec{r} \\&= \int \braket{\vec{r}|\alpha} \ket{\vec{r}+d\vec{r}} d\vec{r} \\&= \int \braket{\vec{r}-d\vec{r}|\alpha} \ket{\vec{r}} d\vec{r} \end{aligned} \tag{0.4.2}

上式是因为积分是对全空间进行的,而r\vec{r}只是个积分变量。这表明,平移态T^(dr)α\hat{\mathcal{T}}(d\vec{r}) \ket{\alpha} 可以通过把rα\braket{\vec{r}|\alpha} 中的r\vec{r}替换成rdr\vec{r}-d\vec{r}得到。

\quad 下面,我们列举一下这种无穷小平移算符 “理所应当” 满足的性质。

  1. 无穷小平移算符T^(dr)\hat{\mathcal{T}}(d\vec{r}) 应当是幺正的。对于一个态矢α\ket{\alpha},我们很自然要求平移前后它都是满足归一条件的:

αα=T^(dr)αT^(dr)α=αT^(dr)T^(dr)α=1T^(dr)T^(dr)=I^(0.4.3)\begin{aligned} \braket{\alpha|\alpha} = \braket{\hat{\mathcal{T}}(d\vec{r})\alpha|\hat{\mathcal{T}}(d\vec{r})\alpha} &= \braket{\alpha|\hat{\mathcal{T}}^\dagger(d\vec{r})\hat{\mathcal{T}}(d\vec{r})\alpha} = 1 \\ &\Downarrow \\ \hat{\mathcal{T}}^\dagger(d\vec{r})&\hat{\mathcal{T}}(d\vec{r}) = \hat{I} \end{aligned} \tag{0.4.3}

  1. 我们要求两个相继无穷小的平移 —— 先平移dr1d\vec{r}_1,再平移dr2d\vec{r}_2,可以等效成一个平移dr1+dr2d\vec{r}_1+d\vec{r}_2 的单一操作:

T^(dr1)T^(dr2)=T^(dr1+dr2)(0.4.4)\hat{\mathcal{T}}(d\vec{r}_1)\hat{\mathcal{T}}(d\vec{r}_2) = \hat{\mathcal{T}}(d\vec{r}_1+d\vec{r}_2) \tag{0.4.4}

  1. 我们要求一个沿相反方向的平移,与沿着原来方向平移的逆是相同的:

T^(dr)=T^1(dr)(0.4.5)\hat{\mathcal{T}}(-d\vec{r}) = \hat{\mathcal{T}}^{-1}(d\vec{r}) \tag{0.4.5}

  1. dr0d\vec{r}\rightarrow0 时,平移算符约化为恒等算符:

limdr=0T^(dr)=I^(0.4.6)\lim_{d\vec{r}=0} \hat{\mathcal{T}}(d\vec{r}) = \hat{I} \tag{0.4.6}

\quad 上面这四条性质,无论是从直觉上,还是实际上都应该是必须满足的。而我们发现,如果将无穷小平移算符取成如下形式:

T^(dr)=I^iK^dr(0.4.7)\hat{\mathcal{T}}(d\vec{r}) = \hat{I} - i \hat{K} \cdot d\vec{r} \tag{0.4.7}

解释一下,上式的K^\hat{K} 完整应该写成K^=K^xex+K^yey+K^zez\hat{K}=\hat{K}_x\vec{e}_x+\hat{K}_y\vec{e}_y+\hat{K}_z\vec{e}_z,并且K^\hat{K} 和它的各个分量K^x,K^y,K^z\hat{K}_x,\hat{K}_y,\hat{K}_z 都是厄密的。那么可以验证这样的形式是满足性质 1 的:

T^(dr)T^(dr)=(I^iK^dr)(I^iK^dr)=(I^+iK^dr)(I^iK^dr)=I^i(K^K^)dr+o[(dr)2]I^(对于无穷小平移,可忽略二阶小量)(0.4.8)\begin{aligned} \hat{\mathcal{T}}^\dagger(d\vec{r})\hat{\mathcal{T}}(d\vec{r}) &= (\hat{I} - i \hat{K} \cdot d\vec{r})^\dagger(\hat{I} - i \hat{K} \cdot d\vec{r}) \\&= (\hat{I} + i \hat{K}^\dagger \cdot d\vec{r})(\hat{I} - i \hat{K} \cdot d\vec{r}) \\&= \hat{I} - i(\hat{K}-\hat{K}^\dagger) \cdot d\vec{r} + o[(d\vec{r})^2] \\&\approxeq \hat{I} \\& \text{\small(对于无穷小平移,可忽略二阶小量)} \end{aligned} \tag{0.4.8}

同样也可以验证它满足性质 2:

T^(dr1)T^(dr2)=(I^iK^dr1)(I^iK^dr2)I^iK^(dr1+dr2)=T^(dr1+dr2)(0.4.9)\begin{aligned} \hat{\mathcal{T}}(d\vec{r}_1)\hat{\mathcal{T}}(d\vec{r}_2) &= (\hat{I} - i \hat{K} \cdot d\vec{r}_1)(\hat{I} - i \hat{K} \cdot d\vec{r}_2) \\&\approxeq \hat{I} -i\hat{K}\cdot(d\vec{r}_1+d\vec{r}_2) \\&= \hat{\mathcal{T}}(d\vec{r}_1+d\vec{r}_2) \end{aligned} \tag{0.4.9}

而性质 3 和 4 也都不难验证是成立的。

# 无穷小平移算符与位置算符的对易关系

\quad 下面,我们再来讨论一下K^\hat{K}r^\hat{r} 直接的对易关系。

x^T^(dx)x=x^x+dx=(x+dx)x+dx(0.4.10)\hat{x} \hat{\mathcal{T}}(dx) \ket{x} = \hat{x} \ket{x+dx} = (x+dx)\ket{x+dx} \tag{0.4.10}

T^(dx)x^x=xT^(dx)x=xx+dx(0.4.11)\hat{\mathcal{T}}(dx) \hat{x} \ket{x} = x \hat{\mathcal{T}}(dx)\ket{x} = x \ket{x+dx} \tag{0.4.11}

上两式相减得:

[x^,T^(dx)]x=dxx+dx=dx(x+o(dx))dxx(0.4.12)[\hat{x},\hat{\mathcal{T}}(dx)] \ket{x} = dx \ket{x+dx} = dx \left( \ket{x} + o(dx) \right) \approxeq dx \ket{x} \tag{0.4.12}

在上式中,我们忽略了近似中产生的二级小量,由此,得:

[x^,T^(dx)]=dx(0.4.13)[\hat{x},\hat{\mathcal{T}}(dx)] = dx \tag{0.4.13}

利用(0.4.7)(0.4.7),可将(0.4.13)(0.4.13) 改写为

[x^,T^(dx)]=[x^,I^iK^xdx]=idx[x^,K^x]=dx[\hat{x},\hat{\mathcal{T}}(dx)] = [\hat{x},\hat{I} - i \hat{K}_x dx] = -idx[\hat{x},\hat{K}_x] = dx

\Downarrow

[x^,K^x]=i(0.4.14)[\hat{x},\hat{K}_x] = i \tag{0.4.14}

同理推广,不难这么有如下对易关系:

[x^j,K^l]=iδkl(0.4.15)[\hat{x}_j,\hat{K}_l] = i \delta_{kl} \tag{0.4.15}

# 动量作为平移的生成元

\quad 我们虽然在(0.4.15)(0.4.15) 式中得到了K^\hat{K} 算符与位置算符r^\hat{r} 的对易关系,但我们从未考虑这个K^\hat{K} 算符究竟是什么。从对易关系来看,它很类似动量算符,只要令:

K^=p^(0.4.16)\hat{K} = \frac{\hat{p}}{\hbar} \tag{0.4.16}

对易关系就完美地回到位置与动量算符的对易关系。因此,我们不妨就直接把无穷小平移算符写成:

T^(dr)=I^ip^dr/(0.4.17)\hat{\mathcal{T}}(d\vec{r}) = \hat{I} - i \hat{p} \cdot d\vec{r} /\hbar \tag{0.4.17}

此时,动量就作为了一个平移的生成元。

\quad 这里我得念叨两句,在这里,我是早就知道了位置算符与动量算符的对易关系,所以直接把K^\hat{K} 变化为动量算符。而实际上,真正的量子力学发展可能并非这样。也许,正确的发展方法是从(0.4.15)(0.4.15) 式中,对应到经典力学的无穷小生成元中。因为在经典力学中,无穷小生成元就是动量,把这一概念借助到量子力学中,才得到了动量算符的对易关系,才定义了动量算符。因此,即使假设我前面完全没学过量子力学,完全也不知道动量算符,在这里,我也可以纯粹地把动量算符理解为,定义为一个量子力学中一个无穷小平移的生成元。

# 有限平移

\quad 到此为止,我们仅限于无穷小平移,而对于一个有限的平移,即一个有限大小的空间平移移动,可以通过相继地组合无穷平移得到。让我们考虑沿xx 方向平移一段Δx\Delta x 距离的有限平移:

T^(Δx)x=x+Δx(0.4.18)\hat{\mathcal{T}}(\Delta x) \ket{x} = \ket{x+\Delta x} \tag{0.4.18}

通过组合NN 次无穷小平移,这样得到表征这一有限平移过程的算符表达式:

T^(Δx)=(T^(ΔxN))N=limN(I^ip^xΔxN)N=exp(ip^xΔx)(0.4.19)\hat{\mathcal{T}}(\Delta x) = \left(\hat{\mathcal{T}}(\frac{\Delta x}{N})\right)^N = \lim_{N \rightarrow \infin} \left( \hat{I}-\frac{i\hat{p}_x\Delta x}{N\hbar} \right)^N = \exp \left(-\frac{i\hat{p}_x\Delta x}{\hbar}\right) \tag{0.4.19}

这里顺带提一下,指数函数上有一个算符,其实代表的是这样的意思:

exp(A^)1+A^+A^22!+(0.4.20)\exp (\hat{A}) \equiv 1 + \hat{A} + \frac{\hat{A}^2}{2!} + \cdots \tag{0.4.20}

\quad 我们再考虑两种平移过程,一种是先沿xx 方向平移Δx\Delta x,再沿yy 方向平移Δy\Delta y;另一种是沿yy 方向平移Δy\Delta y,再沿xx 方向平移Δx\Delta x。那么要问,这两种平移有什么不同呢?答案是:当然没有什么不同!从直觉上我们就可以认识到这一点,如下图。

图5

那么从数学上来看,这两个过程的算符分别是

T^(Δxex)T^(Δyey)=T^(Δxex+Δyey)(0.4.21)\hat{\mathcal{T}}(\Delta x \vec{e}_x) \hat{\mathcal{T}}(\Delta y \vec{e}_y) = \hat{\mathcal{T}}(\Delta x \vec{e}_x + \Delta y \vec{e}_y) \tag{0.4.21}

T^(Δyey)T^(Δxex)=T^(Δyey+Δxex)(0.4.22)\hat{\mathcal{T}}(\Delta y \vec{e}_y) \hat{\mathcal{T}}(\Delta x \vec{e}_x) = \hat{\mathcal{T}}(\Delta y \vec{e}_y + \Delta x \vec{e}_x) \tag{0.4.22}

如果这两个算符没有任何的不同,那么就要让[T^(Δxex),T^(Δyey)]0[\hat{\mathcal{T}}(\Delta x \vec{e}_x),\hat{\mathcal{T}}(\Delta y \vec{e}_y)]\equiv0,我们来利用平移算符的表达式来计算一下这个对易子:

[T^(Δxex),T^(Δyey)]=[exp(ip^xΔx),exp(ip^yΔy)]=[1+(ip^xΔx)+(ip^xΔx)22+,1+(ip^yΔy)+(ip^yΔy)22+]ΔxΔy[p^x,p^y]20(0.4.23)\begin{aligned} [\hat{\mathcal{T}}(\Delta x \vec{e}_x),\hat{\mathcal{T}}(\Delta y \vec{e}_y)] &= \left[\exp \left(-\frac{i\hat{p}_x\Delta x}{\hbar}\right),\exp \left(-\frac{i\hat{p}_y\Delta y}{\hbar}\right)\right] \\&= \left[ 1+(-\frac{i\hat{p}_x\Delta x}{\hbar})+\frac{(-\frac{i\hat{p}_x\Delta x}{\hbar})^2}{2}+\cdots , 1+(-\frac{i\hat{p}_y\Delta y}{\hbar})+\frac{(-\frac{i\hat{p}_y\Delta y}{\hbar})^2}{2}+\cdots \right] \\&\approxeq - \frac{\Delta x \Delta y [\hat{p}_x,\hat{p}_y]}{\hbar^2} \equiv 0 \end{aligned} \tag{0.4.23}

由于Δx,Δy\Delta x,\Delta y 都是任意的,因此要想满足[T^(Δxex),T^(Δyey)]0[\hat{\mathcal{T}}(\Delta x \vec{e}_x),\hat{\mathcal{T}}(\Delta y \vec{e}_y)]\equiv0 的要求,必须导致:

[p^x,p^y]=0(0.4.24)[\hat{p}_x,\hat{p}_y] = 0 \tag{0.4.24}

也不难证明更普遍的形式:

[p^i,p^j]=0(0.4.25)[\hat{p}_i,\hat{p}_j] = 0 \tag{0.4.25}

这个对易关系是沿不同方向平移相互对易的直接后果。在任何情况下,平移的生成元都是互相对易,相应的群被称为 阿贝尔 (Abel) 群。三维平移群就是阿贝尔群。在第三章中,我们将会讨论三维转动群,那个就是非阿贝尔群。

# 位置和动量空间中的波函数

# 位置空间波函数

\quad 很明显,位置算符 x^\hat{x} 的本征谱是一个连续谱,因此我们不难有:

x^x=xx(0.5.1)\hat{x}\ket{x}=x\ket{x} \tag{0.5.1}

xx=δ(xx)(0.5.2)\braket{x'|x}=\delta(x'-x) \tag{0.5.2}

而代表一个物理态的态矢量可以有连续本征谱x\ket{x} 来展开:

α=xαxdx(0.5.3)\ket{\alpha} = \int \braket{x|\alpha} \ket{x} dx \tag{0.5.3}

xα2dx|\braket{x|\alpha}|^2dx 代表了在 xx 附近一个狭窄间隔 dxdx 内找到例子的概率。在我们的形式中,内积 xα\braket{x|\alpha} 就是常说的 α\ket{\alpha} 态在位置表象中的波函数 ψα(x)\psi_\alpha(x)[^3] :

xα=ψα(x)(0.5.4)\braket{x|\alpha} = \psi_\alpha(x) \tag{0.5.4}

在我的很多笔记中,往往贪图方便,把一个态的表示ψk\ket{\psi_k} 写成ψk\psi_k 认真来说这样是不对的,但其实大体上也不会有太大的问题,因为接下来我们可以看到,用波函数描述还是用态矢描述都是可以的。不过以后还是会好好地把\ket{} 写出来的。QAQ
\quad 而两个态矢作内积时,我们也经常会利用位置表象写成下面的形式:

βα=dxβxxα=ψβ(x)ψα(x)dx(0.5.5)\braket{\beta|\alpha} = \int dx \braket{\beta|x} \braket{x|\alpha} = \int \psi_\beta^*(x) \psi_\alpha(x) \ dx \tag{0.5.5}

是不是以前量子力学常用的东西有了更本质的理解了呢?
\quad 还没完,我们接着看,不仅态矢的内积可以换成(位置)波函数的内积来描述,态矢的展开也可以。对于一个态矢 α\ket{\alpha},我们可以用一组本征态 ui\ket{u_i} 来展开,如下:

α=iuiαui(0.5.6)\ket{\alpha} = \sum_i \braket{u_i|\alpha} \ket{u_i} \tag{0.5.6}

如果用波函数的来描述就是:

xα=iuiαxuiψα(x)=iciψui(x)(0.5.7)\braket{x|\alpha} = \sum_i \braket{u_i|\alpha} \braket{x|u_i} \Longrightarrow \psi_\alpha(x) = \sum_i c_i \psi_{u_i}(x) \tag{0.5.7}

\quad 而且,我们计算一个算符的矩阵元 βf^α\bra{\beta}\hat{f}\ket{\alpha} 同样可以用波函数来表示。假设有算符 f^\hat{f} 满足:

f^x=f(x)x(0.5.8)\hat{f}\ket{x} = f(x)\ket{x} \tag{0.5.8}

那么它的矩阵元可以用波函数表示:

βf^α=βxxf^xxαdxdx=ψβ(x)f(x)δ(xx)ψα(x)dxdx=ψβ(x)f(x)ψα(x)dx(0.5.9)\begin{aligned} \bra{\beta}\hat{f}\ket{\alpha} &= \int \int \braket{\beta|x'} \braket{x'|\hat{f}|x''} \braket{x''|\alpha} dx' dx'' \\&= \int \int \psi_\beta^*(x') f(x'') \delta(x'-x'') \psi_\alpha(x'') dx' dx'' \\&= \int \psi^*_\beta(x) f(x) \psi_\alpha(x) dx \end{aligned} \tag{0.5.9}

看!虽然我之前我把波函数和态矢量认为是一个东西,很多时候想表示态矢量ψ\ket{\psi} 时都没加狄拉克 ket,但从上面看来,这样不会导致太大问题,因为态矢量表示的,很多时候也可以等价表示成波函数。不过啊,以后一定要注意,态矢量是态矢量,是一种 “矢量”,而波函数其实是态矢量与态矢量内积后的标量。

# 位置基底中的动量算符

\quad 在以前的量子力学课程中,我的第一认知里,动量算符是如下形式的:

p^x=ix(0.5.10)\hat{p}_x = -i\hbar\frac{\partial}{\partial x} \tag{0.5.10}

但其实这种是作用在波函数上的位置表象中的一种简写。我们现在从更本质的观点出发来看看位置基底中的动量算符。
\quad 我们的出发点是动量的定义是作为无穷小平移生成元,即利用到(0.4.17)(0.4.17) 式,我们将无穷小平移算符作用到态矢α\ket{\alpha} 上:

T^(Δx)α=(1ip^xΔx)α=dxT^(Δx)xxα=dxx+Δxxα=dxxxΔxα=dxx(xαΔxxxα)=αdxx(Δxxxα)(0.5.11)\begin{aligned} \hat{\mathcal{T}}(\Delta x)\ket{\alpha} &= (1-\frac{i\hat{p}_x\Delta x}{\hbar})\ket{\alpha} \\&= \int dx \hat{\mathcal{T}}(\Delta x) \ket{x} \braket{x|\alpha} \\&= \int dx \ket{x+\Delta x} \braket{x|\alpha} \\&= \int dx \ket{x} \braket{x-\Delta x|\alpha} \\&= \int dx \ket{x} (\braket{x|\alpha}-\Delta x \frac{\partial}{\partial x}\braket{x|\alpha}) \\&= \ket{\alpha} - \int dx \ket{x} (\Delta x \frac{\partial}{\partial x}\braket{x|\alpha}) \end{aligned} \tag{0.5.11}

比较上式的左右两边,可得:

p^xα=dxx(ixxα)(0.5.12)\hat{p}_x\ket{\alpha} = \int dx \ket{x} (-i\hbar\frac{\partial}{\partial x}\braket{x|\alpha}) \tag{0.5.12}

或将上式作用一个x\bra{x'},即:

xp^xα=dxxx(ixxα)=ixxα(0.5.13)\braket{x'|\hat{p}_x|\alpha} = \int dx \braket{x'|x} (-i\hbar\frac{\partial}{\partial x}\braket{x|\alpha}) = -i\hbar\frac{\partial}{\partial x'}\braket{x'|\alpha} \tag{0.5.13}

由上式,我们就不难得到动量算符p^x\hat{p}_x 在位置表象下的矩阵元:

xp^xx=ixxx=ixδ(xx)(0.5.14)\braket{x'|\hat{p}_x|x''} = -i\hbar\frac{\partial}{\partial x'}\braket{x'|x''} = -i\hbar\frac{\partial}{\partial x'}\delta(x'-x'') \tag{0.5.14}

上式才是动量算符在位置表象的具体写法,该式在以后会经常用到!
\quad 我们以前往往会这样简写:

p^xψ(x)=ixψ(x)(0.5.15)\hat{p}_x \psi(x) = -i\hbar \frac{\partial}{\partial x}\psi(x) \tag{0.5.15}

但,实际上,更加完整的写法应该如下:

p^xψ(x)=xp^xψ=dxdxxxxp^xxxψ=dxdxδ(xx)(ixδ(xx))ψ(x)=dx(ixδ(xx))ψ(x)=ixψ(x)(0.5.16)\begin{aligned} \hat{p}_x \psi(x) &= \bra{x}\hat{p}_x\ket{\psi} \\&= \iint dx' dx'' \braket{x|x'} \braket{x'|\hat{p}_x|x''} \braket{x''|\psi} \\&= \iint dx' dx'' \delta(x-x') \left(-i\hbar\frac{\partial}{\partial x'}\delta(x'-x'')\right) \psi(x'') \\&= \int dx'' \left(-i\hbar\frac{\partial}{\partial x}\delta(x-x'')\right) \psi(x'') \\&= -i\hbar \frac{\partial}{\partial x} \psi(x) \end{aligned} \tag{0.5.16}

上式的最后一步推导运用到了δ\delta 函数的两个性质:

+f(x)δ(xx0)dx=f(x0)(0.5.17)\int^{+\infin}_{-\infin} f(x) \delta'(x-x_0) dx = -f'(x_0) \tag{0.5.17}

δ(xx0)=δ(x0x)(0.5.18)\delta'(x-x_0) = -\delta'(x_0-x) \tag{0.5.18}

(0.5.16)(0.5.16) 式,波函数表示法和态矢表示法又联系在一起啦!

# 动量空间波函数

\quad 上面我们都是在位置空间 (位置表象) 中进行的,现在我们在换到动量空间 (动量表象)。为了方便讨论,我们在这一小节中仍然只讨论一维空间,下面所写的 p^\hat{p}pp 其实更完整地应该写为 p^x,px\hat{p}_x,p_x
\quad 我们将动量算符 p^\hat{p} 的本征态规定为:

p^p=pp(0.5.19)\hat{p} \ket{p} = p \ket{p} \tag{0.5.19}

且本征态满足正交归一:

pp=δ(pp)(0.5.20)\braket{p'|p''} = \delta(p'-p'') \tag{0.5.20}

任意一个态矢α\ket{\alpha} 可以写成动量算符本质态谱的展开:

α=dppαp(0.5.21)\ket{\alpha} = \int dp \braket{p|\alpha} \ket{p} \tag{0.5.21}

按照惯例,pα\braket{p|\alpha} 就称为动量空间波函数,通常记为ψα(p)\psi_\alpha(p):

ψα(p)=pα(0.5.22)\psi_\alpha(p) = \braket{p|\alpha} \tag{0.5.22}

\quad 下面我们想要求出动量表象和坐标表象波函数的相互转化关系。不过,在这之前,我们先来看看动量本征态在坐标表象下的波函数 xp\braket{x|p} 的表达式是什么。为此,我们利用到 (0.5.13)(0.5.13) 式:

xp^p=ixxp=pxp(0.5.23)\braket{x|\hat{p}|p} = -i\hbar \frac{\partial}{\partial x} \braket{x|p} = p \braket{x|p} \tag{0.5.23}

上式是关于 xp\braket{x|p} 的微分方程,可以解出是:

xp=Nexp(ipx)(0.5.24)\braket{x|p} = N \exp(\frac{ipx}{\hbar} \tag{0.5.24})

其中,NN 是归一化常数。但我们可以看到,上式其实是一个平面波,如果将它规定为全空间都存在,那么很明显这个波函数不存在归一化常数,因为+dxxppx=+dxN2=1\int^{+\infin}_{-\infin}dx\braket{x|p}\braket{p|x}=\int^{+\infin}_{-\infin}dx|N|^2=1 明显要求 N 是无穷小。但如果我们像《固体物理》课程那样,给上周期边界条件,那么 pp 就会变成伪连续的,此时就可以求出归一化常数是 1/2π1/\sqrt{2\pi\hbar}。因此:

xp=12πexp(ipx)(0.5.25)\braket{x|p} = \frac{1}{\sqrt{2\pi\hbar}} \exp(\frac{ipx}{\hbar} \tag{0.5.25})

有了上式,我们就可以写出位置表象和动量表象的波函数相互之间的转化关系了:

xα=dpxppα(0.5.26a)\braket{x|\alpha} = \int dp \braket{x|p} \braket{p|\alpha} \tag{0.5.26a}

pα=dxpxxα(0.5.27b)\braket{p|\alpha} = \int dx \braket{p|x} \braket{x|\alpha} \tag{0.5.27b}

或写出:

傅里叶逆变化:ψα(x)=12πdpexp(ipx)ψα(p)(0.5.28a)\text{傅里叶逆变化:} \psi_\alpha(x) = \frac{1}{\sqrt{2\pi\hbar}} \int dp\ \exp(\frac{ipx}{\hbar}) \psi_\alpha(p) \tag{0.5.28a}

傅里叶变化:ψα(p)=12πdxexp(ipx)ψα(x)(0.5.28b)\text{傅里叶变化:} \psi_\alpha(p) = \frac{1}{\sqrt{2\pi\hbar}} \int dx\ \exp(-\frac{ipx}{\hbar}) \psi_\alpha(x) \tag{0.5.28b}

我们也可以看见,这种变化其实就是傅里叶变换和逆变换!


  1. 这里虽然写了的mix\ket{mix},但严格意义上混合系综不能用态矢来描述,所以mix\ket{mix} 严格来说不能理解成态矢,把它理解成一个态函数吧。 ↩︎

  2. 平移的英文是translationtranslation,记得我在直博资格面试,考察英文文献翻译的时候就出现了这个单词,当时我第一反应是 “翻译” 的意思,但明显结合当时那篇文章的语境是不对的。我想了好久最后才勉强翻译出来是 “平动” 的意思。为了铭记那一时刻,我在这里写平移算符采用了T\mathcal{T},其代表的就是translationtranslation↩︎