Johnson Lilian C, Phelan Frederick R
Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States.
J Phys Chem B. 2023 Aug 10;127(31):7054-7069. doi: 10.1021/acs.jpcb.3c03273. Epub 2023 Jul 31.
In this work, we extend an approach to coarse-grained (CG) modeling for polymer melts in which the conservative potential is parametrized using the iterative Boltzmann inversion (IBI) method and the accelerated dynamics inherent to IBI are corrected using the dissipative Langevin thermostat with a single tunable friction parameter ( 084114). Diffusive measures from picoseconds to nanoseconds are used to determine the Langevin friction factor to apply to the CG model to recover all-atom (AA) dynamics; the resulting friction factors are then compared for consistency. Here, we additionally parametrize the CG dynamics using a material property, the zero-shear viscosity, which we measure using the Green-Kubo (GK) method. Two materials are studied, squalane as a function of temperature and the same polystyrene oligomers previously studied as a function of chain length. For squalane, the friction derived from the long-time diffusive measures and the viscosity all strongly increase with decreasing temperature, showing an Arrhenius-like dependence, and remain consistent with each other over the entire temperature range. In contrast, the friction required for the picosecond diffusive measurement, the Debye-Waller factor, is somewhat lower than the friction from long-time measures and relatively insensitive to temperature. A time-dependent friction would be required to exactly reproduce the AA measurements during the caging transition connecting these two extremes over the entire timespan at this level of coarse-graining. For the polystyrene oligomers for which we previously characterized the diffusive friction, the viscosity-parametrized frictions are consistent with the diffusive measures for the smallest chain length. However, for the longer chains, we find different trends based on measurement method with friction derived from rotational diffusion remaining nearly constant, friction derived from translational diffusion showing a modestly increasing trend, and viscosity-derived friction showing a modest decreasing trend. This seems to indicate that there is some sensitivity of the friction measurement method for systems with increased relaxation times and that in particular, the unsteady dynamics of the individual parametrization schemes plays a role in this. Increased difficulty in applying the GK method with increasing relaxation time of the longer chain systems is also discussed. Overall, we find that when the material is in a high-temperature melt state and the viscosity measurement is reliable, the friction parametrization from the diffusive friction measures is consistent and the lower cost diffusive parametrization is a reliable means for modeling viscosity. Our data give insight into the time-dependent friction one might compute using a non-Markovian approach to enable the recovery of AA dynamics over a wider range of time scales than can be computed using a single friction.
在这项工作中,我们扩展了一种用于聚合物熔体粗粒化(CG)建模的方法,其中保守势通过迭代玻尔兹曼反演(IBI)方法进行参数化,并且使用具有单个可调摩擦参数(084114)的耗散朗之万恒温器来校正IBI固有的加速动力学。使用从皮秒到纳秒的扩散测量来确定应用于CG模型以恢复全原子(AA)动力学的朗之万摩擦因子;然后比较所得的摩擦因子以确保一致性。在这里,我们还使用一种材料属性——零剪切粘度来对CG动力学进行参数化,我们使用格林 - 库博(GK)方法来测量该粘度。我们研究了两种材料,一种是作为温度函数的角鲨烷,以及之前研究过的作为链长函数的相同聚苯乙烯低聚物。对于角鲨烷,从长时间扩散测量得出的摩擦和粘度都随着温度降低而强烈增加,呈现出类似阿累尼乌斯的依赖性,并且在整个温度范围内相互一致。相比之下,皮秒扩散测量所需的摩擦,即德拜 - 瓦勒因子,略低于长时间测量得出的摩擦,并且对温度相对不敏感。在这种粗粒化水平下,需要一个随时间变化的摩擦来在连接这两个极端的笼形转变期间的整个时间跨度内精确再现AA测量结果。对于我们之前表征了扩散摩擦的聚苯乙烯低聚物,粘度参数化的摩擦与最小链长的扩散测量结果一致。然而,对于较长的链,我们根据测量方法发现了不同的趋势,从旋转扩散得出的摩擦几乎保持不变,从平移扩散得出的摩擦呈现适度增加的趋势,而从粘度得出的摩擦呈现适度降低的趋势。这似乎表明对于具有增加的弛豫时间的系统,摩擦测量方法存在一定的敏感性,特别是各个参数化方案的非稳态动力学在其中起了作用。还讨论了随着较长链系统弛豫时间增加应用GK方法的难度增加。总体而言,我们发现当材料处于高温熔体状态且粘度测量可靠时,从扩散摩擦测量得出的摩擦参数化是一致的,并且成本较低的扩散参数化是一种可靠的粘度建模方法。我们的数据深入了解了使用非马尔可夫方法可能计算出的随时间变化的摩擦,以便在比使用单一摩擦计算更广泛的时间尺度上恢复AA动力学。