搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高分子链在分子刷表面吸附的Monte Carlo模拟

王超 周艳丽 吴凡 陈英才

引用本文:
Citation:

高分子链在分子刷表面吸附的Monte Carlo模拟

王超, 周艳丽, 吴凡, 陈英才

Monte Carlo simulation on the adsorption of polymer chains on polymer brushes

Wang Chao, Zhou Yan-Li, Wu Fan, Chen Ying-Cai
PDF
HTML
导出引用
  • 用Monte Carlo方法模拟研究了一条自由高分子链在分子刷表面吸附的静态和动态特性. 结果表明, 随着自由链与分子刷之间吸附作用能(ε)的增大, 自由链出现由脱吸附态到吸附态的相转变, 同时链的扩散由正常模式转为亚扩散模式. 临界吸附能(εC)几乎与自由链长度无关, 但随着分子刷链长度的减小或分子刷链间距的增大而不断增大. 在εC附近, 自由链嵌入分子刷内部, 同时链尺寸达到极小, 而当ε $\gg $ εC时, 自由链处于强吸附态, 链节主要分布于分子刷表层, 同时整个吸附动态过程可分为自由链吸附和分子刷扩散两个阶段.
    The adsorption of polymer on surface is a hot topic in physical, chemical and biological communities, which is influenced by many factors, such as the topological structure and the flexibility of the polymer, the attractive interaction between the polymer and the surface, the detailed structure of the surface, etc. The adsorption of polymers on solid surfaces is extensively studied, while the adsorption behaviors of polymers on soft surfaces are still unclear. In this work, the static and dynamical characters of the adsorption of a free polymer chain on polymer brushes are studied by using Monte Carlo simulation. The brush is formed by grafted polymers with length Nb and distance d. Results indicate that, with increasing the adsorption energy (ε) between the free polymer and the brush, the free polymer shows a phase transition from a desorbed state to an adsorbed state. Based on the dependence of the number of the adsorption segment of the free polymer (mad) on the adsorption energy ε, we defined the critical adsorption point (εC) where the phase transition occurs. εC is nearly independent of the length of the free polymer, but it increases with decreasing the length of the grafted polymer or increasing the distance between the grafted polymers. When ε < εC, the free polymer is desorbed and its size is the same as that in free space. When εεC, the free polymer is sucked into the brush and meanwhile the size is compressed. While when ε $\gg $ εC, the free polymer is strongly adsorbed on the surface of the brush and forms a quasi two-dimensional conformation, and meanwhile the whole adsorption process contains two stages: the adsorption process of the free polymer and the diffusion process of the brush. Moreover, with the increase of ε, the diffusion of the free polymer shows an obvious transition from the normal model to the sub-diffusion model near εC. The transition of the diffusion model maybe useful for separation of polymers with different attractive polymer-brush interactions. For example, one may construct a brush surface and use it as a polymer separation device. Under weak driving force parallel to the surface, polymers with polymer-brush interaction ε < εC can move quickly, while polymers with ε > εC will move slowly or be trapped on the brush.
      通信作者: 王超, chaowang0606@126.com
    • 基金项目: 国家自然科学基金(批准号: 11604232)和浙江省自然科学基金(批准号: LY20A040004, LY16A040004)资助的课题
      Corresponding author: Wang Chao, chaowang0606@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11604232) and the Natural Science Foundation of Zhejiang Province, China (Grant Nos. LY20A040004, LY16A040004)
    [1]

    Liu J, Wu Y, Shen J, Gao Y, Zhang L, Cao D 2011 Phys. Chem. Chem. Phys. 13 13058Google Scholar

    [2]

    Macosko C W, Guegan P, Khandpur A 1996 Macromolecules 29 5590Google Scholar

    [3]

    Bikiaris D, Panayiotou C 1998 J. Appl. Polym. Sci. 70 1503Google Scholar

    [4]

    Díaz M F, Barbosa S E, Capiati N J 2007 Polymer 48 1058Google Scholar

    [5]

    Neyret S, Ouali L, Candau F, Pefferkorn E 1995 J. Colloid Interface Sci. 176 86Google Scholar

    [6]

    Meredith J C, Johnston K P 1998 Macromolecules 31 5518Google Scholar

    [7]

    Teraoka I 1996 Prog. Polym. Sci. 21 89Google Scholar

    [8]

    Jun S, Mulder B 2006 Proc. Natl. Acad. Sci. U.S.A. 103 12388Google Scholar

    [9]

    Williams M C 2007 Proc. Natl. Acad. Sci. U.S.A. 104 11125Google Scholar

    [10]

    Sheng J F, Luo K F 2015 RSC Advances 5 2056Google Scholar

    [11]

    Douah S, Sabeur S A 2018 Macromol. Theory Simul. 27 1700074Google Scholar

    [12]

    Ziebarth J D, Gardiner A A, Wang Y M, Jeong Y, Ahn J, Jin Y, Chang T 2016 Macromolecules 49 8780Google Scholar

    [13]

    Li B, Sun Z Y, An L J 2015 J. Chem. Phys. 143 024908Google Scholar

    [14]

    Li H, Qian C J, Luo M B 2016 J. Chem. Phys. 144 164901Google Scholar

    [15]

    Yang X. Yang Q H, Fu Y, Wu F, Huang J H, Luo M B 2019 Polymer 172 83Google Scholar

    [16]

    Yang X, Wu F, Hu D D, Zhang S, Luo M B 2019 Chin. Phys. Lett. 36 098202Google Scholar

    [17]

    de Carvalho S J, Metzler R, Cherstvy A G 2015 Soft Matter 11 4430Google Scholar

    [18]

    鲁旸, 瞿立建 2017 高分子学报 3 549Google Scholar

    Lu Y, Zhai L J 2017 Acta Polym. Sin. 3 549Google Scholar

    [19]

    王禹, 章林溪 2008 高分子学报 3 216Google Scholar

    Wang Y, Zhang L X 2008 Acta Polym. Sin. 3 216Google Scholar

    [20]

    王禹, 章林溪 2008 物理学报 57 3281Google Scholar

    Wang Y, Zhang L X 2008 Acta. Phys. Sin. 57 3281Google Scholar

    [21]

    李洪, 艾倩雯, 汪鹏君, 高和蓓, 崔毅, 罗孟波 2018 物理学报 67 168201Google Scholar

    Li H, Ai Q W, Wang P J, Gao H B, Cui Y, Luo M B 2018 Acta. Phys. Sin. 67 168201Google Scholar

    [22]

    Descas R, Sommer J, Blumen A 2004 J. Chem. Phys. 120 8831Google Scholar

    [23]

    Luo M B 2008 J. Chem. Phys. 128 044912Google Scholar

    [24]

    Yang Q H, Luo M B 2016 Sci. Rep. 6 37156Google Scholar

    [25]

    Cerda J J, Sintes T 2005 Biophys. Chem. 115 277Google Scholar

    [26]

    Li H, Gong B, Qian C J, Luo M B 2015 Soft Matter 11 3222Google Scholar

    [27]

    Yang Q H, Yang X, Luo M B 2019 Polymer 180 121677Google Scholar

    [28]

    Milchev A, Egorov S A, Binder K 2014 Soft Matter 10 5974Google Scholar

    [29]

    Tessier F, Labrie J, Slater G W 2002 Macromolecules 35 4791Google Scholar

    [30]

    Cherstvy A G, Metzler R 2014 Phys. Rev. E 90 012134Google Scholar

    [31]

    Morrin G T, Schwartz D K 2018 Macromolecules 51 1207Google Scholar

  • 图 1  模型示意图

    Fig. 1.  A sketch of the simulation system.

    图 2  d取不同值时, 自由链吸附链节数mad随吸引作用能ε的变化, 其中Nf = 50, Nb = 10. 插图: Nfd不同时, 比值mad/madsε的变化

    Fig. 2.  The number of the adsorbed segment of the free polymer mad as a function of the adsorption energy ε for different d, where Nf = 50 and Nb = 10. The inset presents the dependence of the ratio mad/mads on ε for different Nf and d.

    图 3  d取不同值时, 自由链均方回转半径$R_{\rm{g}}^{\rm{2}}$ε的变化, 其中Nf = 50, Nb = 10. 插图: $R_{\rm{g}}^{\rm{2}}$的最小值${(R_{\rm{g}}^{\rm{2}})_{\min }}$d的变化

    Fig. 3.  The mean square radius of the free polymer $R_{\rm{g}}^{\rm{2}}$ as a function of ε, where Nf = 50 and Nb = 10. The inset presents the dependence of the minimum of $R_{\rm{g}}^2$, ${(R_{\rm{g}}^{\rm{2}})_{\min }}$, on d.

    图 4  ε取不同值时, 自由链及分子刷链节在z方向上的分布φpφb, 其中Nf = 50, Nb = 10, d = 12

    Fig. 4.  The distribution of the segment of the free polymer and the polymer brush, φpandφb, in the z direction for three different ε, where Nf = 50, Nb = 10 and d = 12.

    图 5  吸引作用能ε不同时自由链质心均方位移(Δr)2随时间t的演化, 其中Nf = 50, Nb = 10, d = 10

    Fig. 5.  The evolution of the mean square displacement of the center of mass of the free polymer (Δr)2 for different ε, where Nf = 50, Nb = 10 and d = 10.

    图 6  d取不同值时, βε的变化, 其中Nf = 50, Nb = 10. 插图: 弱吸附作用下(ε = 1), 自由链扩散系数D随链长度Nf的变化

    Fig. 6.  The dependence of β on ε for different d, where Nf = 50 and Nb = 10. The inset presents the dependence of the diffusion coefficient D on Nf at small ε = 1.

    图 7  自由链吸附链节数mfa和自由链-分子刷链节接触对数mfb随时间t的演化, 其中Nf = 50, Nb = 10, d = 10, ε = 3. 插图: (a)弛豫函数qfa(t)和qfb(t)随时间的演化; (b)自由链吸附时间τfa以及自由链-分子刷链节接触对数弛豫时间τfb随吸引作用能ε的变化, 其中Nf = 50, Nb = 10, d = 10

    Fig. 7.  The evolution of the number of the adsorbed segment of the free polymer (mfa) and that of the number the segment of polymer brush contacting with the free polymer (mfb), where Nf = 50, Nb = 10, d = 10 and ε = 3. The insets: (a) The evolution of the relaxation function qfa(t) and qfb(t); (b) the dependence of the adsorption time τfa of the free polymer and the relaxation time of the number of segment of brush contacting with the free polymer τfb on the adsorption energy ε, where Nf = 50, Nb = 10 and d = 10.

  • [1]

    Liu J, Wu Y, Shen J, Gao Y, Zhang L, Cao D 2011 Phys. Chem. Chem. Phys. 13 13058Google Scholar

    [2]

    Macosko C W, Guegan P, Khandpur A 1996 Macromolecules 29 5590Google Scholar

    [3]

    Bikiaris D, Panayiotou C 1998 J. Appl. Polym. Sci. 70 1503Google Scholar

    [4]

    Díaz M F, Barbosa S E, Capiati N J 2007 Polymer 48 1058Google Scholar

    [5]

    Neyret S, Ouali L, Candau F, Pefferkorn E 1995 J. Colloid Interface Sci. 176 86Google Scholar

    [6]

    Meredith J C, Johnston K P 1998 Macromolecules 31 5518Google Scholar

    [7]

    Teraoka I 1996 Prog. Polym. Sci. 21 89Google Scholar

    [8]

    Jun S, Mulder B 2006 Proc. Natl. Acad. Sci. U.S.A. 103 12388Google Scholar

    [9]

    Williams M C 2007 Proc. Natl. Acad. Sci. U.S.A. 104 11125Google Scholar

    [10]

    Sheng J F, Luo K F 2015 RSC Advances 5 2056Google Scholar

    [11]

    Douah S, Sabeur S A 2018 Macromol. Theory Simul. 27 1700074Google Scholar

    [12]

    Ziebarth J D, Gardiner A A, Wang Y M, Jeong Y, Ahn J, Jin Y, Chang T 2016 Macromolecules 49 8780Google Scholar

    [13]

    Li B, Sun Z Y, An L J 2015 J. Chem. Phys. 143 024908Google Scholar

    [14]

    Li H, Qian C J, Luo M B 2016 J. Chem. Phys. 144 164901Google Scholar

    [15]

    Yang X. Yang Q H, Fu Y, Wu F, Huang J H, Luo M B 2019 Polymer 172 83Google Scholar

    [16]

    Yang X, Wu F, Hu D D, Zhang S, Luo M B 2019 Chin. Phys. Lett. 36 098202Google Scholar

    [17]

    de Carvalho S J, Metzler R, Cherstvy A G 2015 Soft Matter 11 4430Google Scholar

    [18]

    鲁旸, 瞿立建 2017 高分子学报 3 549Google Scholar

    Lu Y, Zhai L J 2017 Acta Polym. Sin. 3 549Google Scholar

    [19]

    王禹, 章林溪 2008 高分子学报 3 216Google Scholar

    Wang Y, Zhang L X 2008 Acta Polym. Sin. 3 216Google Scholar

    [20]

    王禹, 章林溪 2008 物理学报 57 3281Google Scholar

    Wang Y, Zhang L X 2008 Acta. Phys. Sin. 57 3281Google Scholar

    [21]

    李洪, 艾倩雯, 汪鹏君, 高和蓓, 崔毅, 罗孟波 2018 物理学报 67 168201Google Scholar

    Li H, Ai Q W, Wang P J, Gao H B, Cui Y, Luo M B 2018 Acta. Phys. Sin. 67 168201Google Scholar

    [22]

    Descas R, Sommer J, Blumen A 2004 J. Chem. Phys. 120 8831Google Scholar

    [23]

    Luo M B 2008 J. Chem. Phys. 128 044912Google Scholar

    [24]

    Yang Q H, Luo M B 2016 Sci. Rep. 6 37156Google Scholar

    [25]

    Cerda J J, Sintes T 2005 Biophys. Chem. 115 277Google Scholar

    [26]

    Li H, Gong B, Qian C J, Luo M B 2015 Soft Matter 11 3222Google Scholar

    [27]

    Yang Q H, Yang X, Luo M B 2019 Polymer 180 121677Google Scholar

    [28]

    Milchev A, Egorov S A, Binder K 2014 Soft Matter 10 5974Google Scholar

    [29]

    Tessier F, Labrie J, Slater G W 2002 Macromolecules 35 4791Google Scholar

    [30]

    Cherstvy A G, Metzler R 2014 Phys. Rev. E 90 012134Google Scholar

    [31]

    Morrin G T, Schwartz D K 2018 Macromolecules 51 1207Google Scholar

  • [1] 赵明慧, 刘忠军, 姬帅, 刘晨, 敖庆波. 超临界氮气在单壁碳纳米管内吸附行为的GCMC模拟研究. 物理学报, 2022, 71(22): 220201. doi: 10.7498/aps.71.20220765
    [2] 孙立望, 李洪, 汪鹏君, 高和蓓, 罗孟波. 利用神经网络识别高分子链在表面的吸附相变. 物理学报, 2019, 68(20): 200701. doi: 10.7498/aps.68.20190643
    [3] 李洪, 艾倩雯, 汪鹏君, 高和蓓, 崔毅, 罗孟波. 外力驱动作用下高分子链在表面吸附性质的计算机模拟. 物理学报, 2018, 67(16): 168201. doi: 10.7498/aps.67.20180468
    [4] 王超, 陈英才, 周艳丽, 罗孟波. 两嵌段高分子链在周期管道内扩散的Monte Carlo模拟. 物理学报, 2017, 66(1): 018201. doi: 10.7498/aps.66.018201
    [5] 庞宗强, 张悦, 戎舟, 江兵, 刘瑞兰, 唐超. 利用扫描隧道显微镜研究水分子在Cu(110)表面的吸附与分解. 物理学报, 2016, 65(22): 226801. doi: 10.7498/aps.65.226801
    [6] 林文强, 徐斌, 陈亮, 周峰, 陈均朗. 双酚A在氧化石墨烯表面吸附的分子动力学模拟. 物理学报, 2016, 65(13): 133102. doi: 10.7498/aps.65.133102
    [7] 伊丁, 武镇, 杨柳, 戴瑛, 解士杰. 有机分子在铁磁界面处的自旋极化研究. 物理学报, 2015, 64(18): 187305. doi: 10.7498/aps.64.187305
    [8] 郑晖, 张崇宏, 陈波, 杨义涛, 赖新春. 氦离子低温预辐照对不锈钢中氦泡生长抑制作用的Monte Carlo模拟研究. 物理学报, 2014, 63(10): 106102. doi: 10.7498/aps.63.106102
    [9] 刘秀英, 李晓凤, 张丽英, 樊志琴, 马兴科. 甲烷在不同分子筛中吸附的理论对比研究. 物理学报, 2012, 61(14): 146802. doi: 10.7498/aps.61.146802
    [10] 周宇璐, 李仁顺, 张宝玲, 邓爱红, 侯氢. 材料中He深度分布演化的Monte Carlo模拟研究. 物理学报, 2011, 60(6): 060702. doi: 10.7498/aps.60.060702
    [11] 郭宝增, 张锁良, 刘鑫. 钎锌矿相GaN电子高场输运特性的Monte Carlo 模拟研究. 物理学报, 2011, 60(6): 068701. doi: 10.7498/aps.60.068701
    [12] 高茜, 娄晓燕, 祁阳, 单文光. Zn1-xMnxO纳米薄膜磁有序性的Monte Carlo模拟. 物理学报, 2011, 60(3): 036401. doi: 10.7498/aps.60.036401
    [13] 黄平, 杨春. TiO2分子在GaN(0001)表面吸附的理论研究. 物理学报, 2011, 60(10): 106801. doi: 10.7498/aps.60.106801
    [14] 颜超, 段军红, 何兴道. 低能原子沉积在Pt(111)表面的分子动力学模拟. 物理学报, 2010, 59(12): 8807-8813. doi: 10.7498/aps.59.8807
    [15] 姚文静, 王楠. Ni-15%Mo合金熔体热物理性质的Monte Carlo模拟. 物理学报, 2009, 58(6): 4053-4058. doi: 10.7498/aps.58.4053
    [16] 黄朝军, 刘亚锋, 龙姝明, 孙彦清, 吴振森. 烟尘中电磁波传输特性的Monte Carlo模拟. 物理学报, 2009, 58(4): 2397-2404. doi: 10.7498/aps.58.2397
    [17] 肖 沛, 张增明, 孙 霞, 丁泽军. 投影电子束光刻中电子穿透掩膜的Monte Carlo模拟. 物理学报, 2006, 55(11): 5803-5809. doi: 10.7498/aps.55.5803
    [18] 高国良, 钱昌吉, 钟 瑞, 罗孟波, 叶高翔. 非均质基底表面上团簇生长的Monte Carlo模拟. 物理学报, 2006, 55(9): 4460-4465. doi: 10.7498/aps.55.4460
    [19] 张现仁, 沈志刚, 陈建峰, 汪文川. 乙烷在中孔分子筛MCM-41中吸附的计算机分子模拟. 物理学报, 2003, 52(1): 163-168. doi: 10.7498/aps.52.163
    [20] 尚也淳, 张义门, 张玉明. 6H-SiC电子输运的Monte Carlo模拟. 物理学报, 2000, 49(9): 1786-1791. doi: 10.7498/aps.49.1786
