蘇志成，虎 偉

（北京師范大學 物理學系，北京 100875）

球形聚電解質刷的自洽場模擬

蘇志成，虎 偉

（北京師范大學 物理學系，北京 100875）

利用自洽場理論（self-consistent filed theory，SCFT）研究了聚電解質體系鹽濃度、鏈接枝密度及帶電比對球形聚電解質刷相行為的影響，并和標度理論3種刷型的標度關系進行了對照.自洽場理論計算結果不僅給出了3種類型聚電解質刷的行為趨勢，而且可以為不同類型重疊區(qū)的刷行為提供更多細節(jié)信息.離子的密度分布表明體系始終是局部電中性的，加鹽刷與中性刷有相似之處，而滲透壓刷的接枝鏈并不是均勻伸展.鹽濃度下降時，單體自由端分布移向外層并出現(xiàn)排空區(qū).

自洽場理論；聚電解質；鹽濃度；接枝密度；構象

聚電解質刷是由聚電解質鏈通過物理吸附作用［1-3］或化學鍵［4，5］接枝到平面或曲面基底而形成的刷狀結構.聚電解質鏈之間的相互作用使得鏈向接枝面以外的空間伸展.在過去的20年之中，聚電解質由于在制備納米管的陣列［6］、光子晶體［7］等規(guī)整的人工微結構方以及藥物輸送［8］和膠體制備［9］等方面的潛在應用價值而得到了極大的關注.

根據(jù)聚電解質鏈接枝面的幾何形狀的不同，可以分為平面聚電解質刷、柱形聚電解質刷及球形聚電解質刷.球形聚電解質刷從一開始就得到持續(xù)的關注，本文將在理論上對這一體系進行分析，以期能對實驗和模擬研究提供指導.

聚電解質刷相行為的理論理解，對其具體應用是十分必要的.聚電解質刷的相行為受到鏈長度、接枝密度、溶劑的選擇性及接枝面的曲率等眾多因素的影響.Marko和 Witten最先通過理論去檢驗平面聚電解質刷的構象平衡問題，并預測膠束到微管或層狀相的轉變［10，11］.Lai利用蒙特卡羅模擬研究了雙組份刷，研究發(fā)現(xiàn)了少數(shù)鏈從基底向外更加延伸的層狀相［12］.

本文對球形聚電解質刷的自洽場方程組實施實空間數(shù)值求解，研究體系鹽濃度、聚電解質鏈接枝密度、單體帶電比例對聚電解質刷相行為的影響.

1 自洽場理論

在本文所要研究的體系中，np個長度為N的聚電解質鏈接枝均勻地接枝到半徑為Rc的球上，如圖1所示.面接枝密度（單位面積上的鏈數(shù)）為σ，體系共有nM個單體［M指正離子（+），負離子（-）或溶劑（S）］.高分子單體、正負離子的化合價分別用vp、v+和 v-表示.

圖1 球形聚電解質刷示意圖.Rc接枝球半徑，R刷厚度，接枝刷在L處與鹽溶液接觸

這里只給出 SCFT的主要方程，具體推導過程可以參考文獻［13，14］.首先對本文中用到的符號加以說明：N表示接枝鏈的長度，f表示單體帶電比，σ表示接枝密度，b表示單位單體長度，R表示聚電解質刷半徑.ρj（r）和 φj（r）分別表示數(shù)密度和體積比（j=P，S，+，-）.φj（r）=ρj（r）/ρ0，其中 ρ0=（nPN+ns）/V=1/b3.正負離子的本體濃度為 φb±.假設體系是不可壓縮的，且將小離子視為點電荷，于是有：弗洛里-哈金斯（Flory-Huggins）相互作用參數(shù)用χ表示，用來描述聚電解質溶液中不同成分之間的相互作用.kBT是玻爾茲曼常數(shù)與熱力學溫度的乘積.ωj和ψj分別表示場和靜電勢（j=P，S，+，-）.平均場方程通過求體系自由能的最小值獲得.

以kBT為單位的自由能為

其中

其中R（s）表示接枝高分子鏈的構象，即第s個鏈節(jié)在空間可能所在的位置的集合.

令δF/δφj=0，δF/δωj=0和 δF/δψ=0，得到SCFT方程組如下：

