リケラボ論文検索は、全国の大学リポジトリにある学位論文・教授論文を一括検索できる論文検索サービスです。

リケラボ 全国の大学リポジトリにある学位論文・教授論文を一括検索するならリケラボ論文検索大学・研究所にある論文を検索できる

リケラボ 全国の大学リポジトリにある学位論文・教授論文を一括検索するならリケラボ論文検索大学・研究所にある論文を検索できる

大学・研究所にある論文を検索できる 「2次元グランドカノニカルレプリカ交換法による単原子分子の気体、液体、超臨界流体の解析」の論文概要。リケラボ論文検索は、全国の大学リポジトリにある学位論文・教授論文を一括検索できる論文検索サービスです。

コピーが完了しました

URLをコピーしました

論文の公開元へ論文の公開元へ
書き出し

2次元グランドカノニカルレプリカ交換法による単原子分子の気体、液体、超臨界流体の解析

松原, 大貴 名古屋大学

2020.04.02

概要

分子スケールの計算機シミュレーションは、実験環境条件に対応する統計力学的集団(アンサンブル)を再現するように設定された計算を通常は行う。しかし、多自由度系では無数に存在するエネルギー極小状態の一つにシミュレーションが留まってしまい、誤った解を出してしまう問題がある。

この問題を解決する為に、実験環境条件とは異なる人工のアンサンブルに基づく計算も行われてきた。それらの手法の中でも、単一シミュレーション中に平衡状態を保ちながら温度を変化させるレプリカ交換法 (REM) は、その汎用性の高さから広く用いられている。しかし、1 次相転移点を含む系に対しては REM が適合しない事が知られている。そこで、粒子数が変動するグランドカノニカルアンサンブルの Lennard-Jones 粒子系を計算対象とし、手法の有効性を調査した。この系では、少粒子数状態は気体相、多粒子数状態が液体相となり、その境界には 1 次相転移点が存在する。また、それら結んだ線の端点が臨界点と呼ばれる 2 次相転移点となる。さらに、西川他(2000)は CO2、CF3H、H2O 等のX 線小角散乱実験の結果から超臨界流体相においても、臨界数密度よりも低密度側に粒子数揺らぎのピークを持つことを指摘した。

申請者は、多次元レプリカ交換法の枠組みをグランドカノニカルアンサンブルに拡張し、Lennard-Jones 粒子系において温度と化学ポテンシャルを変化させる、2 次元レプリカ交換法を開発し、1 次相転移を含む系の計算を可能にした。これを確かめる為に、気体相・液体相・超臨界流体相の 3 相を含み、さらに 1 次相転移点、2 次相転移点を含む温度、化学ポテンシャル範囲でのシミュレーションを実行した。

シミュレーションの結果、低温・低化学ポテンシャル領域では平均粒子数に不連続性が確認された。この不連続性は同一パラメータ点に二峰性の分布が存在する事を示す為、得られた不連続点は化学ポテンシャル変化に対する 1 次相転移点であることを示唆する。しかし、高温・高化学ポテンシャルになるにつれて、不連続点の段差は徐々に小さくなる事が確認された。すなわち、二峰性の二つの分布が徐々に近づいていき、相図の端点である 2 次相転移点では完全に同化して一つの平坦な分布となる事が確認された。一方、2 次相転移点より高温な領域においても、粒子数揺らぎのピークが相図と繋がったかたちで確認された。この結果により X 線小角散乱実験で示されている超臨界流体相での振舞いを明らかにした。

以上から、設定パラメータ範囲で予想される全ての状態を単一のレプリカ交換シミュレーションから得た。これによりグランドカノニカルアンサンブルへ拡張した多次元レプリカ交換法の相転移系に対する有効性が示された。

参考文献

[1] David E Shaw, Paul Maragakis, Kresten Lindorff-Larsen, Stefano Piana, Ron O Dror, Michael P Eastwood, Joseph A Bank, John M Jumper, John K Salmon, Yibing Shan, et al. Atomic-level characterization of the structural dynamics of proteins. Science, 330(6002):341–346, 2010.

[2] Ayori Mitsutake, Yuji Sugita, and Yuko Okamoto. Generalized-ensemble algo- rithms for molecular simulations of biopolymers. Biopolymers, 60(2):96–123, 2001.

