摂動のソースを表示

{{Otheruses|力学における摂動論|天文学の摂動|摂動 (天文学)}}
{{出典の明記|date=2014年11月}}
'''摂動'''（せつどう、 {{lang-en|perturbation}}）とは、一般に力学系において、主要な力の寄与（'''主要項'''）による運動が、他の副次的な力の寄与（'''摂動項'''）によって乱される現象である。摂動という語は元来、[[古典力学]]において、ある天体の運動が他の天体から受ける引力によって乱れることを指していたが、その類推から[[量子力学]]において、粒子の運動が複数粒子の間に[[相互作用]]が働くことによって乱れることも指すようになった。なお、転じて摂動現象をもたらす副次的な力のことを摂動と呼ぶ場合がある。

==摂動論==
上記のような複数天体間、複数粒子間に相互作用が働くときの運動は数学的に厳密に解くことができないことが知られている（[[多体問題]]）。これらの数学的に厳密に解くことのできない問題の近似解を求める手法の1つに、'''摂動論'''（せつどうろん、 {{lang-en|perturbation theory}}）がある。具体的には、次のような手順で近似解を求める。
* 考えている問題Aを、厳密に解ける問題Bに小さな変更（'''摂動'''）が加えられた問題であるとみなす。
* 問題Aの近似解は、問題Bの厳密解に、摂動が加わったことによって生じる小さな補正（'''摂動項'''）を加えたものであると考える。
*ここで求めるべき摂動項は、問題Bの厳密解の組み合わせ、すなわち[[線型結合|一次結合]]の形で表現出来ると考え、その係数を与えられた条件から順次求める。

== 古典力学における摂動論 ==
天体の運行において、[[月]]と[[地球]]、[[太陽]]と[[地球]]などを扱う[[二体問題]]は厳密に解くことができるが、三体以上の[[多体問題]]を厳密に解くことは（一般的な状況であれば）不可能である。ただし、月と地球、太陽と地球の問題では、他の天体からの引力による相互作用の効果は近似的に非常に小さいとして、これら二体問題に他の天体からの効果を補正項として考慮することによって十分精度の高い近似解を得ることができる<ref>[http://th.nao.ac.jp/MEMBER/tanikawa/airborne/tito-01.html 堀源一郎:''Theory of general perturbations with unspecified canonical variables'', Publ. Astron. Soc. Japan, Vol.18, pp.287-296(1966)]</ref>。

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

=== 時間に依存せず、縮退のない場合 ===
==== 前提 ====
無摂動部分（無摂動項）の[[ハミルトニアン]]を<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''の場合でも''ε<SUB>i</SUB>'' = ''ε<SUB>n</SUB>''、''ε<SUB>m</SUB>'' = ''ε<SUB>n</SUB>''となる場合が存在し、この場合上式二次摂動エネルギーや、一次の摂動波動関数の係数の分母部分が零となり発散してしまう。従って、縮退のある場合には、このような発散を回避する手段を施す必要がある（[[ほとんど自由な電子]]参照）。

摂動は普通、一次の項まで考慮すれば十分であるが、より高次な項を考える必要がある場合も多い（例：[[近藤効果]]は摂動の二次の項まで考慮しないと説明できない）。

=== 縮退のある場合の一次摂動 ===
(摂動のない)シュレディンガー方程式

:<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>個の方程式を解くて、エネルギーの一時摂動及び縮退が解ける様子がわかる。

=== グリーン関数による方法 ===
{{See also|グリーン関数 (多体理論)}}
ここまでに挙げたのは状態ベクトルに対する摂動論であるが、系が時間に依存する場合など、演算子に対する摂動論も便利である。

演算子に対する摂動論として、[[グリーン関数]]を使う方法が知られている。

== 脚注 ==
{{reflist}}


== 学習書 ==
* V. P. マスロフ：「摂動論と漸近的方法」、岩波書店（1976年1月28日）.
* 猪木慶治：「量子力学 II」、講談社、ISBN 978-4-06-153212-0 (1994年3月）※第9章：「摂動論」.
* 柴田正和：「漸近級数と特異摂動法：微分方程式の体系的近似解法」、森北出版、ISBN 978-4-627-07631-0 (2009年1月23日).
* J.J.サクライ：「上級量子力学：第II巻 共変な摂動論」、丸善プラネット、ISBN 978-4-86345-048-6 (2010年4月).
* 数理科学、2016年9月号、No.639、サイエンス社、特集「摂動論を考える」.
* 加藤敏夫:「量子力学の数学理論：摂動論と原子等のハミルトニアン」、近代科学社、ISBN 978-4-7649-0545-0 (2017年11月30日).
* 中野寛之、佐合紀親：「重力波・摂動論」、朝倉書店、ISBN 978-4-254-13531-2 (2022年11月1日）※ 一般相対論の摂動解。
* 日置善郎：「場の量子論：摂動計算の基礎」(第3版)、吉岡書店、ISBN 978-4-8427-0377-0 (2022年11月14日).初版は1999年。改訂版は2005年。

== 外部リンク ==
* {{Wayback|url=http://www.scholarpedia.org/article/Perturbation_methods |title=Perturbation methods |date=20061229092119}} - [[スカラーペディア]]百科事典「摂動法」の項目。
{{authority control}}
{{DEFAULTSORT:せつとう}}
[[Category:力学]]
[[Category:量子力学]]
[[Category:計算物理学]]
[[Category:常微分方程式]]
[[Category:数学に関する記事]]