Nrthugu

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

古典的な調和振動子

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

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

x(t)=Acos⁡ωt+Bsin⁡ωt,{displaystyle x(t)=Acos omega t+Bsin omega t,}

ω=km{displaystyle omega ={sqrt {frac {k}{m}}}} : 調和振動子の角振動数（固有振動数）

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

ハミルトンの運動方程式（正準方程式）から

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

U=12kx2{displaystyle U={frac {1}{2}}kx^{2}}

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

H=12mp2+k2x2{displaystyle H={frac {1}{2m}}p^{2}+{frac {k}{2}}x^{2}}

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

∂x∂t=∂H∂p{displaystyle {frac {partial x}{partial t}}={frac {partial H}{partial p}}}

∂p∂t=−∂H∂x{displaystyle {frac {partial p}{partial t}}=-{frac {partial H}{partial x}}}

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

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

量子的な調和振動子

1次元の調和振動子

量子力学では運動量 p{displaystyle p} を

p=−iℏ∂∂x{displaystyle p=-ihbar {frac {partial }{partial x}}}

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

H=[−ℏ22m∂2∂x2+12kx2]{displaystyle H=left[-{frac {hbar ^{2}}{2m}}{frac {partial ^{2}}{partial x^{2}}}+{frac {1}{2}}kx^{2}right]}

となる。

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

[−ℏ22m∂2∂x2+12kx2]ϕ(x)=Eϕ{displaystyle left[-{frac {hbar ^{2}}{2m}}{frac {partial ^{2}}{partial x^{2}}}+{frac {1}{2}}kx^{2}right]phi (x)=Ephi }

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

ϕn(x)=AHn(ξ)exp⁡(−ξ22){displaystyle phi _{n}(x)=AH_{n}(xi )exp left(-{frac {xi ^{2}}{2}}right)}

ただし、ξ=mωℏx{displaystyle xi ={sqrt {frac {momega }{hbar }}}x}、A{displaystyle A} は規格化定数で次式で与えられる。

A=1n!2nmωπℏ{displaystyle A={sqrt {{frac {1}{n!2^{n}}}{sqrt {frac {momega }{pi hbar }}}}}}

また、エルミート多項式Hn{displaystyle H_{n}}は

