量子力学教程(全)

按顺序阅读。

第一章

波函数 E ( x , t ) = E 0 e i ( k x − ω t ) E(x,t)=E_0e^{i(kx-\omega t)} E(x,t)=E0ei(kxωt),记为 D a w ( A = E 0 , k = k , ω = ω ) Daw(A=E_0,k=k,\omega=\omega) Daw(A=E0,k=k,ω=ω),或 E 0 D k ω ( x , t ) E_0\mathcal{D}_k^\omega(x,t) E0Dkω(x,t)。其中 ω \omega ω为角速度, k k k为波数,即每 2 π 2\pi 2π长度里有多少个波长。

D ω ( t ) = e − i ω t \mathcal{D}^\omega(t)=e^{-i\omega t} Dω(t)=et
D k ( x ) = e i k x \mathcal{D}_k(x)=e^{ikx} Dk(x)=eikx

状态函数 定态波函数

如果 Ψ ( x , t ) = ψ ( x ) f ( t ) \Psi(x,t)=\psi(x)f(t) Ψ(x,t)=ψ(x)f(t),则 Ψ ( x , t ) \Psi(x,t) Ψ(x,t)是定态波函数,否则是状态函数。

f ( t ) f(t) f(t)有确定形式 f ( t ) = D E / ℏ ( t ) = D ω ( t ) f(t)=\mathcal{D}^{E/\hbar}(t)=\mathcal{D}^\omega(t) f(t)=DE/ℏ(t)=Dω(t)

自由粒子

如果是自由粒子或平面波,则 ψ ( x ) = A D k ( x ) \psi(x)=A\mathcal{D}_k(x) ψ(x)=ADk(x),其中 A A A是归一化常数。

于是 Ψ ( x , t ) = A D k ( x ) D ω ( t ) = A D k ω ( x , t ) \Psi(x,t)=A\mathcal{D}_k(x)\mathcal{D}^\omega(t)=A\mathcal{D}_k^\omega(x,t) Ψ(x,t)=ADk(x)Dω(t)=ADkω(x,t)

因此自由粒子和平面波的波函数形式和经典波相同。但根据哥本哈根学派诠释,这个波函数只是概率幅。

狄拉克符号 厄米共轭算符

一个态矢量 ψ ( x ) \psi(x) ψ(x)表示一个粒子的量子态。可以表示为 ∣ ψ ⟩ |\psi\rangle ψ,不再提到 x x x,即与坐标系选取无关。

一个封闭的狄拉克符号表示一个复数,如 C = ⟨ χ ∣ ψ ⟩ C=\langle\chi|\psi\rangle C=χψ。含义是内积,也即 C = ( χ , ψ ) C=(\chi,\psi) C=(χ,ψ)。对于向量可以看做矩阵乘法,对于函数可以看做函数内积。

在中间多乘一个矩阵可以表示为 ⟨ φ ∣ A ∣ ψ ⟩ \langle\varphi|A|\psi\rangle φAψ

投影算符 P φ = ∣ φ ⟩ ⟨ φ ∣ P_\varphi=|\varphi\rangle\langle\varphi| Pφ=φφ

算符 A A A的厄米共轭(即取转置,再每个元素取共轭)记为 A † A^\dagger A。如果 A ∣ α ⟩ = ∣ β ⟩ A|\alpha\rangle=|\beta\rangle Aα=β,那么 ⟨ β ∣ A † = ⟨ α ∣ \langle\beta|A^\dagger=\langle\alpha| βA=α

( ∣ u ⟩ ⟨ v ∣ ) † = ∣ v ⟩ ⟨ u ∣ (|u\rangle\langle v|)^\dagger=|v\rangle\langle u| (uv)=vu

如果 A = A † A=A^\dagger A=A,那么 A A A是厄米算符。

交换性质: ⟨ φ ∣ A ∣ ψ ⟩ = ⟨ ψ ∣ A † ∣ φ ⟩ ∗ \langle\varphi|A|\psi\rangle=\langle\psi|A^\dagger|\varphi\rangle^* φAψ=ψAφ; ⟨ φ ∣ ψ ⟩ = ⟨ ψ ∣ φ ⟩ ∗ \langle \varphi|\psi \rangle=\langle \psi|\varphi\rangle^* φψ=ψφ。其中 ∗ ^* 表示单个数的共轭。

毫无疑问, ∣ ψ ⟩ |\psi\rangle ψ ⟨ ψ ∣ \langle \psi| ψ之间是共轭转置关系,即 ∣ ψ ⟩ = ( ⟨ ψ ∣ ) † |\psi\rangle=(\langle\psi|)^\dagger ψ=(⟨ψ)。Dagger † \dagger 符号和矩阵论里的共轭转置 H H H上标是一回事。

克罗内克符号 正交归一关系

克罗内克符号:Kronecker Delta。

