2017年第9期 物理通报 教育技术应用
教育技术应用
半柔性聚合物中刚性系数K b。nd对
体系影响的Mon te Carlo模拟研究
肖立勇 纪登辉
(六盘水师范学院电气工程学院 贵州六盘水 553004)
(收稿日期:2017—03—24)
摘要:结合自洽场理论的粒子数表象到场表象思想,运用蒙特卡洛模拟方法对半柔性聚合物体系进行计算
模拟.对于一个半柔性聚合物体系,其体系的哈密顿量主要包括键与键之间的键能、体积排斥能和弯曲弹性能,其
中的键能用弹性能来表示,体积排斥能可通过粒子表象到场表象转化的密度场来描述,单链聚合物的弯曲弹性能
可以刚性系数Khena来表述.通过计算模拟发现:1.在稀溶液中的单链聚合物体系中,其回旋半径(R )和末端距
(Ro)随着弹性系数的增加而增长;2.在高浓度的聚合物体系中,其回旋半径和末端距不随弹性系数的变化而变化,
其投影长度(L )随着弹性系数的增加而增加;3.投影长度与单链聚合物的聚合度N的a倍成正比,与弹性系数
K a的b次方成正比,且。≈b≈0.4士0.45;4.投影长度的与体系中聚合物的浓度无关.
关键词:半柔性聚合物 场表象 投影长度 蒙特卡洛模拟
1 引言
半柔性聚合物在自然和生活中领域扮演着重要
角色,半柔性聚合物的性质在一定程度上和一般的
柔性聚合物相似,如能够自组装成各种相态等[1 ],
同时又有自身的性质,如具有一定的刚性,如DNA
分子等.对于半柔性聚合物,在天然和人工合成中,
如多肽、烷基取代纤维素等,具有温度响应性,同时
根据半柔性聚合物的耐热性、加工性优异性,在电
子、汽车部件等领域得到广泛的应用[7 引.在理论
研究方面,C.M.Chen等[9 用Monte Carlo模拟半
柔性聚合物刚性系数对链折叠动力学产生重要的影
响,其发现随着链的刚性系数的增强,其相态出现许
多缺陷,同时Jianhui Song等口。。用格子Monte
Carlo模拟一个半柔性ABA三嵌段共聚物的相行
为,通过研究发现不同刚性的聚合物对聚合物的相
图有很大影响.同时还有Wendi Song等[1q 用一
种新的现实空间数值实现自洽场理论对聚合物的研
究和用自洽场对半柔性线圈的二嵌段共聚物的相行
为熔化进行了研究;Naveen C.Andrewsl等[6 用构
型偏倚蒙特卡洛(CBMC)和非平衡布朗动力学
(NEBD)模拟来研究分析半柔性聚合物的分子动
力学.本文主要用蒙特卡洛(MC)对半柔性聚合物
进行模拟研究,对于本体系中半柔性聚合物的刚性,
我们引入一个刚性系数Kk a,通过刚性系数的引
入,结合自洽场理论对半柔性聚合物的分析和推导,
利用蒙特卡洛模拟方法对体系进行计算模拟,其中
研究和讨论了刚性系数对半柔性聚合的回旋半径
(R )和末端距(R。)的影响.
2 模拟和理论
对于将要模拟的半柔性AB聚合物体系(如图1
所示),其体系中半柔性AB聚合物链为 条,体积为
,且每条半柔性AB聚合物通过粗粒化处理,其聚
合度为N,则半柔性AB聚合物中各组分的单位体
一 r 积片段数密度为』D。一 .设定半柔性AB聚合物键
y
长为b。,其中每条半柔性AB聚合物链中A组分(未
涂灰)占整条链的比例为,^,B组分(涂灰)占整条
链的比例为^,且满足,A+,B一1,A组分与B组分
*国家自然科学基金,项目编号:11504078;贵州省联合基金重点项目,项目编号:黔科合LH字[2014]7449;六盘水师范学院高层次人
才科研启动基金LPSSYKYJJ201404;贵州省教育厅重点项目,项目编号:黔教合KY字[-20151379号
作者简介:肖立勇(1986一 ),男,硕士,讲师,主要从事软凝聚态物理中的纳米粒子粒子复合材料的研究.
通讯作者:纪登辉(1985一 ),男,博士,主要从事磁性纳米功能材料研究.
一89—
2017年第9期 物理通报 教育技术应用
的相互作用势为 N,其中x为Flory—Huggins参
数,半柔性AB聚合物链存在一定的刚性,则其刚性
系数K ,半柔性AB聚合物的回旋半径为R 和末
端距为R。.这样就可以根据自洽场理论把半柔性
AB聚合物体系相应的Hamiltonian量[28]写出来,
其中包括:半柔性AB嵌段聚合物中片段与片段之
间的键能H 和非键能H 以及半柔性AB聚合物
的刚性能H ,即
H[{r, ,0}]一H ]+H [,]+H E03 (1)
图1 半柔性高分子模型
其中第i条半柔性AB嵌段聚合物中片段与片
段之间的键能H 和非键能H 为
H 3k
—
BT ̄ 巡 (2)
H [{,)]一是n丁P。l ,dr[∑ (r) p(r)+
警[1一 (r)一即(r)] ] (3
其中 (,.)和 p(r)表示空间r出A组分和B
组分所对应的密度, 表示体积不可压缩势.
对于半柔性AB聚合物,其聚合物有一定的刚
性,为了能引入其存在的刚性能,我们一个刚性系数
K ,如图1所示,第S片段在第S一2片段到S—l
片段方向上有一定的回复力,对这种回复的强度就
用刚性系数K a来表示,则其单链中所包含的弹性
势能为H
Hrb [s]=kBT Kbendsin2(舞) (4)
其中0 表示第i条链中第S一1个片段到第S个
片段的方向矢量与第s一2个片段到第 一1个片段
的方向矢量的夹角.
则对于半柔性AB聚合物链,其体系的
Hamiltonian量可以表示为
H[{r,s,0)]:==
3k
2
.T  ̄
s=l
+
忌 。l dr[∑ (r) p(r)+
2
(1一 (,)一 (r))。]+
90一
是eT K dsi n2( ) (5)
通过对体系中各对应的参量进行定义和确定,
基本上已经完成了体系Hamiltonian量,但在进行
蒙特卡洛计算模拟中,考虑到片段之的键长存在一
定的伸缩性,则在计算模拟中把键长取在0.8~1.2
之间模拟.
3 模拟结果和讨论
3.1 刚性系数K 对回旋半径(R )和末端距
(R。)的影响
对于半柔性AB聚合物,A组分与B组分存在
排斥势 ABN一30,不可压缩势 一50,其中A组分
占总片段数的比例为0.5,分别对不同聚合物数在
不同刚性系数下得回旋半径和末端距进行研究.
(1)在稀溶液体系中,其中忽略溶剂对半柔性
AB聚合物体系的影响,通过研究模拟发现,如图2
所示,随着刚性系数的增强,其回旋半径(R )和末
端距(R。)也在不断的增长,则说明在稀溶液中刚性
系数对半柔性聚合物的回旋半径和末端距有很大的
影响,且显线性增长.
b
图2 回旋半径R 和末端距Ro随刚性系数K b d的变化情况
2017年第9期 物理通报 教育技术应用
(2)在高浓度体系中,通过对半柔性聚合物回
旋半径和末端距的研究发现,如图3所示,通过研究
发现,对于同一刚性系数的半柔性聚合物链,随着聚
合度的增加其回旋半径和末端距在增加,但对于同
一聚合度的聚合物其回旋半径和末端距不受刚性系
数的影响,则说明在高浓度体系中,半柔性聚合物片
段之间的相互作用势要大于分子本身的刚性作用
势.
图3 不同聚合度F回旋半Rg和末端距R0
随刚性系数Kb 的变化情况
3.2 刚性系数K a对余辉长度(Ls)的影响
对于半柔性聚合物分子,通过上述的研究和模
拟不能发现其末端距和回旋半径的差异,则为了更
好地研究和反应刚性聚合物分子的特点,我们引入
投影长度(L )来研究刚性聚合物分子的特性,则对
于一条聚合物分子的平均余辉长度可以表示为
" N L。一∑∑(, (s)一r (s一1))・
=1 s 1
COS( )/ (6)
其中 为第 条链上的第 一2个片段到第 一
1个片段的矢量与第 一1个片段到第 个片段矢量
的夹角如图1所示.
根据上式模拟发现,对于同一条聚合物,当刚性
系数K 不断加强时,其余辉长度也在不断的增
加,且余辉长度与单条聚合物的聚合度(N)成正比,
如图4所示.
即
L 。C口*N (7)
图4 不I司刚性系数K b。 d F,余辉长度
(Ls)随聚合度(N)的影响
同时为了研究余辉长度与刚性系数K a的关
系,我们研究了K 对余辉长度的影响,通过研究
发现,对聚合物体系刚性系数K 对余辉长度的影
响也存在一定的影响,如图5所示.
图5 不同聚合度(N)下,余辉长度(Ls)
随刚性系数Kb d的变化
通过相应的数据处理和模拟计算,发现余辉长
度与聚合度(N)的关系系数口一0.4±0.045,余辉
长度与刚性系数(K a)的关系系数
b一0.4±0.045
即可以表示为
Ls。C a*N*K (8)
其中a一0.4±0.045,b一0.4±0.045.
3.3 体系浓度对余辉长度(Ls)的影响
为了能充分地了解体系浓度(p)对半柔性聚合
物的相关性质,我们还研究了体系浓度对与余辉长
度的影响,通过研究发现,余辉长度不受体系浓度的
影响如图6所示.
一91—
2017年第9期 物理通报 教育技术应用
图6 不同体系浓度(p)下,余辉长度(Ls)
随刚性系数K d的变化
4 总结与展望
随着高分子研究的不断深入,半柔性聚合物的
研究也在不断的提升,对于半柔性聚合物,其结构比
一般的柔性高分子要复杂许多,则在研究时要考虑
的参数有很多,这样一来对计算模拟带来很大的困
难.本文中主要以简单半柔性聚合物为研究对象,对
其末端距(R。)、回旋半径(R )、余辉长度(Ls)受刚
性系数的影响.通过研究发现刚性系数对余辉长度
的影响很大,则余辉长度L 可以用来反应半柔性聚
合物分子的一个特征.
随着半柔性研究的深入,我们可以逐步来研究
纳米粒子复合材料 、DNA分子、多肽等半柔性
高分子的研究,同时可以结合实验来研究半柔性分
子的一些特性,这对高分子的研究有一定的促进作
用.
参考文献
1 Wendi Song,Ping Tang,Hongdong Zhang,et a1.New
Numerical Impiementation of Self——Consistent Field
Theory for Semiflexible Polymers.Macromolecules,
2009(42):6 300~6 309
2 Wendi Song,Ping Tang,Feng Qiu,et a1.Phase behavior
of semiflexible—coil diblock copolymers:a hybrid
numerical SCFT approach.Soft Matter,2011(7):929~
938
3 Jie Gao,Wendi Song,Ping Tang。et a1.Self——assembly
of semiflexible block copolymers:2D numerical
implementation of self—consistent field theory.Soft
Matter,2011(7):5 208~5 216
4 Jie Gao,Pi ng Tang,Yuliang Yang.Non—lamellae
structures of coil—semifl exible diblock copolymers.
92一
Soft Matter,2013(9):69~81
5 Ying Jiang。Jeff Z.Y.Chen.Influence of Chain Rigidity
on the Phase Behavior of Wormlike Dib1ock
Copolymers.PRL,2013(110):138 305
6 Naveen C.Andrewsl,Anthony J.MeHugh,Jay D.
Schieber.Configuration biased Monte Carlo and
Brownian dynamics simulations of semiflexible
polymers in extensional flows.Macromo1.Theory
Simu1,1998(7):19~26
7 肖立勇,曾云兰,张博凯.掺杂无机性纳米粒子的嵌段聚
合物体系的场理论模拟研究.科学技术与工程,2014,
14(11):1 671~1 815
8 Kostas Ch Daoulas Victor Rtihle,Kurt Kremer.
Simulations of nematic homopolymer melts using
particle based models with interactions expressed
through collective variables.J.Phys.:Condens.
Matter,2O12(24):284 121
9 C M Chen,Paul G Higgs.Monte—Carlo simulations of
polymer crystallization in dilute solution.Journal of
ChemicalPhysics,1998,108(10):4 305~4 314
10 Jianhui Song,Tongfei Shi,Yunqi Li,et a1.Rigidity
effect on phase behavior of sym metric ABAtriblock
copolymers:A Monte Carlo simulation.The Journal Of
Chemical Physics,2008(129):054 906
1 1 Jan Kierfeld,Krzysztof Baczynski,Petra Gutjahr,et a1.
Modelling semiflexible polymers:shape analysis,
b uckling instab ilities, and force generation.Soft
Matter,2010(6):5 764~5 769
12 N Arun Kumar,Venkat Ganesan.Com munication:
Self——assem bly of semiflexible——flexible block
copolymers.The Journal of Chemical Physics,
2012(136):101 101
13 Jan Kierfeld,Krzysztof Baczynski,Petra Gutjahr,et a1.
Semiflexible Polymers and Filaments: From
Variational Prob1ams to Fluctuations.American
Institute of Physics,2008:151~182
14 J N Bright,D R M Williams.Collapse of semi—exible
brushes in poor solvents:Tilting and bending.
Europhys.Lett.1999,48(5):540~546
15 M.Baratlo,H.Fazli.Molecular dynamics simulation
of semiflexible polyampholyte brushes——The effect of
charged monomers sequence.Eur.Phys.J.E,2009(29):
131~138
16 Philip C Nelson.Generalized theory of semiflexible
polymers.Physical Review E,2006(73):031 906
2017年第9期 物理通报 教育技术应用
l7 胡显奎,林少全,刘振兴.聚合物基无机纳米粒子复合
材料的制备技术及应用展望.材料导报,2000,14(10):
62~63
18 王春艳.聚合物/无机纳米复合材料的发展及应用.化
学工程师,2001(6):21~22
19 Matsen M.W,Schick M.Mierophases of a Diblock
Copolymer with Conformational Asym metry.
Marocmolecules,1994(27):4 014~4 015
Mon te Carlo Si m ulation Study on the Effect of Stiffn ess
Coefficient Kbend on the System in Semi Flexible Polymers
Xiao Liyong Ji Denghui
(School of Electrical engineering of Liupanshui Normal University,Liupanshui,Guizhou 553004)
Abstract:The semi—flexible polymer system has been investigated used by the Monte Carlo simulation
method based on the Ideological and theoretical of Particle representation to Density field representation in self—
consistent field theory.The Hamihonian of the system include the bond energy,the volume repulsion energy and
the bending elastic energy,which can be expressed by the elastic energy,the density of polymer segments,the
elasticity Kb d,respectively.
The results of computer simulation indicates that:1.For the dilute solution of single—chain increase of
elasticity Kb d.2。For the high concentration of polymer systems,the Cyclotron polymer system,the Cyclotron
radius(R )and end—to—end distance(R0)increases with the radius and the end—to—end distance does not vary
with the change of elasticity Kb d,the projected length(L。)increases with coefficient of elasticity increasing.3.
The projected length is proportional to aN,and K 。 d,a≈b≈0.4士0.45.4.The projected length has nothing to
do with the concentration of the polymer.
Key words:semi flexible polymer;field representati0n;projection length}monte carlo simulation
(上接第88页)
Measure m en t Device of Refractive In dex
Base on Double Dark Ring Effect
Zhan Gaochao Liu Wei Huang Yan mei Li Sanfeng Liu Zhaohui
(South China Normal University,School of Physic and Telecommunications Engineering,Guangzhou,Guangdong 510006)
Abstract:With the double dark ring effect due to the total reflection,this device analyze the mathematical
relationship between the specific value of the two circular dark rings diameter and the refractive index of the
substance.This is a new method for refractive index measuring.The main measurement devices are laser,glass
container,and smart phone.Proved by experiment,this method is not only easy and simple to handle,but also
with high precision and strong stability.With this method,we designed a household device and developed a
supporting application.Applying to eellphone even make the measure more convenient and fast.People can
measure substances with this device,like the purity of cooking on or the solid content in fruit juice,etc.
Key words:dark ring effect;total reflection;refractive index;non—contact measurement
一93—
本文发布于:2022-12-31 18:30:01,感谢您对本站的认可!
本文链接:http://www.wtabcd.cn/fanwen/fan/90/67094.html
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
留言与评论(共有 0 条评论) |