跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2015年8月実施 II-2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

一次元(座標 xx)量子系の一粒子問題(質量 mm)において、図1のようなポテンシャル V(x)V(x) に束縛された状態を考える。

V(x)=−Ub≤∣x∣≤a+b,(A)V(x)=-U\qquad b\le|x|\le a+b,\tag{A}
V(x)=0∣x∣<b, ∣x∣>a+b.(B)V(x)=0\qquad |x|<b,\ |x|>a+b.\tag{B}

ただし、U>0,a>0,b≥0U>0,a>0,b\ge0 とする。束縛エネルギー E  (<0)E\;(<0) の状態の波動関数(実関数とする)は以下のシュレディンガー方程式を満たす。

−ℏ22md2dx2ϕ(x)+V(x)ϕ(x)=Eϕ(x).(C)-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\phi(x)+V(x)\phi(x)=E\phi(x).\tag{C}

パリティ変換 x→−xx\to-x に対する V(x)V(x) の対称性より、束縛状態は ϕ(x)=ϕ(−x)\phi(x)=\phi(-x) の正パリティ状態もしくは ϕ(x)=−ϕ(−x)\phi(x)=-\phi(-x) の負パリティ状態である。以下では、aa と UU を一定に保ったまま、bb を変えたときに束縛状態がどのように変化するか考える。

図1:対称二重井戸ポテンシャル(原題)

(1) b=0b=0 の場合を考える。このときのポテンシャルを V(0)(x)V^{(0)}(x) と表すと、V(0)(x)V^{(0)}(x) は幅 2a2a、深さ UU の井戸型ポテンシャルである。γ=−2mE/ℏ2, β=2m(E+U)/ℏ2\gamma=\sqrt{-2mE/\hbar^2},\ \beta=\sqrt{2m(E+U)/\hbar^2} として、正および負パリティ状態のそれぞれについて、γ\gamma を β\beta と aa で表せ。

(2) b>0b>0 で bb が十分小さいとき、b=0b=0 の場合の束縛エネルギー E(0)E^{(0)} からのエネルギー変化を摂動論で評価する。

(a) b=0b=0 の場合の波動関数を ϕ(0)(x)\phi^{(0)}(x) と表す。x=0x=0 および x=ax=a での ϕ(0)(x)\phi^{(0)}(x) の値を、それぞれ、ϕ0≡ϕ(0)(0)\phi_0\equiv\phi^{(0)}(0) および ϕa≡ϕ(0)(a)\phi_a\equiv\phi^{(0)}(a) とおいたとき、正および負パリティ状態のそれぞれについて、ϕ0\phi_0 を ϕa,β,a\phi_a,\beta,a を用いて表し、∣ϕ0∣|\phi_0| と ∣ϕa∣|\phi_a| の大小関係を示せ。

(b) bb を0からわずかに増加させた場合のエネルギー変化 ΔE=E−E(0)\Delta E=E-E^{(0)} を考える。bb の値は十分小さく、ポテンシャルの微小変化 ΔV(x)=V(x)−V(0)(x)\Delta V(x)=V(x)-V^{(0)}(x) を摂動項とみなす。エネルギーの一次の補正項を bb の一次までの近似で求めたとき、正および負パリティ状態のそれぞれについての ΔE\Delta E を ϕ0,ϕa,U,b\phi_0,\phi_a,U,b を用いて示せ。また、ΔE\Delta E と0との大小関係を示せ。ただし、ϕ0\phi_0 および ϕa\phi_a は、

∫−∞+∞{ϕ(0)(x)}2 dx=1\int_{-\infty}^{+\infty}\{\phi^{(0)}(x)\}^2\,dx=1

を満たすように規格化された波動関数に対して定義されているものとする。

