変分法を使った計算例
出典: フリー百科事典『ウィキペディア(Wikipedia)』 (2019/02/28 00:15 UTC 版)
「変分法 (解析力学)」の記事における「変分法を使った計算例」の解説
例えば、物性物理学について考えてみよう。多体問題において多体の波動関数を使って固有値問題を解析的かつ厳密に解くことは困難であり、何らかの近似法を用いて解かれる。その近似手法の一つに変分法がある。 ある多体系において、規格化、直交性などの条件の下で任意に選んだ試行関数(変分関数とも言う。ここでは多体の波動関数)を Ψtrial とする。試行関数はいろいろな選び方があるがここでは、Ψtrial は、系を記述する厳密な固有関数(波動関数)Ψi の展開で記述できるとする。 Ψ t r i a l = α 0 Ψ 0 + α 1 Ψ 1 + α 2 Ψ 2 + ⋯ {\displaystyle \Psi _{\mathrm {trial} }=\alpha _{0}\Psi _{0}+\alpha _{1}\Psi _{1}+\alpha _{2}\Psi _{2}+\cdots } ここで、Ψ0 を基底状態の固有関数とする。また、Ψ1, Ψ2, ... は励起状態の固有関数である。系のハミルトニアンを H として、H に対する Ψi に対応する固有値を Ei とすると、試行関数 Ψtrial の固有値 Etrial は、 ⟨ Ψ t r i a l , H Ψ t r i a l ⟩ = E t r i a l {\displaystyle \left\langle \Psi _{\mathrm {trial} },H\Psi _{\mathrm {trial} }\right\rangle =E_{\mathrm {trial} }} であり、 E t r i a l = ∑ i α ¯ i α i ⟨ Ψ i , H Ψ i ⟩ = ∑ i | α i | 2 E i ≥ E 0 {\displaystyle {\begin{aligned}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}\\&\geq E_{0}\end{aligned}}} となる。この時、試行関数の固有値は、必ず基底状態の固有値 E0(これがこの場合の厳密解)に等しいかエネルギー的により高い値となる。そして、展開係数である αi を調節して Etrial の最小値(最適値)Eopt を求める。これが試行関数を使った変分法の手順である。この場合の最適値 Eopt も、真の固有値 Eexact (= E0) に対し、 E o p t ≥ E e x a c t {\displaystyle E_{\mathrm {opt} }\geq E_{\mathrm {exact} }} となる。これが満たされない場合、その変分計算は正しくない。以上では、試行関数は厳密解としての Ψ0 を含むという特殊な場合である。実際の計算では厳密解が得られない場合がほとんどである。尚、以上に出てくる固有値は、系の全エネルギーと置き換えて考えても良い。変分法の結果の良し悪しが、試行関数の選び方に強く依存する場合がある。 試行関数の具体例としては、スレーター行列式を使い、個々の一粒子波動関数を最適化するものや、試行関数にジャストロウ型波動関数を使い量子モンテカルロ法を使って最適値を求めたりする。量子化学的手法やバンド計算も変分法が使われており、様々な場面で利用されている。 試行関数を使用しない変分法も存在する。
※この「変分法を使った計算例」の解説は、「変分法 (解析力学)」の解説の一部です。
「変分法を使った計算例」を含む「変分法 (解析力学)」の記事については、「変分法 (解析力学)」の概要を参照ください。
- 変分法を使った計算例のページへのリンク