预取向半晶态高分子片晶结构形成微观机理及其应力-应变响应特性的分子动力学模拟*

杨俊升 朱子亮 曹启龙

1) (宜宾学院, 计算物理四川省高等学校重点实验室, 宜宾 644000)

2) (潍坊科技学院, 寿光 262700)

3) (青岛大学材料科学与工程学院, 青岛 266003)

通过分子动力学模拟方法对不同预取向聚乙烯醇熔体(polyvinyl alcohol, PVA)形成的半晶态高分子熔体形成核结晶及拉伸过程中的应力−应变响应特性进行了系统地研究.模拟结果显示预取向度高的PVA熔体对应更快的成核动力学.通过追踪全trans伸直链长度(dtt)、成核原子维诺体积(V )和中心对称参数(S)等序参量在不同取向度熔体下的等温成核与结晶演化过程, 给出了PVA熔体成核路径及形成半晶态的分子构象; 通过对形成的半晶态高分子结构进一步分析, 发现随着熔体取向度的增加, 晶体和无定型对应的取向度也会增加, 但是当应变剪切大于5时, 其对应的结晶度、晶体和无定型取向度不再发生变化; 通过对无定型区链结构的定义与分析, 可知取向度越高的熔体对应越高的Tie链数目, 随着熔体取向度的增加, Loop链的数目也会减小; 通过恒速拉伸应力测试可知, 所形成半晶态高分子结构力学响应会随着取向度及Tie链数增加而增加, 当取向及Tie链数目相同时, 应力−应变曲线形状大小也基本保持一致.

1 引 言

高分子材料成型加工过程中, 不可避免地要经历不同维度和强度流场的作用[1−3], 而这些流场的作用对于高分子材料最后的力学、光学和电学等物理特性都会产生一定的影响.归根结底, 影响这些物理特性的本质原因是高分子熔体的结构发生了变化, 这就需要在分子层面对其结构变化的机理进行深入地研究.在高分子制品加工初期对高分子熔体进行预取向是高分子加工的一般步骤, 如流延、挤出或者纺丝, 都是将高分子熔体进行预取向而进行初步加工的, 但是关于高分子预取向以后高分子熔体的成核与结晶特性的研究又非常有限, 因为实验在分子层面也很难解释其中的机理, 所以分子动力学模拟不失为一种有效可行的办法来解决相关的问题[4−15].

实验上对于相关的研究已经发展了更高时间分辨率和更小空间的尺寸探测技术, 如宽角X射线散射 (wide angle X−ray scattering, WAXS)、小角X射线散射(small angle X−ray scattering, SAXS)及小角中子散射(small angle neutron scattering,SANS)等探测技术[2,16,17], 其目的就是为了探索高分子材料在加工过程中结构转变的机理, 以及各结构参数之间的相互关系对加工后高分子材料最终物理特性的影响机制.虽然目前实验技术已经有了很大的改进, 可以根据结构的原位X射线衍射花样来判断其对应的结构与性能的关系, 但是现在的实验技术还是很难从分子层面或者原子层面直接地检测到高分子成核与结晶的相关信息.故而, 分子层面动力学模拟研究成核与结晶演化过程是十分必要的.

目前理论模拟方面也有相关研究, 但是基本上都是基于各向同性的熔体进行的, 对于预取向熔体从无序到有序的成核与结晶机理的研究还十分有限.关于聚合物各向同性的聚合物熔体的研究, 最近Luo和Sommer[18]的研究工作认为高分子的成核与结晶是由高分子最开始的缠结状态决定的.关于预取向的聚合物熔体, Tang等[8]最近的模拟工作显示, 预取向会促进高密度区的形成, 从而促进晶体成核速率增加, 这个过程有别于各向同性聚合物熔体的结晶过程.Tang等[8]研究了温度带来的影响, 但是没有考虑不同取向度对无定型区和力学特性的影响, 这也是进一步研究的重点.关于高分子成核与结晶的研究进展, 中国科学技术大学李良彬教授课题组[19]对于相关的工作也做了展望, 他们提出目前高分子独有的柔性和连接性是高分子有别于小分子的区别, 同时也指出了目前高分子结晶面临的问题：1)柔性链段是如何变成刚性链段的？2)片晶间无定型是如何形成的？3)高分子链是如何到达晶体生长前端的？从而可以看出目前关于高分子的成核与结晶确实需要从分子和理论层面深入分析和解答加工过程中面对的这些问题.