其中q（r，s）表示從球表面開始經(jīng)過s步終止于r位置的子鏈的統(tǒng)計比重；q+（r，s）表示除球面從任意位置開始經(jīng)過s步終止于r位置的子鏈的末端分布函數(shù).它們分別滿足如下的擴散方程：

為了求解以上的擴散方程式（11）、（12）和玻爾茲曼方程式（10），需要有適當?shù)倪吔鐥l件和初值.邊界條件為：

初值條件：

利用文獻［15］的方法，以上方程可以在實空間數(shù)值求解.擴散方程可以用 Crank-Nicolson的求解方法，其他非線性方程用Broyden方法求解［16］.在大多數(shù)情況下，s和 r的離散步長分別為 Δs=0.001N和 Δr=0.25b，L設為2Nb.

末端分布由下式給出：

通過平均計算d=r-Rc，可以得到接枝刷半徑的表達式：

2 標度理論

Zhulina和O.V.Borisov等通過理論分析給出了球形聚電解質刷的標度關系［17］，在此回顧一下結果，并在后面的討論中加以對比.

基于接枝濃度、聚電解質鏈的帶電比和體系中的離子強度，標度理論可以區(qū)別出以下3種主要的類型.

1）Pincus類型.特點：接枝密度比較低、體系中大多反離子不在刷區(qū)域.刷的結構取決于庫倫能和熵能之間的平衡.

2）Osmotic類型（滲透壓型）.特點：接枝密度高、大部分反離子在刷內部.刷的結構取決于反離子的滲透壓和熵能之間的平衡.

3）Salted類型（加鹽型）.特點：體系中加入的鹽濃度超過反離子的平均濃度.刷的結構取決于滲透壓和熵之間的相互作用.

令N表示接枝鏈的長度，f表示單體帶電比，σ表示接枝密度，b表示單位單體長度，R表示聚電解

質刷半徑.lB=e2/εkBT表示 Bjerrum長度，其中ε表示電介常量.

對于Pincus型刷，接枝濃度比較低，帶電比低，則由于庫倫吸引力不足，反離子離開刷區(qū)域.聚電解質刷的平衡半徑由庫倫相互作用能

和構象自由能

的平衡決定，即

對于滲透壓（Osmotic）型刷，接枝密度高，大部分反離子在刷內區(qū)域，屏蔽了刷內庫倫排斥作用，刷內反離子的滲透壓成為接枝鏈伸展的主要動力.刷的半徑取決于滲透壓和支鏈彈力的平衡.

滲透壓力為

彈力可以表示為

二力平衡可得

最值得注意的是，多分枝的球形聚電解質刷半徑與分枝數(shù)（即接枝密度σ）是獨立的.

在溶液中加濃度為ΦS的1∶1鹽的情況下，加鹽型刷的半徑取決于所有自由離子的滲透壓和支鏈彈力的平衡.

自由離子的滲透壓表示為

與彈力平衡，可得

這一結果亦可由下面的方法導出：

體系Debye屏蔽長度κ-1為

由單體間短程相互作用決定的刷半徑為

被屏蔽的庫倫相互作用可以等效為短程體積排除作用，用第二virial系數(shù)表示為將式（23）和式（25）代入式（24），得

與式（22）完全一致.

雖然標度理論可以給出大量有用的信息，但是不能給出聚電解質刷結構的具體細節(jié)，SCFT正好可以做到這點.

3 結果與討論

這一部分給出球形聚電解質刷的SCFT理論的數(shù)值計算結果，并與標度理論進行對比討論.計算過程中，下面這些參數(shù)的值在本文中保持不變.單體長度b=0.5 nm，并在整篇文章中被用為長度單位.室溫下，對于水溶劑而言，Bjerrum長度 lB=0.7 nm.帶電單體只帶一價的負電，所以 v+=-v-=-vP=1，接枝球的半徑為Rc=b.

