Trément Sébastien, Schnell Benoît, Petitjean Laurent, Couty Marc, Rousseau Bernard
Laboratoire de Chimie-Physique, UMR 8000 CNRS, Université Paris-Sud, Orsay, France.
Manufacture Française des Pneumatiques MICHELIN, Centre de Ladoux, 23 place des Carmes, 63000 Clermont-Ferrand, France.
J Chem Phys. 2014 Apr 7;140(13):134113. doi: 10.1063/1.4870394.
We apply operational procedures available in the literature to the construction of coarse-grained conservative and friction forces for use in dissipative particle dynamics (DPD) simulations. The full procedure rely on a bottom-up approach: large molecular dynamics trajectories of n-pentane and n-decane modeled with an anisotropic united atom model serve as input for the force field generation. As a consequence, the coarse-grained model is expected to reproduce at least semi-quantitatively structural and dynamical properties of the underlying atomistic model. Two different coarse-graining levels are studied, corresponding to five and ten carbon atoms per DPD bead. The influence of the coarse-graining level on the generated force fields contributions, namely, the conservative and the friction part, is discussed. It is shown that the coarse-grained model of n-pentane correctly reproduces self-diffusion and viscosity coefficients of real n-pentane, while the fully coarse-grained model for n-decane at ambient temperature over-predicts diffusion by a factor of 2. However, when the n-pentane coarse-grained model is used as a building block for larger molecule (e.g., n-decane as a two blobs model), a much better agreement with experimental data is obtained, suggesting that the force field constructed is transferable to large macro-molecular systems.
我们将文献中可用的操作程序应用于构建粗粒度保守力和摩擦力,以用于耗散粒子动力学(DPD)模拟。完整的程序依赖于自下而上的方法:用各向异性联合原子模型建模的正戊烷和正癸烷的大分子动力学轨迹用作生成力场的输入。因此,预计粗粒度模型至少能半定量地再现基础原子模型的结构和动力学性质。研究了两种不同的粗粒度水平,分别对应每个DPD珠子包含五个和十个碳原子。讨论了粗粒度水平对生成的力场贡献(即保守部分和摩擦部分)的影响。结果表明,正戊烷的粗粒度模型正确地再现了真实正戊烷的自扩散系数和粘度系数,而在环境温度下正癸烷的完全粗粒度模型对扩散的预测高估了两倍。然而,当将正戊烷粗粒度模型用作更大分子(例如,将正癸烷作为双团块模型)的构建块时,与实验数据的吻合度要好得多,这表明构建的力场可转移到大分子系统。