摂動
摂動（せつどう、 英語: perturbation）とは、一般に力学系において、主要な力の寄与（主要項）による運動が、他の副次的な力の寄与（摂動項）によって乱される現象である。摂動という語は元来、古典力学において、ある天体の運動が他の天体から受ける引力によって乱れることを指していたが、その類推から量子力学において、粒子の運動が複数粒子の間に相互作用が働くことによって乱れることも指すようになった。なお、転じて摂動現象をもたらす副次的な力のことを摂動と呼ぶ場合がある。
Contents
摂動論
上記のような複数天体間、複数粒子間に相互作用が働くときの運動は数学的に厳密に解くことができないことが知られている（多体問題）。これらの数学的に厳密に解くことのできない問題の近似解を求める手法の1つに、摂動論（せつどうろん、 英語: perturbation theory）がある。具体的には、次のような手順で近似解を求める。
- 考えている問題Aを、厳密に解ける問題Bに小さな変更（摂動）が加えられた問題であるとみなす。
- 問題Aの近似解は、問題Bの厳密解に、摂動が加わったことによって生じる小さな補正（摂動項）を加えたものであると考える。
- ここで求めるべき摂動項は、問題Bの厳密解の組み合わせ、すなわち一次結合の形で表現出来ると考え、その係数を与えられた条件から順次求める。
古典力学における摂動論
天体の運行において、月と地球、太陽と地球などを扱う二体問題は厳密に解くことができるが、三体以上の多体問題を厳密に解くことは不可能である。ただし、月と地球、太陽と地球の問題では、他の天体からの引力による相互作用の効果は近似的に非常に小さいとして、これら二体問題に他の天体からの効果を補正項として考慮することによって十分精度の高い近似解を得ることができる。
量子力学における摂動論
量子力学における多体問題を解く上においても摂動論は重要な近似解法である。
時間に依存せず、縮退のない場合
前提
無摂動部分（無摂動項）のハミルトニアンを[math]\mathcal{H}_0[/math]とし、摂動部分（摂動項）を[math]\mathcal{H}'[/math]とすると、全体のハミルトニアン[math]\mathcal{H}[/math]は、
- [math] \mathcal{H} = \mathcal{H}_0 + \mathcal{H}' [/math]
となる。この時、ゼロ次（無摂動項）のハミルトニアン[math]\mathcal{H}_0[/math]については、すべての固有値（固有エネルギー）[math]\{\epsilon_n^{(0)}\}[/math]と、対応する固有ベクトル[math]\{|\Psi_n^{(0)}\rangle\}[/math]が完全に分かっているとする。ここで「対応する」とは固有値方程式
- [math] \mathcal{H}_0 |\Psi_n^{(0)}\rangle = \epsilon_n^{(0)} |\Psi_n^{(0)}\rangle [/math]
を満たす関係にあるという意味である。 [math]\mathcal{H}_0[/math]はエルミート演算子である（つまりエネルギーはオブザーバブルである）ので、その固有ベクトル[math]\{|\Psi_n^{(0)}\rangle\}[/math]は完全系を成している。また[math]\{|\Psi_n^{(0)}\rangle\}[/math]は規格直交化されているとする。
ハミルトニアン[math]\mathcal{H}[/math]の固有ベクトル[math]\{|\Psi_n\rangle\} \ [/math]と、対応する固有値[math]\{\epsilon_n\}[/math]を求めたい。ここで[math]\{|\Psi_n\rangle\} \ [/math]と[math]\{\epsilon_n\}[/math]は
- [math] \mathcal{H} |\Psi_n\rangle = \epsilon_n |\Psi_n\rangle [/math]
つまり
- [math] (\mathcal{H}_0 + \mathcal{H}') |\Psi_n\rangle = \epsilon_n |\Psi_n\rangle \cdots (0)[/math]
を満たさなければならない。
摂動論
摂動論では、未知の[math]\mathcal{H}'[/math]、[math]|\Psi_n\rangle \ [/math]、[math]\epsilon_n[/math]を、既知の[math]V \ [/math]、[math]|\Psi_n^{(0)}\rangle[/math]、[math]\epsilon_n^{(0)}[/math]と、未知の[math]\{ |\Psi_n^{(1)}\rangle, |\Psi_n^{(2)}\rangle ,\dotsc \} [/math]、[math]\{\epsilon_n^{(1)}, \epsilon_n^{(2)} , \dotsc \} [/math]、微小係数[math]\lambda \ [/math]を用いて
- [math]\begin{align} \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{align}[/math]
と表す。べき級数の中で既知であるのは、第１項目だけであることに注意。 これで、[math]|\Psi_n\rangle \ [/math]、[math]\epsilon_n[/math]を求める問題は[math]\{ |\Psi_n^{(1)}\rangle, |\Psi_n^{(2)}\rangle ,\dotsc \} [/math]、[math]\{\epsilon_n^{(1)}, \epsilon_n^{(2)} , \dotsc \} [/math]を求める問題に変換された。
これらを(0)式に代入し、任意の[math]\lambda \ [/math]で成立すると仮定すると、
- 未知の[math](|\Psi_n^{(1)}\rangle, \epsilon_n^{(1)})[/math]だけを含む方程式 [math]\cdots (1)[/math]
- 未知の[math](|\Psi_n^{(2)}\rangle, \epsilon_n^{(2)})[/math]と[math](|\Psi_n^{(1)}\rangle, \epsilon_n^{(1)})[/math]だけを含む方程式 [math]\cdots (2)[/math]
- 未知の[math](|\Psi_n^{(3)}\rangle, \epsilon_n^{(3)})[/math]と[math](|\Psi_n^{(2)}\rangle, \epsilon_n^{(2)})[/math]と[math](|\Psi_n^{(1)}\rangle, \epsilon_n^{(1)})[/math]だけを含む方程式 [math]\cdots (3)[/math]
[math]\vdots[/math]
が得られ、未知数を分離することができる。 これらを(1)式、(2)式、・・・の順に解いていくと、[math]\{ |\Psi_n^{(1)}\rangle, |\Psi_n^{(2)}\rangle ,\dotsc \} [/math]、[math]\{\epsilon_n^{(1)}, \epsilon_n^{(2)} , \dotsc \} [/math]が求まる。
これらの式は、未知の[math]\{ |\Psi_n^{(1)}\rangle, |\Psi_n^{(2)}\rangle ,\dotsc \} [/math]を、既知の完全系[math]\{|\Psi_n^{(0)}\rangle\}[/math]の線形結合（重ね合わせ）で展開して、その展開係数[math]c_i \ [/math]を求める問題に変換することで解ける。
- [math] |\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 [/math]
結果
エネルギーの一次の摂動は、[math]|\Psi_n^{(0)}\rangle = |n \rangle[/math]とすると、
- [math] \epsilon_n^{(1)} = \langle n|\mathcal{H}'|n\rangle [/math]
固有ベクトルの一次の摂動の展開係数は、[math]i \neq n[/math]とすると
- [math] c_i = - { \langle i|\mathcal{H}'|n\rangle \over { \epsilon_i^{(0)} - \epsilon_n^{(0)} } } [/math]
二次の摂動エネルギーは、
- [math] \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)} } } [/math]
ここで、[math]\lambda\langle n|V|n \rangle = \langle n|\lambda V|n \rangle = \langle n|\mathcal{H}'|n \rangle[/math]である（他の項も同様）。
縮退のある場合
固有値が縮退している場合は、i ≠ n、m ≠ nの場合でもεi = εn、εm = εnとなる場合が存在し、この場合上式二次摂動エネルギーや、一次の摂動波動関数の係数の分母部分が零となり発散してしまう。従って、縮退のある場合には、このような発散を回避する手段を施す必要がある（ほとんど自由な電子参照）。
摂動は普通、一次の項まで考慮すれば十分であるが、より高次な項を考える必要がある場合も多い（例：近藤効果は摂動の二次の項まで考慮しないと説明できない）。
縮退のある場合の一次摂動
(摂動のない)シュレディンガー方程式
- [math]H|\psi\rangle = E|\psi\rangle[/math]
の固有値[math]E_{n}[/math]がk重縮退していて、その対応する固有状態を[math]|n_{i}\rangle \; (i = 1,2,3,...,k, k\in N)[/math]と表す。
微小な摂動[math]gV[/math]([math]g[/math]は無次元の微小項)を加えた後、エネルギー固有値[math]E_{n}[/math]を持っていた状態に関するシュレディンガー方程式は
- [math](H+gV)|\psi_{n}(g)\rangle=E_{n}(g)|\psi_{n}(g)\rangle[/math]
となる。
ここで
- [math]\begin{align} E_{n}(g) &= E_{n}^{(0)} + gE_{n}^{(1)}+ \dotsb \\ |\psi_{n}(g)\rangle &= |\psi_{n}^{(0)}\rangle + g|\psi_{n}^{(1)}\rangle + \dotsb \end{align}[/math]
と展開できるとして、前述のシュレディンガー方程式の0次項を取り出して、
- [math]H|\psi_{n}^{(0)}\rangle=E_{n}^{(0)}|\psi_{n}^{(0)}\rangle[/math]
を得るが、摂動がない時のシュレディンガー方程式より
- [math]|\psi_{n}^{(0)}\rangle=\sum_{i}c_i|n_{i}\rangle[/math]
とおくことができる。
次に、シュレディンガー方程式の1次項を取り出すと、
- [math]\begin{align} H|\psi_{n}^{(1)}\rangle + V|\psi_{n}^{(0)}\rangle &= E_{n}^{(0)}|\psi_{n}^{(1)}\rangle+E_{n}^{(1)}|\psi_{n}^{(0)}\rangle \\ (V - E_{n}^{(1)}) |\psi_{n}^{(0)}\rangle &= (E_{n}^{(0)} - H) |\psi_{n}^{(1)}\rangle \end{align}[/math]
これに左から[math]\langle n_j|[/math]をかけて
- [math]\langle n_j|(V-E_{n}^{(1)})|\psi_{n}^{(0)}\rangle=0[/math]
よって
- [math]c_jE_{n}^{(1)}=\sum_{i}\langle n_j|V|n_i\rangle c_i[/math]
が成り立つ。
これをすべての[math]j[/math]について出すと、[math]n [/math]個の[math]n+1[/math]元方程式が得られるが、規格化を考えていないため、この[math]n+1[/math]個の方程式を解くて、エネルギーの一時摂動及び縮退が解ける様子がわかる。
グリーン関数による方法
ここまでに挙げたのは状態ベクトルに対する摂動論であるが、系が時間に依存する場合など、演算子に対する摂動論も便利である。
演算子に対する摂動論として、グリーン関数を使う方法が知られている。
外部リンク
- Perturbation methods （英語） - スカラーペディア百科事典「摂動法」の項目。