変分モンテカルロ法のソースを表示

{{翻訳直後|[[:en:Special:Redirect/revision/836870890|en: Variational Monte Carlo]]|date=2018年8月}}

[[計算物理学]]において、'''変分モンテカルロ法'''（へんぶんモンテカルロほう、{{Lang-en-short|variational Monte Carlo method}}, VMC）とは、量子系の[[基底状態]]を近似的に求めるための[[量子モンテカルロ法]]の一つで、[[変分法 (解析力学)|変分法]]を用いる。

その基本的構成要素はなんらかのパラメータ ''<math> a </math>'' に依存する一般化[[波動関数]] <math>| \Psi(a) \rangle </math> である。このパラメータ ''<math> a </math>'' について系のエネルギーを最小化するような最適値を探索する。

具体的には、[[ハミルトニアン]]を <math> \hat{H} </math>、[[多体問題 (量子論)|多体]]配座を <math> X </math> と書くことにすると、エネルギー[[期待値]]は次のように書くことができる。
:<math> E(a) = \frac{\langle \Psi(a) | \hat{H} | \Psi(a) \rangle} {\langle \Psi(a) |  \Psi(a) \rangle } = \frac{\displaystyle \int dX ~ |\Psi(X,a)| ^2 \frac{\hat{H}\Psi(X,a)}{\Psi(X,a)}} { \displaystyle \int dX ~ |\Psi(X,a)|^2}</math>
ここで、[[モンテカルロ法]]により[[積分法|積分]]を評価する際、<math> | \Psi(X,a) | ^2 / \int dX ~ | \Psi(X,a) | ^2  </math>  を[[確率分布]]関数として標本点を選びだせば、エネルギー期待値 <math> E(a) </math> をいわゆる局所エネルギー <math>E_{\textrm{loc}}(X)  = \hat{H}\Psi(X,a) / \Psi(X,a) </math> の平均値として評価することができる。<math> E(a) </math> が所与の変分パラメータ <math> a </math> について得られたならば、エネルギーを最小化するよう変分パラメータを最適化することにより、基底状態波動関数の可能な限り最適な表現を得ることができる。VMC は、多次元積分を数値積分により評価するという点以外は他の変分法と全く変りがない。この問題においては、ありうる全ての配座  から構成される多体ヒルベルト空間の次元が典型的には物理系のサイズに対し指数関数的に大きくなるため、モンテカルロ積分が特に重要となる。エネルギー期待値を数値的に評価するための他の手法は、モンテカルロ法をもちいるものよりもずっと小さい系にしか適用できないことが一般的である。この手法の精度は変分状態の選択に大きく依存する。 最も簡単な選択は典型的には[[平均場近似|平均場]]形式、すなわち <math> \Psi </math> をヒルベルト空間上で因数分解された形で書くことである。この特に単純な形式では、多体効果が無視されているため典型的にはあまり正確ではない。波動関数を分離可能な形で書いた場合に比べた際の精度向上に最も寄与するものの一つは、いわゆる Jastrow 因子の導入である。この場合、波動関数は <math> \Psi(X) = \exp(\sum{u(r_{ij})})</math> のように書き下される。ここで、<math> r_{ij} </math> は粒子対間の距離、<math> u(r) </math> は変分法により決定される関数である。この因子により明示的に粒子・粒子間の相関を考慮することができるが、多体積分は分離不可能となるため、これを効率的に評価できる方法はモンテカルロ法のみとなる。 化学的系においては、30個未満のパラメータで[[電子相関]]エネルギーの80%から90%を得ることができる、わずかだけ洗練されたものが存在する。これに比べて、配置間相互作用法を用いて同等の精度を得るためには、考慮する系によるもののおよそ50,000個のパラメータが必要となる。さらに、VMCは通常粒子数の数乗のオーダーでスケールする。波動関数の関数形によるものの、エネルギー期待値の計算は {{Math|''N''<sup>2–4</sup>}} 程度のオーダーでスケールする。

