We calculate current correlation functions (CCFs) of dissipative particle dy- namics (DPD) and compare them with results of molecular dynamics (MD) and solutions of linearized hydrodynamic equations. In particul...We calculate current correlation functions (CCFs) of dissipative particle dy- namics (DPD) and compare them with results of molecular dynamics (MD) and solutions of linearized hydrodynamic equations. In particular, we consider three versions of DPD, the empirical/classical DPD, coarse-grained (CG) DPD with radial-direction interactions only and full (radial, transversal, and rotational) interactions between particles. To fa- cilitate quantitative discussions, we consider specifically a star-polymer melt system at a moderate density. For bonded molecules, it is straightforward to define the CG variables and to further derive CG force fields for DPD within the framework of the Mori-Zwanzig formalism. For both transversal and longitudinal current correlation functions (TCCFs and LCCFs), we observe that results of MD, DPD, and hydrodynamic solutions agree with each other at the continuum limit. Below the continuum limit to certain length scales, results of MD deviate significantly from hydrodynamic solutions, whereas results of both empirical and CG DPD resemble those of MD. This indicates that the DPD method with Markovian force laws possibly has a larger applicability than the continuum description of a Newtonian fluid. This is worth being explored further to represent gen- eralized hydrodynamics.展开更多
基金funding support of the U.S.Army Research Laboratory with Cooperative Agreement No.W911NF-12-2-0023
文摘We calculate current correlation functions (CCFs) of dissipative particle dy- namics (DPD) and compare them with results of molecular dynamics (MD) and solutions of linearized hydrodynamic equations. In particular, we consider three versions of DPD, the empirical/classical DPD, coarse-grained (CG) DPD with radial-direction interactions only and full (radial, transversal, and rotational) interactions between particles. To fa- cilitate quantitative discussions, we consider specifically a star-polymer melt system at a moderate density. For bonded molecules, it is straightforward to define the CG variables and to further derive CG force fields for DPD within the framework of the Mori-Zwanzig formalism. For both transversal and longitudinal current correlation functions (TCCFs and LCCFs), we observe that results of MD, DPD, and hydrodynamic solutions agree with each other at the continuum limit. Below the continuum limit to certain length scales, results of MD deviate significantly from hydrodynamic solutions, whereas results of both empirical and CG DPD resemble those of MD. This indicates that the DPD method with Markovian force laws possibly has a larger applicability than the continuum description of a Newtonian fluid. This is worth being explored further to represent gen- eralized hydrodynamics.