基于目前理论和实验方面的迫切需求, 本文将对各向同性高分子熔体的静态结晶及不同预取向的高分子熔体采用分子动力学的方法, 通过不同的序参量来表征高分子熔体在分子层面上对应的成核与结晶过程, 并比较分析给出二者之间的成核机理, 从而给出不同高分子熔体对应的成核结构及晶体生长过程对应的分子图像.通过对形成的半晶态高分子结构进一步分析, 给出不同取向度半晶高分子结构的无定型结构特点和应力−应变响应的分子特性.

2 计算与分析方法

2.1 计算方法及模型

本文主要以 Meyer 和 Müller−Plathe[20]发展的粗粒化聚乙烯醇熔体(coarse grain−polyvinyl alcohol, CG−PVA)模型为基础, 采用分子动力学方法研究CG−PVA的成核与结晶过程, 模拟软件采用平行的LAMMPS上的分子动力学进行.对于所研究的体系, 粒子之间的非键相互作用采用Lennard−Jones−9−6 (LJ96)势来描述[21]：

其中, r是两个粒子之间的距离, ε =1.511kBT ,σ0=0.89σ.当 ULJ(r) 最小时, 截断半径为 1.02σ(σ和ε分别为无量纲的距离和能量常数).两单体相邻之间的键长势由下式给出：

其中, b0= 0.5σ, kb= 2704kBT/σ2.CG−PVA 模型的键角势为列表势, 具体参数是由Luo和Sommer[21]给出.

2.2 模型构建

对于研究体系, 本文选取的是具有缠结效应的N = 500的CG−PVA模型结构, 原因是将聚乙烯醇(polyvinyl alcohol, PVA)熔体从熔点以上降温到熔点以下时, 它具有更规整的结构, 这样的过程更有利于研究聚合物的结晶机理, 而且PVA可以代表类似PE的规整高分子结构的结晶过程.为了避免边界效应的影响, 之后将1000根PVA分子采用无归行走的方式添加到周期性盒子里面, 盒子包含了 50 万个粒子, 其大小为 61.2σ × 61.2σ ×61.2σ.为了保证体系的合理性, 在高温T = 1.0下进行107步的时间弛豫.高温弛豫之后, 将体系沿xy方向施加一个 1 ×10-5τ-1的剪切流场(图1(a)).在动力学弛豫和非平衡动力学过程, 本文选取的时间步长都是0.01 τ .图1(b)给出了一个典型的高分子熔体在剪切场下的应力−应变曲线, 里面包含了应力上升、应力屈服、应力下降和应力平台, 结果与实验观察一致[22], 证明预剪切过程是合理的.同时, 以剪切应变 γ = 0—20对应的结构为等温结晶的研究对象, 探索不同预取向对于高分子等温成核与结晶的影响.图1(b)还给出了剪切过程中, 链的末端距Rete及分子链沿着剪切方向的取向参数Pall的变化.曲线参数主要根据链末端矢量与剪切方向夹角θij采用如下所示的经典取向参数定义：

其中θij角对应的是链与链末端矢量间的夹角大小.不难发现, 链的末端距随着剪切应变的增加达到6个应变以后保持不变, 取向达到10个应变之后不再发生变化, 为进一步研究取向对于结晶的影响创造了条件.在等温结晶过程中选取的结晶温度为T = 0.8, 处于结晶温度以下, 高分子熔体可以自发地成核与结晶在这个温度.

