調和振動子

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索

 調和振動子(ちょうわしんどうし、: harmonic oscillator)とは、質点が定点からの距離に比例する引力を受けて運動する系である。調和振動子は定点を中心として振動する系であり、その運動は解析的に解くことができる。

古典的な調和振動子[編集]

ニュートンの運動方程式から[編集]

 一端を壁につないだばね定数  k のばねの他端に質量  m の物体をつなぐ。静止状態から物体を  x だけ手で引っ張り、静かに手を離すと物体は振動を始める。物体に作用するは  -kx である。ニュートンの運動方程式  m\frac{d^2x}{dt^2}=-kx を解くと、一般解は次のようになる。

x(t)=A\cos\omega t+B\sin\omega t,
\omega=\sqrt{\frac{k}{m}} : 調和振動子の角振動数(固有振動数)

 A , B は定数で、初期条件によって決まる。振動数 \omega は、ばね定数と物体の質量にのみ依存する。

 さらに詳しい議論は自由振動を参照。

ハミルトンの運動方程式(正準方程式)から[編集]

 調和振動子のポテンシャル U は次のようになる。

U=\frac{1}{2}kx^2

 ただし xは物体の位置である。ばねが自然長の時の位置を原点とする。ハミルトニアン   H = T + U を求めれば、運動はハミルトンの正準方程式にしたがう。 T運動エネルギー、 p は運動量である。

H=\frac{1}{2m}p^2+\frac{k}{2}x^2

 ハミルトンの正準方程式は

\frac{\partial x}{\partial t}= \frac{\partial H}{\partial p}
\frac{\partial p}{\partial t}=-\frac{\partial H}{\partial x}

である。ハミルトンの正準方程式から連立方程式が得られるが、これを解いても ニュートンの運動方程式  mx''=-kx を得るだけである。したがって、解は古典力学と同じ結果である。 また、ここで用いたハミルトニアンは量子力学でも使用する。


量子的な調和振動子[編集]

1次元の調和振動子[編集]

 量子力学では運動量  p

p=-i\hbar\frac{\partial}{\partial x}

と演算子で書く。 \hbar換算プランク定数、 iは虚数。よってハミルトニアン H

H=\left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}+\frac{1}{2}kx^2\right]

となる。

 1次元の量子的な調和振動子についての時間依存しないシュレーディンガー方程式は、以下のように書ける。

\left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}+\frac{1}{2}kx^2\right]\phi(x)=E\phi

 この方程式は解析的に解くことができ、その解(エネルギー固有状態)はエルミート多項式 H_n を使って以下のように表される。

\phi_n(x)=AH_n(\xi)\exp\left(-\frac{\xi^2}{2}\right)

ただし、\xi=\sqrt{\frac{m\omega}{\hbar}}xA規格化定数で次式で与えられる。

A=\sqrt{\frac{1}{n!2^n}\sqrt\frac{m\omega}{\pi\hbar}}

 また、エルミート多項式 H_n

H_n(x)=(-1)^n\exp\left(x^2\right)\frac{\mathrm{d}^n}{\mathrm{d}x^n}\exp\left(-x^2\right)

で定義される。具体例として n=0,1,2 の場合を示すと

H_0=1
H_1=2x
H_2=4x^2+2

である。  エネルギー固有値は次のようになる。

E_n=\hbar\omega\left(n+\frac{1}{2}\right)  \qquad (n=0,1,2,...)

つまりエネルギー準位は  \hbar\omega という均等な間隔で並ぶ。  n = 0  の状態は零点振動、そのエネルギー固有値  E_n = \frac{1}{2}\hbar\omega  は零点エネルギーと呼ばれる。

より高次元の調和振動子[編集]

 以上は一次元調和振動子の場合であるが、2次元、3次元も同様に解ける。3次元の場合、エネルギー固有値は次のようになる。

E_N=\hbar\omega\left(N+\frac{3}{2}\right)

N は三方向の量子数 ( n_x , n_y , n_z ) の和で、また E_N は、( N +2)( N +1)/2 重に縮退している。これは縮退が見られなかった一次元の場合とは明らかに異なる。

生成消滅演算子[編集]

 調和振動子の扱い方としては他に生成消滅演算子を使用する方法がある。

 以下のような演算子を定義する。

\hat{a}=\sqrt{\frac{\hbar}{2m\omega}}\left(+\frac{\partial}{\partial x}+\frac{m\omega}{\hbar}x\right) : 消滅演算子
\hat{a}^\dagger=\sqrt{\frac{\hbar}{2m\omega}}\left(-\frac{\partial}{\partial x}+\frac{m\omega}{\hbar}x\right) : 生成演算子

 これを使うと、上述のシュレディンガー方程式は次のように書きなおせる。

