This study investigates the efficiency and reliability of the modified Shardlow's (M-Shardlow) method for dissipative particle dynamics (DPD). We show that the M-Shardlow method in which for its construction, the second order velocity Verlet method in the Shardlows method to integrate the Hamiltonian part has been replaced by a symplectic fourth order method, improperly uses some parameters. %In other words, in this paper, it is shown that the initial M-Shardlow method employed some parameters improperly in the fourth order symplectic method that was used for the M-Shardlow method. By numerical experiments and computing, some important configurational quantities such as configurational temperature and radial distribution function (RDF), the M-Shardlow's method is compared with the Shardlow and ABOBA methods. These results indicate that the new method obtained in this way, even with the proper parameters is too costly in the sense of the CPU-time that is required per each step which makes it an inefficient DPD integrator. Besides, by a comparison of the radial distribution function of this method with Shardlow and ABOBA for large time increments, we can observe no considerable improvement in preserving the structure of the system by this new DPD solver.
翻译:暂无翻译