摂動

摂動（せつどう、英語: perturbation）とは、一般に力学系において、主要な力の寄与（主要項）による運動が、他の副次的な力の寄与（摂動項）によって乱される現象である。摂動という語は元来、古典力学において、ある天体の運動が他の天体から受ける引力によって乱れることを指していたが、その類推から量子力学において、粒子の運動が複数粒子の間に相互作用が働くことによって乱れることも指すようになった。なお、転じて摂動現象をもたらす副次的な力のことを摂動と呼ぶ場合がある。

摂動論

上記のような複数天体間、複数粒子間に相互作用が働くときの運動は数学的に厳密に解くことができないことが知られている（多体問題）。これらの数学的に厳密に解くことのできない問題の近似解を求める手法の1つに、摂動論（せつどうろん、英語: perturbation theory）がある。具体的には、次のような手順で近似解を求める。

考えている問題Aを、厳密に解ける問題Bに小さな変更（摂動）が加えられた問題であるとみなす。
問題Aの近似解は、問題Bの厳密解に、摂動が加わったことによって生じる小さな補正（摂動項）を加えたものであると考える。
ここで求めるべき摂動項は、問題Bの厳密解の組み合わせ、すなわち一次結合の形で表現出来ると考え、その係数を与えられた条件から順次求める。

古典力学における摂動論

天体の運行において、月と地球、太陽と地球などを扱う二体問題は厳密に解くことができるが、三体以上の多体問題を厳密に解くことは（一般的な状況であれば）不可能である。ただし、月と地球、太陽と地球の問題では、他の天体からの引力による相互作用の効果は近似的に非常に小さいとして、これら二体問題に他の天体からの効果を補正項として考慮することによって十分精度の高い近似解を得ることができる^[1]。

量子力学における摂動論

量子力学における多体問題を解く上においても摂動論は重要な近似解法である。

時間に依存せず、縮退のない場合

前提

無摂動部分（無摂動項）のハミルトニアンを ${\mathcal {H}}_{0}$ とし、摂動部分（摂動項）を ${\mathcal {H}}'$ とすると、全体のハミルトニアン ${\mathcal {H}}$ は、

{\mathcal {H}}={\mathcal {H}}_{0}+{\mathcal {H}}'

となる。この時、ゼロ次（無摂動項）のハミルトニアン ${\mathcal {H}}_{0}$ については、すべての固有値（固有エネルギー） $\{\epsilon _{n}^{(0)}\}$ と、対応する固有ベクトル $\{|\Psi _{n}^{(0)}\rangle \}$ が完全に分かっているとする。ここで「対応する」とは固有値方程式

{\mathcal {H}}_{0}|\Psi _{n}^{(0)}\rangle =\epsilon _{n}^{(0)}|\Psi _{n}^{(0)}\rangle

を満たす関係にあるという意味である。 ${\mathcal {H}}_{0}$ はエルミート演算子である（つまりエネルギーはオブザーバブルである）ので、その固有ベクトル $\{|\Psi _{n}^{(0)}\rangle \}$ は完全系を成している。また $\{|\Psi _{n}^{(0)}\rangle \}$ は規格直交化されているとする。

ハミルトニアン ${\mathcal {H}}$ の固有ベクトル $\{|\Psi _{n}\rangle \}\$ と、対応する固有値 $\{\epsilon _{n}\}$ を求めたい。ここで $\{|\Psi _{n}\rangle \}\$ と $\{\epsilon _{n}\}$ は

{\mathcal {H}}|\Psi _{n}\rangle =\epsilon _{n}|\Psi _{n}\rangle

つまり

