水溶液中顆粒輸送器的分子動力學模擬

付中玉，孫康，梁棟，徐震

摘 ?要： 由于許多實際納米顆粒的電荷和極性十分微弱甚至可以忽略不計，這阻礙了電場產生直接驅動納米顆粒?；谠谔荻入妶稣T導下的納米顆?？勺园l從電場強度高的位置迅速移動到電場強度低的位置這一原理，本課題使用不同斜率的垂直梯度電場驅動控制水溶液中納米顆粒，實現納米顆粒定向運動。根據分子動力學模擬結果得出當梯度電場斜率合適并且顆粒的大小與輸送顆粒管道直徑相匹配時，納米顆?？梢赃_到較好的輸送效果。該結果對于水溶液中分子馬達、分子水泵、分子篩以及分子開關的設計有著重要的參考價值。

關鍵詞： 水溶液;顆粒輸送器;梯度電場;分子動力學

中圖分類號： TP271+.3 ? ?文獻標識碼： A ? ?DOI：10.3969/j.issn.1003-6970.2019.12.007

本文著錄格式：付中玉，孫康，梁棟，等. 水溶液中顆粒輸送器的分子動力學模擬[J]. 軟件，2019，40（12）：2832

Molecular Dynamics Simulation of Particle Conveyor in Aqueous Solution

FU Zhong-yu， SUN Kang， LIANG Dong， XU Zhen*

（Shanghai University of Engineering Science， Shanghai 201620）

【Abstract】： Since the charge and polarity of many actual nanoparticles are very weak or even negligible， this hinders the electric field from directly driving the nanoparticles. Based on the principle that the nanoparticles induced by the gradient electric field can spontaneously move from the position where the electric field strength is high to the position where the electric field strength is low， this topic uses the vertical gradient electric field with different slopes to drive and control the nanoparticles in the aqueous solution to realize the directional movement of the nanoparticles. According to the results of molecular dynamics simulation， when the slope of the gradient electric field is suitable and the particle size matches the diameter of the transporting particle tube， the nanoparticle can achieve better transport effect. The results have important reference value for the design of molecular motors， molecular water pumps， molecular sieves and molecular switches in aqueous solution.

【Key words】： Aqueous solution; Particle conveyor; Gradient electric field; Molecular dynamics

0 ?引言

納米顆粒的尺寸通常在1～100 nm之間，并被人們認為是宏觀物質與微觀原子、分子之間的紐帶[1]。納米尺度下的納米顆粒表現出很多新穎甚至奇異的物理、化學和生物學性質，因此納米顆粒在材料改性和生物醫療等方面具有重要的應用前景。

納米顆粒可以通過血管、神經突觸和淋巴血管傳輸，而且還能夠有選擇性地積累在不同細胞或者一定的細胞結構中，對進入細胞的分子等細微顆粒進行有效的過濾篩選。納米顆粒的這些特性奠定了其在藥物輸運[2-4]、靶向治療[5-7]及分子篩選[8]等方面的應用潛能。關于納米顆粒的研究近些年來已經取得了一些成果，例如谷紅梅等人[9]研究了納米顆粒與生物膜之間的相互作用對納米顆粒跨膜輸運的影響。該研究表明與生物膜表面存在弱排斥相互作用的納米顆粒更容易進入細胞膜，但從細胞內穿出卻相當困難。郝亮等人[10]利用外加電場驅動管狀顆粒，發現增大電場強度可以增大納米顆粒的跨膜流量并且帶正電的納米顆粒與水溶液的相互作用比較弱，所以優先于帶負電納米顆粒進入納米通道。然而在實際情況下，許多納米顆粒甚至一些生物大分子的電荷和極性十分微弱甚至可以忽略不計。如何驅動控制這些具有極弱電性或不帶電的納米顆粒一直是人們急于解決的問題。由于碳納米管（CNT）的獨特特性[11-18]，Xu等人[19]通過梯度電場（非均勻電場）場強操縱水溶液中的納米顆粒。該研究表明非均勻電場的場強可以改變水溶液能量發生變化并且納米顆粒可以自發從電場強度高的區域到達電場強度低的區域。然而，此研究并未對電場斜率對顆粒運動的影響進行說明。

