調和振動子

調和振動子（ちょうわしんどうし、）とは、質点が定点からの距離に比例する引力を受けて運動する系である。調和振動子は定点を中心として振動する系であり、その運動は解析的に解くことができる。
古典的な調和振動子.
ニュートンの運動方程式から.
一端を壁につないだばね定数 formula_1 のばねの他端に質量 formula_2 の物体をつなぐ。静止状態から物体を formula_3 だけ手で引っ張り、静かに手を離すと物体は振動を始める。物体に作用する力は formula_4 である。ニュートンの運動方程式 formula_5 を解くと、一般解は次のようになる。
"A" , "B" は定数で、初期条件によって決まる。振動数formula_8 は、ばね定数と物体の質量にのみ依存する。
ハミルトンの運動方程式（正準方程式）から.
調和振動子のポテンシャルformula_9 は次のようになる。
ただしformula_3は物体の位置である。ばねが自然長の時の位置を原点とする。ハミルトニアン formula_12を求めれば、運動はハミルトンの正準方程式にしたがう。formula_13 は運動エネルギー、formula_14 は運動量である。
ハミルトンの正準方程式は
である。ハミルトンの正準方程式から連立方程式が得られるが、これを解いても ニュートンの運動方程式 formula_5 を得るだけである。したがって、解は古典力学と同じ結果である。
また、ここで用いたハミルトニアンは量子力学でも使用する。
量子的な調和振動子.
1次元の調和振動子.
量子力学では運動量演算子 formula_19 を
と書く(正準量子化)。formula_21 は換算プランク定数、formula_22は虚数。よってハミルトニアンformula_23は
となる。
1次元の量子的な調和振動子についての時間依存しないシュレーディンガー方程式は、以下のように書ける。
この方程式は解析的に解くことができ、その解（エネルギー固有状態）はエルミート多項式 formula_26 を使って以下のように表される。
ただし、formula_28、formula_29 は規格化定数で次式で与えられる。
また、エルミート多項式formula_26は
で定義される。具体例としてformula_33 の場合を示すと
である。基底状態（formula_37）のエネルギー固有状態はガウス波束であり、formula_38付近に局在している。
エネルギー固有値は次のようになる。
つまりエネルギー準位は formula_40 という均等な間隔で並ぶ。formula_41の状態は零点振動、そのエネルギー固有値formula_42は零点エネルギーと呼ばれる。
より高次元の調和振動子.
以上は一次元調和振動子の場合であるが、2次元、3次元も同様に解ける。3次元の場合、エネルギー固有値は次のようになる。
"N" は三方向の量子数 (formula_44, formula_45, formula_46) の和で、またformula_47 は、("N"+2)("N"+1)/2 重に縮退している。これは縮退が見られなかった一次元の場合とは明らかに異なる。
生成消滅演算子.
調和振動子の扱い方としては、上述の正準変数を用いた方法の他に、生成消滅演算子で書きなおして考える方法がある。
以下のような演算子を定義する。
これを使うと、上述のシュレディンガー方程式は次のように書きなおせる。
1/2の項が出るのは演算子に微分が含まれているためである。エネルギー固有値との比較から、formula_51の固有値は formula_52 に等しいことがわかる。よってformula_51を数演算子と呼びformula_54で表す。
生成・消滅演算子をエネルギー固有状態formula_55に作用させると、formula_54の固有値"n" を増減させる。( formula_52 = formula_58)
つまり formula_52 をなんらかの粒子の数と見なすならば、生成演算子は粒子を一つ作り、消滅演算子は一つ減らす働きをする。また基底状態(粒子数0の状態)に消滅演算子を作用させても、もう粒子は消せない。
この演算子を用いれば、方程式の解を容易に導出できる。
量子場との関係.
場の量子論や量子多体系では、場を量子的な調和振動子に分解することがある。量子的な調和振動子の組があれば、必ずそれをボース粒子の系とみなすことができる。独立な調和振動子からなる系は、エネルギー固有値や平衡状態を議論するかぎり、化学ポテンシャルformula_61の理想ボース気体と数学的に完全に等価である。
ただし全ての場が調和振動子に帰着されるわけではない。調和振動子の集まりと考えることができる場は、双曲線型の微分方程式を満たすものに限られる（詳細は非調和振動子やボゴリューボフ変換を参照）。また粒子像が描けるのは、調和振動子になるような量子場に限られる。たとえばマクスウェルの場の全体が調和振動子の集まりになるわけではなく、遠くのほうに電磁波として伝わっていく成分だけが、調和振動子になる。このとき現れる粒子像が光子である。ただし粒子の数と調和振動子の数には直接的な関係はない。粒子の数が増減すると調和振動子の状態が変化する。
量子的な調和振動子に分解するというのは、量子がもつ粒子性を振幅で解釈し、波動性を振動数で理解しようとする考え方である。この考え方をあえてフェルミ粒子にも適用すると、ボース粒子はいくらでも振幅が大きくなれるが、フェルミ粒子は振幅に制限があるためにあまり大きくなれないと考えることもできる。この量子的な調和振動子の振幅を表すのが生成消滅演算子である。
具体例.
量子力学における1次元の調和振動子の運動をアニメーションで示す（図1）（図2）。青い曲線が粒子の波動関数の実部である。緑の曲線が粒子の存在確率密度である。
量子力学では粒子の運動状態を波動関数で表す。波動関数は一般に複素数で与えられる。波動関数の絶対値の2乗が存在確率密度を表す。図1、図2に示される存在確率密度の変動は古典論での粒子の単振動に対応している。
波動関数は一般に
とかける。ただしformula_63は波束を決定する係数である。初期条件として零点振動の中心をformula_64 だけ変位させた波束
を選ぶ（ただし formula_64 は任意の定数）と、係数formula_63はエルミートの多項式の直交性から
で与えられる(ただし、formula_69 とした)。この場合の粒子の運動が図1、図2である。
図1のアニメーション.
formula_70ではformula_71のformula_72に対してformula_73になる。すなわち波動関数が
となる。波動関数は定常波のように振動する。この振動が零点振動である。存在確率密度が時間変化しない定常状態となる。エネルギー固有値は零点エネルギーformula_75であり、エネルギー状態は基底状態である。基底状態はエネルギーが0の状態ではないので波動関数は運動する。
図2のアニメーション.
formula_76ではformula_63がformula_78でない値を持つformula_52が2つ以上存在する。波動関数はエネルギー状態が基底状態の波動関数と励起状態の波動関数の重ね合わせで表される。波動関数の波形は時間によって変化し、定常状態ではない。波動関数は振動の中心付近で速度が最大になる。ド・ブロイの関係式
により速度が大きくなると波長formula_81が短くなるので波動関数の波長が振動の中心付近では振動の端と比べて短くなっている。
関連項目.