圖2為聚電解質刷半徑R在不同帶電比f時隨ΦS變化的圖像.圖2a中 R以高鹽濃度 ΦS=0.1時的值 Rhigh歸一化，圖2（b）中 R以低鹽濃度 ΦS=1× 10-6時的值 Rlow歸一化.在圖2（a）和（b）中，當鹽濃度 φbs較低時，都可以發(fā)現(xiàn)一個顯著的平臺，曲線平滑，R=Rlow，相當于滲透壓刷的半徑.在相當大的鹽濃度范圍內，刷半徑對帶電比f的變化也并不敏感（圖2（b）），這與單獨的鏈狀聚電解質不同.這說明反離子在分子間屏蔽作用明顯，帶電的聚電解質鏈之間的庫倫相互作用并沒有表現(xiàn)出為對刷半徑的影響.

隨著鹽濃度的增加，體系的整體的滲透壓增大，聚電解質刷開始收縮，以平衡鹽濃度增加導致的體系熵的增加.而且不同帶電比 f的聚電解質刷收縮的程度不同，如圖2（b）所示，帶電比越高，收縮越快，鹽濃度較高時的刷半徑越小.

由圖2（a）可以看出高鹽濃度時，ln R與 ln ΦS有著很強的線性關系.圖中的黑色斜率線-1/5表示R-Φ-1/5S.這與標度理論對加鹽刷的預測式（22）一致.

圖2 聚電解質刷半徑與鹽濃度關系圖（雙對數(shù)坐標）Rhigh歸一化，以Rlow歸一化，鏈長N=100，接技密度σ=0.1

圖3表示聚電解質刷半徑與接枝密度的關系.如圖3所示，鹽濃度很低為 ΦS=1×10-6時，接枝密度較低時，刷半徑隨著接枝密度增大而增加，并呈現(xiàn)出R～σ1/3的關系，這和標度理論對Pincus刷的討論是一致的（式（17））.與中性刷 R-σ1/4相比，庫倫相互作用使得接枝密度增加對擴大刷半徑的貢獻更加明顯.

圖3 聚電解質刷半徑與接枝密度關系圖（雙對數(shù)坐標）鏈長為N=100

接枝密度較高時，這一趨勢受到阻礙，R-σ曲線趨于平緩，聚電解質刷逐漸成為滲透壓型刷，即式（20）：

即刷半徑與接枝密度無顯著關系.考慮到刷內層聚電解質鏈所帶電荷并不完全被反離子所屏蔽，所以刷半徑隨接枝密度 σ增大會緩慢變大，這與 SCFT計算結果一致.

從前面的討論已知，高鹽濃度使聚電解質刷收縮.在高鹽濃度下（ΦS=0.1或 ΦS=0.01），刷半徑隨著接枝密度的增加緩慢增大，斜率約為 1/5，即R-σ1/5，這符合加鹽刷的行為（式（22））.在低鹽濃度下（ΦS=10-5或 ΦS=10-4），接枝密度較低時，ln R與 ln σ有著很顯著的線性關系，符合 Pincus刷標度關系R-σ1/3.刷內區(qū)域反離子的濃度取決于接枝密度σ和帶電比f，而體系鹽濃度和接枝密度較低時，刷內區(qū)域反離子的濃度較小，庫倫相互作用起主導作用，即為 Pincus刷.所以低鹽濃度下，接枝密度σ和帶電比f增大時，聚電解質刷半徑隨鏈長線性擴張.

圖4明顯地展示出3種刷型的聚電解質刷半徑與帶電比關系.接枝密度 σ從小到大變化，其中σ=0.1時，鹽濃度分別取為 ΦS=10-6和 ΦS=10-2兩種，后者是加鹽刷的情況.

圖4 聚電解質刷半徑與帶電比關系圖（雙對數(shù)坐標）鏈長N=100，除特別標明，鹽濃度為ΦS=10-6

如圖4所示，聚電解質刷半徑隨帶電比的增加不斷增大，ln R與ln f呈現(xiàn)出很強的線性關系.在同等條件下，高帶電比的聚電解質刷半徑更大，也更容易向外延伸，這主要是因為大的帶電比例會使聚電解質鏈之間的庫倫相互作用增加，靜電排斥使接枝鏈更加伸展.

不同條件下，ln R與 ln f線性關系比例系數(shù)大小不同.

接枝密度較低，σ=0.01，0.05，0.1時，R-f2/3，且刷半徑隨接枝密度σ增大而增大，正是Pincus刷的標度關系.

接枝密度較高，σ=0.5，0.8時，R-f1/2，但刷半徑隨接枝密度σ增大基本保持不變，這正是滲透壓刷最顯著的特點.

