調和振動子（ちょうわしんどうし、）とは、質点が定点からの距離に比例する引力を受けて運動する系である。調和振動子は定点を中心として振動する系であり、その運動は解析的に解くことができる。一端を壁につないだばね定数 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_18 を得るだけである。したがって、解は古典力学と同じ結果である。また、ここで用いたハミルトニアンは量子力学でも使用する。量子力学では運動量 formula_14 をと演算子で書く。formula_21 は換算プランク定数、formula_22は虚数。よってハミルトニアンformula_23はとなる。1次元の量子的な調和振動子についての時間依存しないシュレーディンガー方程式は、以下のように書ける。この方程式は解析的に解くことができ、その解（エネルギー固有状態）はエルミート多項式 formula_26 を使って以下のように表される。ただし、formula_28、formula_29 は規格化定数で次式で与えられる。また、エルミート多項式formula_26はで定義される。具体例としてformula_33 の場合を示すとである。エネルギー固有値は次のようになる。つまりエネルギー準位は formula_38 という均等な間隔で並ぶ。formula_39の状態は零点振動、そのエネルギー固有値formula_40は零点エネルギーと呼ばれる。以上は一次元調和振動子の場合であるが、2次元、3次元も同様に解ける。3次元の場合、エネルギー固有値は次のようになる。"N" は三方向の量子数 ( formula_42 , formula_43 , formula_44 ) の和で、またformula_45 は、( formula_46 +2)( formula_46 +1)/2 重に縮退している。これは縮退が見られなかった一次元の場合とは明らかに異なる。調和振動子の扱い方としては、上述の正準変数を用いた方法の他に、生成消滅演算子で書きなおして考える方法がある。以下のような演算子を定義する。これを使うと、上述のシュレディンガー方程式は次のように書きなおせる。1/2の項が出るのは演算子に微分が含まれているためである。エネルギー固有値との比較から、formula_51の固有値は formula_52 に等しいことがわかる。よってformula_51を数演算子と呼びformula_54で表す。生成・消滅演算子をエネルギー固有状態formula_55に作用させると、formula_54の固有値"n" を増減させる。( formula_52 = formula_58)つまり formula_52 をなんらかの粒子の数と見なすならば、生成演算子は粒子を一つ作り、消滅演算子は一つ減らす働きをする。また基底状態(粒子数0の状態)に消滅演算子を作用させても、もう粒子は消せない。この演算子を用いれば、方程式の解を容易に導出できる。場の量子論や量子多体系では、場を量子的な調和振動子に分解することがある。ただし全ての場が調和振動子に帰着されるわけではない。調和振動子の集まりと考えることができる場は、双曲線型の微分方程式を満たすものに限られる。量子的な調和振動子に分解するというのは、量子がもつ粒子性を振幅で解釈し、波動性を振動数で理解しようとする考え方である。この考え方だと、ボース粒子はいくらでも振幅が大きくなれるが、フェルミ粒子は振幅に制限があるためにあまり大きくなれない。この振幅を表すのが生成消滅演算子である。量子力学における1次元の調和振動子の運動をアニメーションで示す（図1）（図2）。青い曲線が粒子の波動関数の実部である。緑の曲線が粒子の存在確率密度である。量子力学では粒子の運動状態を波動関数で表す。波動関数は一般に複素数で与えられる。波動関数の絶対値の2乗が存在確率密度を表す。図1、図2に示される存在確率密度の変動は古典論での粒子の単振動に対応している。波動関数は一般にとかける。ただしformula_64は波束を決定する係数である。初期条件として零点振動の中心をformula_65 だけ変位させた波束を選ぶ（ただし formula_65 は任意の定数）と、係数formula_64はエルミートの多項式の直交性からで与えられる(ただし、formula_70 とした)。この場合の粒子の運動が図1、図2である。formula_71ではformula_72のformula_73に対してformula_74になる。すなわち波動関数がとなる。波動関数は定常波のように振動する。この振動が零点振動である。存在確率密度が時間変化しない定常状態となる。エネルギー固有値は零点エネルギーformula_40であり、エネルギー状態は基底状態である。基底状態はエネルギーが0の状態ではないので波動関数は運動する。formula_77ではformula_64がformula_79でない値を持つformula_52が2つ以上存在する。波動関数はエネルギー状態が基底状態の波動関数と励起状態の波動関数の重ね合わせで表される。波動関数の波形は時間によって変化し、定常状態ではない。波動関数は振動の中心付近で速度が最大になる。ド・ブロイの関係式により速度が大きくなると波長formula_82が短くなるので波動関数の波長が振動の中心付近では振動の端と比べて短くなっている。

出典:wikipedia