[3] 岡本祐幸. 拡張アンサンブル法. 「計算科学講座」第9巻「超多自由度系の最適化」第2章. 共立出版, 2013. pp. 110-241.

[4] Koji Hukushima and Koji Nemoto. Exchange Monte Carlo method and ap- plication to spin glass simulations. Journal of the Physical Society of Japan, 65(6):1604–1608, 1996.

[5] Yuji Sugita and Yuko Okamoto. Replica-exchange molecular dynamics method for protein folding. Chemical Physics Letters, 314(1-2):141–151, 1999.

[6] Yukito Iba. Extended ensemble Monte Carlo. International Journal of Modern Physics C, 12(05):623–656, 2001.

[7] Yuji Sugita, Akio Kitao, and Yuko Okamoto. Multidimensional replica- exchange method for free-energy calculations. The Journal of Chemical Physics, 113(15):6042–6051, 2000.

[8] Ayori Mitsutake and Yuko Okamoto. From multidimensional replica-exchange method to multidimensional multicanonical algorithm and simulated temper- ing. Physical Review E, 79(4):047701, 2009.

[9] Ayori Mitsutake and Yuko Okamoto. Multidimensional generalized-ensemble algorithms for complex systems. The Journal of Chemical Physics, 130(21):214105, 2009.

[10] Athanassios Z Panagiotopoulos. Molecular simulation of phase coexistence: Finite-size effects and determination of critical parameters for two-and three- dimensional lennard-jones fluids. International Journal of Thermophysics, 15(6):1057–1072, 1994.

[11] Nigel B Wilding. Simulation studies of fluid critical behaviour. Journal of Physics: Condensed Matter, 9(3):585, 1997.

[12] Keiko Nishikawa and Takeshi Morita. Inhomogeneity of molecular distribution in supercritical fluids. Chemical Physics Letters, 316(3-4):238–242, 2000.

[13] Nicholas Metropolis, Arianna W Rosenbluth, Marshall N Rosenbluth, Au- gusta H Teller, and Edward Teller. Equation of state calculations by fast computing machines. The Journal of Chemical Physics, 21(6):1087–1092, 1953.

[14] Lourdes F Vega, Katherine S Shing, and Luis F Rull. A new algorithm for molecular dynamics simulations in the grand canonical ensemble. Molecular Physics, 82(3):439–453, 1994.

[15] Alan M Ferrenberg and Robert H Swendsen. Phase transitions and critical phenomena. Physical Review Letters, 63:1195–1198, 1989.

[16] Shankar Kumar, John M Rosenberg, Djamal Bouzida, Robert H Swendsen, and Peter A Kollman. The weighted histogram analysis method for free-energy cal- culations on biomolecules. i. the method. Journal of Computational Chemistry, 13(8):1011–1021, 1992.

[17] Rupert G Miller. The jackknife-a review. Biometrika, 61(1):1–15, 1974.

[18] Bernd A Berg. Markov chain Monte Carlo simulations and their statistical analysis: with web-based Fortran code. World Scientific, Singapore, 2004.

[19] Ayori Mitsutake, Yuji Sugita, and Yuko Okamoto. Replica-exchange multi- canonical and multicanonical replica-exchange monte carlo simulations of pep- tides. I. Formulation and benchmark test. The Journal of Chemical Physics, 118(14):6664–6675, 2003.

[20] Ayori Mitsutake. Simulated-tempering replica-exchange method for the multi- dimensional version. The Journal of Chemical Physics, 131(9):094105, 2009.

[21] Michael R Shirts and John D Chodera. Statistically optimal analysis of samples from multiple equilibrium states. The Journal of Chemical Physics, 129(12):124105, 2008.

[22] Nigel B Wilding. Computer simulation of fluid phase transitions. American Journal of Physics, 69(11):1147–1155, 2001.

[23] Kurt Binder. Finite size scaling analysis of ising model block distribution func- tions. Zeitschrift fu¨r Physik B Condensed Matter, 43(2):119–140, 1981.

参考文献をもっと見る

全国の大学の
卒論・修論・学位論文

一発検索!

この論文の関連論文を見る