時間依存しないシュレーディンガー方程式の具体的な解 #
1. 1D TISE #
$$ \begin{align} \left[-\frac{1}{2}\frac{d^2}{dx^2}+V(x)\right]\psi(x)=E\psi(x) \end{align} $$
こちらのウェブツールで試すことができます(FGH法または差分法で解く、1次元時間依存しないシュレーディンガー方程式。
1.1. 無限に深い井戸型ポテンシャル #
- ポテンシャル
$$ \begin{eqnarray} V(x)=\left\{ \begin{aligned} &0 &&( 0\lt x \lt L)\\ &\infty &&( x\le 0\hspace{0.5em} \text{or}\hspace{0.5em} L\le x) \end{aligned} \right. \end{eqnarray} $$
- 固有値
$$ \begin{align} E_n=\frac{(n+1)^2\pi^2}{2 L^2},\hspace{2em}(n=0,1,2,\cdots) \end{align} $$
- 固有関数
$$ \begin{eqnarray} \psi_n(x)=\left\{ \begin{aligned} &\sqrt{\frac{2}{L}}\sin\left[\frac{(n+1)\pi}{L}\left(x+\frac{L}{2}\right)\right] && (0\lt x \lt L)\\ & 0 && ( x\le 0 \hspace{0.5em}\text{or}\hspace{0.5em} L\le x) \end{aligned} \right. \end{eqnarray} $$
- 具体的なエネルギー固有値 ($L=\pi$)
Infinite Square Well: \(E_n = \frac{(n+1)^2\pi^2}{2L^2}\)
-
解法
区間$0<x<L$,ではシュレーディンガー方程式の$V(x)=0$の解、sin,cosの形をしてなければなりません。 また、$\psi(x=0)=0, \psi(x=L)=0$では境界条件より、波動関数の値はゼロです。 よって、解はsin,cosの形で、そのゼロ点がちょうど$x=0,x=L$にくるような関数となります。
1.2. 有限深さの井戸型ポテンシャル #
- ポテンシャル
$$ \begin{eqnarray} V(x)=\left\{ \begin{aligned} &0 &&( -a \lt x \lt a)\\ &V_0 &&( \text{otherwise}) \end{aligned} \right. \end{eqnarray} $$
- 固有値
$$ \begin{align} E_n=\frac{1}{2}\left(\frac{\xi_n}{a}\right)^2 \end{align} $$
対象な解(偶関数である解):
$$ \begin{align} \left\{ \begin{aligned} &\xi \tan{\xi} = \eta \\ &\xi^2+\eta^2=2 V_0a^2 \end{aligned} \right. \hspace{2em} \rightarrow \hspace{2em} \frac{\xi^2}{\cos^2{\xi}}=2 V_0a^2, (\xi>0, \eta>0) \end{align} $$
反対象な解(奇関数である解):
$$ \begin{align} \left\{ \begin{aligned} &\xi \cot{\xi} = -\eta \\ &\xi^2+\eta^2=2 V_0a^2 \end{aligned} \right. \hspace{2em} \rightarrow \hspace{2em} \frac{\xi^2}{\sin^2{\xi}}=2 V_0a^2, (\xi>0, \eta>0) \end{align} $$
- エネルギー固有値
Equation Solver
\( E_n = \dfrac{1}{2}\left(\dfrac{\xi_n}{a}\right)^2 \)
\( \dfrac{\xi^2}{\cos^2\xi} = 2 V_0 a^2 \) and \( \dfrac{\xi^2}{\sin^2\xi} = 2 V_0 a^2 \) with \(\xi, \eta > 0\)
1.3. 調和ポテンシャル #
- ポテンシャル
$$ \begin{align} V(x)=Ax^2 \end{align} $$
- 固有値
$$ \begin{align} E_n=\sqrt{2A}\left(n+\frac{1}{2}\right),~~(n=0,1,2,\cdots) \end{align} $$
- 固有関数
$$ \begin{gather} \psi_n(x)=\frac{1}{\sqrt{\alpha}}\frac{\pi^{-1/4}}{\sqrt{2^n n!}}H_n(x/\alpha)\exp\left(-\sqrt{\frac{A}{2}} x^2\right)\\ \alpha=\left(2A\right)^{-1/4} \end{gather} $$
ここで$H_n(x)$はエルミート多項式を表し、いくつかの具体的な値は下記の通りです。
$$ \begin{align} H_0(z)&= 1 \nonumber\\ H_1(z)&= 2z \nonumber\\ H_2(z)&= 4z^2-2 \nonumber\\ H_3(z)&= 8z^3-12z \nonumber\\ H_4(z)&= 16z^4-48z^2+12 \nonumber\\ \vdots \nonumber \end{align} $$
$A=1/2$の時、
$$ \begin{align} \psi_0(x)&=\pi^{-1/4}\exp(-x^2/2) \\ \psi_1(x)&=\pi^{-1/4}\sqrt{2}x\exp(-x^2/2) \\ \psi_2(x)&=\pi^{-1/4}\frac{1}{\sqrt{2}}(2x^2-1)\exp(-x^2/2) \\ \psi_3(x)&=\pi^{-1/4}\frac{1}{\sqrt{3}}(2x^3-3x)\exp(-x^2/2) \\ \psi_4(x)&=\pi^{-1/4}\frac{1}{2\sqrt{6}}(4x^4-12x^2+3)\exp(-x^2/2) \end{align} $$
- 具体的なエネルギー固有値
Harmonic Oscillator: \(E_n = \sqrt{2A}\,(n+\tfrac{1}{2})\)
-
解法
生成消滅演算子を用いて求める方法や、$x\to\pm\infty$に置ける漸近形を求めることで微分方程式を解く方法などがあります。 有名ですので、ここでは特に述べません。
1.4. 三角量子井戸 #
- ポテンシャル
$$ \begin{eqnarray} V(x)=\left\{ \begin{aligned} &\infty && (x\lt 0) \\ &\alpha x && (0 \lt x) \end{aligned} \right. \end{eqnarray} $$
- 固有値
$$ \begin{align} E_n=-\left(\frac{\alpha^2}{2}\right)^{1/3}a_n \end{align} $$
ここで、$a_n$はAiry関数1のゼロ点です。
- 固有関数
$$ \begin{align} \psi_n(x)=C\cdot \text{Ai}\left(\left(\frac{2}{\alpha^2}\right)^{1/3}(\alpha x-E_n)\right) \end{align} $$
ここで、$\text{Ai}(x)$はAiry関数1,2です。
- 具体的なエネルギー固有値
Triangular Quantum Well: \(E_n = -\left(\tfrac{\alpha^2}{2}\right)^{1/3} a_n\)
-
解法
区間$0<x<\infty$,ではシュレーディンガー方程式の$V(x)=\alpha x$の解、Airy関数の形をしてなければなりません。 また、$\psi(x=0)=0, \psi(x\to \infty)=0$の境界条件を満たします。 よって、解はAiry関数の形で、$x$が漸近で減衰していく方の解($\text{Ai}(x),\text{Bi}(x)$の$\text{Ai}$の方) で、そのゼロ点がちょうど$x=0$にくるような関数となります。
1.5. F.モース3によるポテンシャル #
- ポテンシャル
$$ \begin{align} V(x)=A\left(e^{-2\alpha x}-2e^{-\alpha x}\right) \end{align} $$
- 固有値
$$ \begin{align} E_n=-A\left[1-\frac{\alpha}{\sqrt{2A}}\left(n+\frac{1}{2}\right)\right]^2 \end{align} $$
ここで、$n$は正の整数で、ゼロから始まり、
$$ \begin{align} \frac{\sqrt{2A}}{\alpha}\gt n+\frac{1}{2} \end{align} $$
を満足する最大値$n_\text{max}$まで。
- 固有関数
$$ \begin{align} \psi(x)=e^{-\xi/2} \xi^s w(\xi) \end{align} $$
ここで、$w(\xi), \xi, s$は下記の通りです。
$$ \begin{gather} w(\xi)=M(-n,2s+1,\xi) \\ \xi=\frac{2\sqrt{2A}}{\alpha}e^{-\alpha x} \\ s=\frac{\sqrt{-2E}}{\alpha} \end{gather} $$
ここで、$M(a,b,x)$は合流型超幾何関数です4。
この問題では離散スペクトルは有限個だけ存在します。もしも
$$ \begin{align} \frac{\sqrt{2A}}{\alpha}\gt n+\frac{1}{2} \end{align} $$
であれば、離散スペクトルは一般に存在しません。
- 具体的なエネルギー固有値
Morse Potential: \(E_n = -A\left[1-\dfrac{\alpha}{\sqrt{2A}}\left(n+\tfrac{1}{2}\right)\right]^2\)
- 解法
参考文献3をご覧ください。
1.6. $1/(\cosh^2)$型ポテンシャル #
- ポテンシャル
$$ \begin{align} V(x)=-\frac{V_0}{\cosh^2 \alpha x} \end{align} $$
- 固有値
$$ \begin{align} E_n=-\frac{\alpha^2}{8}\left[-(1+2n)+\sqrt{1+\frac{8V_0}{\alpha^2}}\right]^2 \end{align} $$
で、$n$は
$$ \begin{align} n\lt s = \frac{1}{2}\left(-1+\sqrt{1+\frac{8V_0}{\alpha^2}}\right) \end{align} $$
の条件から、有限個の準位が決まります。
- 固有関数
$$ \begin{gather} \psi(x)=(1-\xi^2)^{\xi/2}F\left(\varepsilon-s, ~\varepsilon+s+1; ~\varepsilon+1; ~\frac{1-\xi}{2}\right)\\ \varepsilon=\frac{\sqrt{-2E}}{\alpha},\hspace{1em}\xi=\tanh \alpha x \end{gather} $$
ここで$F(a,b;c;x)$は超幾何関数です5。
- 具体的なエネルギー固有値
1/cosh² Potential: \(E_n = -\tfrac{\alpha^2}{8}\left[-(1+2n)+\sqrt{1+\tfrac{8V_0}{\alpha^2}}\right]^2\)
-
解法
参考文献3を参照してください。
1.7. クーロンポテンシャル(動径方向) #
- ポテンシャル
$$ \begin{align} V(x)=\frac{l(l+1)}{2x^2}-\frac{Z}{x} \end{align} $$
$l=0,1,\cdots,\hspace{2em} Z>0$
- 固有値
$$ \begin{align} E_n=-\frac{Z^2}{2n^2},\hspace{2em}(n=1,2,\cdots,) \end{align} $$
- 固有関数
$$ \begin{align} \psi(x)=x\cdot \left(\frac{2}{n}\right)^{3/2}\sqrt{\frac{n-l-1}{2n[(n+l)!]}}\cdot e^{-Zx/n}\cdot \left(\frac{2Zx}{n}\right)^l \cdot L_{n-l-1}^{(2l+1)}\left(\frac{2Zx}{n}\right) \end{align} $$
$$ \begin{align} \psi(x=0)=0,\hspace{1em} \psi(x\to\infty)=0 \end{align} $$
$$ \begin{align} \int_{0}^{\infty} |\psi(x)|^2 dx = 1 \end{align} $$
ここで, $C$は規格化定数, $L_n^{(k)}(x)$はラゲール倍多項式6で
$$ \begin{align} L_0^{(k)}(x)&=1 \nonumber\\ L_1^{(k)}(x)&=-x+k+1 \nonumber\\ L_2^{(k)}(x)&=\frac{1}{2}[x^2-2(k+2)x+(k+1)(k+2)] \nonumber\\ L_3^{(k)}(x)&=\frac{1}{6}[-x^3+3(k+3)x^2-3(k+2)(k+3)x+(k+1)(k+2)(k+3)] \nonumber\\ \vdots\nonumber \end{align} $$
を満たします。
Hydrogen-like Atom: \(E_n = -\dfrac{Z^2}{2n^2}\)
-
解法
$x=0,~x\to\infty$の漸近形を考えて、本当の解を
$$ \begin{align} \text{(解) $=$ ($x\to 0$の漸近形) $\times$ ($x\to\infty$の漸近形) $\times$ (未知関数)} \end{align} $$
と仮定します。これをシュレーディンガー方程式に代入して解けば、未知関数がラゲール陪多項式だと分かります。
2. 2D TISE #
3. 参考文献 #
-
Abramowitz and Stegun, HANDBOOK OF MATHEMATICAL FUNCTIONS, p.446,478 http://people.math.sfu.ca/~cbm/aands/page_446.htm ↩︎ ↩︎
-
13 Confluent Hypergeometric Functions §13.2 Definitions and Basic Properties -DLMF ↩︎
-
15 Hypergeometric Function §15.2 Definitions and Analytical Properties -DLMF ↩︎
-
Associated Laguerre Polynomial -wolfram mathworld 式(22)-(25) ↩︎