量子力学Ⅰ/調和振動子 のバックアップ(No.18)

更新


前の単元 <<<              量子力学Ⅰ              >>> 次の単元

目次

量子力学Ⅰ

概要

spring+mass.png

調和振動子とは、理想ばねのように変位に比例する力 f=-Kx により束縛された粒子の系である。*1波数 k と分けるため、ここではバネ定数を大文字 K で表わしている ポテンシャルエネルギーは V(x)=Kx^2/2 となる。 ただしここではばねの自然長位置を x=0 としている。

粒子の質量を m とすれば、古典論では角振動数 \omega=\sqrt{K/m} の単振動

  x=x_0\cos(\omega t+\delta)

が解であり、エネルギーは振幅 x_0 を用いて \varepsilon=Kx_0^2/2 と表せる。

palabora.png

量子力学的なスケールでそのようなばねを思い浮かべることは難しいが、 任意の滑らかで束縛的なポテンシャルの底を2次関数で近似すれば、 そこでの微小振動は一般に調和振動子と見なせる(右図)。

さらに電磁波も、無数の調和振動子の集まりとして表わせることを後に学ぶ。

このように量子力学の重要な問題のあちこちで調和振動子と同等な系が現れるため、 以下ではこの問題を詳細に学んでおく。

演習:1次元の調和振動子

xi.png

(1) 調和振動子のポテンシャルは V(x)=\frac{1}{2}Kx^2 である。 時間に依存しないシュレーディンガー方程式を書け。

(2) このような方程式を解く場合には、変数を無次元化するのが常套手段である。 すなわち、長さの次元を持つ自由変数 x を変数変換して、無次元の量 \xi で記述する。 *2 \xi の書き方の指導法は http://kscalar.kj.yamagata-u.a... を参考にした

ここでは、古典論で得られる角振動数を \omega=\sqrt{\frac{K}{m}} と書いて、

  \xi=\sqrt{\frac{m\omega}{\hbar}}x , \lambda=\frac{2\varepsilon}{\hbar\omega}

と置くと良い。このとき(1)の式が

 &math( \left(-\frac{d^2}{d\xi^2}+\xi^2-\lambda\right)\varphi(\xi)=0 );

となることを示せ。

(3) \xi の大きなところでは \xi^2\gg\lambda となるから、 そこでは \varphi は近似的に次の方程式を満たす。

 &math( \frac{d^2}{d\xi^2}\varphi(\xi)=\xi^2\varphi(\xi) );

ここから予想される

 &math( \varphi(\xi)=X(\xi)e^{\pm\xi^2/2} );

