調和振動子 調和振動子の概要

ウィキペディア

調和振動子

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2021/09/10 01:44 UTC 版)

古典的な調和振動子

ニュートンの運動方程式から

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

$x(t)=A\cos \omega t+B\sin \omega t,$

$\omega ={\sqrt {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}}$

である。ハミルトンの正準方程式から連立方程式が得られるが、これを解いても ニュートンの運動方程式 $m{\ddot {x}}=-kx$ を得るだけである。したがって、解は古典力学と同じ結果である。

また、ここで用いたハミルトニアンは量子力学でも使用する。

量子的な調和振動子

1次元の調和振動子

量子力学では運動量演算子 ${\hat {p}}$ を

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

と書く(正準量子化)。$\hbar$ は換算プランク定数、$i$は虚数。よってハミルトニアン${\hat {H}}$は

${\hat {H}}=-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}+{\frac {1}{2}}m\omega ^{2}x^{2}$

となる。

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

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

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

$\phi _{n}(x)=AH_{n}(\xi )\exp \left(-{\frac {\xi ^{2}}{2}}\right)$

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

$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$

である。基底状態（$n=0$）のエネルギー固有状態はガウス波束であり、$x=0$付近に局在している。 エネルギー固有値は次のようになる。

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

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

より高次元の調和振動子

以上は一次元調和振動子の場合であるが、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,....$)

${\begin{aligned}{\hat {a}}^{\dagger }\phi _{n}(x)&={\sqrt {n+1}}\phi _{n+1}(x)\\{\hat {a}}\phi _{n}(x)&={\sqrt {n}}\phi _{n-1}(x)&(n\geq 1)\\{\hat {a}}\phi _{0}(x)&=0\end{aligned}}$

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

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

量子場との関係

場の量子論や量子多体系では、場を量子的な調和振動子に分解することがある。量子的な調和振動子の組があれば、必ずそれをボース粒子の系とみなすことができる。独立な調和振動子からなる系は、エネルギー固有値や平衡状態を議論するかぎり、化学ポテンシャル$\mu =0$の理想ボース気体と数学的に完全に等価である。

ただし全ての場が調和振動子に帰着されるわけではない。調和振動子の集まりと考えることができる場は、双曲線型の微分方程式を満たすものに限られる（詳細は非調和振動子やボゴリューボフ変換を参照）。また粒子像が描けるのは、調和振動子になるような量子場に限られる。たとえばマクスウェルの場の全体が調和振動子の集まりになるわけではなく、遠くのほうに電磁波として伝わっていく成分だけが、調和振動子になる。このとき現れる粒子像が光子である。ただし粒子の数と調和振動子の数には直接的な関係はない。粒子の数が増減すると調和振動子の状態が変化する。

量子的な調和振動子に分解するというのは、量子がもつ粒子性を振幅で解釈し、波動性を振動数で理解しようとする考え方である。この考え方をあえてフェルミ粒子にも適用すると、ボース粒子はいくらでも振幅が大きくなれるが、フェルミ粒子は振幅に制限があるためにあまり大きくなれないと考えることもできる。この量子的な調和振動子の振幅を表すのが生成消滅演算子である。

例

• 光子（電磁場のフーリエ成分）
• フォノン（格子振動の基準振動）