球座標を用いた変数分離 の履歴(No.7)
更新- 履歴一覧
- 差分 を表示
- 現在との差分 を表示
- ソース を表示
- 量子力学Ⅰ/球座標を用いた変数分離 へ行く。
極座標†
&math( \begin{cases} x=r\sin\theta\cos\phi\\ y=r\sin\theta\sin\phi\\ z=r\cos\theta \end{cases} );
微分の変換†
と との間の変換を考える。
ある関数 について、 が微小量 だけ変化すると、 が だけ変化するとする。
一方、 を 変化させるために をそれぞれ だけ変化させなければならないとする。
このとき、
に、 を代入すれば、
一方、 であるから、
したがって、
&math(\frac{\PD }{\PD x}=\frac{\PD r}{\PD x}\frac{\PD}{\PD r}+\frac{\PD \theta}{\PD x}\frac{\PD}{\PD \theta}+\frac{\PD \phi}{\PD x}\frac{\PD}{\PD \phi}\\ );
同様にして、
&math( \begin{cases} \displaystyle\frac{\PD }{\PD x}=\frac{\PD r}{\PD x}\frac{\PD}{\PD r}+\frac{\PD \theta}{\PD x}\frac{\PD}{\PD \theta}+\frac{\PD \phi}{\PD x}\frac{\PD}{\PD \phi}\\[4mm] \displaystyle\frac{\PD }{\PD y}=\frac{\PD r}{\PD x}\frac{\PD}{\PD r}+\frac{\PD \theta}{\PD x}\frac{\PD}{\PD \theta}+\frac{\PD \phi}{\PD x}\frac{\PD}{\PD \phi}\\[4mm] \displaystyle\frac{\PD }{\PD z}=\frac{\PD r}{\PD x}\frac{\PD}{\PD r}+\frac{\PD \theta}{\PD x}\frac{\PD}{\PD \theta}+\frac{\PD \phi}{\PD x}\frac{\PD}{\PD \phi}\\ \end{cases} );
のように変換される。
計算を進めるには、 などを求める必要がある。
演習:偏微分の計算†
以下、全微分の時と異なり であることに注意せよ。
(1) の関係を用いて、 を で書き表せ。
(2) の関係を用いて、 を で書き表せ。
(3) の関係を用いて、 を で書き表せ。
上記結果を代入すれば、
&math( \begin{cases}
\displaystyle\frac{\PD}{\PD x}= \sin\theta\cos\phi \frac{\PD}{\PD r}
- \frac{1}{r}\cos\theta\cos\phi \frac{\PD}{\PD \theta}
- \frac{\sin\phi}{r\sin\theta} \frac{\PD}{\PD \phi}\\[4mm]
\displaystyle\frac{\PD}{\PD y}= \sin\theta\sin\phi \frac{\PD}{\PD r}
- \frac{1}{r}\cos\theta\sin\phi \frac{\PD}{\PD \theta}
- \frac{\cos\phi}{r\sin\theta} \frac{\PD}{\PD \phi}\\[4mm]
\displaystyle \frac{\PD}{\PD z}= \cos\theta \frac{\PD}{\PD r}
- \frac{1}{r}\sin\theta \frac{\PD}{\PD \theta}\\[4mm]
\end{cases} );
球座標表示のラプラシアン†
に上記を代入するだけで求まる! ・・・ 実際やってみるとえらい大変。→ 計算の詳細
結果だけまとめると、
&math( \nabla^2&=\frac{\PD^2}{\PD r^2}+\frac{2}{r} \frac{\PD}{\PD r}+\frac{1}{r^2}\hat\Lambda\\
&=\frac{1}{r}\frac{\PD^2}{\PD r^2}r+\frac{1}{r^2}\hat\Lambda
);
ただし、
&math(\hat\Lambda=\frac{1}{\sin\theta} \frac{\PD}{\PD \theta} \Big(\sin\theta\frac{\PD}{\PD \theta}\Big)+\frac{1}{\sin^2\theta} \frac{\PD^2}{\PD \phi^2});
球座標の角運動量演算子†
こちらも単に代入すればよいのだけれど、やはり計算は大変 → 詳細はこちら
&math( \begin{cases} \displaystyle \hat l_x=-i\hbar\Big(y\frac{\PD}{\PD z}-z\frac{\PD}{\PD y}\Big) =i\hbar\Big(+\sin\phi\frac{\PD}{\PD\theta}+\frac{\cos\phi}{\tan\theta}\frac{\PD}{\PD\phi}\Big) \\[4mm] \displaystyle \hat l_y=-i\hbar\Big(z\frac{\PD}{\PD x}-x\frac{\PD}{\PD z}\Big) =i\hbar\Big(-\cos\phi\frac{\PD}{\PD\theta}+\frac{\sin\phi}{\tan\theta}\frac{\PD}{\PD\phi}\Big) \\[4mm] \displaystyle \hat l_z=-i\hbar\Big(x\frac{\PD}{\PD y}-y\frac{\PD}{\PD x}\Big) =-i\hbar\frac{\PD}{\PD\phi} \end{cases} );
全角運動量は、
&math( \hat{\bm l}^2&=\hat l_x^2+\hat l_y^2+\hat l_z^2 =-\hbar^2\Big[\frac{1}{\sin\theta}\frac{\PD}{\PD\theta}\Big(\sin\theta\frac{\PD}{\PD\theta}\Big)+\frac{1}{\sin^2\theta}\frac{\PD^2}{\PD\phi^2}\Big] =-\hbar^2\hat\Lambda );
注目すべき重要な結果†
ラプラシアン:
&math(\hat\Lambda=\frac{1}{\sin\theta} \frac{\PD}{\PD \theta} \Big(\sin\theta\frac{\PD}{\PD \theta}\Big)+\frac{1}{\sin^2\theta} \frac{\PD^2}{\PD \phi^2});
軸周りの運動量:
&math( \hat l_z=-i\hbar\frac{\PD}{\PD\phi} );
全角運動量:
&math( \hat{\bm l}^2=-\hbar^2\hat\Lambda );
ラプラシアンの の項の係数は、 全角運動量の演算子と の係数を除いて等しい。
演習:時間に依存しないシュレーディンガー方程式の極座標 変数分離†
球座標表示におけるラプラシアンは、
である。ただし、
&math(\hat\Lambda=\frac{1}{\sin\theta} \frac{\PD}{\PD \theta} \Big(\sin\theta\frac{\PD}{\PD \theta}\Big)+\frac{1}{\sin^2\theta} \frac{\PD^2}{\PD \phi^2});
以下の問いに従って、中心力場中の粒子の運動について考えよ。
(1) を示せ。
(2) 与えられたラプラシアンの表式と (1) の結果を用いて、球座標表示の波動関数 に対する時間に依らないシュレーディンガー方程式を書き下せ。 ただしポテンシャルは動径成分のみに依存するものとし 、 解答には を用いて良い。
(3) 波動関数を と置き、 (2) の方程式を変数分離することにより、以下の方程式を導け。 ただし共通の定数を と置いた。
&math( &-\frac{\hbar^2}{2m}\frac{\PD^2}{\PD r^2}rR(r)+\left\{V(r)+\frac{\hbar^2l(l+1)}{2mr^2}\right\}rR(r)=\varepsilon\,rR(r) );
&math( \hat\Lambda Y(\theta,\phi)=-l(l+1)Y(\theta,\phi) );
(4) (3) の方程式を解いて得られる および が全角運動量 の固有関数であり、その固有値が となることを確かめよ。
(5) 古典論において、質量
の粒子が原点から
の距離を角速度
で回転するときの角運動量は
であり、遠心力は
で与えられる。
ここから遠心力に対するポテンシャルエネルギーが
と書けることを示し、(3) で得た
の方程式に現れる
の項が遠心力の寄与を表わすことを理解せよ。中心力場内では角運動量が保存量となるため、
遠心力とポテンシャルエネルギーとの関係は
一定の元で
であることに注意せよ。
(6) と置いて (3) の第2式を変数分離すると以下の式が得られることを確かめよ。
&math(\left\{\sin\theta \frac{\PD}{\PD \theta}\Big(\sin\theta\frac{\PD}{\PD \theta}\Big)+ l(l+1)\sin^2\theta-m^2\right\}\Theta(\theta)=0);
ただし、共通の定数を と置いた(質量 と紛らわしいが慣例に従った)。
このうち に関する方程式は、
の場合にのみ解を持つことが知られている。
(7) (6), (3) を解いて得られた および が の固有関数であり、その固有値が であることを確かめよ。