计量
  • 文章访问数:  3310
  • PDF下载量:  63
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-18
  • 修回日期:  2020-04-27
  • 上网日期:  2020-05-25
  • 刊出日期:  2020-08-20

高分子链在分子刷表面吸附的Monte Carlo模拟

    基金项目: 国家自然科学基金(批准号: 11604232)和浙江省自然科学基金(批准号: LY20A040004, LY16A040004)资助的课题

摘要: 用Monte Carlo方法模拟研究了一条自由高分子链在分子刷表面吸附的静态和动态特性. 结果表明, 随着自由链与分子刷之间吸附作用能(ε)的增大, 自由链出现由脱吸附态到吸附态的相转变, 同时链的扩散由正常模式转为亚扩散模式. 临界吸附能(εC)几乎与自由链长度无关, 但随着分子刷链长度的减小或分子刷链间距的增大而不断增大. 在εC附近, 自由链嵌入分子刷内部, 同时链尺寸达到极小, 而当ε $\gg $ εC时, 自由链处于强吸附态, 链节主要分布于分子刷表层, 同时整个吸附动态过程可分为自由链吸附和分子刷扩散两个阶段.

English Abstract

    • 高分子链在固液界面上的吸附一直是物理、化学和生物领域的热门课题, 其不仅与许多化工技术相关而且在生物系统中也广泛存在, 如: 高分子纳米粒子复合材料制备、材料表面改性和润滑、胶体稳定、DNA在病毒内的封装等[1-9]. 因此, 深入研究高分子链吸附的静态和动态行为并揭示相应的物理机理, 对化工技术和生物医学技术的研究有重要意义.

      高分子链在界面上的吸附受链的拓扑结构[10-13]、链的组成和柔性[14-16]、界面的几何形状[17,18]、外加作用[19-21]等众多因素影响. 对于均质界面的情形, 高分子链与界面之间的吸引作用能(ε)是一个重要的调控参数. 研究表明, 当吸引作用能超过临界值(εC)时, 高分子链会发生由脱吸附态到吸附态的相转变. 相应地, 高分子链构型由无规线团状转变为薄饼状[22]. Luo[23]曾根据二级相变有限尺寸标度理论得到了高分子链吸附链节数M(序参量)随链长N以及吸引作用能ε的标度关系, 同时发现利用该标度关系可以很好地解释Monte Carlo模拟结果并由此确定了εC. 这也意味着高分子链在平面附近的吸附相转变为二级相变. 另外, 高分子链的动态行为(如: 吸附时间、扩散和弛豫)也与界面吸引作用能有关[24].

      除吸引作用能之外, 界面的细节结构也会对高分子链吸附产生重要影响. Cerda和Sintes[25]曾考虑吸引条纹和中性条纹周期分布的界面, 发现高分子链的吸附存在两个临界吸附点(临界温度)和三种不同的吸附态. 进一步的研究表明, 高分子链在三种吸附态满足不同的扩散规律[26]. 对于吸附点在界面上随机分布的情形, Yang等[27]发现高分子链的吸附时间随吸附点密度的增大而增大, 但在高密度和低密度区域满足不同的吸附机制, 同时还发现高分子链尺寸有可能随吸附点密度的增大而减小, 这与高分子链在均匀界面的吸附不同.

      Milchev等[28]最近考虑更为复杂的界面—分子刷界面, 发现随着自由高分子链与分子刷之间吸引作用的增强, 自由高分子链也会发生脱吸附到吸附的相转变. 但与固体界面不同的是, 根据分子刷密度和吸引强度不同, 自由高分子链在吸附态既有可能嵌入分子刷内部也有可能只吸附于分子刷表层. 这也意味着, 高分子链在分子刷上吸附的动力学过程可能比固液界面情形更为复杂, 然而目前关于这方面的研究还很有限. 受此启发, 本文用Monte Carlo方法模拟研究了一条自由高分子链在分子刷界面上的吸附, 不仅考察自由高分子链由脱吸附态到吸附态的转变和链节空间分布变化, 还研究了自由高分子链的平衡态扩散过程和吸附动态过程, 发现自由链在脱吸附态和吸附态满足不同的扩散规律, 同时还发现整个吸附动态过程可分为自由链吸附和分子刷扩散两个阶段.

    • 整个模型系统建立在三维简立方格子空间(格子常数为长度单位), 如图1所示. 单条长度为Nf的自由高分子链处在分子刷表面附近. 分子刷由nb条长度为Nb的端点吸附链(接枝链)组成. 接枝链固定端均匀地分布在下界面(z = 0位置)的格点上, 近邻固定端间的距离为d. 分子刷疏密程度取决于d, 即d越小, 分子刷越密集, 反之, 分子刷越稀疏. 另外, 为了防止自由链远离分子刷, 假定在z = D位置还存在一个纯排斥上界面, 即自由链限制在分子刷和上界面之间. 系统在x, y方向上采用周期性边界条件.

      图  1  模型示意图

      Figure 1.  A sketch of the simulation system.

      自由链和接枝链均采用八格点键长涨落模型[29], 具体地, 每个链节占据一个立方体的八个格点, 键连链节间的距离可以有$2, \sqrt 5, \sqrt 6, 3, \sqrt {10} $五种取值(高分子链取这些键长可以自动避免键交叉的情况). 自由链链节之间、分子刷链节之间、自由链与分子刷之间、链节与界面之间均满足排除体积作用, 即空间任一格点都不能同时被两个或多个链节占据, 也不能同时被链节和界面占据. 除此之外, 还假定当自由高分子链节中心与分子刷链节中心之间距离为2 (最近距离)时, 自由链与分子刷之间存在吸附作用能 ε.

      模拟开始时, 先在下界面上生成分子刷, 然后在分子刷和上界面之间生成自由高分子链. 系统生成后, 先设定自由链与分子刷之间吸附作用能为零(ε = 0), 并让系统经过长时间的弛豫达到平衡态; 然后开启自由链与分子刷之间的吸附作用, 模拟并考察自由链的吸附和运动. 在模拟中, 分子链的运动通过单链节尝试运动实现. 链节尝试运动采用Monte Carlo重要性抽样算法, 具体如下: 随机选择自由链或分子刷(除吸附端点外)上一个链节并使其随机地向近邻格点移动一个单位长度, 若链节新位置满足模型条件(如: 排除体积条件和键长条件), 则新位置被接收的概率为$\min (1, {{\rm{e}}^{ - \Delta E/kT}})$, 其中ΔE为链节位置移动所引起的能量增量. 定义一个Monte Carlo步长(Monte Carlo step, MCS)为Monte Carlo模拟的时间单位. 在一个时间单位内共有$N_{\rm f} + n_{\rm b}N_{\rm b}$次链节尝试运动.

      为了使自由链能尽快与分子刷接触, 选择分子刷表面到上界面之间的距离略大于自由链的自由回转半径Rg0 (Rg0为自由链在自由空间的回转半径). 在模拟中固定分子刷接枝链总条数nb = 100, 主要考察分子刷链长度Nb、间距d以及吸附作用能ε对自由链的吸附、构型和运动的影响.

    • 首先研究了自由链吸附链节数mad(即与分子刷接触的自由链链节数)随吸附作用能ε的变化, 如图2所示. 可以看出, madε的变化存在三个不同的区域: 1)当ε比较小时, mad → 0, 即自由链处于脱吸附态; 2)当ε比较大时, 自由链处于吸附态, mad趋于饱和值mads; 3)当ε取中等值时, madε的变化存在一个明显的转变区, 在该区域mad由0快速增大为mads. 在吸附态, 自由链可能同时与多条接枝链相吸附, 因此自由链上必然有部分链节介于不同接枝链之间而不被吸附, 称之为桥链节. 随着d的增大, 桥链节数目将不断增大, 从而导致mads不断减小, 如图2所示. mads还与接枝链长度Nb有关. 随着Nb增大, 自由链周围分子刷链节不断增多, 这将使得自由链被吸附链节数mads不断增大(结果未给出). 为了描述madε变化图中转变区的位置, 定义比值mad/mads = 0.5对应的ε为临界吸附能εC, 如图2插图所示. 可以看出, 临界吸附能εC几乎与自由链长度Nf无关, 但随接枝链间距d的增大而增大. 同时可以发现εC随接枝链长度Nb的增大而减小(结果未给出). 随着d增大或Nb减小, 自由链周围分子刷链节不断减少, 因此需要更强的吸附作用才能使自由链吸附在分子刷上, 从而导致εC不断增大, 这与链滴理论(the blob theory)得到的结果定性一致[28]. Luo[23]的研究表明, 高分子链在吸附平面附近的相转变属于二级相变. 虽然分子刷界面与吸附平面在微观结构上有很大不同, 但高分子链在两种界面附近由脱吸附到吸附的宏观转变却非常相似. 因此, 我们认为高分子链在分子刷表面附近的脱吸附-吸附相变类型也为二级相变.

      图  2  d取不同值时, 自由链吸附链节数mad随吸引作用能ε的变化, 其中Nf = 50, Nb = 10. 插图: Nfd不同时, 比值mad/madsε的变化

      Figure 2.  The number of the adsorbed segment of the free polymer mad as a function of the adsorption energy ε for different d, where Nf = 50 and Nb = 10. The inset presents the dependence of the ratio mad/mads on ε for different Nf and d.

      图3给出了自由链均方回转半径$R_{\rm{g}}^{\rm{2}}$随吸附能ε的变化. 由图可以看出$R_{\rm{g}}^{\rm{2}}$ε呈现复杂变化: 当ε < εC时, $R_{\rm{g}}^{\rm{2}} = R_{{\rm{g}}0}^2$(这里$R_{{\rm{g}}0}^2$为自由链在自由空间的均方回转半径); 在εC附近, $R_{\rm{g}}^{\rm{2}}$ε的增大快速减小并达到极小值($ R_{\rm{g}}^{\rm{2}} $)min; 当ε > εC时, $R_{\rm{g}}^{\rm{2}}$ε的增大又逐渐增大. $R_{\rm{g}}^{\rm{2}}$的最小值${(R_{\rm{g}}^{\rm{2}})_{\min }}$与接枝链间距d有关. 随着d的增大, ${(R_{\rm{g}}^{\rm{2}})_{\min }}$先减小然后不断增大, 如图3插图所示.

      图  3  d取不同值时, 自由链均方回转半径$R_{\rm{g}}^{\rm{2}}$ε的变化, 其中Nf = 50, Nb = 10. 插图: $R_{\rm{g}}^{\rm{2}}$的最小值${(R_{\rm{g}}^{\rm{2}})_{\min }}$d的变化

      Figure 3.  The mean square radius of the free polymer $R_{\rm{g}}^{\rm{2}}$ as a function of ε, where Nf = 50 and Nb = 10. The inset presents the dependence of the minimum of $R_{\rm{g}}^2$, ${(R_{\rm{g}}^{\rm{2}})_{\min }}$, on d.

      为了理解自由链尺寸的变化, 分别研究了ε < εC(脱附态), εεC(临界吸附能附近)以及ε > εC(吸附态)时自由链链节在z方向上的分布, 典型结果如图4所示, 其中Nf = 50, Nb = 10, d = 12, 对应的临界吸附能εC ≈ 2.65. 当ε = 1.6时, 自由链处于脱附态, 链节主要分布于分子刷外部(如图4(a)所示); 当ε = 2.6时, 自由链处于临界吸附状态, 链节主要分布于分子刷内部, 即自由链已嵌入到分子刷内部(如图4(b)所示); 当ε = 3.6时, 自由链处于强吸附状态, 此时自由链一旦与分子刷表层接触就很难脱离, 从而导致链节主要分布于分子刷表层(如图4(c)所示), 这与文献[28]的模拟结果一致. 结合自由链链节分布规律, 可以定性理解自由链尺寸$R_{\rm{g}}^{\rm{2}}$ε的非单调变化, 即: 当ε < εC时, 自由链处于分子刷外部, 因此有$R_{\rm{g}}^{\rm{2}} = R_{{\rm{g0}}}^{\rm{2}}$; 当εεC时, 自由链嵌入分子刷内部, 此时在分子刷分子链包围和挤压下自由链尺寸明显减小; 当ε > εC时, 自由链吸附于分子刷表层, 此时为了获得较多的吸附作用能自由链会在分子刷表层形成准二维分布, 从而导致链尺寸增大. 在εC附近, 自由链嵌入分子刷内部的难易程度与接枝链间距d有关, 即d越小, 分子刷越稠密, 自由链就越难嵌入到分子刷内部; 反之, 分子刷越稀疏, 自由链就越容易嵌入到分子刷内部, 但其在分子刷内部受到的包围和挤压就越弱, 从而导致${(R_{\rm{g}}^{\rm{2}})_{\min }}$d出现非单调变化, 如图3插图所示.

      图  4  ε取不同值时, 自由链及分子刷链节在z方向上的分布φpφb, 其中Nf = 50, Nb = 10, d = 12

      Figure 4.  The distribution of the segment of the free polymer and the polymer brush, φpandφb, in the z direction for three different ε, where Nf = 50, Nb = 10 and d = 12.

      接下来, 研究了自由链的扩散. 图5给出了自由链质心均方位移(Δr)2随时间t的演化, 其中Nf = 50, Nb = 10, d = 10, 对应临界吸附能εC ≈ 2.55. 结果表明, 在脱吸附态(ε < εC), 自由链扩散快慢几乎与ε无关; 但在吸附态(ε > εC), 随着ε增大, 自由链扩散进行得越来越慢. 特别地, 当ε比较大时, 自由链处于强吸附状态, 自由链节一旦与分子刷链节吸附就很难脱离, 因此在整个模拟时间内自由链只有微小移动, 如图5所示. 另外发现, 对于任意给定的吸附能ε, (Δr)2随时间t的变化均存在简单的幂次关系${(\Delta r)^2} = 2 D{t^\beta }$. 当ε < εC时, β = 1且扩散系数D与自由链长度Nf之间存在反比关系, 即$D \propto N_{\rm{f}}^{{\rm{ - 1}}}$(如图6插图所示), 这说明自由链的扩散属于正常扩散; 当εεC时, β < 1, 这意味着自由链的扩散属于亚扩散[30]. βε的增大快速减小并最终达到稳定值βe. 在强吸附条件下, 随着d增大, 自由链越来越不易于在接枝链之间转移, 从而导致βed的增大而减小, 如图6所示.

      图  5  吸引作用能ε不同时自由链质心均方位移(Δr)2随时间t的演化, 其中Nf = 50, Nb = 10, d = 10

      Figure 5.  The evolution of the mean square displacement of the center of mass of the free polymer (Δr)2 for different ε, where Nf = 50, Nb = 10 and d = 10.

      图  6  d取不同值时, βε的变化, 其中Nf = 50, Nb = 10. 插图: 弱吸附作用下(ε = 1), 自由链扩散系数D随链长度Nf的变化

      Figure 6.  The dependence of β on ε for different d, where Nf = 50 and Nb = 10. The inset presents the dependence of the diffusion coefficient D on Nf at small ε = 1.

      亚扩散是反常扩散(包括亚扩散(0 < β < 1)和超扩散(β > 1))的一种. 尽管反常扩散的机制可能有很多, 但Cherstvy和Metzler[30]认为反常扩散主要存在于以下三种情形: 1)分形环境中的扩散; 2)黏弹环境中的扩散; 3)粒子在运动过程中存在一系列明显的吸附受限事件. 对于第三种情形, β满足0 < β < 1, 即扩散为亚扩散. 在我们的模型中, 当εεC时, 自由链被分子刷吸附, 自由链上任意链节在运动过程中都会不断地与分子刷链节吸附(受限)和解吸附(脱离受限), 每次吸附到解吸附都需要经历一定的受限时间. 因此, 自由链的扩散属于反常扩散的第三种情形, 即亚扩散. 此时, β的取值与吸附点位的分布密度有关. Morrin和Schwartz[31]曾研究了聚乙二醇(PEG)分子链在疏水三甲基硅烷(TMS)表面的扩散(其中TMS对PEG分子链有吸附作用), 发现目标PEG分子链的扩散为亚扩散且标度指数β与其他PEG在TMS表面的覆盖度有关. PEG浓度越高, 覆盖度越大, TMS对目标PEG分子链的吸附点位就越少. 实验结果表明, β随PEG浓度(覆盖度)呈非单调变化. 在中等浓度(中等覆盖度)区域, β值在0.4—0.7之间且随PEG浓度(覆盖度)的增大而减小. 在我们的模型中, 分子刷链间距d反映分子刷的疏密程度, 即对自由链吸附点位的疏密程度. d越大, 吸附点位就越稀疏. 因此, 指数βd有关, 但没有一个普适值. 结果表明, 当ε $\gg $ εC时, β在0.6附近且随d的增大而略有减小, 如图6所示. 这与实验结果存在一定的相似性.

      最后研究了ε > εC条件下自由链吸附的动态过程. 在吸附过程中, 自由链吸附链节数(mfa)会不断增大, 同时自由链-分子刷链节接触对数(mfb)也会不断增大. 图7给出ε = 3时mfamfb随时间的演化, 其中Nf = 50, Nb = 10, d = 10, 对应临界吸附能εC ≈ 2.55. 由图7可以看出, 自由链的吸附是一个快速的过程, mfamfb均随时间t快速增大并分别趋于饱和值mfaemfbe. 由于自由链的一个链节可以与多个分子刷链节接触, 因此当吸附达到稳定时, mfbe总是大于mfae, 如图7所示.

      图  7  自由链吸附链节数mfa和自由链-分子刷链节接触对数mfb随时间t的演化, 其中Nf = 50, Nb = 10, d = 10, ε = 3. 插图: (a)弛豫函数qfa(t)和qfb(t)随时间的演化; (b)自由链吸附时间τfa以及自由链-分子刷链节接触对数弛豫时间τfb随吸引作用能ε的变化, 其中Nf = 50, Nb = 10, d = 10

      Figure 7.  The evolution of the number of the adsorbed segment of the free polymer (mfa) and that of the number the segment of polymer brush contacting with the free polymer (mfb), where Nf = 50, Nb = 10, d = 10 and ε = 3. The insets: (a) The evolution of the relaxation function qfa(t) and qfb(t); (b) the dependence of the adsorption time τfa of the free polymer and the relaxation time of the number of segment of brush contacting with the free polymer τfb on the adsorption energy ε, where Nf = 50, Nb = 10 and d = 10.

      为了确定自由链的吸附时间, 定义自由链吸附链节数弛豫函数qfa(t):

      ${q_{{\rm{fa}}}}(t) = \frac{{{m_{{\rm{fae}}}} - {m_{{\rm{fa}}}}(t)}}{{{m_{{\rm{fae}}}}}},$

      qfa(t)随时间t呈指数衰减且可表示为[15]

      $ {q_{{\rm{fa}}}}(t) \propto \exp ( - t/{\tau _{{\rm{fa}}}}), $

      其中τfa为自由链吸附时间(或自由链吸附链节数弛豫时间). 模拟结果表明, 在0.4 < qfa < 0.95范围内qfa(t)随t的变化可以很好地用(2)式进行拟合(如图7插图(a)所示), 由此可以得到自由链吸附时间τfa. 类似地, 定义自由链-分子刷链节接触对数弛豫函数qfb(t)

      ${q_{{\rm{fb}}}}(t) = \frac{{{m_{{\rm{fbe}}}} - {m_{{\rm{fb}}}}(t)}}{{{m_{{\rm{fbe}}}}}}.$

      在0.4 < qfb < 0.95范围内, 用公式$ {q_{{\rm{fb}}}}(t) \propto \exp ( - t/{\tau _{{\rm{fb}}}}) $对模拟结果进行拟合(如图7插图(a)所示), 进而得到自由链-分子刷链节接触对数弛豫时间τfb.

      图7插图(b)给出了自由链吸附时间τfa以及自由链-分子刷链节接触对数弛豫时间τfb随吸附作用能ε的变化. 由图可以看出, 随着吸附作用能ε增大, τfaτfb均单调减小, 这说明增大吸附作用能可以有效地加快自由链的吸附过程. 同时, 随着d减小或Nb增大, 自由链周围的分子刷链节增多, 自由链越容易吸附在分子刷上, 从而导致自由链的吸附过程加快(结果未给出). 另外, 我们发现在相同参数条件下τfb总是比τfa大, 这意味着自由链吸附在分子刷上后并没有马上达到稳定吸附, 周围的分子刷链节仍不断地向其靠近并最终形成稳定吸附. 也就是说, 整个吸附过程包含了自由链的吸附和分子刷向自由链的扩散两种运动, 而且可分为两个阶段: 第一阶段, 自由链的吸附占主导, 自由链在分子刷吸附作用下不断靠近并吸附于分子刷上; 第二阶段, 分子刷的扩散占主导, 分子刷链节逐渐向自由链靠近并形成稳定吸附.

    • 用Monte Carlo方法模拟研究了一条自由高分子链在分子刷上的吸附. 模拟结果表明, 随着自由链与分子刷之间吸附能ε的增大, 自由链出现由脱吸附态到吸附态的转变. 临界吸附能εC几乎与自由链长度无关, 但随着分子刷链长Nb的减小或分子刷链间距d的增大而增大. 当ε < εC时, 自由链处于脱附态, 链节主要分布于分子刷之外, 链的构象尺寸与自由空间中的结果一致; 在εC附近, 自由链嵌入分子刷内部, 尺寸达到极小; 当$ \varepsilon \gg \varepsilon_{\rm C} $时, 自由链处于强吸附态, 链节主要分布于分子刷表层, 同时整个吸附过程可分为自由链的吸附和分子刷的扩散两个阶段. 另外, 随着ε的增大, 自由链的扩散在εC附近由正常模式转变为亚扩散模式, 扩散过程迅速减慢. 这种扩散模式(或扩散快慢)的转变可能有助于实验上利用分子刷表面对与分子刷间存在不同吸引作用的高分子多链体系进行分离, 如: 施加外力(平行于基底)使自由高分子链发生定向运动, 处于脱附态的自由链将很快运动, 而处于吸附态的自由链则缓慢运动, 这种运动快慢的差异最终会导致处于脱附态和吸附态的自由链在空间上出现明显分离, 从而达到链分离的目的.

参考文献 (31)

目录

    /

    返回文章
    返回