\hbar\omega\left(\hat{a}^\dagger\hat{a}+\frac{1}{2}\right)\phi=E\phi

1/2の項が出るのは演算子に微分が含まれているためである。エネルギー固有値との比較から、\hat{a}^\dagger\hat{a}の固有値は  n に等しいことがわかる。よって \hat{a}^\dagger\hat{a}数演算子と呼び \hat{n} \ で表す。

 生成・消滅演算子をエネルギー固有状態 \phi_n(x)に作用させると、 \hat{n} \ の固有値n を増減させる。( n = 0,1,2,....)

\hat{a}\phi_n(x)=\sqrt{n}\phi_{n-1}(x)
\hat{a}^\dagger\phi_n(x)=\sqrt{n+1\,}\phi_{n+1}(x)
\hat{a}\phi_0(x)=0

つまり  n をなんらかの粒子の数と見なすならば、生成演算子は粒子を一つ作り、消滅演算子は一つ減らす働きをする。また基底状態(粒子数0の状態)に消滅演算子を作用させても、もう粒子は消せない。

 この演算子を用いれば、方程式の解を容易に導出できる。


具体例[編集]

 量子力学における1次元の調和振動子の運動をアニメーションで示す。(図1)(図2)青い曲線が粒子の波動関数の実部である。緑の曲線が粒子の存在確率密度である。

 量子力学では粒子の運動状態を波動関数で表す。波動関数は一般に複素数で与えられる。波動関数の絶対値の2乗が存在確率密度を表す。図1、図2に示される存在確率密度の変動は古典論での粒子の単振動に対応している。

 波動関数は一般に

\psi(x,t)=\sum^{\infty}_{n=0} C_n\phi_n(x)\exp\left(-i\omega\left(n+\frac{1}{2}\right)t\right)

 とかける。ただし C_nは波束を決定する係数である。初期条件として零点振動の中心を x_0 だけ変位させた波束

\psi(x,0)=A\exp\left(-\frac{m\omega}{2\hbar}(x-x_0)^2\right)

 を選ぶ(ただし  x_0 は任意の定数)と、係数 C_n はエルミートの多項式の直交性から


\begin{align}
C_n&=\int_{-\infty}^{\infty} \phi_n(x)\psi(x,0)\,dx\\[5pt]
&=\frac{1}{n!2^n}(\sqrt\frac{m\omega}{\pi\hbar})(\sqrt\frac{m\omega}{\hbar}x_0)^n\exp(-\frac{m\omega}{4\hbar}x_0^2)\\[5pt]
&=\frac{1}{n!2^n}(\sqrt\frac{m\omega}{\pi\hbar})(\xi_0)^n\exp(-\frac{\xi_0^2}{4})
\end{align}

で与えられる(ただし、\xi_0 = \sqrt\frac{m\omega}{\hbar}x_0 とした)。この場合の粒子の運動が図1,図2である。

図1:\xi_0=0における量子的調和振動子の図
図2:\xi_0=0.45における量子的調和振動子の図

図1のアニメーション[編集]

 \xi_0 = 0.0 では 1 \leqq n  の  n  に対して C_n=0 になる。すなわち波動関数が

\psi(x,t)=C_0\phi_0(x)\exp\left(-\frac{i\omega}{2}t\right)

となる。波動関数は定常波のように振動する。この振動が零点振動である。存在確率密度が時間変化しない定常状態となる。エネルギー固有値は零点エネルギー  E_n = \frac{1}{2}\hbar\omega  であり、エネルギー状態は基底状態である。基底状態はエネルギーが0の状態ではないので波動関数は運動する。

図2のアニメーション[編集]

 \xi_0 = 0.45  では C_nが 0 でない値を持つ n が2つ以上存在する。波動関数はエネルギー状態が基底状態の波動関数と励起状態の波動関数の重ね合わせで表される。波動関数の波形は時間によって変化し、定常状態ではない。波動関数は振動の中心付近で速度が最大になる。ド・ブロイの関係式

p=\frac{\hbar}{\lambda}

により速度が大きくなると波長 \lambda が短くなるので波動関数の波長が振動の中心付近では振動の端と比べて短くなっている。

関連項目[編集]

参考文献[編集]

  • 振動と波 長岡洋介 著 掌華房 1992年
  • 量子力学(I) 小出昭一郎 著 掌華房 1990年
  • 物理学事典(三訂版) 培風館 2005年