接枝密度相同 σ=0.1時，高鹽濃度 ΦS=10-2時，R-f2/5，與加鹽刷的標度關系一致.與低鹽濃度ΦS=10-6相比，加鹽使刷半徑收縮.

由以上討論可知，標度理論給出的聚電解質刷的3種類型是相當準確的，但3種刷沒有明顯的界限，通過調整接枝密度、外加鹽濃度等參數(shù)，可以實現(xiàn)動態(tài)轉變.

上面我們研究了刷的整體性質，下面我們研究聚電解質刷的細致結構.首先研究刷的單體密度分布，如圖5所示.

圖5左下角，直角坐標密度分布圖可以劃分為3個區(qū)域：內核層（r＜r2），線度很小，單體密度最大且基本保持不變直至次外層；中間層（r2＜r＜r1），單體密度較大且以指數(shù)關系減??；外層，（r＞r1），溶液中的離子與聚電解質作用開始顯現(xiàn)，單體密度急劇降低，直至與溶液接觸，即到達聚電解質刷的末端.當溶液中鹽離子濃度升高時，這3個區(qū)域受

到壓縮，內核層、中間層線度更小，單體密度更高，在徑向變化更劇烈.這是由于鹽濃度升高時，球形刷半徑減小，單體分布集中在更小的空間內.

圖5 聚電解質刷單體密度分布（雙對數(shù)坐標）左下角小圖為直角坐標下分布圖象鏈長N=100，接枝密度σ=0.1，帶電比f=0.5

在圖5雙對數(shù)坐標圖中更容易顯示出這3個區(qū)域及其變化趨勢，這都與標度理論的分析結果相類似.

圖5中，由SCFT的數(shù)值結果可知：

內核層聚電解質單體密度保持定值；外層聚電解質單體密度下降，鹽濃度較低（ΦS=10-5）時，如圖上斜率線-4/3表示φp-r-4/3.當鹽濃度升高時，由于聚電解質刷受到壓縮，刷末端單體密度迅速下降.高鹽濃度（ΦS=10-1）的溶液中聚電解質刷末端單體密度下降的斜率小于-4/3，近似為-1.標度關系預測低鹽濃度時 φp-r-2，即聚電解質刷均勻伸展，這與SCFT計算結果不一致.這是由于刷內層帶電聚電解質鏈間庫倫相互作用被反離子屏蔽，接枝鏈未完全伸展.

綜合以上的比較討論可以看出，標度理論對3種類型聚電解質刷的行為趨勢做出了相當準確的預測，但這對于重疊區(qū)的行為描述是不夠的.基于SCFT的數(shù)值計算，可以給出更多細節(jié)，更重要的是，SCFT計算可以展現(xiàn)出整個變化過程，打破3種類型的限制，為研究聚電解質刷行為的本質提供參考.

聚電解質刷與中性聚合物刷的不同在于，前者是帶電的，帶電的接枝鏈之間的庫倫相互作用使聚電解質刷的行為更為復雜.聚電解質刷體系里的帶電離子除聚合物離子外，還有鹽溶液的正負小離子，這些都對聚電解質刷的行為有著重要的影響，因而對于帶電的聚電解質刷，只研究單體密度分布是不夠的.

圖6（見下頁）給出了帶電離子和總電量在不同鹽濃度下的密度分布圖.本文研究的聚合物離子帶負電，因而反離子帶正電.在圖6（a）和（b）中，反離子的密度分布和聚離子基本重疊，除了刷內層的部分區(qū)域.這表明反離子主要被吸引在刷內區(qū)域，局域電中性條件是滿足的.而負離子在刷內的密度很小，相比之下是可以忽略的.在這種情況下，帶電聚電解質單體之間的庫倫排斥作用就基本上被反離子屏蔽了.所以聚電解質刷的結構主要由反離子的滲透壓和聚合鏈的彈性決定.刷內反離子的濃度取決于接枝密度σ和帶電比 f，所以圖6（a）和（b）中雖然鹽濃度相差一個數(shù)量級，但聚電解質刷結構幾乎相同.

