帶電系統(tǒng)中的靜電相互作用-建模、分析與計算.pdf

1、靜電相互作用在很多復雜的帶電系統(tǒng)中起著非常關鍵的作用。如何對這些帶電系統(tǒng)中的靜電相互作用進行有效地建模以及計算是個非常具有挑戰(zhàn)性的問題。其中涉及的困難主要源于系統(tǒng)復雜的幾何構造，靜電相互作用的多體效應以及多尺度現(xiàn)象。經典的泊松-波爾茲曼理論是研究帶電系統(tǒng)中靜電相互作用最為成功的方法之一。它已經被廣泛地被應用于各個領域，如生物分子溶解模型，膠體理論以及電化學等等。但是，泊松-波爾茲曼理論忽略了帶電系統(tǒng)中微粒的體積大小以及微粒之間的相互作用

2、，從而使它的應用范圍有所局限。本篇博士論文探討如何修正泊松-波爾茲曼理論使之能夠考慮系統(tǒng)中微粒的體積大小效應，并把修正后的泊松-波爾茲曼理論應用于計算電解質溶液中帶電表面附近的離子濃度分布。蒙特卡洛模擬結果說明修正的泊松-波爾茲曼理論能夠非常有效地描述系統(tǒng)中離子的體積大小效應。最后，我們還考慮了隱式溶劑變分模型中的靜電相互作用對介電界面運動產生的影響。
　　 在第二章中，我們研究如何修正泊松-波爾茲曼理論，并用它來研究電解質溶液

3、中離子的體積大小效應。首先，我們考慮極小化不含離子大小效應的靜電自由能泛函。由于系統(tǒng)中離子個數(shù)守恒以及整體系統(tǒng)呈電中性，所以系統(tǒng)中每種離子濃度滿足一定的約束條件。應用拉格朗日乘子法把該約束優(yōu)化問題轉換為無約束優(yōu)化問題。從而，我們得到了經典的泊松-波爾茲曼方程的解。其次，我們考慮包含離子大小效應的泊松-波爾茲曼理論?；诶硐霘怏w的格子模型，我們在自由能泛函中熵的部分考慮溶劑分子以及離子的不同體積大小效應。極小化修正后的自由能泛函，我們可以

4、得到一組決定離子濃度分布的非線性方程組。但是，當體系內離子的體積大小以及溶劑分子的體積大小不相同時，我們不能從中得出類似波爾茲曼分布的解析解。由于離子濃度仍需滿足一定的約束，所以我們提出了一種增廣拉格朗日算法來解決相應的約束優(yōu)化問題。數(shù)值例子說明該模型以及相應的算法可以高效地計算體系內的靜電勢和離子濃度分布。在數(shù)值實驗中，我們研究了帶電表面面電荷強度和離子體積大小對體系內離子濃度分布的影響。我們觀察到強帶電表面附近的反離子出現(xiàn)分層現(xiàn)象，

5、并發(fā)現(xiàn)化合價與離子體積比值是決定強帶電表面附近的反離子分層結構的重要因素。本章的工作是和王忠明教授、李波教授共同合作完成。
　　 為了進一步驗證上一章中提出的模型、算法以及結論的正確性，我們在第三章中考慮用蒙特卡洛方法來模擬帶電系統(tǒng)中的離子濃度分布。我們考慮和上一章一樣的帶電體系。我們采用非限制的原始模型(unrestricted primitivemodel)來描述該電解質系統(tǒng)，并用正則系綜蒙特卡洛模擬結合Metropolis

6、準則來研究系統(tǒng)中的離子分布情況。系統(tǒng)中移動的離子用帶電的硬球表示，溶劑的作用通過它的介電系數(shù)來體現(xiàn)。整個蒙特卡洛模擬過程由三個階段組成。第一個階段是增加體系的溫度，從而加速離子運動。第二個階段是在室溫下移動離子使體系達到平衡狀態(tài)。第三個階段是統(tǒng)計離子處于平衡狀態(tài)時的密度分布情況。經過大量的數(shù)值實驗我們發(fā)現(xiàn)離子體積的不同大小效應對電偶層中的離子濃度分布起決定性的作用。在較弱面電荷密度的帶電表面附近，化合價高的反離子在吸附過程中占優(yōu)勢；在較

7、強面電荷密度的帶電表面附近，蒙特卡洛模擬和平均場理論結果均顯示反離子緊密地吸附在帶電表面周圍，并出現(xiàn)分層現(xiàn)象。而且，分層結構受化合價與離子體積比值影響，而不是由化合價唯一決定。具體地說，有著最大化合價與離子體積比值的反離子最先分布在帶電表面周圍，然后再是化合價與離子體積比值較小的反離子。在反離子吸附到強帶電表面的過程中，化合價與離子體積比值很好地體現(xiàn)了自由能泛函中靜電能部分和熵部分之間的相互競爭關系。蒙特卡洛模擬很好地驗證了上一章得出的

8、結論，即化合價與離子體積比值是分析電偶層中反離子濃度分布的關鍵指示量。本章的工作是和聞佳夷、徐振禮教授、李波教授共同合作完成。
　　 在第四章中，我們研究隱式生物分子溶解模型中靜電的相互作用以及區(qū)域的幾何形狀對介電界面(dielectric interface)運動產生的影響。在柱形坐標下，我們考慮界面沿著由表面積能和靜電能組成的自由能泛函的最速下降方向運動。首先，我們求出在柱形坐標下靜電能關于介電界面的形狀導數(shù)。那么，我們可以

9、得到有效的介電界面力的顯式表達式。可以發(fā)現(xiàn)靜電部分產生的力的作用方向永遠是從溶劑區(qū)域指向溶質區(qū)域。當界面外層區(qū)域的介電系數(shù)較大時，界面在運動過程的開始階段逐漸被光滑化，然后朝內部區(qū)域加速運動。當界面內部區(qū)域的介電系數(shù)較大時，表面積能產生的力和靜電能產生的力相互競爭并達到平衡，從而我們可以得到一個穩(wěn)態(tài)界面。線性穩(wěn)定性分析表明該穩(wěn)態(tài)界面只對波數(shù)超過一定閾值的擾動是穩(wěn)定的。我們在數(shù)值實驗中模擬了上述的兩種情況。數(shù)值實驗結果驗證了有效介電界面力