水分子通過(guò)形變碳納米管的運輸行為分析論文

　　水是地球上最常見(jiàn)的物質(zhì)之一，也是生物體的最重要的組成部分，同時(shí)水在生命活動(dòng)中起到重要的作用，在這些活動(dòng)中涉及到水通道蛋白。 由于碳納米管具有獨特的結構和奇特的性質(zhì)[1],故許多學(xué)者[2 -4]將其作為水通道的簡(jiǎn)化模型加以研究。

　　Hummer 等人[5]在 2001 年采用分子動(dòng)力學(xué)模擬的方法，發(fā)現水分子不僅可以在單壁碳納米管中運輸，而且會(huì )加快其運輸速度。 一部分研究人員發(fā)現如果把碳納米管口作部分扭曲后，能很大地加快水分子通過(guò)碳納米管的速率。 2007 年李敬源等人[6]將電荷放置于碳納米管之外，研究了外加電荷對管內水分子的影響，觀(guān)察到了電荷對于水分子運輸的開(kāi)關(guān)方面非常靈敏。 此后，大量的學(xué)者[7 -9]在這些重大成果啟發(fā)之下，對水在碳納米管中動(dòng)力學(xué)性質(zhì)作了許多研究。 如周毅[7]在碳納米管中水的行為的分子動(dòng)力學(xué)模擬中講述： 在納米尺度下的受限空間中，水分子所表現的行為與宏觀(guān)中截然不同。 何俊霞等人[8]通過(guò)改變形變位置探究納米通道中的窄結構對水輸運產(chǎn)生的開(kāi)關(guān)特性。

　　由于碳納米管生長(cháng)和合成實(shí)驗中[10],總有一些拉伸或壓縮的形變碳納米管存在，但形變碳納米管內水分子的運輸特性未見(jiàn)到相關(guān)報道。 故本文運用分子動(dòng)力學(xué)的方法研究碳納米管拉伸或壓縮時(shí)其形變程度對其內部水的運輸行為影響。 即選擇手性矢量（ 7,7） 和（ 8,8） 2 種碳納米管，拉伸或壓縮量 k 分別為： 0. 5,0. 6,0. 8,1. 0,1. 2,1. 5,1. 8 倍的碳納米管進(jìn)行模擬計算，進(jìn)而分析形變碳納米管中水分子的流量、偶極矩概率分布、軸向分布以及徑向分布。

　　1 研究模型與方法

　　本文采用分子動(dòng)力學(xué)軟件 NAMD 和 CHARMM力場(chǎng)文件。 研究模型如圖1 所示，它主要由 1 個(gè)形變的單壁碳納米管和 2 個(gè)石墨烯層構成，并且 2 個(gè)石墨烯成處于碳納米管的 2 個(gè)端口處，在石墨烯層外同時(shí)添加 2 個(gè) TIP3P 水體，其大小為 1. 5 nm ×1. 5 nm × 1. 5 nm,其中形變碳納米管由 NanotubeModeler 軟件生成，其原長(cháng)為 6. 05 nm.

　　模擬時(shí)，設定溫度為 300 K,并采用 Berenden方法進(jìn)行控溫，截斷半徑為 1. 2 nm,利用 Particle -Mesh - Ewald 的方法計算系統的靜電相互作用，網(wǎng)格劃分精度設置為 0. 1 nm,當超過(guò) 1 nm 時(shí)采用過(guò)渡函數修正計算，采用周期性邊界條件其大小為1. 5 nm × 1. 5 nm × 9. 2 nm,模擬步長(cháng) 2 fs,模擬時(shí)長(cháng)均為 10 ns. 為解決錯誤結構的影響本文執行10 000步的能量最小化處理，同時(shí)為了減少運算量和簡(jiǎn)化模擬體系，將碳納米管和石墨烯層進(jìn)行固定。 為使水分子通過(guò)碳納米管，本文沿 z 軸正方向施加 0. 5 v/nm 的電場(chǎng)。 所有的計算任務(wù)均是在高性能服務(wù)器上完成的。

　　2 結果與分析

　　2. 1 水分子通過(guò)碳納米管的流量

　　水分子通過(guò)形變碳納米管的流量如圖 2 所示，可以發(fā)現： 當碳納米管未發(fā)生形變時(shí)，即 k = 1 時(shí)，管徑大的流量大，該結論與王俊等人[11]一致； 當碳納米管發(fā)生較小的形變時(shí)，即 k =0. 8 或 1. 2 時(shí)，管徑大的水分子流量反而減小，但均大于未發(fā)生形變時(shí)的流量； 當形變程度 k≥1. 5 時(shí)，2 種管徑的水分子流量幾乎為零，這說(shuō)明（ 7,7） 和（ 8,8） 型碳納米管被拉伸到 1. 5 倍以上時(shí)，幾乎沒(méi)有水分子能通過(guò)碳納米管。 分析認為碳納米管的軸向形變程度對管內水分子的運輸有重要的影響作用，同時(shí)可以通過(guò)改變碳納米管的形變程度使其起到分子開(kāi)關(guān)的作用。