(3) a,Ua,U を適当にとり、b=0b=0 の場合に 2n+12n+1 個(nn は0か正の整数)の束縛状態が存在する状況を考える。エネルギーの低い方から正パリティ、負パリティの束縛状態が交互に存在している。a,Ua,U を一定に保ったまま、パラメータ bb を変えたとき、bb をある値 bcb_c より大きくとると負パリティの束縛状態が1つ増える。これは、b>bcb>b_c では 2n+22n+2 個の束縛状態が存在し、bb を bcb_c に近づけると下から 2n+22n+2 番目の状態のエネルギー EE はゼロに近づくことに対応する。ここで、bcb_c は小さくなく摂動論は成り立たないものとする。以下では、この 2n+22n+2 番目の状態の波動関数 ϕ2n+2(x)\phi_{2n+2}(x) の振る舞いを考察し、bcb_c の値を求める。

(a) b>bcb>b_c において、領域 x≥a+bx\ge a+b における波動関数の関数形 ϕ2n+2(x)/ϕ2n+2(a+b)\phi_{2n+2}(x)/\phi_{2n+2}(a+b) を γ=−2mE/ℏ2\gamma=\sqrt{-2mE/\hbar^2} を用いて表せ。また、その結果から、b→bcb\to b_c の極限(E→0E\to0)における x≥a+bx\ge a+b での ϕ2n+2(x)/ϕ2n+2(a+b)\phi_{2n+2}(x)/\phi_{2n+2}(a+b) を示せ。

(b) b→bcb\to b_c の極限(E→0E\to0)において、領域 −b≤x≤b-b\le x\le b における波動関数の関数形 ϕ2n+2(x)/ϕ2n+2(b)\phi_{2n+2}(x)/\phi_{2n+2}(b) を E=0E=0 のシュレディンガー方程式を解くことにより求めよ。ただし、ϕ2n+2(x)\phi_{2n+2}(x) は負パリティ状態であることに注意すること。

(c) b>bcb>b_c のとき、ϕ2n+2(x)\phi_{2n+2}(x) は領域 b<x<a+bb<x<a+b に何個のゼロ点(ϕ2n+2(x)=0\phi_{2n+2}(x)=0 となる点)をもつか。(ヒント:ϕ2n+2(x)\phi_{2n+2}(x) は、∣x∣|x| が有限の領域に 2n+12n+1 個のゼロ点をもつ。)

(d) b→bcb\to b_c の極限(E→0E\to0)での波動関数 ϕ2n+2(x)\phi_{2n+2}(x) を、n=1n=1 の場合について、xx の関数として図示せよ。ただし、波動関数は ϕ2n+2(a+b)=1\phi_{2n+2}(a+b)=1 と規格化すること。

(e) β0=2mU/ℏ2\beta_0=\sqrt{2mU/\hbar^2} とし、bcb_c を β0,a\beta_0,a を用いて表せ。

题目描述​

考虑一维(坐标 xx)量子系统中质量为 mm 的单粒子被图1所示势能 V(x)V(x) 束缚的状态:

V(x)=−Ub≤∣x∣≤a+b,(A)V(x)=-U\qquad b\le|x|\le a+b,\tag{A}
V(x)=0∣x∣<b, ∣x∣>a+b.(B)V(x)=0\qquad |x|<b,\ |x|>a+b.\tag{B}

其中 U>0,a>0,b≥0U>0,a>0,b\ge0。束缚能量为 E  (<0)E\;(<0) 的状态的波函数(设为实函数)满足薛定谔方程

−ℏ22md2dx2ϕ(x)+V(x)ϕ(x)=Eϕ(x).(C)-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\phi(x)+V(x)\phi(x)=E\phi(x).\tag{C}

由于 V(x)V(x) 关于宇称变换 x→−xx\to-x 对称,束缚态为满足 ϕ(x)=ϕ(−x)\phi(x)=\phi(-x) 的正宇称态,或满足 ϕ(x)=−ϕ(−x)\phi(x)=-\phi(-x) 的负宇称态。以下固定 a,Ua,U,研究改变 bb 时束缚态如何变化。

(1) 考虑 b=0b=0 的情况。记此时的势为 V(0)(x)V^{(0)}(x),它是宽 2a2a、深 UU 的方势阱。定义 γ=−2mE/ℏ2, β=2m(E+U)/ℏ2\gamma=\sqrt{-2mE/\hbar^2},\ \beta=\sqrt{2m(E+U)/\hbar^2},分别对正、负宇称态,用 β,a\beta,a 表示 γ\gamma。

