logo

無限ポテンシャル井戸における波動関数固有関数とエネルギー固有値の求め方 📂量子力学

無限ポテンシャル井戸における波動関数固有関数とエネルギー固有値の求め方

命題

ポテンシャルが区間[0,a][0, a]上で無限の井戸のような形をしているとき、波動関数のエネルギー(固有値EnE_{n}と波動関数(固有状態ψn\psi_{n}は次のようになる。

En=n2π222ma2ψn(x)=2asin(nπax)n=0,1,2,(0) \begin{align*} E_{n} &=\frac{n^{2}\pi^{2}\hbar^{2}}{2ma^{2}} \\[1em] \psi_{n}{(x)} &= \textstyle \sqrt{\frac{2}{a}}\sin \left( \frac{n\pi}{a}x \right) \end{align*} \qquad\qquad n = 0, 1, 2, \dots \tag{0}

説明

V(x)={,<x<00,0<x<a,a<x< V(x) = \begin{cases} \infty, & -\infty \lt x \lt 0 \\ 0, & 0 \lt x \lt a \\ \infty, & a \lt x \lt \infty \end{cases}

上記のような形式のポテンシャルUU無限ポテンシャル井戸infinite potential wellとと呼ばれる。このシステムは、粒子が特定の区間を絶対に抜け出ることができない状況を描写するモデルであり、箱の中の粒子モデルparticle in a box modelとも呼ばれる。非常に単純なモデルだが、古典力学における結果と大きく異なる重要な例を示している。古典力学では粒子が発見される位置は区間内で全て等しいが、量子力学では位置に応じて粒子が発見される確率が異なる。

波動関数の中でエネルギー準位(=n=n)が最も低い状態を基底状態ground stateという。基底状態ではない状態を励起状態excited stateという。

ground state: ψ1(x)=2asin(πax)first excited state: ψ2(x)=2asin(2πax)second excited state: ψ3(x)=2asin(3πax) \begin{align*} \text{ground state: } & \textstyle \psi_{1}(x)=\sqrt{\frac{2}{a}}\sin\left( \frac{\pi}{a}x \right) \\[1em] \text{first excited state: } & \textstyle \psi_2(x)=\sqrt{\frac{2}{a}}\sin\left( \frac{2\pi}{a}x \right) \\[1em] \text{second excited state: } & \textstyle \psi_{3}(x)=\sqrt{\frac{2}{a}}\sin\left( \frac{3\pi}{a}x \right) \\ \vdots \end{align*}

波動関数は任意のエネルギーを持つことができるのではなく、式(0)(0)に現れる形式のエネルギーのみ持つことができる。エネルギーはn2n^{2}に比例するため、連続的な値ではなく離散的な形式で現れる。これを量子化quantizationされているという。

導出

時間に無関係な一次元シュレーディンガー方程式は次のようである。

(22md2dx2+V(x))ψ(x)=Eψ(x) \left(-\frac{\hbar^{2}}{2m}\frac{ d ^{2} }{ d x^{2} }+V(x)\right)\psi(x) = E\psi(x)

V(x)={,<x<00,0<x<a,a<x< V(x) = \begin{cases} \infty, & -\infty \lt x \lt 0 \\ 0, & 0 \lt x \lt a \\ \infty, & a \lt x \lt \infty \end{cases}

E<0E \lt 0

ポテンシャルは常に00以上であるためエネルギーが負のときは波動関数が存在しない

E>0E \gt 0

[1] x>a|x| \gt aについて

この区間では0<E<V=0 \lt E \lt V = \inftyなので、ψ(x)=0\psi(x) = 0である。

[2] 0<x<a0 \lt x \lt aについて

区間[0,a][0, a]においてポテンシャルはV=0V = 0で、シュレーディンガー方程式は次のようになる。

d2ψdx2+2m2Eψ=0 \frac{d^{2} \psi}{dx^{2}} + \frac{2m}{\hbar^{2}}E\psi=0

EEが正であるため2m2E\dfrac{2m}{\hbar^{2}}Eも正であり、したがってこれをk2k^{2}とする。するとシュレーディンガー方程式は次のようになる。

d2ψdx2+k2ψ=0 \frac{d^{2} \psi}{dx^{2}} + k^{2}\psi = 0

この微分方程式の解は次の通りである。

ψ(x)=Asinkx+Bcoskx \psi(x) = A\sin kx + B\cos kx

[3] 境界条件

波動関数ψ\psiは連続であるため、[1]で求めた関数値と[2]で求めた関数値が境界で一致しなければならない。つまりx=0x = 0のとき、次が成立しなければならない。

0=ψ(0)=Asin0+Bcos0=B    B=0 0 = \psi (0) = A\sin 0 + B\cos 0 = B \implies B=0

境界条件はx=ax = aのときも成立しなければならないため、

0=ψ(a)=Asinka(B=0)(1) 0 = \psi (a) = A\sin ka \qquad (\because B = 0) \tag{1}

サイン関数は整数nnについてsinnπ=0\sin n\pi = 0であるためすべての整数nnについて解が存在し、(2)(2)から次が得られる。

Asinka=0    ka=nπ    k=nπa(2) A \sin ka = 0 \implies ka = n\pi \implies k = \frac{n\pi}{a} \tag{2}

    ψn(x)=Asin(nπax) \implies \psi_{n}(x) = \textstyle A\sin \left( \frac{n\pi}{a}x \right)

波動関数は平方積分されて11でなければならないため

ψnψndx=1=0aA2sin2nπaxdx=A20a12(1cos2nπax)dx=A212[xa2nπsin2nπax]0a=A2a2 \begin{align*} \int\limits_{-\infty}^{\infty} \psi_{n}^{\ast} \psi_{n} dx = 1 &= \int_{0}^{a}\textstyle |A|^{2} \sin^{2} \frac{n\pi}{a}x dx \\ &= |A|^{2} \int_{0}^{a}\textstyle \frac{1}{2}(1-\cos \frac{2n\pi}{a}x)dx \\ &= |A|^{2}\textstyle \frac{1}{2} \left[x-\frac{a}{2n\pi}\sin \frac{2n\pi}{a}x\right]_{0}^a \\ &= |A|^{2}\textstyle \frac{a}{2} \end{align*}

    A2=2a    A=2a(3) \implies |A|^{2} = \frac{2}{a} \implies A = \sqrt{\frac{2}{a}} \tag{3}

第三等号では半角公式が使用された。従って波動関数は次のようである。

ψn(x)=2asin(nπax) \psi_{n}(x)= \textstyle \sqrt{\frac{2}{a}}\sin \left( \frac{n\pi}{a}x \right)

また、上述の通り2m2E=k2\dfrac{2m}{\hbar^{2}}E=k^{2}であることから、波動関数固有関数ψn\psi_{n}に対応するエネルギー固有値EnE_{n}は次のようになる。

En=2k22m=n2π222ma2 E_{n} = \frac{\hbar^{2}k^{2}}{2m} = \frac{n^{2}\pi^{2}\hbar^{2}}{2ma^{2}}