14
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Assessing numerical methods for molecular and particle simulation

      Preprint
      , ,

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          We discuss the design of state-of-the-art numerical methods for molecular dynamics, focusing on the demands of soft matter simulation, where the purposes include sampling and dynamics calculations both in and out of equilibrium. We discuss the characteristics of different algorithms, including their essential conservation properties, the convergence of averages, and the accuracy of numerical discretizations. Formulations of the equations of motion which are suited to both equilibrium and nonequilibrium simulation include Langevin dynamics, dissipative particle dynamics (DPD), and the more recently proposed "pairwise adaptive Langevin" (PAdL) method, which, like DPD but unlike Langevin dynamics, conserves momentum and better matches the relaxation rate of orientational degrees of freedom. PAdL is easy to code and suitable for a variety of problems in nonequilibrium soft matter modeling, our simulations of polymer melts indicate that this method can also provide dramatic improvements in computational efficiency. Moreover we show that PAdL gives excellent control of the relaxation rate to equilibrium. In the nonequilibrium setting, we further demonstrate that while PAdL allows the recovery of accurate shear viscosities at higher shear rates than are possible using the DPD method at identical timestep, it also outperforms Langevin dynamics in terms of stability and accuracy at higher shear rates.

          Related collections

          Most cited references60

          • Record: found
          • Abstract: found
          • Article: found
          Is Open Access

          Canonical sampling through velocity-rescaling

          We present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. We illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent on the thermostat parameter also with regard to the dynamic properties.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: not found

            Langevin dynamics of peptides: the frictional dependence of isomerization rates of N-acetylalanyl-N'-methylamide.

            The rate constant for the transition between the equatorial and axial conformations of N-acetylalanyl-N'-methylamide has been determined from Langevin dynamics (LD) simulations with no explicit solvent. The isomerization rate is maximum at collision frequency gamma = 2 ps-1, shows diffusive character for gamma greater than or equal to 10 ps-1, but does not approach zero even at gamma = 0.01 ps-1. This behavior differs from that found for a one-dimensional bistable potential and indicates that both collisional energy transfer with solvent and vibrational energy transfer between internal modes are important in the dynamics of barrier crossing for this system. It is suggested that conformational searches of peptides be carried out using LD with a collision frequency that maximizes the isomerization rate (i.e., gamma approximately 2 ps-1). This method is expected to be more efficient than either molecular dynamics in vacuo (which corresponds to LD with gamma = 0) or molecular dynamics in solvent (where dynamics is largely diffusive).
              Bookmark
              • Record: found
              • Abstract: not found
              • Article: not found

              The computer study of transport processes under extreme conditions

                Bookmark

                Author and article information

                Journal
                03 November 2017
                Article
                10.1039/C7SM01526G
                1711.01310
                4840129f-0148-4915-adbc-866a94931f4f

                http://arxiv.org/licenses/nonexclusive-distrib/1.0/

                History
                Custom metadata
                Soft Matter (2017)
                physics.comp-ph cond-mat.soft

                Comments

                Comment on this article