基于以上的調研，本研究使用不同斜率的梯度電場驅動納米顆粒（富勒烯小球，C60），并利用分子動力學對所建立的系統進行模擬，觀察梯度電場的斜率對CNT內水流大小、納米顆粒運動的影響。

1 ?模型和數值方法

本研究中所用模型如圖1，圖2所示。在兩個石墨片之間的垂直方向上分別嵌入長度為10 nm的（8，8），（10，10），（12，12），（15，15）的不同管徑的單壁碳納米管（SWCNT），然后將系統放置在尺寸為5.2×5.2×20.4 nm3的立方體盒子中。SWCNT的軸線沿z軸，SWCNT的中心位于整個模擬箱的中心。另外，將密度為998kg/m3的水分子填充在盒子中，但不填充在兩個石墨片之間，并將不同斜率的非均勻電場施加到該系統。非均勻電場函數表達式為：，電場方向沿x軸， 為電場函數的梯度系數，其中z是本研究中電荷原子的Z坐標。這種理想化的非均勻電場物理模型可以由平板電容器構

圖1 ?仿真系統不含納米顆粒的模型圖

Fig.1 ?Model of the simulation system without nanoparticles

圖2 ?仿真系統加入納米顆粒的模型圖

Fig.2 ?Model diagram of the simulation

system adding nanoparticles

成。如果電容器的兩個極板在一個開口端彼此靠近而在另一個開口端彼此遠離，則在改進的極板電容器的內部將形成具有一定梯度的非均勻電場[20]。

此次模擬借助GROMACS 5.1.5軟件[21]，采用的是分子動力學模擬的方法對這個問題進行研究。模擬所用的系綜是NVT，即保持原子數、體積和溫度不變，其中溫度由速度標定法（v-rescale），控制在300K左右。在模擬中所用到的Lennard-Jones（LJ）作用參數為：cc=0.34 nm，cc=0.3664 KJmol–1。碳原子之間的平衡鍵長為0.1418 nm，平衡鍵角為120°，利用harmonic勢能函數進行約束，約束能量系數分別是393960 kJmol–1和527 kJmol–1·rad–2.另外，一個較弱的正則二面角勢能被用來進一步約束碳原子鍵的振動。碳納米管和水分子之間的相互作用采用碳原子和氧原子之間的LJ勢能來描述，其LJ參數為co=0.3275 nm，co=0.4772 KJmol–1。另外，電場對勢能的貢獻項用來描述，其中qi和ri分別表示帶電原子i的電荷量和位置坐標，而E是外加電場。在模擬計算中，LJ勢能采用截斷算法，其截斷半徑為1.2 nm，并將時間步長為1fs，數據采集頻率為每隔0.5 ps記錄一幀。相對于其它算法[22-27]，PME（Particle-Mesh-Ewald）算法在靜電相互作用中應用較廣[28-34]，其中短程作用的截斷半徑為1nm。該模型所用參數在之前的工作中已被廣泛應用[19，32，35-40]。

此次模擬使用 = 0.02，0.05，0.08，0.10和0.12不同斜率的線性梯度電場作用在所建立的系統上，并且對模擬的時間進行設定（不帶顆粒的模型模擬10 ns，帶顆粒的模型模擬20 ns）。在最終得到模型結果進行處理時，本研究借助于UltraEdit編輯軟件對模型處理的軌跡文件和結構文件進行修改，運用1.9.3版本的VMD軟件，對MD軌跡文件進行分析。

2 ?結果與討論

系統在施加斜率為0.08的梯度電場后，（10，10）型CNT管內的水分子變化如圖3所示?？梢钥吹疆斒┘犹荻入妶龊?，系統CNT管內出現一些游離的水分子（如圖3，2 ps時刻的狀態），隨著模擬時間的增長，游離的水分子逐漸聚集到管的一端（如圖3，28 ps時刻的狀態），然后聚集之后的水分子沿著管的一端移動到管的另外一端（如圖3，80-1744 ps的運動過程）。該現象與Xu等人[19]的研究結果相一致。該現象主要是因為，開始施加梯度電場時，作用于系統Z軸方向梯度電場場強小，系統內部原有的水結構體系只是被電場破壞;隨著梯度電場場強沿著Z軸正向增大，系統內部單個游離的水分子增多，甚至出現完全分離現象。相關研究[15，16，41]表明水分子偶極在外加電場的作用下傾向于和電場方向相一致。因此在模擬過程中，當梯度電場沿Z軸增大時，水分子偶極指向被誘導趨向于電場方向相一致，并且在外加電場影響下，管內的水分子偶極分布呈現不對稱性，形成沿Z軸方向的單向水流（如圖3的單向水流）。該結果與Su等人[15]研究結果相吻合，受外加電場的影響單向水流的方向與水分子偶極的方向保持一致。