== VMC における波動関数の最適化 ==
量子モンテカルロ計算においては、試行関数の質が非常に重要となってくる。そのため、波動関数をできる限り基底状態に近くなるよう最適化することが欠かせない。数値シミュレーションの分野において、関数[[数理最適化|最適化]]問題は非常に重要な研究課題である。量子モンテカルロ法では、多次元関数の最小化問題にまつわる通常の困難に加えて、コスト関数（通常はエネルギー）その微分（効率的な最適化のために必要となる）の推定に統計的ノイズが乗るという問題がある。

多体試行関数の最適化には、さまざまなコスト関数とさまざまな戦略が用いられる。量子モンテカルロ最適化においてよく用いられるのはエネルギー、分散、およびそれら二つの線形結合の最適化である。分散最適化法は、波動関数の厳密解の分散が既知である場合に有用である（波動関数の厳密解はハミルトニアンの固有関数であるから、局所エネルギーの分散はゼロである）。したがって、分散最適化は下限があり、正定値で、最低値が既知であるという点で理想的である。しかし、エネルギー最適化のほうが分散最適化よりもより効率的であることを複数の研究者が見い出している。

また、変分モンテカルロ法でも拡散モンテカルロ法でも通常は最低エネルギー状態が興味の対象であって、分散最低状態ではない。さらに、分散最適化は決定的なパラメータの最適化に多数回の反復を要しかつ多数存在する局所最適に囚われてしまう「偽収束」の問題がある。加えて、エネルギー最小化された波動関数は分散最小化された波動関数と比べて、平均的には他の物理量の期待値もより精度良く計算できる。

最適化戦略は3つのカテゴリーに分けられる。1つは相関サンプリングと決定的最適化法に基くものである。このアイデアは第1周期元素については非常に精度のよい結果をもたらすものの、パラメータが波動関数の節に影響する場合は問題が引き起こされ、さらに初期試行関数と現在の試行関数との間の密度比が系の規模に対して指数関数的に増大するという問題がある。2つ目の戦略はコスト関数およびその微分の評価に使うビンを、ノイズが無視できるほど大きくとり、その上で決定的手法を使う方法である。

3つ目のアプローチでは、ノイズ関数を直接扱うために反復的手法を用いる。この種類の手法の最初の例は、確率的勾配近似法 (SGA) である。この手法は構造最適化にも用いられる。近年、確率的再配置 (SR) 法と呼ばれるより進んだ高速なアプローチが提案されている。