图1 (a)高分子熔体剪切示意图; (b)高分子熔体在剪切场温度为T = 1.0下的σxy, Rete及Pall演化曲线Fig.1.(a) Schematic diagram of polymer melts under shear flow; (b) evolutions of σxy, Rete and Pall for polymer melts under shear with temperature of T = 1.0.

2.3 成核结构分析方法

有效地判断聚合物结晶过程中链内和链间的协同作用, 可以采用Luo和Sommer[21]的伸直链段判断晶体的方法, 如图2(b)所示(定义θ大于150°为trans构象, 如果trans连续单体数目大于13, 则其定义为有效伸直刚性链段), 但是并不是所有的伸直刚性链结构都是晶体结构, 故需要一个判定链间结构的方法来进一步判定.本文在定义晶体与成核原子过程中引入了中心对称参数(centro−symmetry parameter, CSP)与伸直链判定相结合的方法来区分晶体[23], CSP的原理公式为

原理图如图2(a)所示.该方法已经被成熟地应用到了没有连接性的粒子体系中, 可以有效地判断这些粒子体系从无序到局部有序的过程.将该方法引入到聚合物体系的成核与结晶过程, 可有效判断半晶态高分子晶区与无定型区的成核原子S值的两个峰值, 如图2(c)所示.将这两相结合, 不难发现通过空间对称分布给出成核与结晶过程, 不只是包含了链内的构象变化, 还考虑了空间链间构象分布的影响.

图2 (a)成核原子判定方法示意图; (b)伸直链判定方法示意图; (c) PVA半晶态高分子结构晶区与无定型区S值的分布Fig.2.Schematic diagram of (a) centro−symmetry paramet−er and (b) stretched chain segment for PVA crystal; (c) S distribution of crystal and melt of PVA semicrystalline polymers.

3 结果与讨论

3.1 不同预取向PVA熔体成核与结晶过程中的链内链间耦合机理

为了探索预取向对于结晶过程中链内构象、链间构象及密度的影响, 图3(a)—(e)分别给出了不同预取向PVA熔体的结晶度 ϕc、伸直链段长度dtt、成核原子平均维诺体积V和平均对称参数S随时间的演化过程.其中结晶度的计算方式与Luo和Sommer[21]定义晶体的方式是一致的, 将二面角大于150°时的结构定义为trans结构, 然后这样的全trans链段必须大于13个单体才能定义有效伸直链段, 将此判断方法与对称参数相结合, 如果S ＜ 0.8, 将此类原子定义为成核原子.