圖3 ?系統施加斜率為0.08的梯度電場后，

（10，10）型CNT內部水分子變化

Fig.3 ?After the system applies a gradient electric field with a slope of 0.08， the water molecules in the

（10，10） type CNTs change

圖4是系統在模擬開始時到系統達到平衡時的一段時間內，對管內的水分子進行統計?？梢钥闯?，受梯度電場斜率的影響，CNT管內都會產生不同大小的單向水流。隨著梯度電場斜率的增大，流經不同管徑內的水分子數基本都呈現先增大后減小的趨勢。因此，當設計水中納米顆粒質量輸運器時，在

圖4 ?在不同斜率梯度電場作用下，

四種不同管徑內水流的變化趨勢

Fig.4 ?Trends of water flow in four different pipe

diameters under different slope gradient electric fields

選好管徑的同時，還應該考慮到斜率對水流的影響，此因素對于顆粒輸運的作用不容忽視。

考慮到C60的直徑與碳納米管的直徑以及碳-碳之間相互作用的關系，只選擇在（10，10）、（12，12）、（15，15）的模擬系統中加入C60顆粒。對于（10，10）型模擬系統，顆粒基本可以在3 ns之前到達管的另一端，并且最后可以在相應的位置處于穩定狀態，實現利用梯度電場輸送納米顆粒的效果。對于（12，12）型模擬系統，顆粒在斜率為0.02時不能實現輸送的目的，而對于斜率為0.05、0.08、0.10、0.12模擬體系可以實現輸送顆粒目的，并且該過程基本可以在3 ns之前到達管的另一端，并且最后可以在相應的位置處于穩定狀態。相對于（10，10）和（12，12）的模擬體系，（15，15）的模擬體系在輸送顆粒過程中，斜率為0.02梯度電場操縱顆粒輸送時表現出較大波動，顆粒抵達指定位置后出現無規律的往復運動。對于斜率為0.05、0.08、0.10、0.12的模擬體系同樣可以達到輸送顆粒的目的。通過圖5、圖6、

圖5 ?在不同斜率電場作用下，（10，10）型碳納米管內

顆粒質心位移隨時間的變化曲線

Fig.5 ?The variation of the centroid displacement of particles in （10，10） type carbon nanotubes with time

under different slope electric fields

圖6 ?在不同斜率電場作用下，（12，12）型碳納米

管內顆粒質心位移隨時間的變化曲線

Fig.6 ?The variation of the centroid displacement of the （12，12） type carbon nanotubes with time under

different slope electric fields

圖7可以看出，斜率為0.05、0.08、0.10的模擬體系，在輸送顆粒過程中，顆粒能夠較快通過管徑抵達指定位置，而且可以穩定在指定位置，輸送效果相對于斜率為0.02和0.12的模擬體系比較好。

圖7 ?在不同斜率電場作用下，（15，15）型碳納米

管內顆粒質心位移隨時間的變化曲線

Fig.7 ?The variation of centroid displacement of

particles in （15，15） type carbon nanotubes with

time under different slope electric fields

圖8是不同管徑內納米顆粒運動的平均變化趨勢。（10，10）型碳納米管內顆粒的運動速度最大為8.94 nm/ns，（12，12）型碳納米管內顆粒的運動速度為7.34 nm/ns，（15，15）型碳納米管內顆粒的運動速度為5.18 nm/ns。可以看出，（10，10）型碳納米管內的納米顆粒運動速度比較快。因為對于（10，10）型模擬系統來講，納米顆粒（C60）的直徑與（10，10）型CNT的管徑相差不多，水分子不能從兩者之間的縫隙穿過，而對于（12，12）和（15，15）型模擬系統，C60與（12，12）和（15，15）型CNT管徑之間的縫隙較大，盡管兩者的水流比（10，10）型模擬系統的水流大，可是卻不能保證所有水分子的偶極全部作用于納米顆粒。因此，當設計水中納米顆粒質量輸

