量子力学教程(全)
按顺序阅读。
第一章
波
波函数E(x,t)=E0ei(kx−ωt)E(x,t)=E_0e^{i(kx-\omega t)}E(x,t)=E0ei(kx−ωt),记为Daw(A=E0,k=k,ω=ω)Daw(A=E_0,k=k,\omega=\omega)Daw(A=E0,k=k,ω=ω),或E0Dkω(x,t)E_0\mathcal{D}_k^\omega(x,t)E0Dkω(x,t)。其中ω\omegaω为角速度,kkk为波数,即每2π2\pi2π长度里有多少个波长。
Dω(t)=e−iωt\mathcal{D}^\omega(t)=e^{-i\omega t}Dω(t)=e−iωt
Dk(x)=eikx\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)=DE/ℏ(t)=Dω(t)f(t)=\mathcal{D}^{E/\hbar}(t)=\mathcal{D}^\omega(t)f(t)=DE/ℏ(t)=Dω(t)。
自由粒子
如果是自由粒子或平面波,则ψ(x)=ADk(x)\psi(x)=A\mathcal{D}_k(x)ψ(x)=ADk(x),其中AAA是归一化常数。
于是Ψ(x,t)=ADk(x)Dω(t)=ADkω(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∣ψ⟩,不再提到xxx,即与坐标系选取无关。
一个封闭的狄拉克符号表示一个复数,如C=⟨χ∣ψ⟩C=\langle\chi|\psi\rangleC=⟨χ∣ψ⟩。含义是内积,也即C=(χ,ψ)C=(\chi,\psi)C=(χ,ψ)。对于向量可以看做矩阵乘法,对于函数可以看做函数内积。
在中间多乘一个矩阵可以表示为⟨φ∣A∣ψ⟩\langle\varphi|A|\psi\rangle⟨φ∣A∣ψ⟩。
投影算符Pφ=∣φ⟩⟨φ∣P_\varphi=|\varphi\rangle\langle\varphi|Pφ=∣φ⟩⟨φ∣。
算符AAA的厄米共轭(即取转置,再每个元素取共轭)记为A†A^\daggerA†。如果A∣α⟩=∣β⟩A|\alpha\rangle=|\beta\rangleA∣α⟩=∣β⟩,那么⟨β∣A†=⟨α∣\langle\beta|A^\dagger=\langle\alpha|⟨β∣A†=⟨α∣。
(∣u⟩⟨v∣)†=∣v⟩⟨u∣(|u\rangle\langle v|)^\dagger=|v\rangle\langle u|(∣u⟩⟨v∣)†=∣v⟩⟨u∣
如果A=A†A=A^\daggerA=A†,那么AAA是厄米算符。
交换性质:⟨φ∣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†符号和矩阵论里的共轭转置HHH上标是一回事。
克罗内克符号 正交归一关系
克罗内克符号:Kronecker Delta。
离散型:δij={1,i=j0,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=0or x→00,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 x→0otherwise
正交基{ui}\{u_i\}{ui}满足⟨ui∣uj⟩=δij\langle u_i|u_j\rangle=\delta_{ij}⟨ui∣uj⟩=δij。
连续基满足(ωα,ωα′)=δ(α−α′)(\omega_\alpha,\omega_{\alpha'})=\delta(\alpha-\alpha')(ωα,ωα′)=δ(α−α′)。
∣ψ⟩|\psi\rangle∣ψ⟩在一组正交基上投影表示为∣ψ⟩=∑iPui∣ψ⟩|\psi\rangle=\sum_iP_{u_i}|\psi\rangle∣ψ⟩=∑iPui∣ψ⟩
必然有∑iPui=I\sum_iP_{u_i}=I∑iPui=I。
连续基对应性质为∫Pωαdα=1\int P_{\omega_\alpha}d\alpha=1∫Pωαdα=1。
如果矩阵S=UTS=UTS=UT,那么S†=(UT)†=T†U†S^\dagger=(UT)^\dagger=T^\dagger U^\daggerS†=(UT)†=T†U†。设eee表示单位正交基。那么Sij=⟨ei∣S∣ej⟩=⟨ei∣U∣T∣ej⟩=⟨u(row=i)∣t(column=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) \rangleSij=⟨ei∣S∣ej⟩=⟨ei∣U∣T∣ej⟩=⟨u(row=i)∣t(column=j)⟩。(S†)ij=⟨ei∣T†∣U†∣ej⟩=⟨t†(row=i)∣u†(column=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=⟨ei∣T†∣U†∣ej⟩=⟨t†(row=i)∣u†(column=j)⟩。如果记⟨ui∣=⟨u(row=i)∣\langle u_i|=\langle u(row=i)|⟨ui∣=⟨u(row=i)∣,那么自然∣ui⟩=∣u†(column=i)⟩|u_i\rangle=|u^\dagger(column=i)\rangle∣ui⟩=∣u†(column=i)⟩。即,若Sij=⟨ui∣tj⟩S_{ij}=\langle u_i|t_j\rangleSij=⟨ui∣tj⟩,则(S†)ij=⟨ti∣uj⟩(S^\dagger)_{ij}=\langle t_i|u_j\rangle(S†)ij=⟨ti∣uj⟩。另一种说明:(S†)ij=(Sji)∗=⟨uj∣ti⟩∗=⟨ti∣uj⟩(S^\dagger)_{ij}=(S_{ji})^*=\langle u_j|t_i\rangle^*=\langle t_i|u_j\rangle(S†)ij=(Sji)∗=⟨uj∣ti⟩∗=⟨ti∣uj⟩。
如果UUU和TTT都是酉矩阵(幺正矩阵),那么SSS也是酉矩阵。我们说{∣ui⟩}\{|u_i\rangle\}{∣ui⟩}或{∣ti⟩}\{|t_i\rangle\}{∣ti⟩}是表象(representation)。SSS称为变换矩阵。
本征值方程 简并度
若A∣ψ⟩=λ∣ψ⟩A|\psi\rangle=\lambda|\psi\rangleA∣ψ⟩=λ∣ψ⟩,则称这个方程为AAA的本征值方程,λ\lambdaλ为其一个本征值(eigenvalue),∣ψ⟩|\psi\rangle∣ψ⟩为其对应的本征矢(eigenvector)。
对于方程A∣ψ⟩=λ∣ψ⟩A|\psi\rangle=\lambda|\psi\rangleA∣ψ⟩=λ∣ψ⟩,其中AAA已知,如果λ\lambdaλ是一个本征值,那么此时方程的解空间称为本征值λ\lambdaλ的本征子空间(eigensubspace)。如果本征子空间里只有一个向量,即是一维的,那么称本征值λ\lambdaλ是非简并的;如果大于一维,那么是简并的(degenerate),子空间维数称为简并度(degeneracy),也就是几何重数(geometric multiplicity)。
如果AAA是厄米的,那么本征值的简并度总是等于重数,即几何重数等于代数重数。
如果AAA是厄米的,那么AAA的本征值都是实数,且不同本征值对应的本征矢都正交。
波函数与归一化 叠加原理
某一时刻,一个粒子的状态用Ψ(q)\Psi(q)Ψ(q)表示,其中qqq是所有可能用到的坐标,这个函数称为波函数。粒子在体积元dqdqdq内的概率为Ψ(q)(Ψ(q))∗dq=∣Ψ(q)∣2dq\Psi(q)(\Psi(q))^*dq=|\Psi(q)|^2dqΨ(q)(Ψ(q))∗dq=∣Ψ(q)∣2dq。以坐标为波函数的参数的量子力学称为坐标表象下的量子力学,或波动力学。
我们自然有归一化条件:∫Ψ(q)(Ψ(q))∗dq=1\int \Psi(q)(\Psi(q))^*dq =1∫Ψ(q)(Ψ(q))∗dq=1。
波函数可以是几个波函数的和。即Ψ=∑iciΨi\Psi=\sum_ic_i\Psi_iΨ=∑iciΨi,其中cic_ici是复数。但这还不能称为叠加原理(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。∣ci∣2|c_i|^2∣ci∣2表示一种无法解释的概率,称为不解释#1概率。也有积分形式Ψ=∫cnΨndn\Psi=\int c_n\Psi_n dnΨ=∫cnΨndn。由于这些状态互斥,因此∑i∣ci∣2=1\sum_i|c_i|^2=1∑i∣ci∣2=1。
对于物理量fff,如果它可能有nnn个取值fif_ifi,称这些值为fff的笨征值(baconvalue)。fff取fif_ifi这个事件用已经归一化的波函数Ψi\Psi_iΨi表示。那么fff的期望fˉ=∑ifi∣ci∣2\bar f=\sum_if_i|c_i|^2fˉ=∑ifi∣ci∣2。(你信吗?)
叠加原理用狄拉克符号可以表示为∣Ψ⟩=∑iai∣Ψi⟩|\Psi\rangle=\sum_i a_i|\Psi_i\rangle∣Ψ⟩=∑iai∣Ψi⟩,将∣Ψi⟩|\Psi_i \rangle∣Ψi⟩记为∣i⟩|i\rangle∣i⟩,则∣Ψ⟩=∑iai∣i⟩|\Psi\rangle=\sum_ia_i|i\rangle∣Ψ⟩=∑iai∣i⟩。
正交条件可以表示为⟨Ψi∣Ψj⟩=δij\langle \Psi_i|\Psi_j\rangle=\delta_{ij}⟨Ψi∣Ψj⟩=δij,即⟨i∣j⟩=δij\langle i|j\rangle=\delta_{ij}⟨i∣j⟩=δij。
一个本征态与波函数的内积就是幅前系数:⟨i∣Ψ⟩=⟨i∣∑jaj∣j⟩=∑jaj⟨i∣j⟩=ai\langle i|\Psi\rangle=\langle i|\sum_ja_j|j\rangle=\sum_ja_j\langle i|j\rangle=a_i⟨i∣Ψ⟩=⟨i∣∑jaj∣j⟩=∑jaj⟨i∣j⟩=ai。也就是说,aia_iai就是波函数到这个本征态的投影长,和投影到一组单位正交基是一个道理。
上式两边取共轭可得⟨Ψ∣i⟩=ai∗\langle \Psi|i\rangle=a_i^*⟨Ψ∣i⟩=ai∗。
那么fff的期望fˉ=∑ifiaiai∗=∑ifiai⟨Ψ∣i⟩=⟨Ψ∣∑ifiai∣i⟩\bar f=\sum_if_ia_ia_i^*=\sum_if_ia_i\langle\Psi|i\rangle=\langle\Psi|\sum_if_ia_i|i\ranglefˉ=∑ifiaiai∗=∑ifiai⟨Ψ∣i⟩=⟨Ψ∣∑ifiai∣i⟩。
物理量对应的算符
物理量fff对应的算符为f^\hat ff^。算符的阶数与fff的可能取值数相同。算符作用于波函数时,会将其每个带系本征态乘上fff的该取值再求和。即f^∣Ψ⟩=∑ifiai∣i⟩\hat f|\Psi\rangle=\sum_if_ia_i|i\ranglef^∣Ψ⟩=∑ifiai∣i⟩。
例如,若Ψ=(2222)\Psi=\begin{pmatrix}\frac{\sqrt 2}{2}\\\frac{\sqrt 2}{2}\end{pmatrix}Ψ=(2222),取∣1⟩=(10)|1\rangle=\begin{pmatrix}1\\0\end{pmatrix}∣1⟩=(10),∣2⟩=(01)|2\rangle=\begin{pmatrix}0\\1\end{pmatrix}∣2⟩=(01)。假设fff可以取333和444,概率用上述波函数表示。那么f^\hat ff^就是(34)\begin{pmatrix}3&\\&4\end{pmatrix}(34)。
那么fˉ\bar ffˉ可以表示为⟨Ψ∣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⟩=fi∣i⟩\hat f|i\rangle=f_i|i\ranglef^∣i⟩=fi∣i⟩,其中fif_ifi表示fff的这个本征态对应的取值。于是我们发现,矩阵表示下,本征态就是该物理量算符的本征矢,对应本征值就是对应笨征值。
本征态也可以是本征函数(eigenfunction),取决于你认为本征态是什么。当然,我们混淆这三个概念也无妨。
双下标等于两边乘写法:fij=⟨i∣f^∣j⟩=⟨i∣fj∣j⟩=fj⟨i∣j⟩=fjδijf_{ij}=\langle i|\hat f|j\rangle=\langle i|f_j|j\rangle=f_j\langle i|j\rangle=f_j\delta_{ij}fij=⟨i∣f^∣j⟩=⟨i∣fj∣j⟩=fj⟨i∣j⟩=fjδij。称双下标等于两边乘写法表示的矩阵[fij][f_{ij}][fij]为跃迁矩阵,是对角矩阵。上面的例子只是刚好f^\hat ff^是跃迁矩阵,若本征态不取单位矩阵的列那么就很可能不是。
物理量与表象 观察算符 得某一结果的概率
物理量fff的本征值、本征态与波函数是确定的,但是这些本征态可以张成多种基矢,这些基矢又张成了同样的波函数,这种现象称为fff有不同的表象。也就是说本征态是fff的其中一个表象。注意,表象的基矢也必须是正交归一的。
设有两组连续表象基矢{∣u⟩}\{|u\rangle\}{∣u⟩},{∣t⟩}\{|t\rangle\}{∣t⟩}。∣Ψ⟩=∫du∣u⟩⟨u∣Ψ⟩|\Psi\rangle=\int du|u\rangle\langle u|\Psi\rangle∣Ψ⟩=∫du∣u⟩⟨u∣Ψ⟩,∣Ψ⟩=∫dt∣t⟩⟨t∣Ψ⟩|\Psi\rangle=\int dt|t\rangle\langle t|\Psi\rangle∣Ψ⟩=∫dt∣t⟩⟨t∣Ψ⟩。显然,利用投影算符的归一性,⟨t∣=∫du⟨t∣u⟩⟨u∣\lang t|=\int du\lang t|u\rang\lang u|⟨t∣=∫du⟨t∣u⟩⟨u∣。
可以记⟨a∣Ψ⟩\lang a|\Psi\rang⟨a∣Ψ⟩为Ψ(a)\Psi(a)Ψ(a)。
如果厄米算符AAA的所有线性无关的本征矢能构成一组态空间(state space)的基,那么算符AAA是观察算符(observation operator)。这组基也是一个表象,一般称为AAA表象。
上述的“物理量对应的算符”,都是观察算符。 也就是说本征子空间肯定与态空间同张,物理量算符肯定是厄米算符。
一个本征值可能对应多个本征矢,即具有简并度,求结果为该本征值的概率时不能忘记简并度。记本征值fif_ifi的简并度为gig_igi,fif_ifi对应的归一化本征矢分别记为∣ψi1⟩,∣ψi2⟩,∣ψi3⟩,…|\psi_i^1\rang,|\psi_i^2\rang,|\psi_i^3\rang,\dots∣ψi1⟩,∣ψi2⟩,∣ψi3⟩,…。则测量fff得结果fif_ifi的概率为pi=∑j=1gi∣⟨ψij∣Ψ⟩∣2p_i=\sum_{j=1}^{g_i}|\lang\psi_i^j|\Psi\rang|^2pi=∑j=1gi∣⟨ψij∣Ψ⟩∣2。
第一章小结
测量一个物理量fff可能得到多个结果。这些结果称为本征值{fi}\{f_i\}{fi}。测量一个物理量的结果用波函数Ψ\PsiΨ表示。测量一个物理量对应一个观察算符,这个观察算符的本征值就是所有可能结果,观察算符的本征态{Ψi}\{\Psi_i\}{Ψi}就是该物理量的本征态。波函数可以由一组单位正交基{ui}\{u_i\}{ui}张成,这组基称为一个表象。表象必然可以由本征子空间张成。波函数与一个本征值对应的所有本征态的内积的模长平方和就是测量得这个值的概率。
第二章
薛定谔方程 自动归一化
一个粒子的波函数为ψ(x,t)\psi(x,t)ψ(x,t)。势能为V(x,t)V(x,t)V(x,t),是实函数。则其波函数满足薛定谔方程
−ℏ22m∂2∂x2ψ+Vψ=iℏ∂∂tψ-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi+V\psi=i\hbar\frac{\partial}{\partial t}\psi−2mℏ2∂x2∂2ψ+Vψ=iℏ∂t∂ψ
众所周知,归一化的波函数才是有意义的。如果我们在t=0t=0t=0时求出了ψ(x)\psi(x)ψ(x)的表达式,必须考虑在之后时间其是否仍然维持归一化。
对薛定谔方程两边取共轭得−ℏ22m∂2∂x2ψ∗+Vψ∗=−iℏ∂∂tψ∗-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi^*+V\psi^*=-i\hbar\frac{\partial}{\partial t}\psi^*−2mℏ2∂x2∂2ψ∗+Vψ∗=−iℏ∂t∂ψ∗
波函数的模长平方对时间的导数∂∂t∣ψ∣2=∂ψ∂tψ∗+∂ψ∗∂tψ=(iℏ2m∇2ψ−iVℏψ)ψ∗+(−iℏ2m∇2ψ∗+iVℏψ∗)ψ=iℏ2m((∇2ψ)ψ∗−(∇2ψ∗)ψ)=iℏ2m∇((∇ψ)ψ∗−ψ(∇ψ∗))\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∂ψ∗ψ=(2miℏ∇2ψ−iℏVψ)ψ∗+(−2miℏ∇2ψ∗+ℏiVψ∗)ψ=2miℏ((∇2ψ)ψ∗−(∇2ψ∗)ψ)=2miℏ∇((∇ψ)ψ∗−ψ(∇ψ∗))
在整个空间上积分有ddt∫dx∣ψ∣2=∫dxiℏ2m∇((∇ψ)ψ∗−ψ(∇ψ∗))=iℏ2m((∇ψ)ψ∗−ψ(∇ψ∗))∣−∞+∞\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}dtd∫dx∣ψ∣2=∫dx2miℏ∇((∇ψ)ψ∗−ψ(∇ψ∗))=2miℏ((∇ψ)ψ∗−ψ(∇ψ∗))∣−∞+∞。
由于波函数在正负无穷处趋于000(否则不能归一化),因此整个空间的概率和随时间不变。即若在t=0t=0t=0归一化,则以后都自动归一化。
(∂2[ψ∣∂x2)∣ψ∗]=∂∂x(∂[ψ∣∂x∣ψ∗])(\frac{\partial^2 [\psi|}{\partial x^2})|\psi^*]=\frac{\partial}{\partial x}(\frac{\partial[\psi|}{\partial x}|\psi^*])(∂x2∂2[ψ∣)∣ψ∗]=∂x∂(∂x∂[ψ∣∣ψ∗])
动量
xxx是一个粒子的一维坐标。其波函数为Ψ(x,t)\Psi(x,t)Ψ(x,t)。则xxx的期望是
∫dx⟨Ψ∣x∣Ψ⟩\int dx\lang\Psi|x|\Psi\rang∫dx⟨Ψ∣x∣Ψ⟩
该期望随时间的导数是(使用分部积分,原部肯定取000,直接写导部)d⟨x⟩dt=∫xdx∂∣ψ∣2∂t=∫xdxiℏ2m∂∂x(∂[ψ∣∂x∣ψ∗])=−iℏ2m∫dx∂[ψ∣∂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^*]dtd⟨x⟩=∫xdx∂t∂∣ψ∣2=∫xdx2miℏ∂x∂(∂x∂[ψ∣∣ψ∗])=−2miℏ∫dx∂x∂[ψ∣∣ψ∗]
再分部积分一次得(明天补充细节)d⟨x⟩dt=−iℏm∫∂ψ∂xψ∗CINVIKdx\frac{d\lang x\rang}{dt}=-\frac{i\hbar}{m}\int \frac{\partial \psi}{\partial x}\psi^*\text{CINVIK}dxdtd⟨x⟩=−miℏ∫∂x∂ψψ∗CINVIKdx
速度的期望就是位置的期望对时间的导数:vˉ=dxˉdt\bar v=\frac{d\bar x}{dt}vˉ=dtdxˉ
所以⟨p⟩=md⟨x⟩dt=−iℏ∫∂ψ∂xψ∗CINVIKdx=∫ψ∗(−iℏ)∂∂xψdx\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⟩=mdtd⟨x⟩=−iℏ∫∂x∂ψψ∗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,−iℏ∂x∂)⟩
自由粒子
如果V≡0V\equiv0V≡0,那么这个粒子是自由粒子(free particle)。
代入薛定谔方程,得知微分方程必有此形式解:ψ(x,t)=ADkω(x,t)\psi(x,t)=A\mathcal{D}_k^\omega(x,t)ψ(x,t)=ADkω(x,t)
由∂Dkω∂t=−iωDkω\frac{\partial\mathcal{D}_k^\omega}{\partial t}=-i\omega\mathcal{D}_k^\omega∂t∂Dkω=−iωDkω∂2Dkω∂x2=−k2Dkω\frac{\partial^2\mathcal{D}_k^\omega}{\partial x^2}=-k^2\mathcal{D}_k^\omega∂x2∂2Dkω=−k2Dkω因此ℏk22m=ω\frac{\hbar k^2}{2m}=\omega2mℏk2=ω
海森堡不确定性原理
德布罗意公式(de Broglie formula)给出,粒子的动量为p=hλ=ℏkp=\frac{h}{\lambda}=\frac{\hbar}{k}p=λh=kℏ
设Δx\Delta xΔx和Δp\Delta pΔp是位置xxx和动量ppp的标准差,那么海森堡不确定性原理(Heisenberg uncertainty principle)给出ΔxΔp≥ℏ2\Delta x\Delta p\ge\frac{\hbar}{2}ΔxΔp≥2ℏ
薛定谔方程的定态解
解形式Ψ(x,t)=ψ(x)φ(t)\Psi(x,t)=\psi(x)\varphi(t)Ψ(x,t)=ψ(x)φ(t)称为薛定谔方程的定态解。
假设势能VVV只与xxx有关。则将定态解代入薛定谔方程得iℏψ(x)dφ(t)dt=V(x)ψ(x)φ(t)−ℏ22mφ(t)d2ψ(x)dx2i\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)−2mℏ2φ(t)dx2d2ψ(x)
两边同时除以波函数,得iℏφ(t)dφ(t)dt=V(x)−ℏ22mψ(x)d2ψ(x)dx2\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)iℏdtdφ(t)=V(x)−2mψ(x)ℏ2dx2d2ψ(x)
该等式左边只与ttt有关,而右边只与xxx有关,这说明等式两边都等于一个常数。让此常数为ℏω\hbar\omegaℏω,代入左边解得φ(t)=ADω(t)\varphi(t)=A\mathcal{D}^\omega(t)φ(t)=ADω(t)
把AAA放入ψ(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\omegaE=ℏω的态,即能量的本征态。这与经典情形相符:势能不随时间变化时,能量应当是定值。
令E=ℏωE=\hbar\omegaE=ℏω,哈密顿算符H^=−ℏ22m∇2+V(r)\hat H=-\frac{\hbar^2}{2m}\nabla^2+V(\bm{r})H^=−2mℏ2∇2+V(r),则薛定谔方程代入定态解两边约去φ(t)\varphi(t)φ(t)给出H^ψ(r)=Eψ(r)\hat H\psi(\bm{r})=E\psi(\bm{r})H^ψ(r)=Eψ(r)哈密顿算符可以看做经典能量在坐标表象下的表示。
哈密顿算符的期望是EEE,即⟨H^⟩=E\lang \hat H\rang=E⟨H^⟩=E。
由H^2ψ=H^(H^ψ)=EH^ψ=E2ψ\hat H^2\psi=\hat H(\hat H\psi)=E\hat H\psi=E^2\psiH^2ψ=H^(H^ψ)=EH^ψ=E2ψ,有HHH的标准差为ΔH=⟨H2⟩−⟨H⟩2=E2−E2=0\Delta H=\sqrt {\lang H^2\rang-\lang H\rang^2}=\sqrt{E^2-E^2}=0ΔH=⟨H2⟩−⟨H⟩2=E2−E2=0。
上述方程就是H^\hat HH^的本征值方程。那么能量的所有可能取值就是H^\hat HH^的本征值。
定态叠加 EEE与VVV的关系
EEE有多个可能取值,对应多个H^\hat HH^的本征函数,记为EnE_nEn和ψn(x)\psi_n(x)ψn(x),n=1,2,3,…n=1,2,3,\dotsn=1,2,3,…。则H^ψn=Enψn\hat H\psi_n=E_n\psi_nH^ψn=Enψn。对应波函数为Ψn(x,t)=ψn(x)Dω(t)=ψn(x)DE/ℏ(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)。
这是线性方程,所以几个波函数叠加的波函数也是方程的解。即Ψ=∑ncnΨn\Psi=\sum_nc_n\Psi_nΨ=∑ncnΨn。特别地,Ψ(x,0)=∑ncnψ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。
能量是动能和势能之和,因此能量必须大于VVV的最小值,否则波函数不能归一化。
一维无限深方势阱(infinite square well)
设V(x)V(x)V(x)在[0,a][0,a][0,a]处为000,其它地方为∞\infty∞。
在势能无限高的地方,粒子不可能到达,此处ψ(x)=0\psi(x)=0ψ(x)=0。
在势阱内,有−ℏ22md2dx2ψ(x)=Eψ(x)-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\psi(x)=E\psi(x)−2mℏ2dx2d2ψ(x)=Eψ(x)令k=2mEℏk=\frac{\sqrt{2mE}}{\hbar}k=ℏ2mE,则d2dx2ψ(x)=−k2ψ(x)\frac{d^2}{dx^2}\psi(x)=-k^2\psi(x)dx2d2ψ(x)=−k2ψ(x)
这是谐振子(harmonic oscillator)方程,其解为ψ(x)=Asin(kx)+Bcos(kx)\psi(x)=A\sin(kx)+B\cos(kx)ψ(x)=Asin(kx)+Bcos(kx),记为Tk(A,B)(x)\mathcal{T}_k(A,B)(x)Tk(A,B)(x),也可以表示为C1eikx+C2e−ikxC_1e^{ikx}+C_2e^{-ikx}C1eikx+C2e−ikx,记为Ck(C1,C2)(x)\mathcal{C}_k(C_1,C_2)(x)Ck(C1,C2)(x)。转换关系为Ck(C1,C2)=C1cos(kx)+iC1sin(kx)+C2cos(−kx)+iC2sin(−kx)=(C1+C2)cos(kx)+(iC1−iC2)sin(kx)=Tk(iC1−iC2,C1+C2)\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)+(iC1−iC2)sin(kx)=Tk(iC1−iC2,C1+C2),反过来是Tk(A,B)=Ck(B−iA2,B+iA2)=iA2Ck(−1,1)+B2Ck(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(2B−iA,2B+iA)=2iACk(−1,1)+2BCk(1,1)。
一般,波函数及其导数是连续的。那么有ψ(0)=ψ(a)=0\psi(0)=\psi(a)=0ψ(0)=ψ(a)=0。于是B=0B=0B=0,即ψ(x)=Tk(A,0)\psi(x)=\mathcal{T}_k(A,0)ψ(x)=Tk(A,0),进而sin(ka)=0\sin(ka)=0sin(ka)=0。于是解出kkk,其中kkk取000是没有意义的:kn=nπa(n=1,2,3,…)k_n=\frac{n\pi}{a}(n=1,2,3,\dots)kn=anπ(n=1,2,3,…)对应的能量为En=ℏ2kn22m=n2π2ℏ22ma2E_n=\frac{\hbar^2k_n^2}{2m}=\frac{n^2\pi^2\hbar^2}{2ma^2}En=2mℏ2kn2=2ma2n2π2ℏ2。
通过归一化可以推得∣A∣=2a|A|=\sqrt{\frac{2}{a}}∣A∣=a2,这里我们直接取A=2aA=\frac{2}{a}A=a2,因为不同的相位也没有意义。于是ψn(x)=2asin(nπxa)\psi_n(x)=\sqrt\frac{2}{a}\sin(\frac{n\pi x}{a})ψn(x)=a2sin(anπx)。
即该解该方程得到无数个解。波函数的图像可以类比波的图像,波的能量与波的频率的平方成正比。因此ψ1\psi_1ψ1能量最低,称为基态(ground state)。ψn\psi_nψn的能量是ψ1\psi_1ψ1的n2n^2n2倍,对于n>1n>1n>1,称为激发态(excited state)。
这些波函数有一些性质:
- 相对于势阱中心,ψ1(x)\psi_1(x)ψ1(x)是偶函数,ψ2(x)\psi_2(x)ψ2(x)是奇函数,ψ3(x)\psi_3(x)ψ3(x)是偶函数,依此类推。
- 称波函数与xxx轴的交点为node(节点,结点)。具体地,ψn(x)\psi_n(x)ψn(x)有n−1n-1n−1个node。
- 它们相互正交(mutually orthogonal):∫(ψm(x))∗ψn(x)dx=δmn\int (\psi_m(x))^*\psi_n(x)dx=\delta_{mn}∫(ψm(x))∗ψn(x)dx=δmn,即⟨ψm∣ψn⟩=δmn\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∞cnψn(x)f(x)=\sum_{n=1}^\infty c_n\psi_n(x)f(x)=∑n=1∞cnψn(x)。参考傅里叶变换。用傅里叶技巧(Fourier’s trick)可以导出cn=∫ψn(x)∗f(x)dxc_n=\int \psi_n(x)^*f(x)dxcn=∫ψn(x)∗f(x)dx,也就是cn=⟨ψn∣f⟩c_n=\lang\psi_n|f\rangcn=⟨ψn∣f⟩。
如果势(potential)是对称的,那么奇偶函数性质成立。节点性质总是成立,不管势是什么形状。正交性质也对很多成立。完备性质也对很多成立。
如果我们知道了Ψ(x,0)\Psi(x,0)Ψ(x,0),那么我们可以先求出cnc_ncn,然后立即得到Ψ(x,t)\Psi(x,t)Ψ(x,t)。对于不是一维无限深方的其它势阱,也可以用此方法,只是ψ\psiψ系不同。
∣cn∣2|c_n|^2∣cn∣2是能量EEE取EnE_nEn的概率,必然满足∑∣cn∣2=1\sum|c_n|^2=1∑∣cn∣2=1。能量的期望值与时间无关,这是能量守恒在量子力学中的体现。
谐振子 升降算符
谐振子的势能是V(x)=12mω2x2=12kx2V(x)=\frac{1}{2}m\omega^2x^2=\frac{1}{2}kx^2V(x)=21mω2x2=21kx2。
代数解法
现在我们看到H=p22m+V=12m[p2+(mωx)2]H=\frac{p^2}{2m}+V=\frac{1}{2m}[ p^2+(m\omega x)^2]H=2mp2+V=2m1[p2+(mωx)2]。
我们希望用复数平方差公式把中括号里的东西展开,即u2+v2=(−iu+v)(iu+v)u^2+v^2=(-iu+v)(iu+v)u2+v2=(−iu+v)(iu+v)。但是现在ppp和xxx不是对易的(commute),即px≠xppx\ne xppx=xp。
令a±=12ℏmω(∓ip+mωx)a_\pm=\frac{1}{\sqrt{2\hbar m\omega}}(\mp ip+m\omega x)a±=2ℏmω1(∓ip+mωx)
那么试算a−a+=12ℏmω[p2+(mωx)2+imωpx−imωxp]=12ℏmω[p2+(mωx)2−imω[x,p]]=12ℏmω[p2+(mωx)2]−i2ℏ[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}a−a+=2ℏmω1[p2+(mωx)2+imωpx−imωxp]=2ℏmω1[p2+(mωx)2−imω[x,p]]=2ℏmω1[p2+(mωx)2]−2ℏi[x,p]其中方括号是对易子(commutator)。
易知[x,p]=iℏ[x,p]=i\hbar[x,p]=iℏ,这被称为canonical commutation relation。
于是a−a+=Hℏω+12a_-a_+=\frac{H}{\hbar \omega}+\frac{1}{2}a−a+=ℏωH+21让HHH到左边,有H=ℏω(a−a+−12)H=\hbar\omega(a_-a_+-\frac{1}{2})H=ℏω(a−a+−21)
计算得a+a−=Hℏω−12a_+a_-=\frac{H}{\hbar \omega}-\frac{1}{2}a+a−=ℏωH−21,于是[a−,a+]=1[a_-,a_+]=1[a−,a+]=1。也有H=ℏω(a+a−+12)H=\hbar\omega(a_+a_-+\frac{1}{2})H=ℏω(a+a−+21)。
我们现在满足Hψ=EψH\psi=E\psiHψ=Eψ。现在计算H(a+ψ)=ℏω(a+a−+12)(a+ψ)=ℏω(a+a−a++12a+)ψ=ℏωa+(a−a++12)ψ=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+a−a++21a+)ψ=ℏωa+(a−a++21)ψ=a+(H+ℏω)ψ=(E+ℏω)(a+ψ)注意升降算符变成H以后就可以到处对易了
同样有H(a−ψ)=(E−ℏω)(a−ψ)H(a_-\psi)=(E-\hbar\omega)(a_-\psi)H(a−ψ)=(E−ℏω)(a−ψ)。
因此我们把a±a_\pma±叫做升降算符(ladder operators),分别是上升算符(raising operator)和下降算符(lowering operator)。
注意,我们可以一直将下降算符作用于波函数,但是一直下去可能会导致波函数不能归一化,实际不存在。必然有一个态ψ0\psi_0ψ0满足a−ψ0=0a_-\psi_0=0a−ψ0=0
尝试解出ψ0\psi_0ψ0。方程是12ℏmω(i(−iℏddx)+mωx)ψ0=0\frac{1}{2\hbar m\omega}(i (-i\hbar\frac{d}{dx})+m\omega x)\psi_0=02ℏmω1(i(−iℏdxd)+mωx)ψ0=0ℏdψ0ψ0=−mωxdxℏlnψ0=−12mωx2+constantψ0=Ae−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ψ0ℏlnψ0ψ0=−mωxdx=−21mωx2+constant=Ae−2ℏmω
归一化给出∣A∣2=mωπℏ|A|^2=\sqrt{\frac{m\omega}{\pi\hbar}}∣A∣2=πℏmω,因此ψ0=(mωπℏ)14e−mω2ℏx2\psi_0=(\frac{m\omega}{\pi\hbar})^\frac{1}{4}e^{-\frac{m\omega}{2\hbar}x^2}ψ0=(πℏmω)41e−2ℏmωx2
代入薛定谔方程得对应能量E0=12ℏωE_0=\frac{1}{2}\hbar\omegaE0=21ℏω于是En=(n+12)ℏωE_n=(n+\frac{1}{2})\hbar\omegaEn=(n+21)ℏω
将ψn\psi_nψn作用于HHH和EnE_nEn右边得ℏω(a+a−+12)ψn=(n+12)ℏωψ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−=na_+a_-=na+a−=n同样可得a−a+=n+1a_-a_+=n+1a−a+=n+1
通过升降关系可以得到其它所有态。注意升降算符作用于邻态后还要归一化才是本态。即ψn+1=An+1a+ψn\psi_{n+1}=A_{n+1}a_+\psi_nψn+1=An+1a+ψn,一般写成cnψn+1=a+ψnc_n\psi_{n+1}=a_+\psi_ncnψn+1=a+ψndnψn−1=a−ψnd_n\psi_{n-1}=a_-\psi_ndnψn−1=a−ψn
对于连续函数fff和ggg,升降算符必然满足性质:(注意升降算符都在函数左边) 证明过两天补充
- ⟨f∣(a+g)⟩=⟨(a−f)∣g⟩\lang f|(a_+g)\rang=\lang(a_-f)|g\rang⟨f∣(a+g)⟩=⟨(a−f)∣g⟩
- ⟨f∣(a−g)⟩=⟨(a+f)∣g⟩\lang f|(a_-g)\rang=\lang(a_+f)|g\rang⟨f∣(a−g)⟩=⟨(a+f)∣g⟩
那么⟨a+ψn∣a+ψn⟩=⟨a−a+ψn∣ψn⟩=(n+1)⟨ψn∣ψn⟩=(n+1)∫∣ψn∣2dx=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+ψn∣a+ψn⟩=⟨a−a+ψn∣ψn⟩=(n+1)⟨ψn∣ψn⟩=(n+1)∫∣ψn∣2dx=n+1,
同样可得⟨a−ψn∣a−ψn⟩=n\lang a_-\psi_n|a_-\psi_n\rang=n⟨a−ψn∣a−ψn⟩=n。
所以cn=n+1c_n=\sqrt{n+1}cn=n+1,dn=nd_n=\sqrt ndn=n。
所以ψ1=0+1a+ψ0=a+ψ0\psi_1=\sqrt{0+1}a_+\psi_0=a_+\psi_0ψ1=0+1a+ψ0=a+ψ0。ψ2=2a+ψ1=2(a+)2ψ0\psi_2=\sqrt 2a_+\psi_1=\sqrt 2(a_+)^2\psi_0ψ2=2a+ψ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)来计算cnc_ncn(见一维无限深方势阱对应部分),即cn=⟨ψn∣f⟩c_n=\lang \psi_n|f\rangcn=⟨ψn∣f⟩。同样,∣cn∣2|c_n|^2∣cn∣2是取到态nnn的概率,即能量取EnE_nEn的概率。
解析解法
令无量纲量ξ=mwℏx\xi=\sqrt\frac{mw}{\hbar}xξ=ℏmwx。令KKK表示能量中半份ℏω\hbar\omegaℏω的份数,即K=2EℏωK=\frac{2E}{\hbar\omega}K=ℏω2E。则薛定谔方程化为d2ψdξ2=(ξ2−K)ψ\frac{d^2\psi}{d\xi^2}=(\xi^2-K)\psidξ2d2ψ=(ξ2−K)ψ
当ξ\xiξ相当大时可以忽略KKK,即ψP2ξ=ξ2ψ\psi P^2\xi=\xi^2\psiψP2ξ=ξ2ψ。这个方程有近似解ψ(x)=Ae−ξ2/2+Beξ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)(ξ)。
xixixi很大的时候B□B\squareB□项显然太大,没法归一化,因此不要了。最后就是ψ(ξ)→(??)e−ξ2/2,ξ很大\psi(\xi)\to (??)e^{-\xi^2/2},\quad\xi很大ψ(ξ)→(??)e−ξ2/2,ξ很大
我们希望(??)(??)(??)填的是h(ξ)h(\xi)h(ξ)。
这内容也太多了。
反正记Hn(ξ)H_n(\xi)Hn(ξ)为厄米多项式(Hermite polynomials),则ψn(x)=(mωπℏ)1/412nn!Hn(ξ)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)=(πℏmω)1/42nn!1Hn(ξ)e−ξ2/2其中二的幂分数和厄米多项式那部分就相当于代数解法里的nnn个上升算符。
束缚态 散射态
束缚态:bound state
散射态:scattering state
德尔塔函数势阱和势垒
这里的德尔塔函数指的是狄拉克德尔塔函数(Dirac delta function)。该函数在0处取无穷大,在非0处取0,并且在实数域上积分是111。事实上这不是函数,可以称作广义函数或者分布。
势能为V(x)=−αδ(x),αis positiveV(x)=-\alpha\delta(x),\alpha\text{ is positive}V(x)=−αδ(x),α is positive。这是势阱情况,后面会讲势垒。
有束缚态E<0E<0E<0和散射态E>0E>0E>0两种解。因为这个势在研究处是0。
束缚态
对于x<0x<0x<0,根据德尔塔函数定义,V=0V=0V=0。则ψP2x=−2mEℏ2ψ≡κ2ψ\psi P^2x=-\frac{2mE}{\hbar^2}\psi\equiv\kappa^2\psiψP2x=−ℏ22mEψ≡κ2ψ其中κ=−2mEℏ\kappa=\frac{\sqrt{-2mE}}{\hbar}κ=ℏ−2mE。
方程的解是ψ=Eκ(A,B)=Ae−κx+Beκx\psi=\mathcal{E}_\kappa(A,B)=Ae^{-\kappa x}+Be^{\kappa x}ψ=Eκ(A,B)=Ae−κx+Beκx
负无穷远处函数应该有界,所以A=0A=0A=0。
对x>0x>0x>0解法类似,解是ψ=Eκ(F,0)。\psi=\mathcal{E}_\kappa(F,0)。ψ=Eκ(F,0)。
标准边界条件(standard boundary conditions)是:
- ψ\psiψ连续
- 在VVV不是无限的地方,ψPx\psi PxψPx连续
在0点处衔接,因此F=BF=BF=B。即ψ(x)=Be−κ∣x∣\psi(x)=Be^{-\kappa |x|}ψ(x)=Be−κ∣x∣。
对薛定谔方程在−ϵ-\epsilon−ϵ到ϵ\epsilonϵ积分,得−ℏ22m(−2κB)+∫−ϵϵV(x)ψ(x)dx=E∫−ϵϵψ(x)dx-\frac{\hbar^2}{2m}(-2\kappa B)+\int_{-\epsilon}^{\epsilon}V(x)\psi(x)dx=E\int_{-\epsilon}^{\epsilon}\psi(x)dx−2mℏ2(−2κB)+∫−ϵϵV(x)ψ(x)dx=E∫−ϵϵψ(x)dxℏ2κBm−αψ(0)=0\frac{\hbar^2\kappa B}{m}-\alpha\psi(0)=0mℏ2κB−αψ(0)=0α=ℏ2κm\alpha=\frac{\hbar^2\kappa}{m}α=mℏ2κκ=mαℏ2\kappa=\frac{m\alpha}{\hbar^2}κ=ℏ2mαE=−mα22ℏE=-\frac{m\alpha^2}{2\hbar}E=−2ℏmα2
其中标准边界条件其实就是Δ(dψdx)=−2mαℏ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<0x<0x<0,ψP2x=−2mEℏ2ψ≡−k2ψ\psi P^2x=-\frac{2mE}{\hbar^2}\psi\equiv-k^2\psiψP2x=−ℏ22mEψ≡−k2ψ其中k=2mEℏk=\frac{\sqrt{2mE}}{\hbar}k=ℏ2mE。
方程的解是Ck(A,B)\mathcal{C}_k(A,B)Ck(A,B)。对x>0x>0x>0,ψ=Ck(F,G)\psi=\mathcal{C}_k(F,G)ψ=Ck(F,G)。零点连续要求A+B=F+GA+B=F+GA+B=F+G。
Δ(ψPx)=ik(F−G−A+B)\Delta(\psi Px)=ik(F-G-A+B)Δ(ψPx)=ik(F−G−A+B)。ψ(0)=A+B\psi(0)=A+Bψ(0)=A+B。代入标准边界条件化简为F−G=(1+2iβ)A−(1−2iβ)Bwhere β=mαℏ2kF-G=(1+2i\beta)A-(1-2i\beta)B\quad\text{where }\beta=\frac{m\alpha}{\hbar^2k}F−G=(1+2iβ)A−(1−2iβ)Bwhere β=ℏ2kmα
可以用波的观点来看这四个系数。我们考虑一般散射实验,粒子是从一边射进来的,比如说左边,那么从右边射进来的就是0,即G=0G=0G=0。
解之前的方程有B=iβ1−iβAB=\frac{i\beta}{1-i\beta}AB=1−iβiβA,F=11−iβAF=\frac{1}{1-i\beta}AF=1−iβ1A。
反射系数(reflection coefficient)R=∣B∣2∣A∣2=β2/(1+β2)R=\frac{|B|^2}{|A|^2}=\beta^2/(1+\beta^2)R=∣A∣2∣B∣2=β2/(1+β2)
透射系数(transmission coefficient)T=∣F∣2∣A∣2=1/(1+β2)T=\frac{|F|^2}{|A|^2}=1/(1+\beta^2)T=∣A∣2∣F∣2=1/(1+β2)
R+T=1R+T=1R+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<−ax<-ax<−a,V=0V=0V=0,又见到了熟悉的ψP2x=κ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\squareA□会爆炸,因此ψ=Eκ(0,B)\psi=\mathcal{E}_\kappa(0,B)ψ=Eκ(0,B)。
对x∈[−a,a]x\in[-a,a]x∈[−a,a],V=−V0V=-V_0V=−V0。根据束缚态归一化条件E>Vmin=−V0E>V_{\min}=-V_0E>Vmin=−V0,因此E+V0>0E+V_0>0E+V0>0。有ψP2x=−l2ψ\psi P^2x=-l^2\psiψP2x=−l2ψ其中l=2m(E+V0)ℏl=\frac{\sqrt{2m(E+V_0)}}{\hbar}l=ℏ2m(E+V0)。
通解是ψ=Tl(C,D)\psi=\mathcal{T}_l(C,D)ψ=Tl(C,D)。对于x>ax>ax>a,无须赘述,ψ=Eκ(F,0)\psi=\mathcal{E}_\kappa(F,0)ψ=Eκ(F,0)。
由于势能函数是偶函数,因此波函数必须要么是奇函数要么是偶函数(证明待补充)。这样做的好处是我们只用考虑x≥0x\ge 0x≥0或x≤0x\le0x≤0。原来刚才用T\mathcal{T}T符号而不是C\mathcal{C}C符号表示的原因也在于利用奇偶性质。现在先假设是偶函数解一下。2022-10-30 01:12:23
加入sin\sinsin成分会让函数变得不偶,因此C=0C=0C=0。在aaa处,由标准边界条件1,有Fe−κa=Dcos(la)Fe^{-\kappa a}=D\cos(la)Fe−κa=Dcos(la)。由标准边界条件2,−κFe−κa=−Dlsin(la)-\kappa Fe^{-\kappa a}=-Dl\sin(la)−κFe−κa=−Dlsin(la)。用上式除以上上式,得κ=ltan(la)\kappa=l\tan(la)κ=ltan(la)
κ\kappaκ和lll都是正比于EEE的函数,因此这个方程实际上是解出EEE的方程。
令z=laz=laz=la,则tan(z)=κl=−2mE2m(E+V0)=−EE+V0=V0E+V0−1=2mV0(lℏ)2−1=2ma2V0ℏ2z2−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+V0−E=E+V0V0−1=(lℏ)22mV0−1=z2ℏ22ma2V0−1
令z0=aℏ2mV0z_0=\frac{a}{\hbar}\sqrt{2mV_0}z0=ℏa2mV0,则tan(z)=(z0z)2−1\tan(z)=\sqrt{(\frac{z_0}{z})^2-1}tan(z)=(zz0)2−1
解这个超越(transcendental)方程就可以求出能量,用计算机解决。画图表示,右边是一个在0处取无穷大并且逐渐下降,在x>0某点与x轴相交的函数。
有很多个解,分别记为zn,n=1,2,3,…z_n,n=1,2,3,\dotszn,n=1,2,3,…,对应的能量为EnE_nEn。其关系为En=ln2ℏ22m−V0=zn2ℏ22ma2−V0E_n=\frac{l_n^2\hbar^2}{2m}-V_0=\frac{z_n^2\hbar^2}{2ma^2}-V_0En=2mln2ℏ2−V0=2ma2zn2ℏ2−V0
我们考虑两个有趣的特殊情况。
再考虑奇函数情况:Fe−κa=Csin(la)Fe^{-\kappa a}=C\sin(la)Fe−κa=Csin(la);−κFe−κa=lCcos(la)-\kappa Fe^{-\kappa a}=lC\cos(la)−κFe−κa=lCcos(la)。上式除以上上式得κ=−lcot(la)\kappa=-l\cot(la)κ=−lcot(la)
令z=laz=laz=la,则cot(z)=−κl=−(z0z)2−1\cot(z)=-\frac{\kappa}{l}=-\sqrt{(\frac{z_0}{z})^2-1}cot(z)=−lκ=−(zz0)2−1
分析两个特殊情况:
散射态
在左边的零势能处,有ψ=Ck(A,B)\psi=\mathcal{C}_k(A,B)ψ=Ck(A,B),kkk的定义同前。
在势阱里面,有ψ=Tl(C,D)\psi=\mathcal{T}_l(C,D)ψ=Tl(C,D),lll的定义同前。
类似之前,假设粒子从左边射入而不从右边射入,即对于右边的零势能,ψ=Ck(F,0)\psi=\mathcal{C}_k(F,0)ψ=Ck(F,0)。那么AAA是入射幅(incident amplitude),BBB是反射幅(reflected amplitude),FFF是透射幅(transmitted amplitude)。
在−a-a−a和aaa处分别有2个边界条件,一共4个。最后可以消去CCC和DDD,最后用AAA表示FFF,FFF表示BBB。
可以给出透射系数T=1T=1T=1(即完全透射,完全透明)时的条件En+V0=n2π2ℏ28ma2E_n+V_0=\frac{n^2\pi^2\hbar^2}{8ma^2}En+V0=8ma2n2π2ℏ2这个形式是惊人的。也就是说,束缚态里的能量允许取值可以让散射态完全透射。
总结
- 上一篇: 9个最好的PowerPoint模板201
- 下一篇: 幻灯片轮播图(含自动播放和手动点击播放,