Scalability and Performance of a Program that Uses Domain Decomposition for Monte Carlo Simulation of Molecular Liquids

被引：0
作者：
Teplukhin A.V. [1 ]
机构：
[1] Institute of Mathematical Problems of Biology of RAS - the Branch of Keldysh Institute of Applied Mathematics of RAS, Pushchino
来源：
Supercomputing Frontiers and Innovations | 2022年 / 9卷 / 03期
关键词：
biopolymers; Monte Carlo; MPI; parallel calculation;
D O I：
10.14529/jsfi220303
中图分类号：
学科分类号：
摘要：
The main factors hindering the development of supercomputer programs for molecular simulation by the Monte Carlo method within the framework of classical physics are considered, and possible ways to eliminate the problems that arise in this case are discussed. Thus, the use of molecular models with moderate stiffness of covalent bonds between fragments makes it possible not only to increase the efficiency of scanning the configuration space of the model, but also to abandon the complex apparatus of kinematics with rigid links, which significantly limits the possibilities of domain decomposition. Based on the domain decomposition strategy and a simplified treatment of the deformation energy of covalent bonds and angles, an original parallel algorithm for calculating the properties of large all-atomic models of aqueous solutions of biopolymers by the Monte Carlo method was developed. To speed up computations within the framework of this approach, each domain is assigned its own group of processors/cores using local data replication and splitting the loop over the interacting partners. The article discusses the logical scheme of the computational algorithm and the main components of the software package (fortran77, MPI 1.2). Test calculations performed for water and n-hexane demonstrated the high performance and scalability of the program in which the proposed algorithm was implemented © The Author 2022. This paper is published with open access at SuperFri.org
引用
收藏
页码：51 / 64
页数：13
相关论文
共 35 条
  • [1] Allen F.H., Kennard O., Watson D.G., Et al., Tables of bond lengths determined by X-ray and neutron diffraction. Part 1. Bond lengths in organic compounds, J. Chem. Soc., Perkin Trans, 2, pp. S1-S19, (1987)
  • [2] Allen M.P., Tildesley D.J., Computer Simulation of Liquids, (1987)
  • [3] Anderson J.A., Jankowski E., Grubbb T.L., Et al., Massively parallel Monte Carlo for many-particle simulations on GPUs, J. Comput. Phys, 254, pp. 27-38, (2013)
  • [4] Berendsen H.J.C., Grigera J.R., Straatsma T.P., The missing term in effective pair potentials, J. Phys. Chem, 91, pp. 6269-6271, (1987)
  • [5] Cornell W.D., Cieplak P., Bayly C.I., Et al., A second generation force field for the simulation of proteins, nucleic acids, and organic molecules, J. Am. Chem. Soc, 117, pp. 5179-5197, (1995)
  • [6] Dubbeldam D., Torres-Knoop A., Walton K.S., On the inner workings of Monte Carlo codes, Mol. Simul, 39, pp. 1253-1292, (2013)
  • [7] Durstenfeld R., Algorithm 235: Random permutation, Commun. ACM, 7, (1964)
  • [8] Go N., Scheraga H.A., Analysis of the contribution of internal vibrations to the statistical weights of equilibrium conformations of macromolecules, J. Chem. Phys, 51, pp. 4751-4767, (1969)
  • [9] Hammersley J.M., Handscomb D.C., Monte Carlo Methods, (1964)
  • [10] Heffelfinger G.S., Lewitt M.E., A comparison between two massively parallel algorithms for Monte Carlo computer simulation: An investigation in the grand canonical ensemble, J. Comput. Chem, 17, pp. 250-265, (1996)
1234