Simulation of ring polymer melts with GPU acceleration

Abstract We implemented the elastic lattice polymer model on the GPU (Graphics Processing Unit), and show that the GPU is very efficient for polymer simulations of dense polymer melts. The implementation is able to perform up to 4.1 ⋅ 10 9 Monte Carlo moves per second. Compared to our standard CPU implementation, we find an effective speed-up of a factor 92. Using this GPU implementation we studied the equilibrium properties and the dynamics of non-concatenated ring polymers in a melt of such polymers, using Rouse modes. With increasing polymer length, we found a very slow transition to compactness with a growth exponent ν ≈ 1 / 3 . Numerically we find that the longest internal time scale of the polymer scales as N 3.1 , with N the molecular weight of the ring polymer.

[1]  Raoul D. Schram,et al.  Reaction-diffusion model Monte Carlo simulations on the GPU , 2013, J. Comput. Phys..

[2]  Debabrata Panja,et al.  Passage times for unbiased polymer translocation through a narrow pore. , 2005, Physical review letters.

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

[4]  Ralf Everaers,et al.  Viscoelasticity and primitive path analysis of entangled polymer liquids: from F-actin to polyethylene. , 2007, The Journal of chemical physics.

[5]  G. Barkema,et al.  Rouse modes of self-avoiding flexible polymers. , 2009, The Journal of chemical physics.

[6]  T. Sakaue.  Ring polymers in melts and solutions: scaling and crossover. , 2011, Physical review letters.

[7]  Ye Zhao,et al.  Lattice Boltzmann based PDE solver on the GPU , 2008, The Visual Computer.

[8]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[9]  Kurt Kremer,et al.  Molecular dynamics simulation study of nonconcatenated ring polymers in a melt. I. Statics. , 2011, The Journal of chemical physics.

[10]  A. Grosberg.  Annealed lattice animal model and Flory theory for the melt of non-concatenated rings: towards the physics of crumpling. , 2013, Soft matter.

[11]  Y. Matsushita,et al.  Topological effect in ring polymers investigated with Monte Carlo simulation. , 2008, The Journal of chemical physics.

[12]  Rubinstein.  Discretized model of entangled-polymer dynamics. , 1987, Physical review letters.

[13]  Reaching large lengths and long times in polymer dynamics simulations , 2003, cond-mat/0306280.

[14]  Marko,et al.  Electrophoresis of charged polymers: Simulation and scaling in a lattice model of reptation. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[16]  G. Barkema,et al.  Structural modes of a polymer in the repton model. , 2011, The Journal of chemical physics.

[17]  S. Edwards,et al.  The Theory of Polymer Dynamics , 1986 .

[18]  G. Barkema,et al.  Unwinding dynamics of double-stranded polymers. , 2010, The Journal of chemical physics.

[19]  Kurt Kremer,et al.  Statistics of polymer rings in the melt: a numerical simulation study , 2009, Physical biology.

[20]  Michael Lang,et al.  GPU implementations of the bond fluctuation model , 2012, J. Comput. Phys..