蒙特卡罗模拟中基于双向链表的元胞链表方法
PDF下载 (230)王逸梅,王少云,童朝晖.蒙特卡罗模拟中基于双向链表的元胞链表方法[J].宁波大学学报(理工版),2021,34(4):86-92.DOI:
WANG Yimei,WANG Shaoyun,TONG Chaohui.Cell lists method based on doubly linked lists for Monte Carlo simulation[J].Journal of Ningbo University(Natural Science & Engineering Edition),2021,34(4):86-92.DOI:
| Title: | Cell lists method based on doubly linked lists for Monte Carlo simulation |
| 作者: | 王逸梅, 王少云, 童朝晖 |
| Author(s): | WANG Yimei, WANG Shaoyun, TONG Chaohui |
| 关键词: | 蒙特卡罗方法; 非局域移动; 双向链表; 元胞链表方法 |
| Keywords: | Monte Carlo method; nonlocal move; cell lists method; doubly linked lists |
| 分类号: | O411.3 |
| 文献标识码: | A |
| 摘要: | 较之分子动力学, 蒙特卡罗能够实现非局域的粒子移动, 从而解决一些分子动力学不容易模拟的问题. 非局域的粒子移动主要包括模拟化学反应时粒子产生和消失的过程, 高分子模拟时的扭折-跳跃、绕枢轴转动和蠕动以及位形偏倚蒙特卡罗中链的回溯和再生. 然而在蒙特卡罗方法处理非局域移动时, 并不存在一种计算短程作用的计算复杂度为 的算法, 从而限制了蒙特卡罗方法的应用. 本文基于双向链表的数据结构, 发展了蒙特卡罗模拟中因粒子删除和插入而引起的短程势能变化的计算复杂度为 的元胞链表方法. 所有非局域的粒子移动可以转化为粒子的删除和插入, 因此该方法适用于上述所有情形. 此外, 由于Metropolis算法中给某粒子一个随机位移的过程可以看成旧位置粒子的删除以及新位置粒子的插入, 因此该方法也适用于Metropolis算法中粒子的随机移动. |
| Abstract: | Compared with molecular dynamics, Monte Carlo method involves some nonlocal moves which accelerate the simulation. These nonlocal moves include kink-jump, pivot, reptation move in polymer simulation, the retrace and regrowth of chains in configurational biased Monte Carlo, and the breaking and generating of bonds in chemical reaction and acid/base ionization. However, there lacks an algorithm with computational complexity in calculation of the short-range potential because Verlet lists method fails for nonlocal moves. In this paper, a cell lists method based on doubly linked lists and with complexity is developed for particle deletion and insertion in Monte Carlo simulation. Because the random moves in Metropolis algorithm can be reduced to particle deletion at old position and particle insertion at new position, this method can also be used in Metropolis algorithm. In addition, the method is verified by comparing the simulation results obtained using the presented method with the one using Verlet lists method. Finally, the space and time complexity of the proposed method is also discussed. |
| 参考文献 /References: | [1] Frenkel D, Smit B. Understanding Molecular Simulation: From Algorithms to Applications[M]. 2nd ed. San Diego: Academic, 2002. [2] Allen M P, Tildesley D J. Computer Simulation of Liquids[M]. 2nd ed. Oxford: Oxford University Press, 2017. [3] Hockney R W, Eastwood J W. Computer Simulation Using Particles[M]. New York: McGraw-Hill, 1981. [4] Heath Turner C, Brennan J K, Lísal M, et al. Simulation of chemical reaction equilibria by the reaction ensemble Monte Carlo method: A review[J]. Molecular Simulation, 2008, 34:119-146. [5] Reed C E, Reed W F. Monte Carlo study of titration of linear polyelectrolytes[J]. The Journal of Chemical Physics, 1992, 96(2):1609-1620. [6] Landsgesell J, Holm C, Smiatek J. Simulation of weak polyelectrolytes: A comparison between the constant pH and the reaction ensemble method[J]. The European Physical Journal Special Topics, 2017, 226:725-736. [7] Carmesin I, Kremer K. The bond fluctuation method: A new effective algorithm for the dynamics of polymers in all spatial dimensions[J]. Macromolecules, 1988, 21:2819-2823. [8] Kremer K, Binder K. Monte Carlo simulation of lattice models for macromolecules[J]. Computer Physics Reports, 1988, 7:259-310. [9] Sadus R J. Molecular Simulation of Liquids: Theory, Algorithms and Objection-Orientation[M]. Amsterdam: Elsevier, 2002. [10] Siepmann J I, Frenkel D. Configurational bias Monte Carlo: A new sampling scheme for flexible chains[J]. Molecular Physics, 1992, 75:59-70. [11] Mazzeo M D, Ricci M, Zannoni C. The linked neighbour list (LNL) method for fast off-lattice Monte Carlo simulations of fluids[J]. Computer Physics Communications, 2010, 181:569-581. [12] Drozdek A. Data Structures and Algorithms in C++[M]. 4th ed. Boston: Cengage Learning, 2013. [13] Welling U, Germano G. Efficiency of linked cell algorithms[J]. Computer Physics Communications, 2011, 182:611-615. [14] Heinz T N, Hünenberger P H. A fast pairlist-construction algorithm for molecular simulations under periodic boundary conditions[J]. Journal of Computational Chemistry, 2004, 25:1474-1486. [15] Gonnet P. A simple algorithm to accelerate the computation of non-bonded interactions in cell-based molecular dynamics simulations[J]. Journal of Computational Chemistry, 2007, 28:570-573. [16] Mattson W, Rice B M. Near-neighbor calculations using a modified cell-linked list method[J]. Computer Physics Communication, 1999, 119:135-148. |
| 备注/Memo: | 收稿日期: 2020－12－08. 宁波大学学报（理工版）网址: http://journallg.nbu.edu.cn/ 基金项目: 国家自然科学基金（21774067）. 第一作者: 王逸梅（1994－）, 女, 安徽淮北人, 在读硕士研究生, 主要研究方向: 高分子物理. E-mail: wangyimei0817@163.com *通信作者: 童朝晖（1968－）, 男, 湖南常德人, 教授, 主要研究方向: 聚电解质理论. E-mail: tongchaohui@nbu.edu.cn 宁波大学学报（理工版）网址：http://journallg.nbu.edu.cn/ |
