非調和性

古典力学における非調和性（ひちょうわせい、英語: anharmonicity）とは、系の調和振動子からのずれのこと。単振動で振動しない振動子は非調和振動子（英語: anharmonic oscillator）と呼ばれ、系は調和振動子に近似することができ、摂動理論を用いて非調和性を計算することができる。非同調性が大きい場合は、他の数値解析を使用する必要がある。 その結果、 ω {\displaystyle \omega } は振動子の基本周波数とすると、 2 ω {\displaystyle 2\omega } や 3 ω {\displaystyle 3\omega } などの振動数をもつ振動子が現れる。 さらに、振動数 ω {\displaystyle \omega } は調和振動子の振動数 ω 0 {\displaystyle \omega _{0}} からずれる。 第一近似では、振動数のシフト Δ ω = ω − ω 0 {\displaystyle \Delta \omega =\omega -\omega _{0}} は振動子の振幅 A {\displaystyle A} の二乗に比例する。 Δ ω ∝ A 2 {\displaystyle \Delta \omega \propto A^{2}} ω α {\displaystyle \omega _{\alpha }} , ω β {\displaystyle \omega _{\beta }} , ... の固有振動数をもつ振動子の系では、非調和性により振動数 ω α ± ω β {\displaystyle \omega _{\alpha }\pm \omega _{\beta }} をもつ振動子が得られる。 例として次のようなハミルトニアンで表される非調和振動子を考える。 H = p 2 2 m + m ω 2 2 x 2 + λ x 4 {\displaystyle H={\frac {p^{2}}{2m}}+{\frac {m\omega ^{2}}{2}}x^{2}+\lambda x^{4}} 非調和項 λ x 4 {\displaystyle \lambda x^{4}} が十分に小さい（ λ ≪ ℏ ω {\displaystyle \lambda \ll \hbar \omega } ）として1次の摂動まで考えると、非調和振動子のエネルギー準位は次のように調和振動子のエネルギー準位からずれる。 E n ≈ ℏ ω ( n + 1 2 ) + ⟨ n | λ x 4 | n ⟩ {\displaystyle E_{n}\approx \hbar \omega \left(n+{\frac {1}{2}}\right)+\langle n|\lambda x^{4}|n\rangle } ここで | n ⟩ {\displaystyle |n\rangle } は、調和振動子の数演算子の固有状態である。ここで ⟨ n | λ x 4 | n ⟩ {\displaystyle \langle n|\lambda x^{4}|n\rangle } に x = ℏ 2 m ω ( a + a † ) {\displaystyle x={\sqrt {\frac {\hbar }{2m\omega }}}(a+a^{\dagger })} を代入すると、生成消滅演算子についての16個の項が得られる。生成消滅演算子の昇降性により、ゼロでない期待値を与えるのは２個の a {\displaystyle a} と２個の a † {\displaystyle a^{\dagger }} を含む項のみである。よってこの項のみを計算すると、次のようになる。よって調和振動子のように等間隔なエネルギー準位ではないことがわかる。 E n ≈ ℏ ω ( n + 1 2 ) + 3 λ ( ℏ 2 m ω ) 2 ( 2 n 2 + 2 n + 1 ) {\displaystyle E_{n}\approx \hbar \omega \left(n+{\frac {1}{2}}\right)+3\lambda \left({\frac {\hbar }{2m\omega }}\right)^{2}(2n^{2}+2n+1)} 全ハミルトニアン H {\displaystyle H} が自由状態 H 0 {\displaystyle H_{0}} と非調和相互作用 H ′ {\displaystyle H'} の和で表され、それらが2種類のフェルミ粒子の生成消滅演算子 c k , d k {\displaystyle c_{k},d_{k}} で表される場合を考える。 H = H 0 + H ′ {\displaystyle H=H_{0}+H'} H 0 = ∑ k ϵ k ( c k † c k + d k † d k ) {\displaystyle H_{0}=\sum _{k}\epsilon _{k}(c_{k}^{\dagger }c_{k}+d_{k}^{\dagger }d_{k})} H ′ = ∑ k f k ( c k d − k + d − k † d − k † ) {\displaystyle H'=\sum _{k}f_{k}(c_{k}d_{-k}+d_{-k}^{\dagger }d_{-k}^{\dagger })} この全ハミルトニアン H {\displaystyle H} は、ボゴリューボフ変換 C k = c k cos ⁡ θ k − d − k † sin ⁡ θ k {\displaystyle C_{k}=c_{k}\cos \theta _{k}-d_{-k}^{\dagger }\sin \theta _{k}} D − k = d − k cos ⁡ θ k + c k † sin ⁡ θ k {\displaystyle D_{-k}=d_{-k}\cos \theta _{k}+c_{k}^{\dagger }\sin \theta _{k}} によって次のような対角形になり、固有値を求める事ができる。 H = ∑ k E k ( C k † C k + D k † D k ) + W 0 {\displaystyle H=\sum _{k}E_{k}(C_{k}^{\dagger }C_{k}+D_{k}^{\dagger }D_{k})+W_{0}} ここで E k = ϵ k 2 + f k 2 {\displaystyle E_{k}={\sqrt {\epsilon _{k}^{2}+f_{k}^{2}}}} は各量子のエネルギー、 W 0 = ∑ k ( ϵ k − E k ) {\displaystyle W_{0}=\sum _{k}(\epsilon _{k}-E_{k})} は系全体のエネルギーの自由状態からのずれである。よって相互作用ハミルトニアンに現れる関数 f k {\displaystyle f_{k}} の大きさに関わらず量子像は保存される。 全ハミルトニアンが2種類のボース粒子の生成消滅演算子 a k , b k {\displaystyle a_{k},b_{k}} で表される場合を考える。 H = H 0 + H ′ {\displaystyle H=H_{0}+H'} H 0 = ∑ k ϵ k ( a k † a k + b k † b k ) {\displaystyle H_{0}=\sum _{k}\epsilon _{k}(a_{k}^{\dagger }a_{k}+b_{k}^{\dagger }b_{k})} H ′ = ∑ k f k ( a k b − k + b − k † b − k † ) {\displaystyle H'=\sum _{k}f_{k}(a_{k}b_{-k}+b_{-k}^{\dagger }b_{-k}^{\dagger })} このとき全ハミルトニアンの非対角項が消えるような変換ができるのは、 ϵ k 2 > f k 2 {\displaystyle \epsilon _{k}^{2}>f_{k}^{2}} のときだけである。つまりボース粒子でボゴリューボフ変換が使えるのは相互作用が小さいときのみである。相互作用が大きいときには量子像が壊れるのみならず、エネルギーに下限が無くなり、物理的解釈が困難になる。