　　2. 2 水分子通過(guò)碳納米管的徑向分布

　　碳納米管中水分子的徑向分布函數如圖 3 所示，其中圖 3（ a） 、（ b） 分別表示（ 7,7） 和（ 8,8） 型碳納米管中水分子的徑向分布圖。 由圖 3（ a） 可以發(fā)現不同形變程度下，碳納米管內水分子的徑向分布函數均只出現 1 個(gè)峰值并位于 r = 0. 16 nm 處，說(shuō)明水分子在該處出現的概率最大，這表明水分子在管內是以水柱面的形式通過(guò)碳納米管； 同時(shí)還可以發(fā)現，當 k 越小時(shí)其峰越尖銳，這說(shuō)明水分子僅僅從柱面處通過(guò)，當 k 增加時(shí)其水柱內也有水分子通過(guò)。 由圖 3（ b） 同樣可以發(fā)現徑向分布函數僅存在1 個(gè)位于 r = 0. 23 nm 處的峰，同樣也說(shuō)明水分子是以水柱面的形式通過(guò)碳納米管。 同時(shí)可以發(fā)現管徑大的其水柱面的半徑大。

　　2. 3 偶極矩概率分布

　　為了認識水分子通過(guò)碳納米管時(shí)其取向分布，圖 4 給出了水分子偶極矩方向與 z 軸正方向的夾角的概率分布。 由圖4（ a） 可以發(fā)現對于（ 7,7） 型碳納米管中水分子的偶極矩方向的概率分布僅出現1 個(gè)的峰，說(shuō)明水分子是以較規則的方式通過(guò)碳納米管。 同時(shí)還可以發(fā)現，出現峰的位置與碳納米管的形變程度 k 有關(guān)，如當 k<1 k="" z="">1 時(shí)，其偶極矩方向的概率分布的峰位置隨著(zhù) k 值的增加而增加，這說(shuō)明當碳納米管拉伸越大時(shí)水分子的偶極矩取向越垂直于 z軸。 由于 k 為 1. 5 和 1. 8 時(shí)碳納米管中水分子的流量較小，從而采樣較小，所以水分子偶極矩方向分布沒(méi)有出現尖銳的峰值。 由圖 4（ b） 可知對于（ 8,8） 型碳納米管，當 k < 1. 5 時(shí)，其水分子偶極矩方向的概率分布均只出現 1 個(gè)峰值，并且峰的位置與碳納米管的形變程度無(wú)關(guān)，這說(shuō)明管徑大的其內部水分子的取向與管本身的形變程度無(wú)關(guān)。 同樣由于 k為 1. 5 和 1. 8 時(shí)，管內水分子的流量較小，所以水分子概率分布未出現尖銳的峰值。

　　3 總結

　　本文通過(guò)分子動(dòng)力學(xué)的方法研究了水分子通過(guò)形變碳納米管的運輸行為，即探究碳納米管在拉伸或壓縮時(shí)其形變程度 k 與管內水分子的流量、徑向分布函數和偶極矩概率分布的關(guān)系。 結果表明：碳納米管的軸向形變程度對管內水分子的運輸有重要的影響作用； 水分子在管內是以水柱面的形式通過(guò)碳納米管，當 k 越小時(shí)其僅能從水柱的表面通過(guò)，并且發(fā)現管徑大的其水柱面的半徑大； 對于（ 7,7） 型碳納米管中水分子的偶極矩方向的概率分布僅出現 1 個(gè)的峰，并且峰的位置與碳納米管的形變程度 k 有關(guān)； 對于（ 8,8） 型碳納米管其水分子偶極矩方向的概率分布也均只出現 1 個(gè)峰值，但峰的位置與碳納米管的形變程度無(wú)關(guān)。 分析認為碳納米管的形變程度對管內水分子的運輸有重要的影響作用。 本文的研究對認識生命的新陳代謝活動(dòng)、受限水的行為和納米運輸器件有著(zhù)重要的參考作用。

　　參考文獻：

　　[1]王普寶。 碳納米管的相關(guān)力學(xué)問(wèn)題研究[D]. 大連： 大連理工大學(xué)，2007.

　　[2]王禹。 水分子在碳納米管中的分子動(dòng)力學(xué)模擬[D]. 上海： 復旦大學(xué)，2012.

　　[3]曾明穎。 多壁碳納米管彎曲角度的研究[D]. 廈門(mén)： 廈門(mén)大學(xué)，2005.

　　[4]呂小彬。 受限于碳納米管中水的分子動(dòng)力學(xué)模擬[D].太原： 中北大學(xué)，2014.

【水分子通過(guò)形變碳納米管的運輸行為分析論文】相關(guān)文章：

案例分析論文11-27

案例分析論文07-15

案例分析論文[優(yōu)選]07-17

【實(shí)用】案例分析論文07-17

行為金融投資管理論文11-07

《圖形變變變》中班教案07-30

論文的問(wèn)卷分析報告范文11-05

案例分析論文15篇10-12

案例分析論文15篇(經(jīng)典)07-17

笛子的演奏特點(diǎn)分析的論文09-23