Massively parallelized replica-exchange simulations of polymers on GPUs

Article history: We discuss the advantages of parallelization by multithreading on graphics processing units (GPUs) for parallel tempering Monte Carlo computer simulations of an exemplified bead-spring model for homopolymers. Since the sampling of a large ensemble of conformations is a prerequisite for the precise estimation of statistical quantities such as typical indicators for conformational transitions like the peak structure of the specific heat, the advantage of a strong increase in performance of Monte Carlo simulations cannot be overestimated. Employing multithreading and utilizing the massive power of the large number of cores on GPUs, being available in modern but standard graphics cards, we find a rapid increase in efficiency when porting parts of the code from the central processing unit (CPU) to the GPU. © 2011 Elsevier B.V. All rights reserved.

[1]  Tom R. Halfhill.  NVIDIA's Next-Generation CUDA Compute and Graphics Architecture, Code-Named Fermi, Adds Muscle for Parallel Processing , 2009 .

[2]  Jeffrey S. Vetter,et al.  Accuracy and performance of graphics processors: A Quantum Monte Carlo application case study , 2009, Parallel Comput..

[3]  Wolfhard Janke,et al.  Thermodynamics of lattice heteropolymers. , 2003, The Journal of chemical physics.

[4]  Junqi Yin,et al.  Phase diagram and critical behavior of the square-lattice Ising model with competing nearest-neighbor and next-nearest-neighbor interactions. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  A. Arnold,et al.  Harvesting graphics power for MD simulations , 2007, 0709.3225.

[6]  Tobias Preis,et al.  Fluctuation patterns in high-frequency financial asset returns , 2008 .

[8]  Wolfhard Janke,et al.  Multicanonical study of coarse-grained off-lattice models for folding heteropolymers. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  Steve B. Jiang,et al.  GPU-based ultrafast IMRT plan optimization , 2009, Physics in medicine and biology.

[10]  Wolfhard Janke,et al.  Elastic Lennard-Jones polymers meet clusters: differences and similarities. , 2009, The Journal of chemical physics.

[11]  Martin Weigel,et al.  Simulating spin models on GPU , 2010, Comput. Phys. Commun..

[12]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[13]  Jie Cheng,et al.  Programming Massively Parallel Processors. A Hands-on Approach , 2010, Scalable Comput. Pract. Exp..

[14]  Steve B Jiang,et al.  GPU-based ultra-fast dose calculation using a finite size pencil beam model. , 2009, Physics in medicine and biology.

[15]  Vijay S. Pande,et al.  Accelerating molecular dynamic simulation on graphics processing units , 2009, J. Comput. Chem..

[16]  Sanford Weisberg,et al.  Computing science and statistics : proceedings of the 30th Symposium on the Interface, Minneapolis, Minnesota, May 13-16, 1998 : dimension reduction, computational complexity and information , 1998 .

[17]  D. Landau,et al.  Collapse transitions in a flexible homopolymer chain: application of the Wang-Landau algorithm. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[18]  Wolfgang Paul,et al.  GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model , 2009, J. Comput. Phys..

[19]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[20]  Lister Staveley-Smith,et al.  GPU accelerated radio astronomy signal convolution , 2008 .

[21]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[22]  Wolfhard Janke,et al.  Multicanonical chain-growth algorithm. , 2003, Physical review letters.

[23]  Wolfhard Janke,et al.  Surface effects in the crystallization process of elastic flexible polymers , 2009, 1002.2118.

[24]  Peter J. Cameron,et al.  Galois groups of multivariate Tutte polynomials , 2010, 1006.3869.

[25]  Eric B. Ford.  Parallel algorithm for solving Kepler’s equation on Graphics Processing Units: Application to analysis of Doppler exoplanet searches , 2009 .