(2) 当 b>0b>0 且足够小时,用微扰论估计相对于 b=0b=0 时束缚能量 E(0)E^{(0)} 的能量变化。

(a) 记 b=0b=0 时的波函数为 ϕ(0)(x)\phi^{(0)}(x),在 x=0x=0 和 x=ax=a 处的值分别为 ϕ0≡ϕ(0)(0)\phi_0\equiv\phi^{(0)}(0)、ϕa≡ϕ(0)(a)\phi_a\equiv\phi^{(0)}(a)。分别对正、负宇称态,用 ϕa,β,a\phi_a,\beta,a 表示 ϕ0\phi_0,并比较 ∣ϕ0∣|\phi_0| 与 ∣ϕa∣|\phi_a| 的大小。

(b) 考虑 bb 从0略微增大时的能量变化 ΔE=E−E(0)\Delta E=E-E^{(0)}。设 bb 足够小,将势能的微小变化 ΔV(x)=V(x)−V(0)(x)\Delta V(x)=V(x)-V^{(0)}(x) 视为微扰项。当一阶能量修正保留至 bb 的一阶近似时,分别对正、负宇称态,用 ϕ0,ϕa,U,b\phi_0,\phi_a,U,b 表示 ΔE\Delta E,并比较 ΔE\Delta E 与0的大小。这里 ϕ0,ϕa\phi_0,\phi_a 是对满足归一化条件

∫−∞+∞{ϕ(0)(x)}2 dx=1\int_{-\infty}^{+\infty}\{\phi^{(0)}(x)\}^2\,dx=1

的波函数定义的。

(3) 适当选取 a,Ua,U,使 b=0b=0 时存在 2n+12n+1 个束缚态,其中 nn 为0或正整数。从低能到高能依次为正、负宇称态交替排列。保持 a,Ua,U 不变而改变参数 bb,当 bb 大于某个 bcb_c 时,会新增一个负宇称束缚态。也就是说,b>bcb>b_c 时存在 2n+22n+2 个束缚态,而当 bb 趋于 bcb_c 时,自低能起第 2n+22n+2 个状态的能量 EE 趋于零。设 bcb_c 并不小,微扰论不再适用。以下考察第 2n+22n+2 个状态的波函数 ϕ2n+2(x)\phi_{2n+2}(x),并求 bcb_c。

(a) 当 b>bcb>b_c 时,用 γ=−2mE/ℏ2\gamma=\sqrt{-2mE/\hbar^2} 表示区域 x≥a+bx\ge a+b 中的波函数比值 ϕ2n+2(x)/ϕ2n+2(a+b)\phi_{2n+2}(x)/\phi_{2n+2}(a+b)。利用该结果,给出 b→bcb\to b_c(E→0E\to0)时该区域中这一比值的极限。

(b) 在 b→bcb\to b_c(E→0E\to0)的极限下,通过求解 E=0E=0 的薛定谔方程,求区域 −b≤x≤b-b\le x\le b 内的波函数比值 ϕ2n+2(x)/ϕ2n+2(b)\phi_{2n+2}(x)/\phi_{2n+2}(b)。注意该波函数是负宇称态。

(c) 当 b>bcb>b_c 时,ϕ2n+2(x)\phi_{2n+2}(x) 在区域 b<x<a+bb<x<a+b 内有几个零点,即满足 ϕ2n+2(x)=0\phi_{2n+2}(x)=0 的点?提示:该波函数在 ∣x∣|x| 有限的区域内共有 2n+12n+1 个零点。

(d) 对 n=1n=1,画出 b→bcb\to b_c(E→0E\to0)极限下波函数 ϕ2n+2(x)\phi_{2n+2}(x) 随 xx 的图像。取归一化条件 ϕ2n+2(a+b)=1\phi_{2n+2}(a+b)=1。

(e) 定义 β0=2mU/ℏ2\beta_0=\sqrt{2mU/\hbar^2},用 β0,a\beta_0,a 表示 bcb_c。

Kai​

(1)​

井戸内の解は正パリティで Acos⁡βxA\cos\beta x、負パリティで Asin⁡βxA\sin\beta x。x>ax>a では減衰解 Be−γ(x−a)Be^{-\gamma(x-a)} であり、x=ax=a における対数微分の一致から