== 参照文献 ==
* {{Cite journal|last=McMillan|first=W. L.|date=1965-04-19|title=<nowiki>Ground State of Liquid ${\mathrm{He}}^{4}$</nowiki>|url=https://link.aps.org/doi/10.1103/PhysRev.138.A442|journal=Physical Review|volume=138|issue=2A|pages=A442–A451|doi=10.1103/PhysRev.138.A442}}
* {{Cite journal|last=Ceperley|first=D.|last2=Chester|first2=G. V.|last3=Kalos|first3=M. H.|date=1977-10-01|title=Monte Carlo simulation of a many-fermion study|url=https://link.aps.org/doi/10.1103/PhysRevB.16.3081|journal=Physical Review B|volume=16|issue=7|pages=3081–3099|doi=10.1103/PhysRevB.16.3081}}
* 変分モンテカルロ法による波動関数の最適化
** {{Cite journal|last=Snajdr|first=Martin|last2=Rothstein|first2=Stuart M.|date=2000-03-15|title=Are properties derived from variance-optimized wave functions generally more accurate? Monte Carlo study of non-energy-related properties of H2, He, and LiH|url=https://aip.scitation.org/doi/10.1063/1.481047|journal=The Journal of Chemical Physics|volume=112|issue=11|pages=4935–4941|language=en|doi=10.1063/1.481047|issn=0021-9606}}
** {{Cite journal|last=Bressanini|first=Dario|last2=Morosi|first2=Gabriele|last3=Mella|first3=Massimo|date=2002-04|title=Robust wave function optimization procedures in quantum Monte Carlo methods|url=https://aip.scitation.org/doi/10.1063/1.1455618|journal=The Journal of Chemical Physics|volume=116|issue=13|pages=5345–5350|language=en|doi=10.1063/1.1455618|issn=0021-9606|arxiv=physics/0110003}}
** {{Cite journal|last=Umrigar|first=C. J.|last2=Wilson|first2=K. G.|last3=Wilkins|first3=J. W.|date=1988-04-25|title=Optimized trial wave functions for quantum Monte Carlo calculations|url=https://link.aps.org/doi/10.1103/PhysRevLett.60.1719|journal=Physical Review Letters|volume=60|issue=17|pages=1719–1722|doi=10.1103/PhysRevLett.60.1719}}
** {{Cite journal|last=Kent|first=P. R. C.|last2=Needs|first2=R. J.|last3=Rajagopal|first3=G.|date=1999-05-15|title=Monte Carlo energy and variance-minimization techniques for optimizing many-body wave functions|url=https://link.aps.org/doi/10.1103/PhysRevB.59.12344|journal=Physical Review B|volume=59|issue=19|pages=12344–12351|doi=10.1103/PhysRevB.59.12344}}
** {{Cite journal|last=Lin|first=Xi|last2=Zhang|first2=Hongkai|last3=Rappe|first3=Andrew M.|date=2000-02-08|title=Optimization of quantum Monte Carlo wave functions using analytical energy derivatives|url=https://aip.scitation.org/doi/10.1063/1.480839|journal=The Journal of Chemical Physics|volume=112|issue=6|pages=2650–2654|language=en|doi=10.1063/1.480839|issn=0021-9606|arxiv=physics/9911005}}
** {{Cite journal|last=Harju|first=A.|last2=Barbiellini|first2=B.|last3=Siljamäki|first3=S.|last4=Nieminen|first4=R. M.|last5=Ortiz|first5=G.|date=1997-08-18|title=Stochastic Gradient Approximation: An Efficient Method to Optimize Many-Body Wave Functions|url=https://link.aps.org/doi/10.1103/PhysRevLett.79.1173|journal=Physical Review Letters|volume=79|issue=7|pages=1173–1177|doi=10.1103/PhysRevLett.79.1173}}
** {{Cite journal|last=Tanaka|first=Shigenori|date=1994-05-15|title=Structural optimization in variational quantum Monte Carlo|url=https://aip.scitation.org/doi/10.1063/1.466885|journal=The Journal of Chemical Physics|volume=100|issue=10|pages=7416–7420|language=en|doi=10.1063/1.466885|issn=0021-9606}}
** {{Cite journal|last=Casula|first=Michele|last2=Attaccalite|first2=Claudio|last3=Sorella|first3=Sandro|date=2004-10-15|title=Correlated geminal wave function for molecules: An efficient resonating valence bond approach|url=https://aip.scitation.org/doi/10.1063/1.1794632|journal=The Journal of Chemical Physics|volume=121|issue=15|pages=7110–7126|language=en|doi=10.1063/1.1794632|issn=0021-9606|arxiv=cond-mat/0409644}}
** {{Cite journal|last=Drummond|first=N. D.|last2=Needs|first2=R. J.|date=2005-08-18|title=Variance-minimization scheme for optimizing Jastrow factors|url=https://link.aps.org/doi/10.1103/PhysRevB.72.085124|journal=Physical Review B|volume=72|issue=8|pages=085124|doi=10.1103/PhysRevB.72.085124}}

== 関連項目 ==
* {{仮リンク|時間依存変分モンテカルロ法|en|Time-dependent_variational_Monte_Carlo}}：変分モンテカルロ法を[[量子状態|純粋量子状態]]の動力学を扱えるように拡張したもの。

{{デフォルトソート:へんふんもんてかるろほう}}
[[Category:最適化]]
[[Category:量子化学]]