在圖6（c）和（d）中，顯著數(shù)量的負離子和反離子分布在刷內層區(qū)域.當鹽濃度低時，如圖6（a）和（b）中右下角局部放大圖所示，反離子與負離子的徑向密度分布有一個很長的尾跡.而在圖6（c）和（d）中，離子濃度在刷末端降至溶液鹽濃度.而分布在刷內層區(qū)域的反離子與負離子濃度分別大于和小于溶液鹽濃度，這是電中性約束條件的結果.

在圖6（d）中，刷內刷外所有離子的濃度處在同一數(shù)量級.在這種情況下，聚電解質刷接枝鏈伸展的驅動力主要來自于所有自由離子的滲透壓.圖6（d）和圖6（c）對比，鹽濃度的升高，滲透壓降低，因而使聚電解質刷半徑收縮.

在圖6的4幅圖中，總帶電荷量的分布具有相同的特點.分布曲線從峰值降到谷底，然后緩慢上升，最后穩(wěn)定在0的水平線.峰值代表這一區(qū)域沒有單體，被反離子占據(jù).最小值在接枝表面附近，其絕對值隨著鹽濃度的增加而減小.這些現(xiàn)象與柱狀聚電解質刷的SCFT結果類似［18］.與之不同的是，球狀聚電解質刷的總帶電量分布曲線在穩(wěn)定之前還要經(jīng)歷小的峰值，如圖6（a）和（b）中右下角局部放大圖所示，這與平面聚電解質刷的SCFT結果類似［13］.而柱狀的曲線從谷底緩慢上升至0的水平線，沒有觀測到峰值.

圖6 離子密度分布圖右下角（a2）、（b2）分別為（a）、（b）的局部放大補充圖

圖7表示聚電解質刷自由端密度分布圖.從圖中可以看出，當鹽濃度很高的時候，自由端分布符合高斯型.隨著鹽濃度的降低，自由端密度最大值開始向鏈邊緣移動，即鏈進一步伸展.當鹽濃度進一步下降時，自由端分布逐漸變寬且變對稱.

圖7 球形聚電解質刷自由端密度分布鏈長N=100，接技密度σ=0.1，帶電比f=0.5

另一個重要的現(xiàn)象是，鹽濃度下降時，對于自由端，聚電解質刷內層逐漸出現(xiàn)明顯的排空區(qū)域.這是由于隨著鹽濃度降低，Debye屏蔽長度開始增加，單體之間的庫倫相互作用更加強，鏈更容易向外伸展.因而自由端出現(xiàn)在內層的概率變小.

4 結論

本文通過對球形聚電解質刷自洽場方程在實空間的數(shù)值求解，探究了體系鹽濃度、鏈接枝密度及帶電比對聚電解質刷相行為的影響，并和標度理論進行了對照.標度理論對3種類型聚電解質刷的行為趨勢做出了相當準確的預測，而 SCFT計算結果可以給出更多細節(jié)信息，尤其是各相過渡區(qū)的刷行為.

溶液未加鹽的情況下，接枝密度較低時，刷半徑隨著接枝密度增大而增加；高接枝密度下，刷半徑與接枝密度不再有相關關系.與中性刷相比，庫倫相互作用使得接枝密度增加對擴大刷半徑的貢獻更加明顯.在同等條件下，帶電比高的聚電解質刷半徑更大一些.

增加體系鹽濃度會使聚電解質刷半徑收縮，末端分布曲線的最大值點向球內層方向移動，同時自

由端密度分布的對稱性也遭到破壞.

鹽濃度較低時，聚電解質刷結構主要由反離子的滲透壓和聚合鏈的彈性決定.反之，接枝鏈伸展的力來自所有自由離子的滲透壓.

由于影響聚電解質體系穩(wěn)定性的因素非常多，本文在討論的過程中也做了一些取舍.在處理靜電相互作用時利用了泊松-玻爾茲曼方程，導致對于體系中小離子的描述不足.

致謝：

本文得到了北京師范大學物理學系嚴大東教授的悉心指導，其課題組的瞿立建博士亦給予很大幫助，謹致謝意.

本文的研究工作受到了北京師范大學“國家級大學生創(chuàng)新創(chuàng)業(yè)訓練計劃”（105602）的資助.

［1］ Kent M S.A quantitative study of tethered chains in various solution conditions using Langmuir diblock copolymer monolayers［J］.Macromol Rapid Commun，2000，21（6）：243-270.

