Dissipative particle dynamics thermostat: a novel thermostat for molecular dynamics simulation of liquid crystals with Gay–Berne potential

作     者：Yuting Ouyang Liang Hao Yanping Ma Hongxia Guo 

作者机构：Beijing National Laboratory for Molecular Sciences (BNLMS) State Key Laboratory for Polymer Physics and ChemistryInstitute of Chemistry Chinese Academy of Sciences 

出 版 物：《Science China Chemistry》 (中国科学（化学英文版）)

年 卷 期：2015年第58卷第4期

页      面：694-707页

核心收录：

学科分类：07[理学] 070205[理学-凝聚态物理] 08[工学] 080501[工学-材料物理与化学] 0805[工学-材料科学与工程（可授工学、理学学位）] 0703[理学-化学] 0702[理学-物理学] 

基　　金：support of the National Natural Science Foundation of China(20674093) 

主　　题：molecular dynamics simulation liquid crystal thermostat dissipative particle dynamics Gay-Berne potential 

摘      要：The Gay-Berne （GB） model has been proved to be highly successful in the simulation of liquid crystal phases via both molec- ular dynamics （MD） and nonequilibrium molecular dynamics （NEMD）. However, the conventional thermostats used in the simulations of GB systems, such as Nose-Hoover and Langevin thermostats, have serious shortcomings especially in NEMD simulations. Recently, dissipative particle dynamics （DPD） has established itself as a useful thermostat for soft matter simulations, whereas the application of DPD thermostat in （NE）MD simulations is limited to the spherically isotropic potential models, such as the Lennard-Jones model. Considering the virtues of the DPD thermostat, that is, local, momentum conserved, and Galilean invariant, we extend the DPD thermostat to the non-spherical GB model. It is interesting to find that the translational DPD and rotational DPD thermostats can be used in the GB system independently and both can achieve the thermostatting effects. Also, we compared the performance of the DPD thermostat with other commonly used thermostats in NEMD simulations by investigating the streaming velocity profiles and the dynamics of phase separation in a typical but simple binary GB mixture under shear field. It is revealed that the known virtues of DPD thermostats, such as Galilean invariant, shear velocity profile-unbiased, and unscreened hydrodynamic interactions, are still intact when applying to GB systems. Finally, the appro- priate parameters for the DPD thermostat in the GB system are identified for future investigations.