离散型: δ i j = { 1 , i = j 0 , i ≠ j \delta_{ij}=\left\{\begin{aligned}1,i=j\\0,i\ne j\end{aligned}\right. δij={1,i=j0,i=j

连续型: δ ( x ) = { 1 , x = 0  or  x → 0 0 , otherwise \delta(x)=\left\{\begin{aligned}1,&x=0\text{ or }x\to 0\\0,&\text{otherwise}\end{aligned}\right. δ(x)={1,0,x=0 or x0otherwise

正交基 { u i } \{u_i\} {ui}满足 ⟨ u i ∣ u j ⟩ = δ i j \langle u_i|u_j\rangle=\delta_{ij} uiuj=δij

连续基满足 ( ω α , ω α ′ ) = δ ( α − α ′ ) (\omega_\alpha,\omega_{\alpha'})=\delta(\alpha-\alpha') (ωα,ωα)=δ(αα)

∣ ψ ⟩ |\psi\rangle ψ在一组正交基上投影表示为 ∣ ψ ⟩ = ∑ i P u i ∣ ψ ⟩ |\psi\rangle=\sum_iP_{u_i}|\psi\rangle ψ=iPuiψ

必然有 ∑ i P u i = I \sum_iP_{u_i}=I iPui=I

连续基对应性质为 ∫ P ω α d α = 1 \int P_{\omega_\alpha}d\alpha=1 Pωαdα=1

如果矩阵 S = U T S=UT S=UT,那么 S † = ( U T ) † = T † U † S^\dagger=(UT)^\dagger=T^\dagger U^\dagger S=(UT)=TU。设 e e e表示单位正交基。那么 S i j = ⟨ e i ∣ S ∣ e j ⟩ = ⟨ e i ∣ U ∣ T ∣ e j ⟩ = ⟨ u ( r o w = i ) ∣ t ( c o l u m n = j ) ⟩ S_{ij}=\langle e_i|S|e_j\rangle=\langle e_i|U|T|e_j\rangle=\langle u(row=i)|t(column=j) \rangle Sij=eiSej=eiUTej=u(row=i)t(column=j)⟩ ( S † ) i j = ⟨ e i ∣ T † ∣ U † ∣ e j ⟩ = ⟨ t † ( r o w = i ) ∣ u † ( c o l u m n = j ) ⟩ (S^\dagger)_{ij}=\langle e_i|T^\dagger |U^\dagger|e_j\rangle=\langle t^\dagger(row=i)|u^\dagger(column=j)\rangle (S)ij=eiTUej=t(row=i)u(column=j)⟩。如果记 ⟨ u i ∣ = ⟨ u ( r o w = i ) ∣ \langle u_i|=\langle u(row=i)| ui=u(row=i),那么自然 ∣ u i ⟩ = ∣ u † ( c o l u m n = i ) ⟩ |u_i\rangle=|u^\dagger(column=i)\rangle ui=u(column=i)⟩。即,若 S i j = ⟨ u i ∣ t j ⟩ S_{ij}=\langle u_i|t_j\rangle Sij=uitj,则 ( S † ) i j = ⟨ t i ∣ u j ⟩ (S^\dagger)_{ij}=\langle t_i|u_j\rangle (S)ij=tiuj。另一种说明: ( S † ) i j = ( S j i ) ∗ = ⟨ u j ∣ t i ⟩ ∗ = ⟨ t i ∣ u j ⟩ (S^\dagger)_{ij}=(S_{ji})^*=\langle u_j|t_i\rangle^*=\langle t_i|u_j\rangle (S)ij=(Sji)=ujti=tiuj

如果 U U U T T T都是酉矩阵(幺正矩阵),那么 S S S也是酉矩阵。我们说 { ∣ u i ⟩ } \{|u_i\rangle\} {ui⟩} { ∣ t i ⟩ } \{|t_i\rangle\} {ti⟩}表象(representation) S S S称为变换矩阵

本征值方程 简并度

A ∣ ψ ⟩ = λ ∣ ψ ⟩ A|\psi\rangle=\lambda|\psi\rangle Aψ=λψ,则称这个方程为 A A A的本征值方程, λ \lambda λ为其一个本征值(eigenvalue), ∣ ψ ⟩ |\psi\rangle ψ为其对应的本征矢(eigenvector)。

对于方程 A ∣ ψ ⟩ = λ ∣ ψ ⟩ A|\psi\rangle=\lambda|\psi\rangle Aψ=λψ,其中 A A A已知,如果 λ \lambda λ是一个本征值,那么此时方程的解空间称为本征值 λ \lambda λ的本征子空间(eigensubspace)。如果本征子空间里只有一个向量,即是一维的,那么称本征值 λ \lambda λ非简并的;如果大于一维,那么是简并的(degenerate),子空间维数称为简并度(degeneracy),也就是几何重数(geometric multiplicity)。

如果 A A A是厄米的,那么本征值的简并度总是等于重数,即几何重数等于代数重数。

如果 A A A是厄米的,那么 A A A的本征值都是实数,且不同本征值对应的本征矢都正交。

波函数与归一化 叠加原理

某一时刻,一个粒子的状态用 Ψ ( q ) \Psi(q) Ψ(q)表示,其中 q q q是所有可能用到的坐标,这个函数称为波函数。粒子在体积元 d q dq dq内的概率为 Ψ ( q ) ( Ψ ( q ) ) ∗ d q = ∣ Ψ ( q ) ∣ 2 d q \Psi(q)(\Psi(q))^*dq=|\Psi(q)|^2dq Ψ(q)(Ψ(q))dq=∣Ψ(q)2dq。以坐标为波函数的参数的量子力学称为坐标表象下的量子力学,或波动力学。

我们自然有归一化条件: ∫ Ψ ( q ) ( Ψ ( q ) ) ∗ d q = 1 \int \Psi(q)(\Psi(q))^*dq =1 Ψ(q)(Ψ(q))dq=1

波函数可以是几个波函数的和。即 Ψ = ∑ i c i Ψ i \Psi=\sum_ic_i\Psi_i Ψ=iciΨi,其中 c i c_i ci是复数。但这还不能称为叠加原理(superposition principle)。波函数可以表示一种事件,比如猫是死的,或者骰子掷到6点(怎么理解?见之后的章节)。如果 Ψ i \Psi_i Ψi表示的事件之间是互斥的,例如 Ψ 1 \Psi_1 Ψ1表示一只猫是死的, Ψ 2 \Psi_2 Ψ2表示这只猫是活的(猫不能同时死和活),那么这才称为状态叠加原理。这些态称为本征态(eigenstate)。互斥也就是说这两个波函数是正交的,即 ( Ψ i , Ψ j ) = 0 (\Psi_i,\Psi_j)=0 (Ψi,Ψj)=0 ∣ c i ∣ 2 |c_i|^2 ci2表示一种无法解释的概率,称为不解释#1概率。也有积分形式 Ψ = ∫ c n Ψ n d n \Psi=\int c_n\Psi_n dn Ψ=cnΨndn。由于这些状态互斥,因此 ∑ i ∣ c i ∣ 2 = 1 \sum_i|c_i|^2=1 ici2=1

对于物理量 f f f,如果它可能有 n n n个取值 f i f_i fi,称这些值为 f f f的笨征值(baconvalue)。 f f f f i f_i fi这个事件用已经归一化的波函数 Ψ i \Psi_i Ψi表示。那么 f f f的期望 f ˉ = ∑ i f i ∣ c i ∣ 2 \bar f=\sum_if_i|c_i|^2 fˉ=ifici2。(你信吗?)

叠加原理用狄拉克符号可以表示为 ∣ Ψ ⟩ = ∑ i a i ∣ Ψ i ⟩ |\Psi\rangle=\sum_i a_i|\Psi_i\rangle ∣Ψ=iaiΨi ∣ Ψ i ⟩ |\Psi_i \rangle Ψi记为 ∣ i ⟩ |i\rangle i,则 ∣ Ψ ⟩ = ∑ i a i ∣ i ⟩ |\Psi\rangle=\sum_ia_i|i\rangle ∣Ψ=iaii

正交条件可以表示为 ⟨ Ψ i ∣ Ψ j ⟩ = δ i j \langle \Psi_i|\Psi_j\rangle=\delta_{ij} ΨiΨj=δij,即 ⟨ i ∣ j ⟩ = δ i j \langle i|j\rangle=\delta_{ij} ij=δij

一个本征态与波函数的内积就是幅前系数: ⟨ i ∣ Ψ ⟩ = ⟨ i ∣ ∑ j a j ∣ j ⟩ = ∑ j a j ⟨ i ∣ j ⟩ = a i \langle i|\Psi\rangle=\langle i|\sum_ja_j|j\rangle=\sum_ja_j\langle i|j\rangle=a_i i∣Ψ=ijajj=jajij=ai。也就是说, a i a_i ai就是波函数到这个本征态的投影长,和投影到一组单位正交基是一个道理。

上式两边取共轭可得 ⟨ Ψ ∣ i ⟩ = a i ∗ \langle \Psi|i\rangle=a_i^* Ψ∣i=ai

那么 f f f的期望 f ˉ = ∑ i f i a i a i ∗ = ∑ i f i a i ⟨ Ψ ∣ i ⟩ = ⟨ Ψ ∣ ∑ i f i a i ∣ i ⟩ \bar f=\sum_if_ia_ia_i^*=\sum_if_ia_i\langle\Psi|i\rangle=\langle\Psi|\sum_if_ia_i|i\rangle fˉ=ifiaiai=ifiaiΨ∣i=Ψ∣ifiaii

物理量对应的算符

物理量 f f f对应的算符为 f ^ \hat f f^。算符的阶数与 f f f的可能取值数相同。算符作用于波函数时,会将其每个带系本征态乘上 f f f的该取值再求和。即 f ^ ∣ Ψ ⟩ = ∑ i f i a i ∣ i ⟩ \hat f|\Psi\rangle=\sum_if_ia_i|i\rangle f^∣Ψ=ifiaii

例如,若 Ψ = ( 2 2 2 2 ) \Psi=\begin{pmatrix}\frac{\sqrt 2}{2}\\\frac{\sqrt 2}{2}\end{pmatrix} Ψ=(22 22 ),取 ∣ 1 ⟩ = ( 1 0 ) |1\rangle=\begin{pmatrix}1\\0\end{pmatrix} ∣1=(10) ∣ 2 ⟩ = ( 0 1 ) |2\rangle=\begin{pmatrix}0\\1\end{pmatrix} ∣2=(01)。假设 f f f可以取 3 3 3 4 4 4,概率用上述波函数表示。那么 f ^ \hat f f^就是 ( 3 4 ) \begin{pmatrix}3&\\&4\end{pmatrix} (34)

那么 f ˉ \bar f fˉ可以表示为 ⟨ Ψ ∣ f ^ ∣ Ψ ⟩ \langle\Psi|\hat f|\Psi\rangle Ψ∣f^∣Ψ。可以简写为 ⟨ f ⟩ \langle f \rangle f

另外,算符作用于 ∣ i ⟩ |i\rangle i时,由于 ∣ i ⟩ |i\rangle i只有一个本征态,即自己, ∣ i ⟩ = 1 × ∣ i ⟩ |i\rangle=1\times|i\rangle i=1×i,那么 f ^ ∣ i ⟩ = f i ∣ i ⟩ \hat f|i\rangle=f_i|i\rangle f^i=fii,其中 f i f_i fi表示 f f f的这个本征态对应的取值。于是我们发现,矩阵表示下,本征态就是该物理量算符的本征矢,对应本征值就是对应笨征值。

本征态也可以是本征函数(eigenfunction),取决于你认为本征态是什么。当然,我们混淆这三个概念也无妨。

双下标等于两边乘写法: f i j = ⟨ i ∣ f ^ ∣ j ⟩ = ⟨ i ∣ f j ∣ j ⟩ = f j ⟨ i ∣ j ⟩ = f j δ i j f_{ij}=\langle i|\hat f|j\rangle=\langle i|f_j|j\rangle=f_j\langle i|j\rangle=f_j\delta_{ij} fij=if^j=ifjj=fjij=fjδij。称双下标等于两边乘写法表示的矩阵 [ f i j ] [f_{ij}] [fij]跃迁矩阵,是对角矩阵。上面的例子只是刚好 f ^ \hat f f^是跃迁矩阵,若本征态不取单位矩阵的列那么就很可能不是。

物理量与表象 观察算符 得某一结果的概率

物理量 f f f的本征值、本征态与波函数是确定的,但是这些本征态可以张成多种基矢,这些基矢又张成了同样的波函数,这种现象称为 f f f有不同的表象。也就是说本征态是 f f f的其中一个表象。注意,表象的基矢也必须是正交归一的。

设有两组连续表象基矢 { ∣ u ⟩ } \{|u\rangle\} {u⟩} { ∣ t ⟩ } \{|t\rangle\} {t⟩} ∣ Ψ ⟩ = ∫ d u ∣ u ⟩ ⟨ u ∣ Ψ ⟩ |\Psi\rangle=\int du|u\rangle\langle u|\Psi\rangle ∣Ψ=duuu∣Ψ ∣ Ψ ⟩ = ∫ d t ∣ t ⟩ ⟨ t ∣ Ψ ⟩ |\Psi\rangle=\int dt|t\rangle\langle t|\Psi\rangle ∣Ψ=dttt∣Ψ。显然,利用投影算符的归一性, ⟨ t ∣ = ∫ d u ⟨ t ∣ u ⟩ ⟨ u ∣ \lang t|=\int du\lang t|u\rang\lang u| t=dutuu

可以记 ⟨ a ∣ Ψ ⟩ \lang a|\Psi\rang a∣Ψ Ψ ( a ) \Psi(a) Ψ(a)

如果厄米算符 A A A的所有线性无关的本征矢能构成一组态空间(state space)的基,那么算符 A A A观察算符(observation operator)。这组基也是一个表象,一般称为 A A A表象。

上述的“物理量对应的算符”,都是观察算符。 也就是说本征子空间肯定与态空间同张,物理量算符肯定是厄米算符。

一个本征值可能对应多个本征矢,即具有简并度,求结果为该本征值的概率时不能忘记简并度。记本征值 f i f_i fi的简并度为 g i g_i gi f i f_i fi对应的归一化本征矢分别记为 ∣ ψ i 1 ⟩ , ∣ ψ i 2 ⟩ , ∣ ψ i 3 ⟩ , … |\psi_i^1\rang,|\psi_i^2\rang,|\psi_i^3\rang,\dots ψi1,ψi2,ψi3,。则测量 f f f得结果 f i f_i fi的概率为 p i = ∑ j = 1 g i ∣ ⟨ ψ i j ∣ Ψ ⟩ ∣ 2 p_i=\sum_{j=1}^{g_i}|\lang\psi_i^j|\Psi\rang|^2 pi=j=1giψij∣Ψ2

第一章小结

测量一个物理量 f f f可能得到多个结果。这些结果称为本征值 { f i } \{f_i\} {fi}。测量一个物理量的结果用波函数 Ψ \Psi Ψ表示。测量一个物理量对应一个观察算符,这个观察算符的本征值就是所有可能结果,观察算符的本征态 { Ψ i } \{\Psi_i\} {Ψi}就是该物理量的本征态。波函数可以由一组单位正交基 { u i } \{u_i\} {ui}张成,这组基称为一个表象。表象必然可以由本征子空间张成。波函数与一个本征值对应的所有本征态的内积的模长平方和就是测量得这个值的概率。

第二章

薛定谔方程 自动归一化

一个粒子的波函数为 ψ ( x , t ) \psi(x,t) ψ(x,t)。势能为 V ( x , t ) V(x,t) V(x,t),是实函数。则其波函数满足薛定谔方程
− ℏ 2 2 m ∂ 2 ∂ x 2 ψ + V ψ = i ℏ ∂ ∂ t ψ -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi+V\psi=i\hbar\frac{\partial}{\partial t}\psi 2m2x22ψ+Vψ=itψ

众所周知,归一化的波函数才是有意义的。如果我们在 t = 0 t=0 t=0时求出了 ψ ( x ) \psi(x) ψ(x)的表达式,必须考虑在之后时间其是否仍然维持归一化。

对薛定谔方程两边取共轭得 − ℏ 2 2 m ∂ 2 ∂ x 2 ψ ∗ + V ψ ∗ = − i ℏ ∂ ∂ t ψ ∗ -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi^*+V\psi^*=-i\hbar\frac{\partial}{\partial t}\psi^* 2m2x22ψ+Vψ=itψ

波函数的模长平方对时间的导数 ∂ ∂ t ∣ ψ ∣ 2 = ∂ ψ ∂ t ψ ∗ + ∂ ψ ∗ ∂ t ψ = ( i ℏ 2 m ∇ 2 ψ − i V ℏ ψ ) ψ ∗ + ( − i ℏ 2 m ∇ 2 ψ ∗ + i V ℏ ψ ∗ ) ψ = i ℏ 2 m ( ( ∇ 2 ψ ) ψ ∗ − ( ∇ 2 ψ ∗ ) ψ ) = i ℏ 2 m ∇ ( ( ∇ ψ ) ψ ∗ − ψ ( ∇ ψ ∗ ) ) \frac{\partial}{\partial t}|\psi|^2=\frac{\partial \psi}{\partial t}\psi^*+\frac{\partial \psi^*}{\partial t}\psi=(\frac{i\hbar}{2m}\nabla^2\psi-i\frac{V}{\hbar}\psi)\psi^*+(-\frac{i\hbar}{2m}\nabla^2\psi^*+\frac{iV}{\hbar}\psi^*)\psi=\frac{i\hbar}{2m}((\nabla^2\psi)\psi^*-(\nabla^2\psi^*)\psi)=\frac{i\hbar}{2m}\nabla((\nabla \psi)\psi^*-\psi(\nabla \psi^*)) tψ2=tψψ+tψψ=(2mi2ψiVψ)ψ+(2mi2ψ+iVψ)ψ=2mi((2ψ)ψ(2ψ)ψ)=2mi((ψ)ψψ(ψ))

在整个空间上积分有 d d t ∫ d x ∣ ψ ∣ 2 = ∫ d x i ℏ 2 m ∇ ( ( ∇ ψ ) ψ ∗ − ψ ( ∇ ψ ∗ ) ) = i ℏ 2 m ( ( ∇ ψ ) ψ ∗ − ψ ( ∇ ψ ∗ ) ) ∣ − ∞ + ∞ \frac{d}{dt}\int dx|\psi|^2=\int dx\frac{i\hbar}{2m}\nabla((\nabla \psi)\psi^*-\psi(\nabla \psi^*))=\frac{i\hbar}{2m}((\nabla \psi)\psi^*-\psi(\nabla \psi^*))|^{+\infty}_{-\infty} dtddxψ2=dx2mi((ψ)ψψ(ψ))=2mi((ψ)ψψ(ψ))+

由于波函数在正负无穷处趋于 0 0 0(否则不能归一化),因此整个空间的概率和随时间不变。即若在 t = 0 t=0 t=0归一化,则以后都自动归一化。

( ∂ 2 [ ψ ∣ ∂ x 2 ) ∣ ψ ∗ ] = ∂ ∂ x ( ∂ [ ψ ∣ ∂ x ∣ ψ ∗ ] ) (\frac{\partial^2 [\psi|}{\partial x^2})|\psi^*]=\frac{\partial}{\partial x}(\frac{\partial[\psi|}{\partial x}|\psi^*]) (x22[ψ)ψ]=x(x[ψψ])

动量

x x x是一个粒子的一维坐标。其波函数为 Ψ ( x , t ) \Psi(x,t) Ψ(x,t)。则 x x x的期望是
∫ d x ⟨ Ψ ∣ x ∣ Ψ ⟩ \int dx\lang\Psi|x|\Psi\rang dxΨ∣x∣Ψ

该期望随时间的导数是(使用分部积分,原部肯定取 0 0 0,直接写导部) d ⟨ x ⟩ d t = ∫ x d x ∂ ∣ ψ ∣ 2 ∂ t = ∫ x d x i ℏ 2 m ∂ ∂ x ( ∂ [ ψ ∣ ∂ x ∣ ψ ∗ ] ) = − i ℏ 2 m ∫ d x ∂ [ ψ ∣ ∂ x ∣ ψ ∗ ] \frac{d\lang x\rang}{dt}=\int xdx\frac{\partial |\psi|^2}{\partial t}=\int xdx \frac{i\hbar}{2m}\frac{\partial}{\partial x}(\frac{\partial [\psi|}{\partial x}|\psi^*])=-\frac{i\hbar}{2m}\int dx\frac{\partial [\psi|}{\partial x}|\psi^*] dtdx=xdxtψ2=xdx2mix(x[ψψ])=2midxx[ψψ]

再分部积分一次得(明天补充细节) d ⟨ x ⟩ d t = − i ℏ m ∫ ∂ ψ ∂ x ψ ∗ CINVIK d x \frac{d\lang x\rang}{dt}=-\frac{i\hbar}{m}\int \frac{\partial \psi}{\partial x}\psi^*\text{CINVIK}dx dtdx=mixψψCINVIKdx

速度的期望就是位置的期望对时间的导数: v ˉ = d x ˉ d t \bar v=\frac{d\bar x}{dt} vˉ=dtdxˉ

所以 ⟨ p ⟩ = m d ⟨ x ⟩ d t = − i ℏ ∫ ∂ ψ ∂ x ψ ∗ CINVIK d x = ∫ ψ ∗ ( − i ℏ ) ∂ ∂ x ψ d x \lang p\rang=m\frac{d\lang x\rang}{dt}=-i\hbar\int \frac{\partial \psi}{\partial x}\psi^*\text{CINVIK}dx=\int\psi^*(-i\hbar)\frac{\partial}{\partial x}\psi dx p=mdtdx=ixψψCINVIKdx=ψ(i)xψdx

那么对于物理量 Q ( x , p ) Q(x,p) Q(x,p),有 ⟨ Q ( x , p ) ⟩ = ⟨ Q ( x , − i ℏ ∂ ∂ x ) ⟩ \lang Q(x,p)\rang=\lang Q(x,-i\hbar\frac{\partial}{\partial x})\rang Q(x,p)⟩=Q(x,ix)⟩

自由粒子

如果 V ≡ 0 V\equiv0 V0,那么这个粒子是自由粒子(free particle)。

代入薛定谔方程,得知微分方程必有此形式解: ψ ( x , t ) = A D k ω ( x , t ) \psi(x,t)=A\mathcal{D}_k^\omega(x,t) ψ(x,t)=ADkω(x,t)

∂ D k ω ∂ t = − i ω D k ω \frac{\partial\mathcal{D}_k^\omega}{\partial t}=-i\omega\mathcal{D}_k^\omega tDkω=Dkω ∂ 2 D k ω ∂ x 2 = − k 2 D k ω \frac{\partial^2\mathcal{D}_k^\omega}{\partial x^2}=-k^2\mathcal{D}_k^\omega x22Dkω=k2Dkω因此 ℏ k 2 2 m = ω \frac{\hbar k^2}{2m}=\omega 2mk2=ω

海森堡不确定性原理

德布罗意公式(de Broglie formula)给出,粒子的动量为 p = h λ = ℏ k p=\frac{h}{\lambda}=\frac{\hbar}{k} p=λh=k

Δ x \Delta x Δx Δ p \Delta p Δp是位置 x x x和动量 p p p的标准差,那么海森堡不确定性原理(Heisenberg uncertainty principle)给出 Δ x Δ p ≥ ℏ 2 \Delta x\Delta p\ge\frac{\hbar}{2} ΔxΔp2

薛定谔方程的定态解

解形式 Ψ ( x , t ) = ψ ( x ) φ ( t ) \Psi(x,t)=\psi(x)\varphi(t) Ψ(x,t)=ψ(x)φ(t)称为薛定谔方程的定态解。

假设势能 V V V只与 x x x有关。则将定态解代入薛定谔方程得 i ℏ ψ ( x ) d φ ( t ) d t = V ( x ) ψ ( x ) φ ( t ) − ℏ 2 2 m φ ( t ) d 2 ψ ( x ) d x 2 i\hbar\psi(x)\frac{d \varphi(t)}{dt}=V(x)\psi(x)\varphi(t)-\frac{\hbar^2}{2m}\varphi(t)\frac{d^2\psi(x)}{dx^2} iψ(x)dtdφ(t)=V(x)ψ(x)φ(t)2m2φ(t)dx2d2ψ(x)

两边同时除以波函数,得 i ℏ φ ( t ) d φ ( t ) d t = V ( x ) − ℏ 2 2 m ψ ( x ) d 2 ψ ( x ) d x 2 \frac{i\hbar}{\varphi(t)}\frac{d \varphi(t)}{dt}=V(x)-\frac{\hbar^2}{2m\psi(x)}\frac{d^2\psi(x)}{dx^2} φ(t)idtdφ(t)=V(x)2mψ(x)2dx2d2ψ(x)

该等式左边只与 t t t有关,而右边只与 x x x有关,这说明等式两边都等于一个常数。让此常数为 ℏ ω \hbar\omega ω,代入左边解得 φ ( t ) = A D ω ( t ) \varphi(t)=A\mathcal{D}^\omega(t) φ(t)=ADω(t)

A A A放入 ψ ( x ) \psi(x) ψ(x)中,得定态解形式 Ψ ( x , t ) = ψ ( x ) D ω ( t ) \Psi(x,t)=\psi(x)\mathcal{D}^\omega(t) Ψ(x,t)=ψ(x)Dω(t)

普朗克-爱因斯坦关系说明,定态就是能量为确定值 E = ℏ ω E=\hbar\omega E=ω的态,即能量的本征态。这与经典情形相符:势能不随时间变化时,能量应当是定值。

E = ℏ ω E=\hbar\omega E=ω,哈密顿算符 H ^ = − ℏ 2 2 m ∇ 2 + V ( r ) \hat H=-\frac{\hbar^2}{2m}\nabla^2+V(\bm{r}) H^=2m22+V(r),则薛定谔方程代入定态解两边约去 φ ( t ) \varphi(t) φ(t)给出 H ^ ψ ( r ) = E ψ ( r ) \hat H\psi(\bm{r})=E\psi(\bm{r}) H^ψ(r)=Eψ(r)哈密顿算符可以看做经典能量在坐标表象下的表示。

哈密顿算符的期望是 E E E,即 ⟨ H ^ ⟩ = E \lang \hat H\rang=E H^=E

H ^ 2 ψ = H ^ ( H ^ ψ ) = E H ^ ψ = E 2 ψ \hat H^2\psi=\hat H(\hat H\psi)=E\hat H\psi=E^2\psi H^2ψ=H^(H^ψ)=EH^ψ=E2ψ,有 H H H的标准差为 Δ H = ⟨ H 2 ⟩ − ⟨ H ⟩ 2 = E 2 − E 2 = 0 \Delta H=\sqrt {\lang H^2\rang-\lang H\rang^2}=\sqrt{E^2-E^2}=0 ΔH=H2H2 =E2E2 =0

上述方程就是 H ^ \hat H H^的本征值方程。那么能量的所有可能取值就是 H ^ \hat H H^的本征值。

定态叠加 E E E V V V的关系

E E E有多个可能取值,对应多个 H ^ \hat H H^的本征函数,记为 E n E_n En ψ n ( x ) \psi_n(x) ψn(x) n = 1 , 2 , 3 , … n=1,2,3,\dots n=1,2,3,。则 H ^ ψ n = E n ψ n \hat H\psi_n=E_n\psi_n H^ψn=Enψn。对应波函数为 Ψ n ( x , t ) = ψ n ( x ) D ω ( t ) = ψ n ( x ) D E / ℏ ( t ) \Psi_n(x,t)=\psi_n(x)\mathcal{D}^\omega(t)=\psi_n(x)\mathcal{D}^{E/\hbar}(t) Ψn(x,t)=ψn(x)Dω(t)=ψn(x)DE/ℏ(t)。特别地, Ψ n ( x , 0 ) = ψ n ( x ) \Psi_n(x,0)=\psi_n(x) Ψn(x,0)=ψn(x)

这是线性方程,所以几个波函数叠加的波函数也是方程的解。即 Ψ = ∑ n c n Ψ n \Psi=\sum_nc_n\Psi_n Ψ=ncnΨn。特别地, Ψ ( x , 0 ) = ∑ n c n ψ n ( x ) \Psi(x,0)=\sum_nc_n\psi_n(x) Ψ(x,0)=ncnψn(x)

定态下, ⟨ x ⟩ \lang x\rang x是定值,因此 ⟨ p ⟩ = 0 \lang p\rang=0 p=0

能量是动能和势能之和,因此能量必须大于 V V V的最小值,否则波函数不能归一化。

一维无限深方势阱(infinite square well)

V ( x ) V(x) V(x) [ 0 , a ] [0,a] [0,a]处为 0 0 0,其它地方为 ∞ \infty

在势能无限高的地方,粒子不可能到达,此处 ψ ( x ) = 0 \psi(x)=0 ψ(x)=0

在势阱内,有 − ℏ 2 2 m d 2 d x 2 ψ ( x ) = E ψ ( x ) -\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\psi(x)=E\psi(x) 2m2dx2d2ψ(x)=Eψ(x) k = 2 m E ℏ k=\frac{\sqrt{2mE}}{\hbar} k=2mE ,则 d 2 d x 2 ψ ( x ) = − k 2 ψ ( x ) \frac{d^2}{dx^2}\psi(x)=-k^2\psi(x) dx2d2ψ(x)=k2ψ(x)

这是谐振子(harmonic oscillator)方程,其解为 ψ ( x ) = A sin ⁡ ( k x ) + B cos ⁡ ( k x ) \psi(x)=A\sin(kx)+B\cos(kx) ψ(x)=Asin(kx)+Bcos(kx),记为 T k ( A , B ) ( x ) \mathcal{T}_k(A,B)(x) Tk(A,B)(x),也可以表示为 C 1 e i k x + C 2 e − i k x C_1e^{ikx}+C_2e^{-ikx} C1eikx+C2eikx,记为 C k ( C 1 , C 2 ) ( x ) \mathcal{C}_k(C_1,C_2)(x) Ck(C1,C2)(x)。转换关系为 C k ( C 1 , C 2 ) = C 1 cos ⁡ ( k x ) + i C 1 sin ⁡ ( k x ) + C 2 cos ⁡ ( − k x ) + i C 2 sin ⁡ ( − k x ) = ( C 1 + C 2 ) cos ⁡ ( k x ) + ( i C 1 − i C 2 ) sin ⁡ ( k x ) = T k ( i C 1 − i C 2 , C 1 + C 2 ) \mathcal{C}_k(C_1,C_2)=C_1\cos(kx)+iC_1\sin(kx)+C_2\cos(-kx)+iC_2\sin(-kx)=(C_1+C_2)\cos(kx)+(iC_1-iC_2)\sin(kx)=\mathcal{T}_k(iC_1-iC_2,C_1+C_2) Ck(C1,C2)=C1cos(kx)+iC1sin(kx)+C2cos(kx)+iC2sin(kx)=(C1+C2)cos(kx)+(iC1iC2)sin(kx)=Tk(iC1iC2,C1+C2),反过来是 T k ( A , B ) = C k ( B − i A 2 , B + i A 2 ) = i A 2 C k ( − 1 , 1 ) + B 2 C k ( 1 , 1 ) \mathcal{T}_k(A,B)=\mathcal{C}_k(\frac{B-iA}{2},\frac{B+iA}{2})=\frac{iA}{2}\mathcal{C}_k(-1,1)+\frac{B}{2}\mathcal{C}_k(1,1) Tk(A,B)=Ck(2BiA,2B+iA)=2iACk(1,1)+2BCk(1,1)

一般,波函数及其导数是连续的。那么有 ψ ( 0 ) = ψ ( a ) = 0 \psi(0)=\psi(a)=0 ψ(0)=ψ(a)=0。于是 B = 0 B=0 B=0,即 ψ ( x ) = T k ( A , 0 ) \psi(x)=\mathcal{T}_k(A,0) ψ(x)=Tk(A,0),进而 sin ⁡ ( k a ) = 0 \sin(ka)=0 sin(ka)=0。于是解出 k k k,其中 k k k 0 0 0是没有意义的: k n = n π a ( n = 1 , 2 , 3 , …   ) k_n=\frac{n\pi}{a}(n=1,2,3,\dots) kn=a(n=1,2,3,)对应的能量为 E n = ℏ 2 k n 2 2 m = n 2 π 2 ℏ 2 2 m a 2 E_n=\frac{\hbar^2k_n^2}{2m}=\frac{n^2\pi^2\hbar^2}{2ma^2} En=2m2kn2=2ma2n2π22

通过归一化可以推得 ∣ A ∣ = 2 a |A|=\sqrt{\frac{2}{a}} A=a2 ,这里我们直接取 A = 2 a A=\frac{2}{a} A=a2,因为不同的相位也没有意义。于是 ψ n ( x ) = 2 a sin ⁡ ( n π x a ) \psi_n(x)=\sqrt\frac{2}{a}\sin(\frac{n\pi x}{a}) ψn(x)=a2 sin(ax)

即该解该方程得到无数个解。波函数的图像可以类比波的图像,波的能量与波的频率的平方成正比。因此 ψ 1 \psi_1 ψ1能量最低,称为基态(ground state) ψ n \psi_n ψn的能量是 ψ 1 \psi_1 ψ1 n 2 n^2 n2倍,对于 n > 1 n>1 n>1,称为激发态(excited state)

这些波函数有一些性质:

  • 相对于势阱中心 ψ 1 ( x ) \psi_1(x) ψ1(x)是偶函数, ψ 2 ( x ) \psi_2(x) ψ2(x)是奇函数, ψ 3 ( x ) \psi_3(x) ψ3(x)是偶函数,依此类推。
  • 称波函数与 x x x轴的交点为node(节点,结点)。具体地, ψ n ( x ) \psi_n(x) ψn(x) n − 1 n-1 n1个node。
  • 它们相互正交(mutually orthogonal): ∫ ( ψ m ( x ) ) ∗ ψ n ( x ) d x = δ m n \int (\psi_m(x))^*\psi_n(x)dx=\delta_{mn} (ψm(x))ψn(x)dx=δmn,即 ⟨ ψ m ∣ ψ n ⟩ = δ m n \lang \psi_m|\psi_n\rang=\delta_{mn} ψmψn=δmn。我们说 ψ \psi ψ系( ψ ’s \psi\text{'s} ψ’s)是正交归一的(orthonormal)
  • 它们是完备的(complete),意思是任何一个函数 f ( x ) f(x) f(x)都可以用它们的线性组合表示, f ( x ) = ∑ n = 1 ∞ c n ψ n ( x ) f(x)=\sum_{n=1}^\infty c_n\psi_n(x) f(x)=n=1cnψn(x)。参考傅里叶变换。用傅里叶技巧(Fourier’s trick)可以导出 c n = ∫ ψ n ( x ) ∗ f ( x ) d x c_n=\int \psi_n(x)^*f(x)dx cn=ψn(x)f(x)dx,也就是 c n = ⟨ ψ n ∣ f ⟩ c_n=\lang\psi_n|f\rang cn=ψnf

如果势(potential)是对称的,那么奇偶函数性质成立。节点性质总是成立,不管势是什么形状。正交性质也对很多成立。完备性质也对很多成立。

如果我们知道了 Ψ ( x , 0 ) \Psi(x,0) Ψ(x,0),那么我们可以先求出 c n c_n cn,然后立即得到 Ψ ( x , t ) \Psi(x,t) Ψ(x,t)。对于不是一维无限深方的其它势阱,也可以用此方法,只是 ψ \psi ψ系不同。

∣ c n ∣ 2 |c_n|^2 cn2是能量 E E E E n E_n En的概率,必然满足 ∑ ∣ c n ∣ 2 = 1 \sum|c_n|^2=1 cn2=1。能量的期望值与时间无关,这是能量守恒在量子力学中的体现。

谐振子 升降算符

谐振子的势能是 V ( x ) = 1 2 m ω 2 x 2 = 1 2 k x 2 V(x)=\frac{1}{2}m\omega^2x^2=\frac{1}{2}kx^2 V(x)=21mω2x2=21kx2

代数解法

现在我们看到 H = p 2 2 m + V = 1 2 m [ p 2 + ( m ω x ) 2 ] H=\frac{p^2}{2m}+V=\frac{1}{2m}[ p^2+(m\omega x)^2] H=2mp2+V=2m1[p2+(x)2]

我们希望用复数平方差公式把中括号里的东西展开,即 u 2 + v 2 = ( − i u + v ) ( i u + v ) u^2+v^2=(-iu+v)(iu+v) u2+v2=(iu+v)(iu+v)。但是现在 p p p x x x不是对易的(commute),即 p x ≠ x p px\ne xp px=xp

a ± = 1 2 ℏ m ω ( ∓ i p + m ω x ) a_\pm=\frac{1}{\sqrt{2\hbar m\omega}}(\mp ip+m\omega x) a±=2ℏ 1(ip+x)

那么试算 a − a + = 1 2 ℏ m ω [ p 2 + ( m ω x ) 2 + i m ω p x − i m ω x p ] = 1 2 ℏ m ω [ p 2 + ( m ω x ) 2 − i m ω [ x , p ] ] = 1 2 ℏ m ω [ p 2 + ( m ω x ) 2 ] − i 2 ℏ [ x , p ] \begin{aligned}a_-a_+&=\frac{1}{2\hbar m\omega}[p^2+(m\omega x)^2+im\omega px-im\omega xp]\\&=\frac{1}{2\hbar m \omega}[p^2+(m\omega x)^2-im\omega[x,p]]\\&=\frac{1}{2\hbar m\omega}[p^2+(m\omega x)^2]-\frac{i}{2\hbar}[x,p]\end{aligned} aa+=2ℏ1[p2+(x)2+imωpximωxp]=2ℏ1[p2+(x)2imω[x,p]]=2ℏ1[p2+(x)2]2ℏi[x,p]其中方括号是对易子(commutator)。

易知 [ x , p ] = i ℏ [x,p]=i\hbar [x,p]=i,这被称为canonical commutation relation。

于是 a − a + = H ℏ ω + 1 2 a_-a_+=\frac{H}{\hbar \omega}+\frac{1}{2} aa+=ωH+21 H H H到左边,有 H = ℏ ω ( a − a + − 1 2 ) H=\hbar\omega(a_-a_+-\frac{1}{2}) H=ω(aa+21)

计算得 a + a − = H ℏ ω − 1 2 a_+a_-=\frac{H}{\hbar \omega}-\frac{1}{2} a+a=ωH21,于是 [ a − , a + ] = 1 [a_-,a_+]=1 [a,a+]=1。也有 H = ℏ ω ( a + a − + 1 2 ) H=\hbar\omega(a_+a_-+\frac{1}{2}) H=ω(a+a+21)

我们现在满足 H ψ = E ψ H\psi=E\psi Hψ=Eψ。现在计算 H ( a + ψ ) = ℏ ω ( a + a − + 1 2 ) ( a + ψ ) = ℏ ω ( a + a − a + + 1 2 a + ) ψ = ℏ ω a + ( a − a + + 1 2 ) ψ = a + ( H + ℏ ω ) ψ 注意升降算符变成 H 以后就可以到处对易了 = ( E + ℏ ω ) ( a + ψ ) \begin{aligned}H(a_+\psi)&=\hbar\omega(a_+a_-+\frac{1}{2})(a_+\psi)\\&=\hbar\omega(a_+a_-a_++\frac{1}{2}a_+)\psi\\&=\hbar\omega a_+(a_-a_++\frac{1}{2})\psi \\&= a_+(H+\hbar\omega)\psi &注意升降算符变成H以后就可以到处对易了\\&=(E+\hbar\omega)(a_+\psi)\end{aligned} H(a+ψ)=ω(a+a+21)(a+ψ)=ω(a+aa++21a+)ψ=ωa+(aa++21)ψ=a+(H+ω)ψ=(E+ω)(a+ψ)注意升降算符变成H以后就可以到处对易了

同样有 H ( a − ψ ) = ( E − ℏ ω ) ( a − ψ ) H(a_-\psi)=(E-\hbar\omega)(a_-\psi) H(aψ)=(Eω)(aψ)

因此我们把 a ± a_\pm a±叫做升降算符(ladder operators),分别是上升算符(raising operator)和下降算符(lowering operator)。

注意,我们可以一直将下降算符作用于波函数,但是一直下去可能会导致波函数不能归一化,实际不存在。必然有一个态 ψ 0 \psi_0 ψ0满足 a − ψ 0 = 0 a_-\psi_0=0 aψ0=0

尝试解出 ψ 0 \psi_0 ψ0。方程是 1 2 ℏ m ω ( i ( − i ℏ d d x ) + m ω x ) ψ 0 = 0 \frac{1}{2\hbar m\omega}(i (-i\hbar\frac{d}{dx})+m\omega x)\psi_0=0 2ℏ1(i(idxd)+x)ψ0=0 ℏ d ψ 0 ψ 0 = − m ω x d x ℏ ln ⁡ ψ 0 = − 1 2 m ω x 2 + constant ψ 0 = A e − m ω 2 ℏ \begin{aligned}\hbar\frac{d\psi_0}{\psi_0}&=-m\omega xdx\\\hbar\ln\psi_0&=-\frac{1}{2}m\omega x^2+\text{constant}\\\psi_0&=Ae^{-\frac{m\omega}{2\hbar}}\end{aligned} ψ0dψ0lnψ0ψ0=xdx=21x2+constant=Ae2ℏ

归一化给出 ∣ A ∣ 2 = m ω π ℏ |A|^2=\sqrt{\frac{m\omega}{\pi\hbar}} A2=π ,因此 ψ 0 = ( m ω π ℏ ) 1 4 e − m ω 2 ℏ x 2 \psi_0=(\frac{m\omega}{\pi\hbar})^\frac{1}{4}e^{-\frac{m\omega}{2\hbar}x^2} ψ0=(π)41e2ℏx2

代入薛定谔方程得对应能量 E 0 = 1 2 ℏ ω E_0=\frac{1}{2}\hbar\omega E0=21ω于是 E n = ( n + 1 2 ) ℏ ω E_n=(n+\frac{1}{2})\hbar\omega En=(n+21)ω

ψ n \psi_n ψn作用于 H H H E n E_n En右边得 ℏ ω ( a + a − + 1 2 ) ψ n = ( n + 1 2 ) ℏ ω ψ n \hbar\omega(a_+a_-+\frac{1}{2})\psi_n=(n+\frac{1}{2})\hbar\omega\psi_n ω(a+a+21)ψn=(n+21)ωψn所以 a + a − = n a_+a_-=n a+a=n同样可得 a − a + = n + 1 a_-a_+=n+1 aa+=n+1

通过升降关系可以得到其它所有态。注意升降算符作用于邻态后还要归一化才是本态。即 ψ n + 1 = A n + 1 a + ψ n \psi_{n+1}=A_{n+1}a_+\psi_n ψn+1=An+1a+ψn,一般写成 c n ψ n + 1 = a + ψ n c_n\psi_{n+1}=a_+\psi_n cnψn+1=a+ψn d n ψ n − 1 = a − ψ n d_n\psi_{n-1}=a_-\psi_n dnψn1=aψn

对于连续函数 f f f g g g,升降算符必然满足性质:(注意升降算符都在函数左边) 证明过两天补充

  • ⟨ f ∣ ( a + g ) ⟩ = ⟨ ( a − f ) ∣ g ⟩ \lang f|(a_+g)\rang=\lang(a_-f)|g\rang f(a+g)⟩=⟨(af)g
  • ⟨ f ∣ ( a − g ) ⟩ = ⟨ ( a + f ) ∣ g ⟩ \lang f|(a_-g)\rang=\lang(a_+f)|g\rang f(ag)⟩=⟨(a+f)g

那么 ⟨ a + ψ n ∣ a + ψ n ⟩ = ⟨ a − a + ψ n ∣ ψ n ⟩ = ( n + 1 ) ⟨ ψ n ∣ ψ n ⟩ = ( n + 1 ) ∫ ∣ ψ n ∣ 2 d x = n + 1 \lang a_+\psi_n|a_+\psi_n\rang=\lang a_-a_+\psi_n|\psi_n\rang=(n+1)\lang \psi_n|\psi_n\rang=(n+1)\int |\psi_n|^2dx=n+1 a+ψna+ψn=aa+ψnψn=(n+1)ψnψn=(n+1)ψn2dx=n+1
同样可得 ⟨ a − ψ n ∣ a − ψ n ⟩ = n \lang a_-\psi_n|a_-\psi_n\rang=n aψnaψn=n

所以 c n = n + 1 c_n=\sqrt{n+1} cn=n+1 d n = n d_n=\sqrt n dn=n

所以 ψ 1 = 0 + 1 a + ψ 0 = a + ψ 0 \psi_1=\sqrt{0+1}a_+\psi_0=a_+\psi_0 ψ1=0+1 a+ψ0=a+ψ0 ψ 2 = 2 a + ψ 1 = 2 ( a + ) 2 ψ 0 \psi_2=\sqrt 2a_+\psi_1=\sqrt 2(a_+)^2\psi_0 ψ2=2 a+ψ1=2 (a+)2ψ0 ψ 3 = 6 ( a + ) 3 ψ 0 \psi_3=\sqrt 6(a_+)^3\psi_0 ψ3=6 (a+)3ψ0

正交归一性:可以证明 ψ \psi ψ系是正交归一的。这意味着我们同样可以用傅里叶技巧(Fourier’s trick)来计算 c n c_n cn(见一维无限深方势阱对应部分),即 c n = ⟨ ψ n ∣ f ⟩ c_n=\lang \psi_n|f\rang cn=ψnf。同样, ∣ c n ∣ 2 |c_n|^2 cn2是取到态 n n n的概率,即能量取 E n E_n En的概率。

解析解法

令无量纲量 ξ = m w ℏ x \xi=\sqrt\frac{mw}{\hbar}x ξ=mw x。令 K K K表示能量中半份 ℏ ω \hbar\omega ω的份数,即 K = 2 E ℏ ω K=\frac{2E}{\hbar\omega} K=ω2E。则薛定谔方程化为 d 2 ψ d ξ 2 = ( ξ 2 − K ) ψ \frac{d^2\psi}{d\xi^2}=(\xi^2-K)\psi dξ2d2ψ=(ξ2K)ψ

ξ \xi ξ相当大时可以忽略 K K K,即 ψ P 2 ξ = ξ 2 ψ \psi P^2\xi=\xi^2\psi ψP2ξ=ξ2ψ。这个方程有近似 ψ ( x ) = A e − ξ 2 / 2 + B e ξ 2 / 2 \psi(x)=Ae^{-\xi^2/2}+Be^{\xi^2/2} ψ(x)=Aeξ2/2+Beξ2/2,记作 Q ( A , B ) ( ξ ) \mathcal{Q}(A,B)(\xi) Q(A,B)(ξ)

x i xi xi很大的时候 B □ B\square B项显然太大,没法归一化,因此不要了。最后就是 ψ ( ξ ) → ( ? ? ) e − ξ 2 / 2 , ξ 很大 \psi(\xi)\to (??)e^{-\xi^2/2},\quad\xi很大 ψ(ξ)(??)eξ2/2,ξ很大

我们希望 ( ? ? ) (??) (??)填的是 h ( ξ ) h(\xi) h(ξ)

这内容也太多了。

反正记 H n ( ξ ) H_n(\xi) Hn(ξ)厄米多项式(Hermite polynomials),则 ψ n ( x ) = ( m ω π ℏ ) 1 / 4 1 2 n n ! H n ( ξ ) e − ξ 2 / 2 \psi_n(x)=(\frac{m\omega}{\pi\hbar})^{1/4}\frac{1}{\sqrt{2^nn!}}H_n(\xi)e^{-\xi^2/2} ψn(x)=(π)1/42nn! 1Hn(ξ)eξ2/2其中二的幂分数和厄米多项式那部分就相当于代数解法里的 n n n个上升算符。

束缚态 散射态

束缚态:bound state

散射态:scattering state

德尔塔函数势阱和势垒

这里的德尔塔函数指的是狄拉克德尔塔函数(Dirac delta function)。该函数在0处取无穷大,在非0处取0,并且在实数域上积分是 1 1 1。事实上这不是函数,可以称作广义函数或者分布。

势能为 V ( x ) = − α δ ( x ) , α  is positive V(x)=-\alpha\delta(x),\alpha\text{ is positive} V(x)=αδ(x),α is positive。这是势阱情况,后面会讲势垒。

有束缚态 E < 0 E<0 E<0和散射态 E > 0 E>0 E>0两种解。因为这个势在研究处是0。

束缚态

对于 x < 0 x<0 x<0,根据德尔塔函数定义, V = 0 V=0 V=0。则 ψ P 2 x = − 2 m E ℏ 2 ψ ≡ κ 2 ψ \psi P^2x=-\frac{2mE}{\hbar^2}\psi\equiv\kappa^2\psi ψP2x=22mEψκ2ψ其中 κ = − 2 m E ℏ \kappa=\frac{\sqrt{-2mE}}{\hbar} κ=2mE

方程的解是 ψ = E κ ( A , B ) = A e − κ x + B e κ x \psi=\mathcal{E}_\kappa(A,B)=Ae^{-\kappa x}+Be^{\kappa x} ψ=Eκ(A,B)=Aeκx+Beκx

负无穷远处函数应该有界,所以 A = 0 A=0 A=0

x > 0 x>0 x>0解法类似,解是 ψ = E κ ( F , 0 ) 。 \psi=\mathcal{E}_\kappa(F,0)。 ψ=Eκ(F,0)

标准边界条件(standard boundary conditions)是:

  • ψ \psi ψ连续
  • V V V不是无限的地方, ψ P x \psi Px ψPx连续

在0点处衔接,因此 F = B F=B F=B。即 ψ ( x ) = B e − κ ∣ x ∣ \psi(x)=Be^{-\kappa |x|} ψ(x)=Beκx

对薛定谔方程在 − ϵ -\epsilon ϵ ϵ \epsilon ϵ积分,得 − ℏ 2 2 m ( − 2 κ B ) + ∫ − ϵ ϵ V ( x ) ψ ( x ) d x = E ∫ − ϵ ϵ ψ ( x ) d x -\frac{\hbar^2}{2m}(-2\kappa B)+\int_{-\epsilon}^{\epsilon}V(x)\psi(x)dx=E\int_{-\epsilon}^{\epsilon}\psi(x)dx 2m2(2κB)+ϵϵV(x)ψ(x)dx=Eϵϵψ(x)dx ℏ 2 κ B m − α ψ ( 0 ) = 0 \frac{\hbar^2\kappa B}{m}-\alpha\psi(0)=0 m2κBαψ(0)=0 α = ℏ 2 κ m \alpha=\frac{\hbar^2\kappa}{m} α=m2κ κ = m α ℏ 2 \kappa=\frac{m\alpha}{\hbar^2} κ=2mα E = − m α 2 2 ℏ E=-\frac{m\alpha^2}{2\hbar} E=2ℏmα2

其中标准边界条件其实就是 Δ ( d ψ d x ) = − 2 m α ℏ 2 ψ ( 0 ) \Delta(\frac{d\psi}{dx})=-\frac{2m\alpha}{\hbar^2}\psi(0) Δ(dxdψ)=22mαψ(0)

归一化条件给出 B = κ = m α ℏ B=\sqrt \kappa=\frac{\sqrt{m\alpha}}{\hbar} B=κ =mα 。有唯一一个束缚态解。

散射态

x < 0 x<0 x<0 ψ P 2 x = − 2 m E ℏ 2 ψ ≡ − k 2 ψ \psi P^2x=-\frac{2mE}{\hbar^2}\psi\equiv-k^2\psi ψP2x=22mEψk2ψ其中 k = 2 m E ℏ k=\frac{\sqrt{2mE}}{\hbar} k=2mE

方程的解是 C k ( A , B ) \mathcal{C}_k(A,B) Ck(A,B)。对 x > 0 x>0 x>0 ψ = C k ( F , G ) \psi=\mathcal{C}_k(F,G) ψ=Ck(F,G)。零点连续要求 A + B = F + G A+B=F+G A+B=F+G

Δ ( ψ P x ) = i k ( F − G − A + B ) \Delta(\psi Px)=ik(F-G-A+B) Δ(ψPx)=ik(FGA+B) ψ ( 0 ) = A + B \psi(0)=A+B ψ(0)=A+B。代入标准边界条件化简为 F − G = ( 1 + 2 i β ) A − ( 1 − 2 i β ) B where  β = m α ℏ 2 k F-G=(1+2i\beta)A-(1-2i\beta)B\quad\text{where }\beta=\frac{m\alpha}{\hbar^2k} FG=(1+2iβ)A(12iβ)Bwhere β=2kmα

可以用波的观点来看这四个系数。我们考虑一般散射实验,粒子是从一边射进来的,比如说左边,那么从右边射进来的就是0,即 G = 0 G=0 G=0

解之前的方程有 B = i β 1 − i β A B=\frac{i\beta}{1-i\beta}A B=1iβiβA F = 1 1 − i β A F=\frac{1}{1-i\beta}A F=1iβ1A

反射系数(reflection coefficient) R = ∣ B ∣ 2 ∣ A ∣ 2 = β 2 / ( 1 + β 2 ) R=\frac{|B|^2}{|A|^2}=\beta^2/(1+\beta^2) R=A2B2=β2/(1+β2)
透射系数(transmission coefficient) T = ∣ F ∣ 2 ∣ A ∣ 2 = 1 / ( 1 + β 2 ) T=\frac{|F|^2}{|A|^2}=1/(1+\beta^2) T=A2F2=1/(1+β2)
R + T = 1 R+T=1 R+T=1

势垒情况

由于Vmin=0,因此束缚态不存在,因为E<0不能归一化。观察散射态,反射、透射系数与 β 2 \beta^2 β2,也就是 α 2 \alpha^2 α2有关,也就是说势垒变成势阱反射、透射情况竟然不变。即使对于无限高势垒,也可能穿过。这称为量子隧穿(tunneling)

有限深方势阱

捏麻麻地,第二张终于要写完了。Lapuok it mimard koi Jee’manana, bojp si dai.

在-a~a,V=-V_0。其它地方V=0。也就是说这里V最大是0,而无限深方势阱V最小是0。

这个同样有束缚态和散射态两种情况。

束缚态

对于 x < − a x<-a x<a V = 0 V=0 V=0,又见到了熟悉的 ψ P 2 x = κ 2 ψ \psi P^2x=\kappa^2\psi ψP2x=κ2ψ κ \kappa κ的定义同前。

通解(general solution)是 ψ = E κ ( A , B ) \psi=\mathcal{E}_\kappa(A,B) ψ=Eκ(A,B),但是 A □ A\square A会爆炸,因此 ψ = E κ ( 0 , B ) \psi=\mathcal{E}_\kappa(0,B) ψ=Eκ(0,B)

x ∈ [ − a , a ] x\in[-a,a] x[a,a] V = − V 0 V=-V_0 V=V0。根据束缚态归一化条件 E > V min ⁡ = − V 0 E>V_{\min}=-V_0 E>Vmin=V0,因此 E + V 0 > 0 E+V_0>0 E+V0>0。有 ψ P 2 x = − l 2 ψ \psi P^2x=-l^2\psi ψP2x=l2ψ其中 l = 2 m ( E + V 0 ) ℏ l=\frac{\sqrt{2m(E+V_0)}}{\hbar} l=2m(E+V0)

通解是 ψ = T l ( C , D ) \psi=\mathcal{T}_l(C,D) ψ=Tl(C,D)。对于 x > a x>a x>a,无须赘述, ψ = E κ ( F , 0 ) \psi=\mathcal{E}_\kappa(F,0) ψ=Eκ(F,0)

由于势能函数是偶函数,因此波函数必须要么是奇函数要么是偶函数(证明待补充)。这样做的好处是我们只用考虑 x ≥ 0 x\ge 0 x0 x ≤ 0 x\le0 x0。原来刚才用 T \mathcal{T} T符号而不是 C \mathcal{C} C符号表示的原因也在于利用奇偶性质。现在先假设是偶函数解一下。2022-10-30 01:12:23

加入 sin ⁡ \sin sin成分会让函数变得不偶,因此 C = 0 C=0 C=0。在 a a a处,由标准边界条件1,有 F e − κ a = D cos ⁡ ( l a ) Fe^{-\kappa a}=D\cos(la) Feκa=Dcos(la)。由标准边界条件2, − κ F e − κ a = − D l sin ⁡ ( l a ) -\kappa Fe^{-\kappa a}=-Dl\sin(la) κFeκa=Dlsin(la)。用上式除以上上式,得 κ = l tan ⁡ ( l a ) \kappa=l\tan(la) κ=ltan(la)

κ \kappa κ l l l都是正比于 E E E的函数,因此这个方程实际上是解出 E E E的方程。

z = l a z=la z=la,则 tan ⁡ ( z ) = κ l = − 2 m E 2 m ( E + V 0 ) = − E E + V 0 = V 0 E + V 0 − 1 = 2 m V 0 ( l ℏ ) 2 − 1 = 2 m a 2 V 0 ℏ 2 z 2 − 1 \tan(z)=\frac{\kappa}{l}=\frac{\sqrt{-2mE}}{\sqrt{2m(E+V_0)}}=\sqrt{\frac{-E}{E+V_0}}=\sqrt{\frac{V_0}{E+V_0}-1}=\sqrt{\frac{2mV_0}{(l\hbar)^2}-1}=\sqrt{\frac{\frac{2ma^2V_0}{\hbar^2}} {z^2}-1} tan(z)=lκ=2m(E+V0) 2mE =E+V0E =E+V0V01 =(l)22mV01 =z222ma2V01

z 0 = a ℏ 2 m V 0 z_0=\frac{a}{\hbar}\sqrt{2mV_0} z0=a2mV0 ,则 tan ⁡ ( z ) = ( z 0 z ) 2 − 1 \tan(z)=\sqrt{(\frac{z_0}{z})^2-1} tan(z)=(zz0)21

解这个超越(transcendental)方程就可以求出能量,用计算机解决。画图表示,右边是一个在0处取无穷大并且逐渐下降,在x>0某点与x轴相交的函数。

有很多个解,分别记为 z n , n = 1 , 2 , 3 , … z_n,n=1,2,3,\dots zn,n=1,2,3,,对应的能量为 E n E_n En。其关系为 E n = l n 2 ℏ 2 2 m − V 0 = z n 2 ℏ 2 2 m a 2 − V 0 E_n=\frac{l_n^2\hbar^2}{2m}-V_0=\frac{z_n^2\hbar^2}{2ma^2}-V_0 En=2mln22V0=2ma2zn22V0

我们考虑两个有趣的特殊情况。

  1. 很宽很深的势阱。那么右手边函数很高,与 tan ⁡ z \tan z tanz总是在接近 n π / 2 = z n n\pi/2= z_n /2=zn的地方相交。即 E n + V 0 ≈ n 2 π 2 ℏ 2 8 m a 2 E_n+V_0\approx\frac{n^2\pi^2\hbar^2}{8ma^2} En+V08ma2n2π22右边实际上就是对应的无限深势阱的能量(对应章节里宽度是 a a a,而本章节宽度是 2 a 2a 2a)。对于任意有限的 V V V态的数量都是有限的。
  2. 很浅很窄的势阱。那么右手边函数很贴近0,并且零点越来越左移。最后只剩下一个态(但始终至少有1个态)。

再考虑奇函数情况 F e − κ a = C sin ⁡ ( l a ) Fe^{-\kappa a}=C\sin(la) Feκa=Csin(la) − κ F e − κ a = l C cos ⁡ ( l a ) -\kappa Fe^{-\kappa a}=lC\cos(la) κFeκa=lCcos(la)。上式除以上上式得 κ = − l cot ⁡ ( l a ) \kappa=-l\cot(la) κ=

你可能感兴趣的:(物理,线性代数,机器学习,矩阵)