［2］ Cosgrove T，Heath T G，Ryan K，et al.Neutron scattering from adsorbed polymer layers［J］.Macromolecules，1987，20（11）：2879-2882.

［3］ Parsonage E，Ryan T G，Crowley T L.Adsorption of poly（2-vinylpyridine）-poly（styrene）block copolymers from toluene solutions［J］.Macromolecules，1991，24（8）：1987-1995.

［4］ Tirrell M，Patel S，Hadziioannou G.Polymeric amphiphiles at solid-fluid interfaces：Forces between layers of adsorbed block copolymers［J］.Proc Natl Acad Sci USA，1987，84（14）：4725-4728.

［5］ Zhao B，Brittain W J.Polymer brushes：surface-immobilized macromolecules［J］.Prog Polym Sci，2000，25（5）：677-710.

［6］ Park M，Harrison C，Chaikin P M，et al.Block copolymer lithography：periodic arrays of similar to 1011holes in 1 square centimeter［J］.Science，1997，276（5317）：1401-1404.

［7］ Morkved T L，Wiltzius P，Jaeger H M，et al.Mesoscopic self-assembly of gold islands and diblock-copolymer films［J］.Appl Phys Lett，1994，64（4）：442-444.

［8］ Zhao B，Zhu L.Mixed Polymer Brush-Grafted Particles：A New Class of Environmentally Responsive Nanostructured Materials［J］.Macromolecules，2009，42（42）：9369-9383.

［9］ Luzinov I，Minko S，Tsukruk V V.Designing Patterned Surfaces by Grafting Y-Shaped Copolymers［J］.Prog Polym Sci，1996，29（4）：2667-2676.

［10］ Marko J F，Witten T A.Phase separation in a grafted polymer layer［J］.Phys Rev Lett，1991，66（34）：1541-1692.

［11］ Marko J F，Witten T A.Correlations in grafted polymer layers［J］.Macromolecules 1992，25（1）：296-307.

［12］ Lai P Y.Binary mixture of grafted polymer chains：A Monte Carlo simulation［J］.J Chem Phys，1994，100（4）：3351-3357.

［13］ Seki H，Suzuki Y Y，Orland H.Self-Consistent Field Study of Polyelectrolyte Brushes［J］.J Phys Soc Jap，2007，76（10）：104601.

［14］ Shi A C，Noolandi J.Theory of inhomogeneous weakly charged polyelectrolytes［J］.Macromol Theory Simul，1999，8（3）：214-229.

［15］ Wang Q.Charge Inversion by Flexible Polyelectrolytes on Flat Surfaces from Self-Consistent Field Calculations［J］.Macromolecules，2005，38（21）：8911-8922.

［16］ Press W H，Teukolsky S A，Vetterling W T，et al.P.Numerical Recipes［M］.3rd Ed.New York：Cambridge University Press，2002.

［17］ Borisov O V，Zhulina E B.Effects of ionic strength and charge annealing in star-branched polyelectrolytes［J］.Eur Phys J B，1998，4（2）：205-217.

［18］ Qu L J，Jin X G，Liao Q.Numerical Self-Consistent Field Theory of Cylindrical Polyelectrolyte Brushes［J］.Macromol Theory Simul，2009，18（3）：162-170.

Numerical simulation of spherical polyelectrolyte brushes by self-consistent filed theory

SU Zhi-cheng，HU Wei
（Department of physics，Beijing Normal University，Beijing 100875，China）

A spherical polyelectrolyte brush is studied by self-consistent field theory（SCFT）and the results are compared with predictions from scaling theory.It is shown that the SCFT results give the general trends as well as insight into the crossover regions between different regimes.The density profiles of the polyions and small ions indicate that the systems are locally electroneutral.The salted brush bears characteristics similar to those of a neutral brush.Counter-intuitively，the chains are not uniformly stretched in the osmotic regime.The free-end monomers shift to the outer region and an exclusion zone appears and grows with decreasing of salt concentration.

SCFT；polyelectrolyte；salt concentration；grafting density；morphology

O 552.4+2

A

1000-0712（2016）11-0045-07

2013-03-13；

2016-03-09

蘇志成（1989- ），男，甘肅平?jīng)鋈?，北京師范大學物理學系2009級本科生.