《武汉工程大学学报》 2014年05期
5-8
出版日期：2014-05-31
ISSN:1674-2869
CN:42-1779/TQ
环碳酸酯型聚氨酯的玻璃化温度分子模拟
0引言聚氨酯由于其良好的机械性能、加工性能、耐腐蚀性能、易于改性,广泛应用于黏合剂、涂料、塑料、橡胶等工业领域.玻璃化温度是衡量高分子材料性能的重要指标之一,通过预测玻璃化温度对于新材料的合成具有应用价值.随着计算机技术的发展,分子模拟的精确程度也在逐渐增加,当前的分子模拟已经可以部分替代实验,比如药物的设计,详细的分子模拟可以提供许多通常手段难以观察到的信息,从长远的角度,分子模拟终将成为化学实验的必备工具之一.分子动力学模拟方法已广泛应用于聚合物体系,它被用来预测玻璃化温度［14］.本研究提出一个简单可行方法以预测线型非异氰酸酯聚氨酯的玻璃化温度.J M Barton研究了一个线型环氧模型［5］,将所获得的信息作为建模的思路［6］.本研究首先构建模拟体系的分子模型；然后对模型选择系综动力学进行模拟.最后比较模拟与实验结果来检验模拟方法的准确性以及观察由于结构的不同对非异氰酸酯聚氨酯的玻璃化温度的影响.1线型非异氰酸酯聚氨酯的分子结构建模基于高分子材料模拟软件Materials Studio 4.2,本实验的模拟过程不考虑反应物的交联效应.所选的组分为双酚A型环碳酸酯(BisA5CC)(图1)、丁二醇二缩水环碳酸酯(C45CC)(图2)与1,6己二胺(HDA）(图3),乙二烯三胺(DETA)(图4).此处假设组分之间充分反应形成了无规共聚物,本研究的方法适合制备高转化率的非异氰酸酯聚氨酯［7］.图1双酚A环碳酸酯Fig.12,2bis［4(1,3dioxolane2one4ylmethoxy)phenyl］propane图2丁二醇二缩水甘油环碳酸酯Fig.21,4bis（1,3dioxolane2one4ylmethoxy）butaneNH2CH26NH2 图31,6己二胺Fig.31.6hexanediamineNH2CH2CH2NHCH2CH2NH2 图4二乙烯三胺Fig.4Diethylenetriamine1.1分子模型通过MS软件建立环碳酸酯预聚物(图1)与胺组分(图3)的小分子结构,对于图1分子进行开环处理,并手动添加未饱和的氧原子上的氢原子使其变成羟基.如图5所示.图5分子构建方式Fig.5Molecular way to build通过MS软件定义无规共聚物单元的头和尾,并以头尾的连接方式形成一条聚合度为10的分子链.然后通过forcite模块进行能量小化处理,结合Qeq电荷平衡对分子链分配电荷.优化方法使用smart优化,因为此时的分子链所含的原子348个,可以直接优化.然后以Monte Carlo方法产生无定型聚合物的代表构象,这种方法基于旋转异构态理论(RIS)，认为稳定构象的贡献大大高于其他构象,这样大大减少了计算量.以单链构象为基础,结合统计力学原理进行相关推导,得出配分函数与宏观热力学量之间的关系.利用MS软件的Amorphous cell模块,产生一个无定型元胞,其中含有10个代表构象.对产生的元胞经行能量最小化和初步的动力学松弛,系综选择NPT,进行小步长的松弛(1 000步),此时的系统原子数已达到3 480个,需要进行对结构进行一定的松弛.这个模型被期望反映开环聚合后的体系.建模结果如图6所示.图6模型示意图Fig.6Schematic model第5期樊庆春，等：环碳酸酯型聚氨酯的玻璃化温度分子模拟武汉工程大学学报第36卷1.2玻璃化温度的分子动力学模拟初始模型建立以后,对系统能量优化直至收敛到稳定值.然后进行NVT系综模拟,以1 fs(1 fs=10-15 s)为时间步长,温度600 K,模拟3 000个时间步,不使用截断距离,然后再经过NPT系综模拟,同样是1 fs的时间步长,执行5 000个时间步.重复以上的模拟步骤,温度设置为550 K,交叉进行NVT,NPT模拟,在450 K时,温差缩短为每25 K进行一次实验模拟,以此类推直到温度降为200 K.图1分子与图4分子,图2分子与图3分子，图2分子与图4分子进行同样的模拟处理.记录模拟系统的体积,通过多个温度点,得出温度体积曲线,由曲线的拟合处理得到玻璃化温度.整个动力学模拟过程使用DREIDING力场,体系总势能分为成键作用和非键作用.成键作用包括键的伸缩能、弯曲能、扭转能、振动能等等；非键作用包括van der Waals相互作用和Coulomb静电作用.其他参数设置为默认值.2结果与讨论对模拟体系的体积与温度作图观察拐点的出现.因为玻璃化温度的附近会出现广度性质和强度性质的变化,通过对拐点两边的点进行分段线性拟合.结果如图7~10所示.图7HDA+C45CC模拟体系的VT图Fig.7HDA+C45CC simulation system VT diagram 图8HDA+BisA5CC模拟体系的VT图Fig.8HDA+BisA5CC simulation system VT diagram图9DETA+C45CC模拟体系VT图Fig.9DETA+C45CC simulation system VT diagram图10DETA+BisA45CC模拟体系的VT图Fig.10DETA+BisA45CC simulation system VT diagram表1不同体系玻璃化温度的实验值和模拟值对比Table 1Experimental values of different systems andglass transition temperature constant analog values模拟体系模拟值/K实验值/KHDA+C45CC251247HDA+BisA5CC340320DETA+C45CC255246DETA+BisA5CC370322结果与实验值［7］比较,对于含有线型柔性链结构的非异氰酸酯聚氨酯,模拟结果能较好预测玻璃化温度如图6和图8的结果,对于含有刚性苯环结构的生成的非异氰酸酯聚氨酯，模拟值对实验值存在较大偏差如图7与图9.产生偏差的原因可能是本实验的MD分子动力学模拟考虑的是典型分子构象的运动,实验测得的玻璃化温度是通过DSC、DMA仪器测得,它们依据的是Tg附近产生的性质的突变效应,反映到热效应上,通过传感器测出变化.是一种宏观热效应的表现,而模拟是基于分子层次的表现，包含了更原始的信息.分子模拟考察了分子内和分子间的详细的相互作用.对于柔性分子的预测具有一定的指导作用.具有体型结构的非异氰酸酯聚氨酯的玻璃化温度的预测还需要提出更接近实际情况的分子模型\[89\].3结语采用Monte Carlo法选取代表性构象,通过分子动力学模拟预测了两种线型模型的玻璃化温度.柔性结构的聚氨酯模拟值与实验值接近,可用于预测柔性结构的非异氰酸酯聚氨酯的玻璃化温度;刚性结构聚氨酯的模拟值与实验值偏差较大,说明对于刚性结构的模型还需要进一步研究.致谢感谢武汉工程大学化工与制药学院分子模拟研究室的帮助.