Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
J Phys Chem B. 2012 Feb 23;116(7):2075-89. doi: 10.1021/jp210764h. Epub 2012 Feb 8.
In this paper, with the aim to establish a rational coarse-grained (CG) model for the 4-cyano-4'-pentylbiphenyl (5CB) molecule, we construct three possible CG models (5P, 6P, and 7P) and then determine the bonded and nonbonded interaction parameters separately. For the intramolecular bonded interactions, the bond and angle distributions of the 5CB bulk phase are used as the target properties. For the nonbonded interactions between CG particles, we combine the structure-based and thermodynamic quantities-based methods for the parametrization of CG interaction potentials and attempt to use several fragment molecular systems to derive the CG nonbonded interaction parameters in order to maintain the transferability of our CG models to some extent. Finally, we fix the optimal nonbonded LJ parameters between CG bead pairs such that the results from CG simulations not only correctly reproduce the experimental density and the nematic LC state at 300 K and 1 atm but also reasonably approximate the local structural properties calculated from the underlying atomistic model. Through comparison of the resulting CG data with target properties, the 6P model is found to be the best one among the three, and then we use this model to investigate the phase behavior and dynamic properties. Our results show that the phase transition temperature from nematic to isotropic phase and the diffusion coefficients are reproduced very well, demonstrating the rationality of the 6P model. Our coarse-grained process should have implications for constructing CG models for nCB series or molecules with similar architectures.
在本文中,我们旨在为 4-氰基-4'-戊基联苯(5CB)分子建立一个合理的粗粒(CG)模型,为此我们构建了三个可能的 CG 模型(5P、6P 和 7P),然后分别确定了键合和非键合相互作用参数。对于分子内的键合相互作用,我们使用 5CB 本体相的键和角度分布作为目标性质。对于 CG 粒子之间的非键相互作用,我们结合基于结构和基于热力学量的方法来参数化 CG 相互作用势能,并尝试使用几个片段分子系统来推导 CG 非键相互作用参数,以在一定程度上保持我们 CG 模型的可转移性。最后,我们固定 CG 珠粒之间的最佳非键 LJ 参数,使得 CG 模拟的结果不仅正确再现了实验密度和 300 K 和 1 atm 下的向列液晶态,而且还合理地近似了从基础原子模型计算得出的局部结构性质。通过将得到的 CG 数据与目标性质进行比较,发现 6P 模型在这三个模型中是最好的,然后我们使用这个模型来研究相行为和动态性质。我们的结果表明,从向列相向各向同性相的转变温度和扩散系数得到了很好的再现,证明了 6P 模型的合理性。我们的粗粒化过程应该对构建 nCB 系列或具有类似结构的分子的 CG 模型具有启示意义。