圖8 ?表示在不同斜率電場作用下，

不同管徑內納米顆粒平均速度的變化趨勢

Fig. 8 ?It shows the change trend of the average

velocity of nanoparticles in different diameters

under different slope electric fields

運器時，需要根據顆粒大小選取合適的輸送管道，從而實現最佳的輸送效果。

3 ?結論

本文研究了不同斜率梯度電場作用下納米顆粒輸送效果。根據分子動力學模擬的結果發現，納米顆粒自身的大小以及梯度電場的斜率都與顆粒的輸送效果相關。當梯度電場斜率合適并且顆粒的大小與輸送顆粒管道直徑相匹配時（輸送管道中的水分子偶極可全部作用于納米顆粒），納米顆粒會在短時間內到達指定位置，并處于穩定狀態。此研究成果對于存儲納米顆粒、輸送生物分子、制作納米流體器件、海水淡化等有著重要的參考價值。

參考文獻

[1]白春禮. 納米科技及其發展前景[J]. 世界科學， 2001， 52（4）： 37-37.

[2]于海洋. 順鉑高分子納米抗腫瘤藥物輸送體系的設計與應用[J]. 2015.

[3]李艷， 羅成， 周婕. 納米沉淀法及其在納米藥物載體領域的應用研究進展[J]. 中國生物醫學工程學報， 2017， （4）.

[4]康垚， 王素真， 樊江莉， 等. 無機納米藥物載體在腫瘤診療中的研究進展[J]. 化工學報， 2018.

[5]莊四良. 載順鉑mPEG-PGA納米微球靶向偶聯物的制備研究[D]. 2015.

[6]楊雪垠， 金武龍. 納米藥物載體在口腔癌靶向治療中的研究進展[J]. 中國醫藥導報， 2016， （36）.

[7]吳江. 基于介孔多功能納米材料靶向乳腺癌的影像診斷和治療一體化研究[D]. 2015.

[8]唐亞林， 向俊鋒， 周秋菊， 等. 基于生物靶分子識別的天然抗腫瘤分子篩選鑒定一體化技術[C]. 中國分析測試協會科學技術獎發展回顧， 2015.

[9]谷紅梅， 梁清， 陳慶虎. 相互作用對納米顆?？缟锬ぽ斶\的影響[J]. 生物物理學報， 2011， （5）.

[10]郝亮， 蘇加葉， 郭洪霞. 水溶液中管狀帶電納米顆粒跨膜輸運的分子動力學模擬[J]. 高分子學報， 2016， （12）.

[11]Yeh I C， Hummer G. Nucleic acid transport through carbon nanotube membranes[J]. Proceedings of the National Academy of Sciences， 2004， 101（33）： 12177-12182.

[12]Li J， Gong X， Lu H， et al. Electrostatic gating of a nanometer water channel[J]. Proceedings of the National Academy of Sciences of the United States of America， 2007， 104（10）： 3687-3692.

[13]Lu H， Li J， Gong X， et al. Water permeation and wavelike density distributions inside narrow nanochannels[J]. Physical Review B， 2008， 77（17）.

[14]Joseph S， Aluru N R. Pumping of confined water in carbon nanotubes by rotation-translation coupling[J]. Physical Review Letters， 2008， 101（6）： 064502.

[15]Su J， Guo H. Control of Unidirectional Transport of Single- File Water Molecules through Carbon Nanotubes in an Electric Field[J]. ACS Nano， 2011， 5（1）： 351-359.

[16]Wang Y， Zhao Y J， Huang J P. Giant Pumping of Single-File Water Molecules in a Carbon Nanotube[J]. The Journal of Physical Chemistry B， 2011， 115（45）： 13275-13279.

[17]姜麗麗， 李傳通， 于海濤， 等. 石墨烯/碳納米管復合材料的制備方法及應用進展[J]. 化學與生物工程， 2017， （2）.

[18]王浩然. 碳納米管及其應用[J]. 山東化工， 2017， （01）： 68- 70.

[19]Xu Z， Wang C， Sheng N， et al. Manipulation of a neutral and nonpolar nanoparticle in water using a nonuniform electric field[J]. The Journal of Chemical Physics， 2016， 144014302.