Hn(x)=(−1)nexp⁡(x2)dndxnexp⁡(−x2){displaystyle 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{displaystyle n=0,1,2} の場合を示すと

H0=1{displaystyle H_{0}=1}

H1=2x{displaystyle H_{1}=2x}

H2=4x2+2{displaystyle H_{2}=4x^{2}+2}

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

En=ℏω(n+12)(n=0,1,2,...){displaystyle E_{n}=hbar omega left(n+{frac {1}{2}}right)qquad (n=0,1,2,...)}

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

より高次元の調和振動子

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

EN=ℏω(N+32){displaystyle E_{N}=hbar omega left(N+{frac {3}{2}}right)}

N は三方向の量子数 ( nx{displaystyle n_{x}} , ny{displaystyle n_{y}} , nz{displaystyle n_{z}} ) の和で、またEN{displaystyle E_{N}} は、( N{displaystyle N} +2)( N{displaystyle N} +1)/2 重に縮退している。これは縮退が見られなかった一次元の場合とは明らかに異なる。

生成消滅演算子

調和振動子の扱い方としては、上述の正準変数を用いた方法の他に、生成消滅演算子で書きなおして考える方法がある。

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

a^=ℏ2mω(+∂∂x+mωℏx){displaystyle {hat {a}}={sqrt {frac {hbar }{2momega }}}left(+{frac {partial }{partial x}}+{frac {momega }{hbar }}xright)} : 消滅演算子

a^†=ℏ2mω(−∂∂x+mωℏx){displaystyle {hat {a}}^{dagger }={sqrt {frac {hbar }{2momega }}}left(-{frac {partial }{partial x}}+{frac {momega }{hbar }}xright)} : 生成演算子

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

ℏω(a^†a^+12)ϕ=Eϕ{displaystyle hbar omega left({hat {a}}^{dagger }{hat {a}}+{frac {1}{2}}right)phi =Ephi }

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

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

a^ϕn(x)=nϕn−1(x){displaystyle {hat {a}}phi _{n}(x)={sqrt {n}}phi _{n-1}(x)}

a^†ϕn(x)=n+1ϕn+1(x){displaystyle {hat {a}}^{dagger }phi _{n}(x)={sqrt {n+1,}}phi _{n+1}(x)}

a^ϕ0(x)=0{displaystyle {hat {a}}phi _{0}(x)=0}

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

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

量子場との関係

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

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

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

例

• 
  光子（電磁場のフーリエ成分）


• 
  フォノン（格子振動の基準振動）

具体例

図1:ξ0=0{displaystyle xi _{0}=0}における量子的調和振動子の図

図2:ξ0=0.45{displaystyle xi _{0}=0.45}における量子的調和振動子の図

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

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

波動関数は一般に

ψ(x,t)=∑n=0∞Cnϕn(x)exp⁡(−iω(n+12)t){displaystyle psi (x,t)=sum _{n=0}^{infty }C_{n}phi _{n}(x)exp left(-iomega left(n+{frac {1}{2}}right)tright)}

とかける。ただしCn{displaystyle C_{n}}は波束を決定する係数である。初期条件として零点振動の中心をx0{displaystyle x_{0}} だけ変位させた波束

ψ(x,0)=Aexp⁡(−mω2ℏ(x−x0)2){displaystyle psi (x,0)=Aexp left(-{frac {momega }{2hbar }}(x-x_{0})^{2}right)}

を選ぶ（ただし x0{displaystyle x_{0}} は任意の定数）と、係数Cn{displaystyle C_{n}}はエルミートの多項式の直交性から

Cn=∫−∞∞ϕn(x)ψ(x,0)dx=1n!2n(mωπℏ)(mωℏx0)nexp⁡(−mω4ℏx02)=1n!2n(mωπℏ)(ξ0)nexp⁡(−ξ024){displaystyle {begin{aligned}C_{n}&=int _{-infty }^{infty }phi _{n}(x)psi (x,0),dx\[5pt]&={frac {1}{n!2^{n}}}({sqrt {frac {momega }{pi hbar }}})({sqrt {frac {momega }{hbar }}}x_{0})^{n}exp(-{frac {momega }{4hbar }}x_{0}^{2})\[5pt]&={frac {1}{n!2^{n}}}({sqrt {frac {momega }{pi hbar }}})(xi _{0})^{n}exp(-{frac {xi _{0}^{2}}{4}})end{aligned}}}

で与えられる(ただし、ξ0=mωℏx0{displaystyle xi _{0}={sqrt {frac {momega }{hbar }}}x_{0}} とした)。この場合の粒子の運動が図1、図2である。

図1のアニメーション

ξ0=0.0{displaystyle xi _{0}=0.0}では1≦n{displaystyle 1leqq n}のn{displaystyle n}に対してCn=0{displaystyle C_{n}=0}になる。すなわち波動関数が

ψ(x,t)=C0ϕ0(x)exp⁡(−iω2t){displaystyle psi (x,t)=C_{0}phi _{0}(x)exp left(-{frac {iomega }{2}}tright)}

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

図2のアニメーション

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

p=ℏλ{displaystyle p={frac {hbar }{lambda }}}

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

引用

1. 
   ^ 田崎, 晴明 『統計力学 II』 培風館、2008年12月5日。ISBN 978-4563024383。
   


2. ^ ^ab高橋康 『物理数学ノート<2>力学I』 講談社、1993年12月。ISBN 978-4061532083。
   


3. 
   ^ 鈴木博「場の量子論の考え方」数理科学, No.41, p.6-11 (2008).
   

参考文献

• 振動と波長岡洋介著掌華房1992年


• 量子力学（I）小出昭一郎著掌華房1990年


• 物理学事典（三訂版）培風館2005年

関連項目

• 共振


• 格子振動


• 分子内振動


• 原子核振動


• 非調和振動子


• RLC回路