レプリカ交換法

ウィキペディア

レプリカ交換法

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2024/08/02 04:08 UTC 版)

レプリカ交換法（レプリカこうかんほう、英: replica exchange method、レプリカ交換MCMCサンプリング）はパラレルテンパリング（parallel tempering、並列焼きもどし）法としても知られ、モンテカルロシミュレーションやマルコフ連鎖モンテカルロ法(MCMC)のサンプリング効率を改善するための方法である。SwendsenとWang^[1]によって開発され、Geyer^[2]によって拡張され、その後、特に、福島・根本^[3]およびジョルジョ・パリージ^[4][5]によって発展した。杉田と岡本はパラレルテンパリングの分子動力学法版を考案した^[6]。これはレプリカ交換分子動力学（replica-exchange molecular dynamics、REMD）として知られている。

手法としては、始めに異なる温度でランダムに初期化された N 個の系のコピーを走らせ、メトロポリス法の基準でそれぞれ温度間で系の状態を交換するものである。

この方法の重要な点は、低温のシミュレーションで高温の設定が（またはその逆も）できることである。低エネルギー配置と高エネルギー配置の両方をサンプリングできるため、とても安定にかつ正確なシミュレーションを行うことができる。このようにして、正準集団では一般にうまく計算されない比熱といった熱力学特性がかなり正確に計算できる。

詳細

メトロポリス・ヘイスティングス法を用いるモンテカルロ法は1つの系の確率過程によって成り立っており、採択と棄却に影響する温度 T は1つである。温度が高い場合、よりエネルギーの差分が大きい更新でも採択される確率が比較的高い。系の相関が高い場合は棄却される確率が高く、致命的なシミュレーション速度の低下が起きると言われる。

温度が ΔT 離れた2つの系のシミュレーションが走っているとして、ΔT が十分小さければ、それぞれ N 回のモンテカルロステップを行って得られたエネルギーをヒストグラムにして得られるサンプル分布はいくらか重なる。分布の重なりは正規化されたサンプル分布のヒストグラムが重なる面積で量ることができる。ΔT → 0 において重なりは 1 に近づく。

他の定義の仕方として、T₁ の系におけるサンプルが T₂ においてどれだけ出現しうるかということもできる。

マルコフ連鎖は過去の記憶を持たないために、マルコフ連鎖が2つの温度の系で成り立っていてもよい。

モンテカルロ法においてそれぞれの系が持つ温度をメトロポリスヘイスティングス基準で交換する。交換確率を p とすると以下のようになり、

${\displaystyle p=\min \left(1,{\frac {\exp \left(-{\frac {E_{j}}{kT_{i}}}-{\frac {E_{i}}{kT_{j}}}\right)}{\exp \left(-{\frac {E_{i}}{kT_{i}}}-{\frac {E_{j}}{kT_{j}}}\right)}}\right)=\min \left(1,\exp \left\{{(E_{i}-E_{j})\left({\frac {1}{kT_{i}}}-{\frac {1}{kT_{j}}}\right)}\right\}\right)}$

1 − p で棄却する。

詳細釣り合い条件は満たされている^[7]。

この方法は2つ以上の系にも一般化できる。温度と系の数を適切に決定することで、追加でかかるモンテカルロ法の計算コストを上回る性能向上を達成することができる。レプリカ交換法は、緩和が遅くかつ極小値に陥いりやすい焼きなまし法を改善することができる。

参考文献

1. ^ Swendsen RH, Wang JS (1986). “Replica Monte Carlo simulation of spin glasses”. Phys. Rev. Lett. 57: 2607-2609. doi:10.1103/PhysRevLett.57.2607. PMID 10033814. 
2. ^ Geyer, C. J. (1991). “Markov Chain Monte Carlo Maximum Likelihood”. In Keramidas, E. M.. Computing Science and Statistics. Proceedings of the 23rd Symposium on the Interface. New York: American Statistical Association. pp. 156-163. http://www.stat.umn.edu/geyer/f05/8931/c.pdf 
3. ^ Hukushima, K. & Nemoto, K. (1996). “Exchange Monte Carlo method and application to spin glass simulations”. J. Phys. Soc. Jan. 65 (6): 1604-1608. doi:10.1143/JPSJ.65.1604. 
4. ^ Falcioni, M.; Deem, M. W. (1999). “A Biased Monte Carlo Scheme for Zeolite Structure Solution”. J. Chem. Phys. 110 (3): 1754. arXiv:cond-mat/9809085. Bibcode: 1999JChPh.110.1754F. doi:10.1063/1.477812. 
5. ^ Earl, D. J.; Deem, M. W. (2005). “Parallel tempering: Theory, applications, and new perspectives”. Phys. Chem. Chem. Phys. 7: 3910. doi:10.1039/b509983h. 
6. ^ Sugita, Y.; Okamoto, Y. (1999). “Replica-exchange molecular dynamics method for protein folding”. Chemical Physics Letters 314: 141–151. doi:10.1016/S0009-2614(99)01123-9. 
7. ^ Radford M. Neal (1996). “Sampling from multimodal distributions using tempered transitions”. Statistics and Computing 6 (4): 353–366. doi:10.1007/BF00143556.