[20]Rinne K F， Gekle S， Bonthuis D J， et al. Nanoscale pumping of water by AC electric fields[J]. Nano Letters， 2012， 12（4）： 1780-1783.

[21]Hess B， Kutzner C， Spoel D V D， et al. GROMACS 4： Algorithms for highly efficient， load-balanced， and scalable molecular simulation[J]. Joural of Chemical Theory and Computation， 2008， 4（3）： 435-447.

[22]李斌， 李蓉， 周蕾. 分布式K-means聚類算法研究與實現[J]. 軟件， 2018， （1）： 35-38.

[23]吳江， 孫劍偉. 一種基于云模型的數據預測算法[J]. 軟件， 2015， 36（12）： 212-215.

[24]陳向華， 劉可昂. 基于FP-Tree的最大頻繁項目集挖掘算法[J]. 軟件， 2015， 36（12）： 98-102.

[25]曹盟盟， 姚文斌. 基于改進粒子群算法的虛擬機放置算法[J]. 軟件， 2015， 36（12）： 89-92.

[26]唐漢香， 劉宇暢， 唐一超. 基于局部搜索的Memetic差分演化算法[J]. 軟件， 2015， 36（10）： 13-16.

[27]趙乃剛. 一種新的基于模擬退火的粒子群算法[J]. 軟件， （7）： 1-4.

[28]Darden T， York D， Pedersen L. Particle mesh ewald： an N?log（N） method for ewald sums in large systems[J]. The Journal of Chemical Physics， 1993， 98（12）： 10089-10092.

[29]Rahman M M， Chowdhury M M， Alam M K. Rotating-electric-field-induced carbon-nanotube-based nanomotor in water： a molecular dynamics study[J]. Small， 2017， 13（19）： 1603978.

[30]Xu Z， Wang C， Sheng N， et al. Manipulation of a neutral and nonpolar nanoparticle in water using a nonuniform electric field[J]. The Journal of Chemical Physics， 2016， 144（1）： 014302.

[31]Winarto， Takaiwa D， Yamamoto E， et al. Separation of water-ethanol solutions with carbon nanotubes and electric fields[J]. Physical Chemistry Chemical Physics， 2016， 18（48）： 33310-33319.

[32]Liu L， Patey G N. Simulated conduction rates of water through a （6，6） carbon nanotube strongly depend on bulk properties of the model employed[J]. The Journal of Chemical Physics， 2016， 144（18）： 184502.

[33]Gong X， Li J， Lu H， et al. A charge-driven molecular water pump[J]. Nat Nanotechnol， 2007， 2（11）： 709-712.

[34]Rahman M M， Chowdhury M M， Alam M K. Dynamics of fullerene self-insertion into carbon nanotubes in water[C]. International Conference on Electrical & Computer Engineering. IEEE， 2015.

[35]Ostler D， Kannam S K， Daivis P J， et al. Electropumping of water in functionalised carbon nanotubes using rotating electric fields[J]. The Journal of Physical Chemistry C， 2017， 12128158-28165.

[36]Sisan， Thomas B， Lichter S. Solitons transport water through narrow carbon nanotubes[J]. Physical Review Letters 2014， 112（4）： 044501.

[37]Liu L， Patey G N. Simulations of water transport through carbon nanotubes： How different water models influence the conduction rate[J]. The Journal of Chemical Physics 2014， 141（18）： 18C518.

[38]Secchi E， Marbach S， Niguès， et al. Massive radius- dependent flow slippage in carbon nanotubes[J]. Nature， 2016， 537（7619）： 210-213.

[39]Da?i? M， Stankovi? I， Gkagkas K. Molecular dynamics investigation of the influence of the shape of the cation on the structure and lubrication properties of ionic liquids[J]. Physical Chemistry Chemical Physics， 2019， 21.

[40]Sam A， Kannam S K， Hartkamp R， et al. Water flow in carbon nanotubes： The effect of tube flexibility and thermostat[J]. The Journal of Chemical Physics， 2017， 146234701.

[41]Joseph S， Aluru N R. Pumping of Confined Water in Carbon Nanotubes by Rotation-Translation Coupling[J]. Physical Review Letters， 2008， 101（6）： 064502.