Mapping continuous potentials to discrete forms.

The optimal conversion of a continuous inter-particle potential to a discrete equivalent is considered here. Existing and novel algorithms are evaluated to determine the best technique for creating accurate discrete forms using the minimum number of discontinuities. This allows the event-driven molecular dynamics technique to be efficiently applied to the wide range of continuous force models available in the literature, and facilitates a direct comparison of event-driven and time-driven molecular dynamics. The performance of the proposed conversion techniques are evaluated through application to the Lennard-Jones model. A surprising linear dependence of the computational cost on the number of discontinuities is found, allowing accuracy to be traded for speed in a controlled manner. Excellent agreement is found for static and dynamic properties using a relatively low number of discontinuities. For the Lennard-Jones potential, the optimized discrete form outperforms the original continuous form at gas densities but is significantly slower at higher densities.

[1]  Athanassios Z. Panagiotopoulos,et al.  Phase behavior of the restricted primitive model and square-well fluids from Monte Carlo simulations in the grand canonical ensemble , 1999 .

[2]  Marcus N. Bannerman,et al.  DynamO: a free ${\cal O}$(N) general event‐driven molecular dynamics simulator , 2010, Journal of computational chemistry.

[3]  Mookyung Cheon,et al.  Spontaneous formation of twisted Aβ(16-22) fibrils in large-scale molecular-dynamics simulations. , 2011, Biophysical journal.

[4]  B. Smit,et al.  Phase diagrams of Lennard‐Jones fluids , 1992 .

[5]  Molecular dynamics for discontinuous potential. IV. Lennard-Jonesium , 1989 .

[6]  Jeremy Schofield,et al.  Event-driven dynamics of rigid bodies interacting via discretized potentials. , 2007, The Journal of chemical physics.

[7]  Stephan Kabelac,et al.  Transport coefficients of the Lennard-Jones model fluid. II Self-diffusion. , 2004, The Journal of chemical physics.

[8]  Hans-Jörg Limbach,et al.  ESPResSo - an extensible simulation package for research on soft matter systems , 2006, Comput. Phys. Commun..

[9]  Dennis J. Miller,et al.  Application of the Step Potential for Equilibria and Dynamics (SPEAD) Method to Bioderived Esters and Acetals , 2012 .

[10]  T. Poeschel,et al.  Event-driven molecular dynamics of soft particles , 2013, 1301.3543.

[11]  F. Ding,et al.  Scaling behavior and structure of denatured proteins. , 2005, Structure.

[12]  Wilding Critical-point and coexistence-curve properties of the Lennard-Jones fluid: A finite-size scaling study. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[13]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[14]  J. Richard Elliott,et al.  Optimization of Transferable Site–Site Potentials Using a Combination of Stochastic and Gradient Search Algorithms , 2012 .

[15]  S. Buldyrev,et al.  Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model. , 2004, Biophysical journal.

[16]  C. Hall,et al.  Phase diagrams describing fibrillization by polyalanine peptides. , 2004, Biophysical journal.

[17]  J. Richard Elliott,et al.  Correlation of Mixture Vapor-Liquid Equilibria with the SPEADMD Model , 2008 .

[18]  B. Alder,et al.  Phase Transition for a Hard Sphere System , 1957 .

[19]  Stephan Kabelac,et al.  Transport coefficients of the Lennard-Jones model fluid. I. Viscosity. , 2004, The Journal of chemical physics.

[20]  Nikolay V Dokholyan,et al.  Studies of folding and misfolding using simplified models. , 2006, Current opinion in structural biology.

[21]  T. G. Cowling,et al.  The mathematical theory of non-uniform gases , 1939 .

[22]  J. Richard Elliott,et al.  Phase diagrams for a multistep potential model of n-alkanes by discontinuous molecular dynamics and thermodynamic perturbation theory , 2002 .

[23]  J. Barker,et al.  Perturbation Theory and Equation of State for Fluids: The Square‐Well Potential , 1967 .

[24]  Leo Lue,et al.  Exact on-event expressions for discrete potential systems. , 2010, The Journal of chemical physics.

[25]  Carol K. Hall,et al.  Molecular dynamics simulations of DPPC bilayers using "LIME", a new coarse-grained model. , 2013, The journal of physical chemistry. B.

[26]  Transferable step potentials for amines, amides, acetates, and ketones , 2005 .

[27]  G. A. Chapela,et al.  Discrete perturbation theory for the hard-core attractive and repulsive Yukawa potentials. , 2010, The Journal of chemical physics.

[28]  J. Barker,et al.  Perturbation Theory and Equation of State for Fluids. II. A Successful Theory of Liquids , 1967 .

[29]  Thorsten Pöschel,et al.  Computational Granular Dynamics , 2005 .

[30]  Paolo Valentini,et al.  A combined Event-Driven/Time-Driven molecular dynamics algorithm for the simulation of shock waves in rarefied gases , 2009, J. Comput. Phys..

[31]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[32]  J. Elliott,et al.  Transferable step potentials for perfluorinated hydrocarbons , 2008 .

[33]  J. Elliott,et al.  Transferable Step Potentials for the Straight-Chain Alkanes, Alkenes, Alkynes, Ethers, and Alcohols , 2004 .

[34]  G. Galliéro,et al.  Thermal conductivity of the Lennard-Jones fluid: An empirical correlation , 2008 .

[35]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[36]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[37]  S. Miller,et al.  Event-driven molecular dynamics in parallel , 2004 .

[38]  J. M. Haile,et al.  Molecular dynamics simulation : elementary methods / J.M. Haile , 1992 .

[39]  N. B. Wilding Computer simulation of fluid phase transitions , 2001 .

[40]  A. Galindo,et al.  Application of the generalised SAFT-VR approach for long-ranged square-well potentials to model the phase behaviour of real fluids , 2009 .

[41]  J. Alejandre,et al.  Liquid-vapor phase diagram and surface properties in oppositely charged colloids represented by a mixture of attractive and repulsive Yukawa potentials. , 2013, The Journal of chemical physics.

[42]  Ana Laura Benavides,et al.  Discrete perturbation theory applied to Lennard-Jones and Yukawa potentials. , 2010, The Journal of chemical physics.