の形の式を (2) の式に代入し、 X(\xi) に対する条件

 &math( X''(\xi)=2\xi X'(\xi)+(1-\lambda) X(\xi) );

を導出せよ。(系が x=0 付近に束縛されていることから、複号は負を取れ)

(4) X(\xi)=\sum_{l=0}^\infty c_l\xi^l と置いて (3) の式に代入し、 係数 c_l の漸化式

  c_{l+2}=\frac{2l+1-\lambda}{(l+2)(l+1)}c_l

を導け。


さて、(4) により、 c_0 を決めると c_2,c_4,c_6,\dots が、 c_1 を決めると c_3,c_5,c_7,\dots が、 それぞれすべて決まる。

c_l が途中でゼロにならない場合、漸化式より l\to \infty において

  \frac{c_{l+2}}{c_l}=\frac{2l+1-\lambda}{(l+2)(l+1)}\to \frac{\,2\,}{l}

が成り立ち、この係数の比は

  \xi^2 e^{\xi^2}=\sum_{n=1}^{\infty}\frac{2^{n-1}}{(2n-2)!!}\xi^{2n}  あるいは、   \xi e^{\xi^2}=\sum_{n=1}^{\infty}\frac{2^{n-1}}{(2n-2)!!}\xi^{2n-1}

n\to\infty の時の係数の比と同じである*3!! は二重階乗(http://ja.wikipedia.org/wiki/%...)を表わす。このようになっていては \xi\to\pm\infty において、 X(\xi)e^{-\xi^2/2}\sim \xi^? e^{\xi^2/2} のように発散してしまう。 求める解は原点付近に束縛された状態であるから、 \xi\to\pm\infty においてゼロに収束するはずであり、そのような解は境界条件を満たさない。

かといって c_0=c_1=0 では \varphi(x)=0 になってしまうから、 c_0\ne 0 あるいは c_1\ne 0 である。

境界条件を満たすのは、ある n\ge 0 において c_n\ne 0 であるものの 2n+1-\lambda=0 が成立することにより、 それ以降の l>n ですべて c_l=0 となる場合である。 このとき、 c_0 あるいは c_1 の片方はゼロでなくて構わないが、 もう一方はゼロでなければならない。


(5) n\ge 0 を量子数として \lambda_n,\varepsilon_n n,\hbar,\omega で表わせ。

(6) X_4(\xi) を求めよ。規格化はしなくて良い。

解答:1次元の調和振動子

>>> 解答はこちら

エネルギー固有値

\varepsilon_0=\hbar\omega/2 がゼロ点エネルギーとなる。

ここでは示さないが、ポテンシャルエネルギーの期待値は \varepsilon_n/2 であり、
運動エネルギーの期待値も \varepsilon_n/2 であり、両者は等しくなる。

すなわち、調和振動子は基底状態においても運動エネルギー、ポテンシャルエネルギーがゼロではない。

\varepsilon_{n+1}-\varepsilon_n=\hbar\omega より、エネルギー固有値は等間隔に並ぶ。

固有関数

  • n=0 のとき X_0(\xi)=c_0
  • n=1 のとき X_1(\xi)=c_1\xi
  • n=2 のとき X_2(\xi)=c_0\left(1-2\xi^2\right)
  • n=3 のとき X_3(\xi)=c_1\left(\xi-\frac{2}{3}\xi^3\right)
  • n=4 のとき X_4(\xi)=c_0\left(1-4\xi^2+\frac{4}{3}\xi^4\right)
  • ・・・

ここでもポテンシャルの対称性を反映して X_{2n}(\xi)e^{-\xi^2/2} は偶関数であり、 X_{2n+1}(\xi)e^{-\xi^2/2} は奇関数である。

固有関数の形

\varphi_n(x) および |\varphi_n(x)|^2 を下から n=0,1,\dots,10 の範囲でプロットした。それぞれ上方向に \varepsilon_n 分だけオフセットしてある。

harmonic2.png   harmonic1.png

二次曲線は y=Kx^2/2 であり、古典的な調和振動子ではこの外には出られない。 有限な箱形ポテンシャルの場合と同様、量子力学的な解は外側にも少しはみ出しており、 \varepsilon>V の領域では振動し、 \varepsilon<V では指数関数的に減衰している。

\varphi_n(x) n+1 個の腹と n 個の節を持ち、 基底状態では偶関数、その後では奇関数と偶関数が交互に現れる。

k を波数とすれば、

  \frac{\hbar^2k^2}{2m}+V(x)=\varepsilon_n  つまり、 \frac{\hbar^2k^2}{2m}=\varepsilon_n-V(x)

の関係が成り立つから、 n が大きいところの解では x=0 の近辺で波長が短く、端の方で波長が長くなっているのが分かる。

同様に、振幅は中央で小さく両端で大きくなっている。 上図の n=10 の付近を拡大して下に示す。

harmonic-density.png

網を掛けたカーブは古典粒子の存在確率を表わしている。 古典論では中央部で速度が速く、両端で速度が遅いため、比較的中央部の存在確率が低くなる。*4古典的な場合の存在確率は P(x)=\frac{2}{v(x)T}=\frac{1}{\pi\sqrt{a^2-x^2}} で与えられる。ただし T は振動周期、 v(x) は速度である。

細かく振動していること、古典振幅を越えて染み出すこと、 の2点を除いて量子力学的な解が古典力学の解と類似していることが分かる。 中央部で \hbar k すなわち運動量が大きいことも古典論と一致している。

下図に示すように、 X_n(\xi) (以下に示すようにエルミート多項式 H_n(\xi) の定数倍になる) は |\xi| の増加に伴い急速に振幅の大きくなる関数であるが、 e^{-\xi^2/2} が外側の振幅を抑えることで上記のような中央部で一定の振幅を持つ波形が得られている。

h10.png

さらに詳しく学ぶ

調和振動子のエネルギーが等間隔に量子化している様子は、 電磁波を光子の集まりと見なせることと数学的に同じ形をしており、 量子力学において非常に重要な例題の一つとなっている。

そこで、上で解いた調和振動子の問題についてさらに詳しく学んでおく。

エルミート多項式

X_n(\xi) エルミート多項式と呼ばれる一連の多項式 H_n(\xi) の定数倍になっている。

エルミート多項式 H_n(\xi) は、母関数 となる次の S(\xi,t) をテイラー展開した際に現れる係数として定義される。

 &math( S(\xi,t)=e^{-t^2+2\xi t}=\sum_{n=0}^\infty\frac{1}{n!}H_n(\xi)t^n );

H_n(\xi) が次の方程式を満たすことを、 母関数を \xi t で微分した形を調べることにより証明できる。 → 詳しくはこちら

  \Big(\frac{d^2}{d\xi^2}-2\xi\frac{d}{d\xi}+\underbrace{2n}_{\lambda-1}\Big)H_n(\xi)=0

この式は X_n(\xi) に対する方程式の \lambda n で書き換えたものと一致する。

また関数系 \big\{H_n(\xi)e^{-\xi^2/2}\big\} に関する、

漸化式:

 &math( H_{n+1}(\xi)e^{-\xi^2/2}=\left(\xi-\frac{d}{d\xi}\right)\left(H_n(\xi)e^{-\xi^2/2}\right) );

 &math( 2nH_{n-1}(\xi)e^{-\xi^2/2}=\left(\xi+\frac{d}{d\xi}\right)\left(H_n(\xi)e^{-\xi^2/2}\right) );

直交性:

 &math(\int_{-\infty}^\infty H_n(\xi)e^{-\xi^2/2}H_m(\xi)e^{-\xi^2/2}d\xi= \left\{\begin{array}{ll} 0&n\ne m\\ 2^nn!\sqrt\pi & n=m \end{array}\right. );

も母関数を使って証明できる。 → 詳しくはこちら

H_0(\xi)=1 および漸化式を用いれば、

  • H_0(\xi)=1
  • H_1(\xi)=2\xi
  • H_2(\xi)=-2+4\xi^2
  • H_3(\xi)=-12\xi+8\xi^3
  • H_4(\xi)=12-48xi^2+16\xi^4
  • H_5(\xi)=120\xi-160\xi^3+32\xi^5
  • H_6(\xi)=-120+720\xi^2-480\xi^4+64\xi^6
  • ・・・

が得られる。

固有関数の性質

上記エルミート多項式の直交性を表わす式を用いることで \varphi_n(x) を正規化できて、

 &math( \varphi_n(x)=\sqrt{ \frac{1}{2^nn!\sqrt\pi} \sqrt{\frac{m\omega}{\hbar}}} H_n\Big(\sqrt{\frac{m\omega}{\hbar}}x\Big) \exp\Big(-\frac{1}{2}\frac{m\omega}{\hbar}x^2\Big) );

このとき正規直交性は、

  \int_{-\infty}^\infty\varphi_n(x)\varphi_m(x)dx=\delta_{nm}

漸化式は、

 &math( \sqrt{\frac{m\omega}{2\hbar}}\left(x-\frac{\hbar}{m\omega}\frac{d}{dx}\right)\varphi_n(x)=\hat a^\dagger\varphi_n(x)=\sqrt{n+1}\,\varphi_{n+1}(x) );

 &math( \sqrt{\frac{m\omega}{2\hbar}}\left(x+\frac{\hbar}{m\omega}\frac{d}{dx}\right)\varphi_n(x)=\hat a\varphi_n(x)=\sqrt{n}\,\varphi_{n-1}(x) );

となる。

\hat a^\dagger,\hat a \hat p=\frac{\hbar}{i}\frac{d}{dx} を用いて、

  \hat a=\sqrt{\frac{m\omega}{2\hbar}}\left(x+\frac{i}{m\omega}\hat p\right)

  \hat a^\dagger=\sqrt{\frac{m\omega}{2\hbar}}\left(x-\frac{i}{m\omega}\hat p\right)

と書ける*5教科書では2つ目の式が \hat a^* として書かれているが、 \hat p^*\ne \hat p なので正しくない。 \hat p^\dagger= \hat p i^\dagger=-i なので \hat a^\dagger と書くのが正しい。←ちょっと自信が無くなってきた。あとでよく調べる。

ここから、

  \hat a^\dagger\hat a\varphi_n(x)=n\varphi_n(x)

  \hat a\hat a^\dagger\varphi_n(x)=(n+1)\varphi_n(x)

  (\hat a\hat a^\dagger-\hat a^\dagger\hat a)\varphi_n(x)=\varphi_n(x)

  \hbar\omega(\hat a^\dagger\hat a+1/2)\varphi_n(x)=\frac{\hbar\omega}{2}(\hat a\hat a^\dagger+\hat a^\dagger\hat a)\varphi_n(x)=\hbar\omega(n+1/2)\varphi_n(x)=\varepsilon_n\varphi_n(x) *6左辺が \hat H\varphi_n(x) に等しいことを 量子力学Ⅰ/調和振動子/メモ#pc1eedc1(https://dora.bk.tsukuba.ac.jp:...) のように示せる

が得られる。これらを

  \hat a^\dagger\hat a=\hat n

  \hat a\hat a^\dagger-\hat a^\dagger\hat a=\hat 1

  \hbar\omega(\hat a^\dagger\hat a+1/2)=\frac{\hbar\omega}{2}(\hat a\hat a^\dagger+\hat a^\dagger\hat a)=\hat H

のように書くこともできる。

ここで学んだ生成・消滅演算子( \hat a^\dagger,\hat a )、 数演算子( \hat n=\hat a^\dagger\hat a )の考え方は、 第二量子化後の波動関数においてフェルミ粒子・ボーズ粒子を扱う際に多用する。 また、上記の交換関係 \hat a\hat a^\dagger-\hat a^\dagger\hat a=\hat 1 はボーズ粒子のそれに相当する。 応用範囲の広い内容であるため、しっかり復習しておくこと。

3次元の調和振動子

粒子が位置 \bm r にあるとき、原点との距離は |\bm r|=\sqrt{x^2+y^2+z^2} であるから、ポテンシャルを V(\bm r)=\frac{K}{2}|\bm r|^2=\frac{K}{2}(x^2+y^2+z^2) と置いて、

  \underbrace{\left(-\frac{\hbar^2}{2m}\nabla^2+V(\bm r)\right)}_{=\hat H(\bm r)}\varphi(\bm r)=\varepsilon\varphi(\bm r)

を解くことになる。

演習

(1) \hat H_x(x)=-\frac{\hbar^2}{2m}\frac{\PD^2}{\PD x^2}+\frac{K}{2}x^2 などと置けば、 \hat H(\bm r)=\hat H_x(x)+\hat H_y(y)+\hat H_z(z) と表せることを確かめよ。

(2) \varphi(\bm r)=X(x)Y(y)Z(z) と置くことにより 3次元調和振動子に対する時間に依らないシュレーディンガー方程式を変数分離し、 それが1次元調和振動子のハミルトニアン \hat H_x(x),\hat H_y(y),\hat H_z(z) の固有値問題に帰着することを示せ。 またそのとき、 \hat H(\bm r) のエネルギー固有値が \hat H_x(x),\hat H_y(y),\hat H_z(z) それぞれの固有値の和として、 \varepsilon=\varepsilon_x+\varepsilon_y+\varepsilon_z のように表せることを示せ。

(3) 1次元調和振動子について得られた結果を用いて、3次元調和振動子の基底状態、 第1励起状態、第2励起状態、第3励起状態のエネルギーを求めよ( \omega=\sqrt{K/m} を用いてよい)。また、それらがそれぞれ何重に縮退しているか答えよ。

(4) 第 n 励起状態が何重に縮退しているか答えよ。


>>> 解答はこちら

形状

(0,0,0), (1,0,0), (1,1,0), (2,0,0), (1,1,1), (2,1,0), (3,0,0), (2,1,1), (2,2,0), (3,1,0), (4,0,0) について |\varphi(\bm r)|^2=0.001 となる面を m\omega/\hbar=1 として描いた。

harmonic000.jpg harmonic100.jpg harmonic110.jpg harmonic200.jpg harmonic111.jpg harmonic210.jpg harmonic300.jpg harmonic211.jpg harmonic220.jpg harmonic310.jpg harmonic400.jpg

状態密度

単位エネルギーあたりどれだけの状態があるか、という量を「状態密度」と言う。

3次元調和振動子では (4) よりほぼ二次関数となる。 横軸を \varepsilon/\hbar\omega=n+3/2 、縦軸を縮退する状態の数としてプロットすれば、

harmonic-dos.jpg



前の単元 <<<              量子力学Ⅰ              >>> 次の単元

質問・コメント





*1 波数 k と分けるため、ここではバネ定数を大文字 K で表わしている
*2 \xi の書き方の指導法は http://kscalar.kj.yamagata-u.ac.jp/~endo/greek/zeta_xi.html を参考にした
*3 !! は二重階乗を表わす
*4 古典的な場合の存在確率は P(x)=\frac{2}{v(x)T}=\frac{1}{\pi\sqrt{a^2-x^2}} で与えられる。ただし T は振動周期、 v(x) は速度である。
*5 教科書では2つ目の式が \hat a^* として書かれているが、 \hat p^*\ne \hat p なので正しくない。 \hat p^\dagger= \hat p i^\dagger=-i なので \hat a^\dagger と書くのが正しい。←ちょっと自信が無くなってきた。あとでよく調べる。
*6 左辺が \hat H\varphi_n(x) に等しいことを 量子力学Ⅰ/調和振動子/メモ#pc1eedc1 のように示せる

Counter: 164263 (from 2010/06/03), today: 1, yesterday: 0