γ=βtan⁡βa(正),γ=−βcot⁡βa(負).\boxed{\gamma=\beta\tan\beta a\quad(\text{正}),\qquad \gamma=-\beta\cot\beta a\quad(\text{負})}.

(2a)​

正パリティでは

ϕ0=ϕacos⁡βa,∣ϕ0∣>∣ϕa∣.\boxed{\phi_0=\frac{\phi_a}{\cos\beta a},\qquad |\phi_0|>|\phi_a|}.

束縛状態では γ>0\gamma>0 なので ∣cos⁡βa∣<1|\cos\beta a|<1。負パリティでは

ϕ0=0,∣ϕ0∣<∣ϕa∣.\boxed{\phi_0=0,\qquad |\phi_0|<|\phi_a|}.

(2b)​

∣x∣<b|x|<b に +U+U、a<∣x∣<a+ba<|x|<a+b に −U-U のポテンシャル変化がある。したがって

ΔE(1)=2U[∫0b(ϕ(0))2dx−∫aa+b(ϕ(0))2dx]=2Ub(ϕ02−ϕa2)+O(b2).\begin{gathered} \Delta E^{(1)}=2U\left[\int_0^b(\phi^{(0)})^2dx-\int_a^{a+b}(\phi^{(0)})^2dx\right]\\ =\boxed{2Ub(\phi_0^2-\phi_a^2)+O(b^2)}. \end{gathered}

正パリティでは増加、負パリティでは減少する。

(3a)​

ϕ(x)ϕ(a+b)=e−γ[x−(a+b)]⟶1\boxed{\frac{\phi(x)}{\phi(a+b)}=e^{-\gamma[x-(a+b)]}\longrightarrow1}

である。これは固定された有限の xx での極限であり、閾値解自体は全空間で規格化できない。

(3b)​

E=V=0E=V=0 の中央区で ϕ′′=0\phi''=0、負パリティから定数項はゼロ。

ϕ(x)ϕ(b)=xb.\boxed{\frac{\phi(x)}{\phi(b)}=\frac xb}.

(3c)​

中央区の負パリティ解は E<0E<0 で sinh⁡γx\sinh\gamma x に比例し、節は原点だけ。外側の指数関数に節はない。残る 2n2n 個は左右の井戸に対称に分布するから、右井戸内には n 個\boxed{n\text{ 個}}。

(3d)​

a+bca+b_c で値 1、微分ゼロとなるようにすると、正の半軸上の閾値解は

ϕ(x)={(x/bc)cos⁡(β0a)0≤x≤bc,cos⁡[β0(x−a−bc)]bc≤x≤a+bc,1x≥a+bc,ϕ(−x)=−ϕ(x).\boxed{\begin{gathered} \phi(x)=\begin{cases} (x/b_c)\cos(\beta_0a)&0\le x\le b_c,\\ \cos[\beta_0(x-a-b_c)]&b_c\le x\le a+b_c,\\ 1&x\ge a+b_c, \end{cases}\\ \phi(-x)=-\phi(x) \end{gathered}}.

n=1n=1 では右井戸内に節が一つあり、中央で負の傾き、左右の外側でそれぞれ −1,+1-1,+1 に平坦になる。例として β0a=5π/4\beta_0a=5\pi/4 を選ぶと、右井戸内で −1-1 の極小を一度とり、その後ゼロを横切って +1+1 へ達する。

(3e)​

x=bcx=b_c の微分の連続性より

cos⁡(β0a)bc=β0sin⁡(β0a),bc=cot⁡(β0a)β0.\frac{\cos(\beta_0a)}{b_c}=\beta_0\sin(\beta_0a), \qquad\boxed{b_c=\frac{\cot(\beta_0a)}{\beta_0}}.

b=0b=0 で 2n+12n+1 個の束縛状態を持つ範囲では nπ<β0a≤(n+1/2)πn\pi<\beta_0a\le(n+1/2)\pi、したがって bc≥0b_c\ge0。端点では bc=0b_c=0 として極限をとる。

奇パリティ状態の束縛閾値

Reference​