半刚性链高分子体系相行为的自洽场理论研究

摘要

经过50多年的不断发展和完善,自洽平均场理论(Self-Consistent Field Theory,SCFT)已经成为高分子凝聚态物理理论中十分重要且有效的理论方法,由于可以考虑不同的链拓扑结构和给出链构象信息,因此被广泛应用在非均相多组分高分子体系的热力学平衡态结构研究中。然而,关于蠕虫状链模型的自洽场理论,其数值算法远远没有高斯链自洽场理论那样完善。一方面,由于蠕虫状链模型中引入单位矢量u描述链的弯曲刚性,高分子链的构象概率不仅与空间位置r有关,还与链节取向方向u有关。对应的链节概率函数(传播因子)q(r,u,s)的扩散方程,在三维空间中是一个高达六维的偏微分方程;另一方面,蠕虫状链扩散方程中包含对取向角度u(θ,ψ)的二阶梯度算符▽2u,与平面Laplace算符▽2r相比,其数值计算方法还不完善。除此之外,在半刚性链高分子体系中,由于链的刚性效应往往存在各向异性取向相互作用,其描述形式比各向同性相互作用更加复杂。因此,探索运用蠕虫状链自洽场理论及其有效的数值求解算法研究半刚性高分子的相行为,一直是高分子物理领域理论研究者长期以来致力于解决的难题。
　　至今为止,关于蠕虫状链模型的自洽场理论,其数值算法还处于起步和探索阶段。其中较为具有代表性的是基于球谐函数Y1,m(u)展开的思想,对体系中函数的取向依赖性(u-dependence)进行变量分离,在球谐空间处理取向问题和完成角度Laplace算符的求解,再通过“合成”方法获得实空间下关于取向依赖u(θ,ψ)的信息。由于数值计算量较大,球谐函数展开法在研究半刚性链体系时大多对取向进行简化处理,仅考虑对角度θ的依赖性,即令m=0,因此只得到nematic、smectic-A等液晶结构,却无法观察到smectic-C相。虽然借助免费的SPHEREPACK软件包,可以完成在实空间(θ,ψ)和球谐空间(l,m)之间的快速变换,在一定程度上提高了计算效率,同时也可以考虑链节取向对角度ψ的依赖性(m≠0)。然而,由于目前借助于SPHEREPACK软件包的球谐变换只适合于描述在球面上分布连续光滑的函数,因此该法不适用于强取向相互作用的体系。因此,本论文的思路是探索有效的数值算法,求解基于蠕虫状链模型的自洽场理论,并初步探讨其在实际半刚性链高分子体系中的应用。
　　本论文的第一章首先从统计物理的基本原理出发,介绍不同柔顺性高分子链的相关模型(从高斯链到蠕虫状链直至完全刚性链),并以rod-coil嵌段共聚物为代表,介绍近四十多年来关于半刚性链高分子体系的实验和理论模拟研究进展。
　　第二章将介绍基于蠕虫状链模型自洽场理论的实空间算法。其基本思想是:通过三维笛卡尔直角坐标系下,三角网格化的单位球面描述体系的取向变量u,并结合曲面Laplace-Beltrami算子的有限体积算法求解角度Laplace算符。该实空间算法的特点是,不需要对取向变量的描述空间进行变换(球谐分解与合成),表达具有直观性并且操作简单。更重要的是,由于对取向变量u的描述是真实的三维空间表达,不需要对体系的取向对称性进行预先假定,因此能够得到液晶取向与界面法向不在同一方向的有序结构(如smectic-C),这是目前大多研究中采用球谐函数展开法无法得到的。本章模拟了rod-coil两嵌段共聚物的相分离,计算了体系在Onsager排斥体积作用下,由于熵效应可以形成具有部分双层结构的smectic-A液晶相。在完全相同的参数条件下,与球谐函数展开法得到结果完全相同。
　　Rod-coil嵌段共聚物是一类典型的半刚性链高分子体系,近年来在功能材料领域的应用前景就得到了人们越来越广泛的关注。尤其是以共轭高分子(conjugatedpolymer)作为刚性嵌段的一类聚合物,由于共轭结构的半导体和光学特性,其在有机光电器件领域的开发得到了深入研究和长足发展。然而,rod-coil的自组装行为与传统coil-coil柔性两嵌段共聚物有显著不同,究其原因主要包括以下两点:一方面,体系中不仅由于化学组成的差异性产生各向同性相互作用,而且由于链节的刚性效应产生各向异性取向相互作用;另一方面,rod和coil在分子链尺寸的描述上也有显著不同:对于柔性高斯链,通常使用无扰均方回转半径标度,Rg～N1/2;而对于半刚性或刚性链,则通常使用链的围线长度描述,L～N。因此,rod-coil体系的相图具有明显的不对称性,而且层状相稳定区域明显变大。特别的是,在有序层状结构中,rod-coil又可以排列成各种不同smectic微观结构:根据rod液晶取向与层状法线方向之间的夹角θt,可以区分为smectic-A(θt=0)和smectic-C(θt≠0),根据rod的排列方式可以区分为bilayersmectic、monolayersmectic和foldingsmectic。综上所述,rod-coil嵌段共聚物体系具有更为复杂的参数空间,既存在微相分离与取向相互作用的耦合,又存在rod与coil之间的尺寸差异效应。这些因素会造成rod-coil体系中具有特殊的分子链排列规律,从而表现出特殊的液晶相行为和微观聚集形态。因此,从理论角度研究rod-coil嵌段共聚物的相行为和组装规律,对指导实验合成、自组装和材料性能优化,具有重要的指导意义和实际应用价值。鉴于此目的,本章基于蠕虫状链模型的自洽场理论,系统研究了rod-coil嵌段共聚物在稀溶液和本体中的相行为和有序微观结构。
　　第三章使用蠕虫状链模型描述刚性rod和柔性coil嵌段,忽略两种嵌段之间的各向同性Flory-Huggins相互作用,仅通过弯曲刚性模量(ζR=10和ζC=0.1)来区分rod与coil。使用溶致液晶理论中的Onsager排斥体积相互作用描述所有链节(包括rod-rod、rod-coil、coil-coil)之间的取向相互作用,以聚合物浓度G和rod嵌段体积分数fR为变量,计算得到了体系在一维空间下的相图(G-fR)。其中包括典型的isotropic、nematic及smectic相,它们对应的相区分布与球谐函数展开法得到的结果基本一致。特别的是,由于实空间算法对取向描述具有直观和完整性,无需对体系的液晶取向进行假定和简化,因此在一维空间相图中进一步区分了smectic-C与smectic-A结构。
　　第四章对rod-coil理论模型和相互作用描述进一步完善,同时也改进了数值算法,将体系中函数的空间依赖性使用基函数展开(Matsen谱方法),从而提高了对空间算符▽r求解的精确度和稳定性。该章节重点讨论微相分离与取向相互作用比值μ/X,和rod与coil间的尺寸不对称效应β=L/Rg,对rod-coil嵌段共聚物本体相行为和smectic微观结构的影响。结果发现,在相同的μ/χ参数下,β较小有利于双层结构的稳定;β较大则倾向于形成单层smectic,进一步增大β有利于smectic-C的稳定。另外,当β确定而体系微相分离作用明显时,rod-coil倾向于发生微相分离,同时nematic相区显著减小;而当取向相互作用占主导时,体系具有明显的液晶相行为,nematic稳定相区增加。以上结论与已有的实验结果和理论分析基本一致。另外还发现,体系仅在Maier-Saupe取向作用下就能形成isotropic、nematic与smectic液晶结构,这与第三章中使用Onsager排斥体积模型得到的结果基本类似。特别的,基于Maier-Saupe平均场理论的取向相互作用(μN)对rod/coil两相界面面积也有显著影响:随着μN增大,刚性rod间的取向程度增大,排列更加紧密,导致界面面积减小。此时,由于coil的构象熵效应,体系会向界面面积增大的微观结构转变。因此,在相图中观察到了有趣的foldingsmectic-A(fA)和foldingsmectic-C(fC)结构。
　　在有机光电高分子材料的应用研究中发现,激子(电子-空穴对)的有效扩散距离仅为10nm,因此寻取合理途径来控制纳米结构的聚集尺度,并获得优越的两相界面和微区取向,成为优化材料性能的一项关键技术。实验研究发现,向rod-coil嵌段共聚物中共混化学组分相同的rod或coil均聚物,即能获得rod与coil组分的合适比例,是一种高效且简便的调控体系自组装形态的手段。
　　有鉴于此,第五章沿用第四章中的理论模型和混合数值算法,针对rod-coil嵌段共聚物与rod或coil均聚物共混体系的相行为展开研究。结果发现,rod-coil/rod和rod-coil/coil共混物在相同的作用条件(χN和μN)和coil体积分数(fC)下具有更宽的近晶相(smectic)稳定区间。另外,根据自洽场理论得到的嵌段浓度分布和链节分布信息,还重点考察了有序结构中微区尺寸(包括层状周期、rod富集区尺寸、coil富集区尺寸等)和界面形态的变化规律,发现rod和coil均聚物在rod-coil嵌段共聚物本体中的相容机理有显著不同。Rod均聚物会与rod嵌段相互穿插、并排列成有序的液晶取向微区,此时rod/coil的界面面积变化会诱导体系在单层smectic和双层smectic之间发生转变。Coil均聚物则由于分子量不同,与coil嵌段具有不同的相容性。当分子量较大时,coll均聚物倾向于形成独立的浓度富集区,并显著增大coil组分的微区尺寸,但rod/coil界面和rod微区形态基本不发生变化;当分子量较小时,coil均聚物会溶胀在coil嵌段的微区内,并向界面附近渗透,造成rod/coil两相界面变宽。