图3 不同应变下 PVA 熔体的 ϕc , dtt, V 和 S 随时间的演化过程 (a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4; (f) γ = 8Fig.3.Evolutions of ϕc , dtt, V and S of nucleation atoms under different stains：(a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4;(f) γ = 8.

通过图3(a)—(f)中 ϕc的变化不难看出, 预取向的确可以促进高分子熔体的结晶过程, 而且随着剪切应变越大, 对应的成核诱导期也会越短.dtt代表链内构象的变化, 对于各向同性的聚合物熔体(图3(a)), 比较dtt与 ϕc的变化不难发现, 伸直链长度开始明显增长的时间明显要比结晶度开始快速增长的时间要提前, 说明在均相成核与结晶过程中, 伸直链段的形成是成核前提, 而且随着伸直链段的形成, V先于伸直链段和中心对称参数降低,密度降低之后, S也会进一步降低, 伸直链长度会发生一个突然的增加, 进一步说明在稳定核开始形成过程中, 成核链段需要通过成核原子密度的波动形成高密度区, 进而不断地调整自己的构象(S值的降低), 以促进稳定晶核的形成.稳定晶核形成以后伴随的就是晶体稳定的生长过程(图3(a)).对于预取向PVA熔体, 不难发现拉伸会产生更多的高密度区, 而这些高密度区会促进伸直链的进一步形成, 从而促进稳定核的形成.与各向同性的熔体比较, 预取向主要是促进了链间的密度有序, 密度有序会进一步促进晶核的形成.通过链内构象、链间密度及局部序参量, 可以构建出不同取向度高分子熔体的成核与结晶生长机理.这个结论恰好验证了实验上通过SAXS或者WAXS研究结晶的一些猜想.Wang等[24]通过SAXS和WAXS对于高密度聚乙烯的研究表明, 预有序体可以促使晶核的形成, 而且预有序体属于高密度相, 这个结论与本文的模拟数据也是相符的.

3.2 不同预取向半晶态高分子晶区与无定型区的结构特性

图4 (a)半晶态高分子结构中无定型对应的链结构分类模型; (b), (c) PVA半晶态高分子结构中晶体和无定型的原子结构Fig.4.(a) Schematic diagram of crystalline and amorphous chain structure for the semicrystalline polymers; (b) and (c) atomic structure of crystal and melt of PVA semicrystalline polymers.

图5 (a)不同剪切应下PVA半晶态高分子对应的 ϕc 、晶体取向参数Pc和无定型结构取向参数Pa; (b)不同剪切应变下PVA半晶态高分子无定型链结构数目的演化Fig.5.(a) ϕc , crystalline order parameter Pc, and amorphous order parameter Pa for PVA semicrystalline polymers with different shear strains; (b) the evolution of the numbers of amorphous chains for PVA semicrystalline polymers with different shear strains.

高分子在熔体阶段预取向的不同会导致结晶后力学性能和加工性能存在差异, 然而, 这些加工特性和物理特性的区别主要归因于产生的半晶态高分子晶区与无定型区的结构在预剪切过程中对于高分子成核的影响.具体发生了那些改变？实验很难准确地检测到晶体和无定型结构的变化, 动力学结果的统计与分析恰好可以弥补实验上在这方面的缺陷和不足.为了研究无定型结构的区别, 根据实验理解与可能存在的结构, 我们定义了三种类型的无定型结构(图4(a)).Tie 链是连接两个片晶之间的无定型链段, Loop链是折叠链结构, Tail链是链端结构.从图4(c)可以看出, 本文判定方法可以准确有效地判定出晶体和无定型链的结构分布,为进一步的研究提供了可靠的理论依据.

根据图4所示的结构定义, 图5(a)给出了不同剪切应变下半晶态高分子结构晶区与无定型结构晶体取向和无定型取向的变化.从图5(a)中不难发现, 应变越大的熔体形成的半晶态高分子结构对应的结晶度越大, 同时晶体和无定型的取向也越大.分析不同的无定型链结构在半晶态高分子结构中的数量发现, 应变越大的熔体的Tie链数目也越多, 而Loop链的数目会随着应变的增加而减小.在小应变下, Tie链数目越小, Loop链数目越多表明, 没有预取向的半晶态高分子晶区主要是以链内的折叠为主; 在高应变下, 主要是以链间的有序为主, 通过Tie链连接片晶与片晶.Tie链数目的改变及Loop链的数目进而会影响到力学性能的改变.

3.3 不同预取向半晶态高分子结构在恒应变场下的应力-应变响应曲线

为了进一步探索不同取向度半晶态高分子结构的力学特性, 对形成的半晶态高分子结构施加了一个拉伸速率为 1 × 10-5τ-1的恒应变拉伸场.相同的拉伸速率下, 不难发现不同取向度半晶态高分子结构对应的力学性能会有很大的差别.图6为不同初始剪切应变下体系的应力−应变曲线汇总.拉伸的应力−应变曲线与最近实验的结果一致[25−28],进一步从侧面反映了模拟结果的合理性, 但是实验目前只能从拉伸过程中X射线衍射花样来判断晶体的熔融重结晶过程.相比于实验动力学的轨迹,分子动力学可以直接给出结构的演化, 从而进一步验证实验的猜想和推断.分析应力−应变曲线, 发现随着初始剪切应变的增加, 体系的屈服应变逐渐增加, 且初始模量也有所增长, 但是在初始应变达到一定数值后, 屈服应变和初始模量的变化都较小.各体系的初始模量取决于各晶体的取向程度, 取向度越大, 体系初始模量越大; 通过进一步的构象分析可知, 在早期拉伸过程中, 晶体的旋转和破坏是重要的应力来源, 随着晶体与拉伸方向取向一致,晶体破坏减少, 此时主要发生取向链段结晶和晶体间的无定形链拉伸, 因此在大应变初始预剪切体系中, 晶体取向较快, 力学曲线进而与无定形链末端距变化近乎一致; 从微观的原子受力角度来看, 晶区原子沿取向方向的受力大于非晶区原子, 整体上随着链段取向的增加, 原子沿取向方向受力也随之增加; 预剪切对缠结密度影响不大, 但是对于拉伸过程的缠结度变化率有较大影响, 整体上缠结度随着拉伸过程下降, 且初始预剪切应变越大体系最终缠结度越小.

图6 不同预取向半晶态高分子结构在恒定速率为1×10-5τ-1拉伸场下对应的应力−应变曲线Fig.6.Stress−strain curves of PVA semicrystalline polymer with different orientation degree under the stretched rate of 1×10-5τ-1.

进一步分析单链构象的变化.图7给出了应变等于6时形成的半静态高分子单链在不用拉伸应变下的结构.力学曲线结合构象的演化可以得到以下结论.1)预取向可以提高沿取向方向的力学强度, 大应变预剪切和小应变预剪切的力学曲线形状存在差异.2) 拉伸前期, 应力来源于晶体的旋转和破坏, 之后是链段的拉伸和取向结晶.3)链段沿拉伸方向力与取向相关性很强, 取向越大的链, 原子受力越大.因为晶区有着更大的取向, 在拉伸方向力更大.4)对整个拉伸过程的理解：首先是键角尺度上的变化引起的弹性阶段, 然后是晶体旋转和破坏同时伴随无定形链拉伸取向结晶, 之后是晶体基本取向并且在拉伸下继续破坏, 随着链的取向进一步加强取向结晶数目增多, 反映为总结晶度开始上升.