({\mathcal {H}}_{0}+{\mathcal {H}}')|\Psi _{n}\rangle =\epsilon _{n}|\Psi _{n}\rangle \cdots (0)

を満たさなければならない。

摂動論

摂動論では、未知の ${\mathcal {H}}'$ 、 $|\Psi _{n}\rangle \$ 、 $\epsilon _{n}$ を、既知の $V\$ 、 $|\Psi _{n}^{(0)}\rangle$ 、 $\epsilon _{n}^{(0)}$ と、未知の $\{|\Psi _{n}^{(1)}\rangle ,|\Psi _{n}^{(2)}\rangle ,\dotsc \}$ 、 $\{\epsilon _{n}^{(1)},\epsilon _{n}^{(2)},\dotsc \}$ 、微小係数 $\lambda \$ を用いて

{\begin{aligned}{\mathcal {H}}'&=\lambda V\\|\Psi _{n}\rangle &=|\Psi _{n}^{(0)}\rangle +\lambda |\Psi _{n}^{(1)}\rangle +\lambda ^{2}|\Psi _{n}^{(2)}\rangle +\dotsb \\\epsilon _{n}&=\epsilon _{n}^{(0)}+\lambda \epsilon _{n}^{(1)}+\lambda ^{2}\epsilon _{n}^{(2)}+\dotsb \end{aligned}}

と表す。べき級数の中で既知であるのは、第１項目だけであることに注意。これで、 $|\Psi _{n}\rangle \$ 、 $\epsilon _{n}$ を求める問題は $\{|\Psi _{n}^{(1)}\rangle ,|\Psi _{n}^{(2)}\rangle ,\dotsc \}$ 、 $\{\epsilon _{n}^{(1)},\epsilon _{n}^{(2)},\dotsc \}$ を求める問題に変換された。

これらを(0)式に代入し、任意の $\lambda \$ で成立すると仮定すると、

未知の $(|\Psi _{n}^{(1)}\rangle ,\epsilon _{n}^{(1)})$ だけを含む方程式　　 $\cdots (1)$
未知の $(|\Psi _{n}^{(2)}\rangle ,\epsilon _{n}^{(2)})$ と $(|\Psi _{n}^{(1)}\rangle ,\epsilon _{n}^{(1)})$ だけを含む方程式　　 $\cdots (2)$
未知の $(|\Psi _{n}^{(3)}\rangle ,\epsilon _{n}^{(3)})$ と $(|\Psi _{n}^{(2)}\rangle ,\epsilon _{n}^{(2)})$ と $(|\Psi _{n}^{(1)}\rangle ,\epsilon _{n}^{(1)})$ だけを含む方程式　　 $\cdots (3)$

　　　 $\vdots$

が得られ、未知数を分離することができる。これらを(1)式、(2)式、・・・の順に解いていくと、 $\{|\Psi _{n}^{(1)}\rangle ,|\Psi _{n}^{(2)}\rangle ,\dotsc \}$ 、 $\{\epsilon _{n}^{(1)},\epsilon _{n}^{(2)},\dotsc \}$ が求まる。

これらの式は、未知の $\{|\Psi _{n}^{(1)}\rangle ,|\Psi _{n}^{(2)}\rangle ,\dotsc \}$ を、既知の完全系 $\{|\Psi _{n}^{(0)}\rangle \}$ の線形結合（重ね合わせ）で展開して、その展開係数 $c_{i}\$ を求める問題に変換することで解ける。

|\Psi _{n}^{1}\rangle =c_{1}|\Psi _{1}^{(0)}\rangle +c_{2}|\Psi _{2}^{(0)}\rangle +c_{3}|\Psi _{3}^{(0)}\rangle +\dotsb =\sum _{i}c_{i}|\Psi _{i}^{(0)}\rangle

結果

エネルギーの一次の摂動は、 $|\Psi _{n}^{(0)}\rangle =|n\rangle$ とすると、

\epsilon _{n}^{(1)}=\langle n|{\mathcal {H}}'|n\rangle

固有ベクトルの一次の摂動の展開係数は、 $i\neq n$ とすると

c_{i}=-{\langle i|{\mathcal {H}}'|n\rangle  \over {\epsilon _{i}^{(0)}-\epsilon _{n}^{(0)}}}

二次の摂動エネルギーは、

\epsilon _{n}^{(2)}=-\sum _{\{m|m\neq n\}}{\langle n|{\mathcal {H}}'|m\rangle \langle m|{\mathcal {H}}'|n\rangle  \over {\epsilon _{m}^{(0)}-\epsilon _{n}^{(0)}}}

ここで、 $\lambda \langle n|V|n\rangle =\langle n|\lambda V|n\rangle =\langle n|{\mathcal {H}}'|n\rangle$ である（他の項も同様）。

縮退のある場合

固有値が縮退している場合は、i ≠ n、m ≠ nの場合でもε_i = ε_n、ε_m = ε_nとなる場合が存在し、この場合上式二次摂動エネルギーや、一次の摂動波動関数の係数の分母部分が零となり発散してしまう。従って、縮退のある場合には、このような発散を回避する手段を施す必要がある（ほとんど自由な電子参照）。