Nicholson David A, Rutledge Gregory C
Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA.
J Chem Phys. 2016 Dec 28;145(24):244903. doi: 10.1063/1.4972894.
Non-equilibrium molecular dynamics is used to study crystal nucleation of n-eicosane under planar shear and, for the first time, uniaxial extension. A method of analysis based on the mean first-passage time is applied to the simulation results in order to determine the effect of the applied flow field type and strain rate on the steady-state nucleation rate and a characteristic growth rate, as well as the effects on kinetic parameters associated with nucleation: the free energy barrier, critical nucleus size, and monomer attachment pre-factor. The onset of flow-enhanced nucleation (FEN) occurs at a smaller critical strain rate in extension as compared to shear. For strain rates larger than the critical rate, a rapid increase in the nucleation rate is accompanied by decreases in the free energy barrier and critical nucleus size, as well as an increase in chain extension. These observations accord with a mechanism in which FEN is caused by an increase in the driving force for crystallization due to flow-induced entropy reduction. At high applied strain rates, the free energy barrier, critical nucleus size, and degree of stretching saturate, while the monomer attachment pre-factor and degree of orientational order increase steadily. This trend is indicative of a significant diffusive contribution to the nucleation rate under intense flows that is correlated with the degree of global orientational order in a nucleating system. Both flow fields give similar results for all kinetic quantities with respect to the reduced strain rate, which we define as the ratio of the applied strain rate to the critical rate. The characteristic growth rate increases with increasing strain rate, and shows a correspondence with the nucleation rate that does not depend on the type of flow field applied. Additionally, a structural analysis of the crystalline clusters indicates that the flow field suppresses the compaction and crystalline ordering of clusters, leading to the formation of large articulated clusters under strong flow fields, and compact well-ordered clusters under weak flow fields.
非平衡分子动力学被用于研究正二十烷在平面剪切以及首次在单轴拉伸下的晶体成核过程。一种基于平均首次通过时间的分析方法被应用于模拟结果,以确定所施加的流场类型和应变率对稳态成核率、特征生长率的影响,以及对与成核相关的动力学参数的影响:自由能垒、临界核尺寸和单体附着前置因子。与剪切相比,流动增强成核(FEN)的起始在拉伸时发生在较小的临界应变率下。对于大于临界应变率的应变率,成核率的快速增加伴随着自由能垒和临界核尺寸的减小,以及链伸展的增加。这些观察结果符合一种机制,即FEN是由流动诱导的熵减少导致的结晶驱动力增加所引起的。在高施加应变率下,自由能垒、临界核尺寸和拉伸程度达到饱和,而单体附着前置因子和取向有序度稳步增加。这种趋势表明在强流动下对成核率有显著的扩散贡献,这与成核系统中的全局取向有序度相关。对于所有动力学量,相对于我们定义为施加应变率与临界应变率之比的约化应变率,两种流场给出了相似的结果。特征生长率随应变率增加而增加,并且显示出与成核率的对应关系,该关系不依赖于所施加的流场类型。此外,对晶体团簇的结构分析表明,流场抑制了团簇的压实和晶体有序化,导致在强流场下形成大的铰接团簇,在弱流场下形成紧凑的有序团簇。