Gao Peiyuan, Guo Hongxia
Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
Phys Chem Chem Phys. 2015 Dec 21;17(47):31693-706. doi: 10.1039/c5cp05017k.
In this paper, we evaluate the transferability of the coarse-grained (CG) force field (FF) of trans-1,4-polybutadiene which was built via a combined structure-based and thermodynamic quantity-based CG method at 413 K and 1 atm by systematically examining CG simulated structural and thermodynamic properties against the underlying atomistic simulation results at different temperatures. Interestingly, the derived CG force field exhibits good "state-point transferability" to some extent. For example, when applying this CG FF to the nearby state point (e.g., amorphous phase at 500 K), the resulting local conformation statistics, chain size, and local packing properties as well as density values for the CG models closely match the atomistic simulated data. When further applying this CG force field to the crystalline state at 300 K, the structural and thermodynamic properties of the crystalline phase formed for these CG and atomistic MD simulations still match within a certain level of accuracy. Furthermore, the CG torsion potential has a dual effect: for the amorphous state, the presence of an intramolecular energy barrier against rotation improves the capability of CG models to more precisely reproduce the structural properties, while in the crystalline state this CG torsion barrier suppresses the formation of the more fully stretched chain with a higher trans content. As a result, in the crystalline phase the CG model chains without torsion potentials possess a more stretched chain conformation, pack more efficiently, and have a higher crystallinity degree than its counterpart with CG torsion potentials as well as its underlying atomistic model. However, the dual effect of CG torsion potentials does not mean that we have to use different dihedral parameters to describe different state points. Both CG FFs, one with and another without torsion potentials, are able to represent the melt and the crystalline states. Overall, the phase and its structural consistency between the CG and atomistic models over other state points (e.g., crystalline phase) for which CG FFs were not explicitly parameterized very encouraging such that the combined structure-based and thermodynamic quantity-based CG method can be used to derive an optimized CG FF for multi-scale simulation of polymer systems under different thermodynamic conditions.
在本文中,我们评估了通过基于结构和基于热力学量的组合粗粒化(CG)方法在413 K和1 atm下构建的反式1,4 - 聚丁二烯的粗粒化力场(FF)的可转移性,通过系统地将CG模拟的结构和热力学性质与不同温度下的基础原子模拟结果进行对比。有趣的是,所推导的CG力场在一定程度上表现出良好的“状态点可转移性”。例如,当将此CG FF应用于附近的状态点(如500 K的非晶相)时,CG模型的局部构象统计、链尺寸、局部堆积性质以及密度值与原子模拟数据紧密匹配。当进一步将此CG力场应用于300 K的晶态时,这些CG和原子分子动力学模拟形成的晶相的结构和热力学性质在一定精度水平内仍然匹配。此外,CG扭转势具有双重作用:对于非晶态,分子内旋转能垒的存在提高了CG模型更精确再现结构性质的能力,而在晶态中,此CG扭转势垒抑制了具有更高反式含量的更完全伸展链的形成。结果,在晶相中,没有扭转势的CG模型链具有更伸展的链构象,堆积更有效,并且比具有CG扭转势的对应物及其基础原子模型具有更高的结晶度。然而,CG扭转势的双重作用并不意味着我们必须使用不同的二面角参数来描述不同的状态点。有无扭转势的两种CG FF都能够代表熔体和晶态。总体而言,对于未明确参数化CG FF的其他状态点(如晶相),CG和原子模型之间的相及其结构一致性非常令人鼓舞,使得基于结构和基于热力学量的组合CG方法可用于推导用于不同热力学条件下聚合物系统多尺度模拟的优化CG FF。