変分法 (解析力学)
変分法 (へんぶんほう、英: Variational method, Calculus of variations)とは、関数を取り値を返す対応である汎関数についての微分にあたる手法を言う。オイラー及びラグランジュらによって導入された[1]。
解析力学における重要な方程式は最小作用の原理を元に変分法を用いて導出される。
Contents
変分法を使った原理
変分法を使った計算例
例えば、物性物理学について考えてみよう。多体問題において多体の波動関数を使って固有値問題を解析的かつ厳密に解くことは困難であり、何らかの近似法を用いて解かれる。その近似手法の一つに変分法がある。
ある多体系において、規格化、直交性などの条件の下で任意に選んだ試行関数(変分関数とも言う。ここでは多体の波動関数)を Ψtrial とする。試行関数はいろいろな選び方があるがここでは、Ψtrial は、系を記述する厳密な固有関数(波動関数)Ψi の展開で記述できるとする。
- [math] \Psi_\mathrm{trial} = \alpha_0 \Psi_0 + \alpha_1 \Psi_1 + \alpha_2 \Psi_2 + \cdots [/math]
ここで、Ψ0 を基底状態の固有関数とする。また、Ψ1, Ψ2, ... は励起状態の固有関数である。系のハミルトニアンを H として、H に対する Ψi に対応する固有値を Ei とすると、試行関数 Ψtrial の固有値 Etrial は、
- [math] \left\langle \Psi_\mathrm{trial}, H \Psi_\mathrm{trial} \right\rangle = E_\mathrm{trial} [/math]
であり、
- [math] \begin{align} E_\mathrm{trial} & = \sum_i \overline{\alpha}_i\alpha_i \left\langle \Psi_i, H \Psi_i\right\rangle \\ & = \sum_i | \alpha_i |^2 E_i \\ & \ge E_0 \end{align} [/math]
となる。この時、試行関数の固有値は、必ず基底状態の固有値 E0(これがこの場合の厳密解)に等しいかエネルギー的により高い値となる。そして、展開係数である αi を調節して Etrial の最小値(最適値)Eopt を求める。これが試行関数を使った変分法の手順である。この場合の最適値 Eopt も、真の固有値 Eexact (= E0) に対し、
- [math] E_\mathrm{opt} \ge E_\mathrm{exact} [/math]
となる。これが満たされない場合、その変分計算は正しくない。以上では、試行関数は厳密解としての Ψ0 を含むという特殊な場合である。実際の計算では厳密解が得られない場合がほとんどである。尚、以上に出てくる固有値は、系の全エネルギーと置き換えて考えても良い。変分法の結果の良し悪しが、試行関数の選び方に強く依存する場合がある。
試行関数の具体例としては、スレーター行列式を使い、個々の一粒子波動関数を最適化するものや、試行関数にジャストロウ型波動関数を使い量子モンテカルロ法を使って最適値を求めたりする。量子化学的手法やバンド計算も変分法が使われており、様々な場面で利用されている。
試行関数を使用しない変分法も存在する。
脚注
参考文献
- 日本数学会 『岩波数学辞典』 岩波書店、1985、第 3 版。ISBN 4000800167。
- 戸川, 隼人 『変分法と有限要素法』 日本評論社、1994。
- R.クーラン、D.ヒルベルト 『数理物理学の方法1』 斎藤 利弥(監訳)、丸山 滋弥(訳)、東京図書株式会社、1959年。
- エルンスト・マッハ 『マッハ力学史 - 力学の発展と批判-(下)』 岩野 秀明(訳)、筑摩書房、2006年。