图7 当γ = 6时对应的单链构象随着应变ε的演化过程Fig.7.Structural evolution of single PVA chain as a func−tion of strain ε when γ = 6.

4 总 结

通过分子动力学模拟方法探究了各向同性与不同预取向PVA熔体的成核路径.通过比较 ϕc,dtt, V和S随着时间的演化过程, 发现取向度越高的PVA熔体对应越快的成核动力学.并且得到以下更详细的成核动力学路径和力学响应特性：1)对于各向同性的PVA熔体, 链段结构由柔性向刚性转变先于链间有序的发生, 进而是局部有序的发生, 最后才是成核或密度的转变发生, 这个分子图像被定义成链内有序到局部有序结构波动诱导的成核.2)当PVA熔体在剪切场下取向后, 发现在成核结晶前发生了由于构象有序链段的聚集、共同取向而形成的高密度区, 及密度波动先于成核结晶的发生; 这样的图像被总结为流动下密度波动诱导的成核结晶.3)预取向可以改变半晶态高分子结构的晶体和无定型的取向.对于无定型区, 剪切应变不同, 对应的Tie链、Loop链的数目也会发生明显的变化, 进而影响其结晶加工过程的力学特性.4)预取向可以提高沿取向方向的力学强度, 大应变预剪切和小应变预剪切的力学曲线形状存在差异.拉伸前期, 应力来源有晶体的旋转和破坏,之后是链段的拉伸和取向结晶.链段沿拉伸方向力与取向相关性很强, 取向越大的链, 原子受力越大.