石墨烯狹縫受限孔道中水分子的分子動力學模擬

趙夢堯 楊雪平 楊曉寧

(南京工業大學化學化工學院,材料化學工程國家重點實驗室,南京 210009)

石墨烯狹縫受限孔道中水分子的分子動力學模擬

趙夢堯 楊雪平 楊曉寧*

(南京工業大學化學化工學院,材料化學工程國家重點實驗室,南京 210009)

石墨烯是一種具有廣泛應用前景的納米材料,特別是由石墨烯片層自組裝形成的二維納米通道能夠應用于物質的過濾分離.本文采用分子動力學模擬方法研究了原態石墨烯/羥基改性石墨烯狹縫孔道中水分子的微觀行為,模擬計算了水的界面結構性質和擴散動力學性質,所研究的石墨烯孔寬為0.6–1.5 nm.模擬結果表明,在石墨烯狹縫孔道中,水分子受限結構呈現層狀分布,在超微石墨烯孔道(0.6–0.8 nm)中水分子可形成特殊的環狀有序結構,石墨烯表面可誘導產生特殊的水分子界面取向.在石墨烯孔道中,水分子的擴散運動低于主體相水分子的擴散運動,羥基化石墨烯孔道可以促使水分子的擴散能力降低.對于改性石墨烯狹縫孔道,由于羥基的作用,水分子可以自發滲入0.6 nm的石墨烯孔道內.模擬所得到的受限水分子的動力學性質與水分子在石墨烯孔道內的有序結構有關.本文研究結果將有助于分析理解水分子通過石墨烯納米通道的滲透機理,為設計基于石墨烯的納米膜提供理論指導.

分子模擬; 石墨烯; 水; 受限結構

Key Words:Molecular simulation; Graphene; Water; Confined structure

1 引 言

水在納米尺寸環境(例如生物分子腔、沸石膜、二氧化硅納米孔)1–3的行為已引起廣泛關注.實驗和理論4–11報道受限水具有不同于體態水的流動行為和特殊結構.X射線衍射實驗12發現水在溫度為148到303 K,直徑為1.13 nm的碳納米管下可以形成像固體一樣的結構.Koga等13發現在壓力50–500 MPa下,水分子在直徑為1.11–1.26 nm的椅型碳納米管中可以形成n-多邊形環狀結構.

石墨烯具有大的比表面積、優越的電子性能和機械性能等特點,可用于超級電容器、光催化、納米過濾和鹽水處理等.14–18石墨烯層狀膜由單層石墨烯堆砌而成,具有納米尺度的層間距,19,20這個特殊的層間結構可以允許水分子通過,而截留其他分子.Han等21實驗設計了帶負電的超薄石墨烯納米過濾膜用于水的凈化,發現對有機染料具有高達99%的截留率,對于鹽離子的截留率為20%–60%,同時發現純水滲透通過該石墨烯膜時具有很高的流率.

氧化石墨烯膜是石墨烯表面經過氧化的一種特殊的石墨烯層狀材料,其主要制備方法有:化學氧化法和電化學氧化法.氧化石墨烯表面存在多種官能團,如羥基、環氧基、羧基等,而不同的氧化制備工藝將會影響氧化石墨烯表面官能團的分布結構及氧化程度.由于氧化石墨烯制造工藝簡單、價格低廉,且層間距具有可調節性等優點,氧化石墨烯引起廣泛的關注.在干燥條件下,氧化石墨烯膜層間距為～0.3 nm,22在潮濕的空氣或浸在水中時,由于水化作用,狹縫寬度可增加至～1.0 nm.23Nair等24實驗報道了水滲透通過氧化石墨烯膜未受阻礙,而氦不能通過氧化石墨烯膜.Joshi等23實驗研究了利用氧化石墨烯作為分離膜,結果發現小離子可快速滲透通過該膜.近期研究25–27表明,使用不同官能團對氧化石墨烯膜進行修飾,可以調節層間距,達到分離廢水中特定有機分子及離子的效果.

為了進一步開發上述石墨烯膜在實際過程中的應用,需更深入地了解水溶劑在原態石墨烯/親水性石墨烯層狀狹縫納米孔道中的相關性質.由于現有實驗方法很難直接觀測到如此小的受限空間內分子行為,分子模擬成為研究受限性質的有效方法.28–31之前已有文獻4,21,23,24,32報道采用分子模擬方法考察水分子通過石墨烯納米通道時流動行為,發現水分子可以快速通過石墨烯孔道,研究認為該現象主要由于石墨烯表面原子與水分子非常弱的相互作用能產生的,但這些研究主要針對納米尺度在大于2.0 nm的范圍.對于受限孔寬在0.6–1.5 nm的超微孔中,水分子結構和動力學行為將發生明顯改變,33其微觀行為的研究仍然是一個前沿熱點科學問題.Bocquet等33發現當扶手椅式碳納米管半徑小于0.5 nm時,管內水分子的結構性質發生改變可導致水分子發生快速滑移現象,這種小孔徑石墨烯納米孔道非常適于水的滲透分離.此外,目前對兩種不同親水特性的石墨烯孔道對受限水的作用缺乏深入的研究,因此,考察水分子在0.6–1.5 nm的原態石墨烯/羥基化石墨烯納米通道內的結構特征與動力學性質十分重要.此外,以前的模擬研究28,34–36往往采用將水分子直接嵌入孔道的周期性邊界方法,未考慮主體相中水分子與孔內分子交互運動的影響,從而可能導致孔內水分子密度不是相應模擬體系下的自然密度,使模擬結果存在一定的差異.

本文采用分子動力學模擬的方法,研究孔寬為0.6–1.5 nm的原態石墨烯/親水性石墨烯狹縫孔中,水分子結構及動力學行為.由于氧化石墨烯結構的多樣性和復雜性,為了突出研究的關鍵問題,在本研究中,通過直接的羥基修飾來考察親水性表面對受限水分子的結構及動力學行為的影響.在初始構型中,我們采用將石墨烯納米狹縫孔道的兩邊與外部主體水相(water reservoir)相連,這樣可以保持外部水分子在納米孔道內的水分子存在相互作用,更為接近真實體系下水分子穿越石墨烯狹縫孔道的結構和動力學性質,較為合理地反映實際膜孔道的特征.通過本文的研究,可以進一步幫助認識受限環境對水分子結構和運動的影響,更為具體地了解水分子在原態石墨烯/羥基化石墨烯狹縫孔道內的結構性質及動力學行為的差異.進而有助于對石墨烯納米通道的設計、改進,為石墨烯膜在分離過程中的應用提供理論指導.

2 模擬方法

2.1 模擬模型

模擬模型如圖1所示,從圖中可以看出,納米通道由兩層構成,每層是由3塊石墨烯組成,石墨烯之間采用Bernal(ABA)模型排布,37層與層之間的距離為0.335 nm.38,39石墨烯板x、y方向的尺寸分別為2.3332和5.1048 nm,狹縫寬度H分別為0.6、0.7、0.8、0.9、1.0、1.5 nm,模擬盒子大小:Lx=2.475 nm,Ly=15.84 nm,Lz=6.912 nm.

圖1 模擬體系的構型Fig.1 Configuration of the simulation system (a) the pristine graphene slits;(b) the hydroxyl modified graphene slits;(c) the hydroxyl arrangement on the graphene surface;H:slit pore width

本研究中水分子選取SPC/E40模型,其中O―H鍵的鍵長為0.1 nm,H―O―H的鍵角是109.47°,模擬過程中石墨烯的碳原子電荷采用中性處理,其中C―C的Lennard-Jones(LJ)相互作用是采用Cheng和Steele41報道的參數.羥基化石墨烯中的羥基直接連接在石墨烯表面.Xu等42研究表明:當羥基分布在石墨烯平面時,羥基官能團呈隨機分布狀態與在表面呈“晶格型”有序分布狀態對計算所得的水分子潤濕角變化影響較小,而在較為復雜的排布情況下,如:當羥基在石墨烯表面呈局部團簇排布時,水分子的潤濕角可發生一定的變化,該潤濕行為的變化可能影響水分子動力學行為.為了減少復雜因素和突出表面親水特征,我們采用文獻43中的處理方法來構建一種親水石墨烯模型,即修飾的羥基官能團以矩形(0.491 nm × 0.425 nm)點陣結構在板面有序排列,表面覆蓋率近似為13%,羥基中的氧原子和氫原子之間的距離γOH=0.0945 nm,氧原子與被連接的碳原子的距離γCO=0.143 nm,OH鍵與石墨烯表面法向量成71.5°夾角,該OH取向角對應著普遍化OPLS力場中的平衡位置的烷基醇C―O―H間的平衡夾角.44在模擬過程中,采取將修飾官能團整體保持固定不動的簡化模型,并以此考察親水性表面對受限水分子動力學性質的影響.盡管這種親水性石墨烯不能完全代表氧化石墨烯,但是其結果可以用來比較兩種不同石墨烯表面對受限水的作用,類似的表面官能團處理方式在以前文獻45,46中也有報道.不同原子之間的LJ相互作用參數采用的是Lorentz-Berthelot混合規則,具體模型參數如表1所示.

表1 模擬使用的力場參數Table 1 Force field parameters used in the simulation

2.2 模擬細節

本文模擬使用LAMMPS軟件,47選取NVT系綜,在溫度為298 K下考察原態石墨烯/羥基化石墨烯孔道對水分子動力學性質的影響.整個模擬過程中,原態石墨烯/羥基化石墨烯浸在水主體相中,并保持整體結構固定.運動方程采用velocity-Verlet算法求解,庫侖相互作用采用particle-particle-particle-mesh(PPPM)方法計算,精確度為10–4,非鍵范德華力采用12-6 L-J勢函數計算,其L-J相互作用勢能和靜電相互作用的截斷半徑都為1.2 nm.主體相中水分子個數為8204,密度控制為1.0 g·cm–3,原態石墨烯/羥基化石墨烯模擬體系中總原子個數分別為27492、27732.模擬采用三維周期性邊界條件,時間步長采用1 fs,總的模擬時間長達3 ns.

3 結果與討論

3.1 受限水的結構性質

3.1.1 密度分布

首先我們考察了水分子滲入原態石墨烯/羥基化石墨烯狹縫孔道中動力學行為.圖2所示為受限孔道水分子數目隨模擬時間的變化,在原態石墨烯狹縫孔道中,當模擬時間達到100 ps時,孔道內水分子數目達到平衡,當納米孔寬為0.6 nm時,水分子不能自發地進入石墨烯狹縫孔中,這與文獻47報道相一致,這是由于憎水性微孔毛細管壓力阻礙了水的進入.這也體現了水分子進入原態石墨烯狹縫孔道內獲得的范德華能量不足以補償其進入孔道內氫鍵變化而損耗的能量.48,49隨著孔寬的升高,狹縫內水分子數目增加;對于羥基化石墨烯狹縫孔道,在羥基的阻礙作用下,水分子數目需更長的時間才能達到平衡,并且由于羥基的空間效應,孔道內水分子的數目相比于同樣原態石墨烯狹縫孔道要小.不過,對于羥基化的狹縫,在0.6 nm的孔道內仍有水分子自發滲入.這是由該親水性孔道毛細管推動力所導致的,我們的結果與Boukhvalov等50的報道相一致.

圖2 狹縫孔中水分子數目隨時間的變化Fig.2 Changes of number of water molecules over time inside slit pores(a) the pristine graphene slit pores;(b) the hydroxyl modified graphene slit pores

水分子在原態石墨烯和羥基修飾石墨烯狹縫孔中的密度分布見圖3.對于原態石墨烯狹縫孔道,如圖3(a)所示,水的密度分布關于狹縫中心對稱并呈現出明顯的層狀結構,在靠近碳板界面處存在明顯的流體減少,同時在距離石墨烯表面～0.32 nm處,呈現水密度最大峰.8,38并且由圖可知,在0.7 nm孔中形成單層的結構,0.8 nm中是介于單層與雙層之間的過渡層,0.9和1.0 nm都是出現雙層結構,1.5 nm出現了四層結構.出現不同分層結構不僅與狹縫寬度有關同時還與σCO有關,通過我們研究發現層狀結構與狹縫寬度存在如下近似關系:

式中σCO=0.3283 nm,σOO=0.3166 nm.其中H是狹縫寬度,N是在狹縫中水形成的層數.由此近似方程可估算形成單層和雙層狹縫寬度分別是0.66和0.97 nm,這與模擬的結果相一致,同時可推測狹縫寬度為1.5 nm時能夠形成四層結構,這與1.5 nm時的密度分布圖呈現四層相吻合.對于0.7 nm原態石墨烯狹縫孔道,Najafi等34與Hirunsit和Balbuena36采用與本文相同的力場參數,而水分子采取直接嵌入的周期邊界孔模型方法時,得出水分子的層狀數目為2,而Deshmukh等51采用與本文類似的由主體相自發進入狹縫孔道的模型方法時,水分子同樣形成單層分布狀態.因此,石墨烯狹縫孔道的不同方式對水分子的受限結構存在一定的影響,這些不一致的結果反映了不同孔道模型導致了受限孔中流體有效密度差異,進而影響受限結構.

對于羥基化石墨烯納米孔道,從圖3(b)可看出水的密度分布同樣關于狹縫中心對稱并呈現出明顯的層狀結構.進一步對比分析可知,當羥基化石墨烯狹縫寬度為0.8、1.0、1.5 nm時,形成的水層數目與原態石墨烯狹縫有所不同,這主要由于羥基的親水性及狹縫有效寬度協同作用的結果.值得注意的是,在1.0 nm的羥基化石墨烯狹縫孔道內,水分子呈現三層分布狀態.而在相同孔寬條件下,Xu等52采用直接嵌入的方式模擬得出水分子在羥基化石墨烯孔道內的層狀數目卻為2.該模擬結果的差異可能同樣由于孔道模擬不同所致.

從水分子的結構模擬構型圖3(c–f)可看出:0.7 nm原態石墨烯狹縫孔中水分子呈現單層分布,并形成四元環結構.5,36與原態石墨烯孔道不同,在圖3(d,e)中,羥基化狹縫孔道內水分子形成單層分布,通過構型圖可以發現:在0.6 nm改性的石墨烯狹縫孔道中大多數水分子可形成四元環結構,而在0.7 nm改性孔道中形成了六元環結構.5對于0.8 nm的羥基改性石墨烯狹縫孔,水分子形成雙層結構,在構型圖3(f)中展示了第一層水分子結構,盡管部分第二層水分子顯露在第一層水分子的上方,但整體上,我們可以清楚觀察到水分子在狹縫孔道內呈線性分布方式.對于其他尺度狹縫孔道沒有出現特殊的結構.

3.1.2 徑向分布函數

圖3 受限水分子沿z軸的密度分布圖及在狹縫孔道內的構型圖Fig.3 Density profiles along the z-axis and the configurations inside the slit pores for confined water molecules

圖4 不同狹縫孔道內靠近石墨烯板面水分子的徑向分布函數(gO-O(r))Fig.4 Lateral radial distribution function(gO-O(r)) for the water molecules adjacent to the graphene wall in different slit pores

徑向分布函數(RDF)是表征體系在原子水平上結構性質的一種重要性質,本文采用RDF分析固-液界面層水的結構.從圖4(a)和4(b)中可以看出受限體系和主體相中gO-O(r)的第一個峰大約在0.27 nm處,與文獻32,34報道的第一個峰位置相一致,該峰的位置對應于水分子之間第一配位圈距離,此距離不受石墨烯表面及其基團性質的影響.根據圖形可以看出體態水的第二個比較寬的峰大約在0.45 nm處,其它文獻40也得到了類似的結果,而當H=0.8,0.9,1.0,1.5 nm時第二個比較寬的峰大約在0.55 nm處.尤其可以看出,對于H=0.7 nm,發現接近界面水層的RDF表現為長程有序的特征,這表明水在狹縫為0.7 nm的受限狹縫孔中形成一種類固態的結構.34從圖4(b)可以看出0.6、0.7、0.8 nm孔道內,水分子的RDF也表現為長程有序的特征,表明其在改性狹縫孔中同樣形成類似的固態結構,這與圖3(d–f)中所展現的構型圖是一致的.當羥基修飾的狹縫孔大于0.8 nm時,RDF沒有表現出長程有序的特征.

3.1.3 取向分布

為了進一步研究水分子在孔道內的分布狀態,我們分析水的OH鍵與表面法向量的夾角φ,以及偶極方向與表面法向量的夾角ψ來判定水分子在孔道內的取向特征,結果如圖5所示.取向分布顯示,在原態狹縫孔道中,水的OH鍵及水分子的偶極與z軸的夾角在約90°處達到最大值,說明水分子以平行于該石墨烯板的形式存在,這將有利于受限孔中氫鍵的形成,這與之前36,53研究結論相符.從圖中還可以看出,當石墨烯狹縫寬度大于0.9 nm時,由于水分子的運動加劇,取向分布的范圍更廣.在羥基改性的狹縫孔道中,當孔寬為0.6、0.7 nm時,與原態石墨烯狹縫孔道內水分子的羥基取向分布對比可以發現,水分子的OH鍵與z軸的夾角在35°、145°角度處出現明顯地凸起,同時,在對應的羥基化石墨烯狹縫孔道中,水分子的偶極與z軸的夾角在90°處出現最大值,表明水分子的偶極方向平行于板面,兩種取向結合可以說明:當狹縫孔寬為0.6–0.7 nm時,相比于原態石墨烯狹縫孔道,在羥基功能化的石墨烯狹縫孔道中,更多水分子的―OH鍵傾向于指向改性石墨烯表面,主要由于部分水分子與羥基化石墨烯表面的羥基可能形成氫鍵作用.由此可知羥基的存在影響水分子在狹縫孔道內的取向分布.當狹縫寬度超過0.7 nm時,φ和ψ都主要集中在60°–120°之間,取向角度分布變得相對無序.

3.2 受限水的動力學性質

3.2.1 均方位移

受限體系中水分子的動力學性質與體態水不同.因此,我們進一步計算了平行于納米通道方向水的均方位移(MSDxy),如圖6(a,b)所示.水的MSDxy要比主體相要小,表明微孔對其擴散施加一定的限定作用.基于MSDxy的斜率,可計算出主體相的擴散系數是2.63 × 10–9m2·s–1,這與文獻10報道的SPC/E水模型擴散系數(2.67 × 10–9m2·s–1)有很好的吻合.水分子在原態石墨烯狹縫孔道中,擴散系數表現出非單調行為.相比于其它狹縫孔道,在H=0.8 nm的水分子呈現較大的擴散系數,這可能由于該孔道水分子受限密度較低所致.在羥基化石墨烯狹縫孔道中,水分子的均方位移隨著孔寬的增大呈上升趨勢.由于羥基的親水性及孔道內有效空間的變化,水分子的擴散能力比相同寬度原態石墨烯狹縫孔道低1–2個數量級.

圖5 受限水分子的取向角(φ,ψ)分布Fig.5 Orientation angle(φ,ψ) distributions of confined water molecules(a,b) the pristine graphene slit pores;(c,d) the hydroxyl modified graphene slit pores

圖6 (a,b)不同受限環境及主體相中多種類型水分子xy平面的MSDxy;(c,d)受限水分子沿z軸的運動軌跡Fig.6 (a,b) Mean square displacement of xy plane(MSDxy) for several typical water molecules in different confinement environments and bulk phase;(c,d) trajectories of the confined water molecules along the z-axis (a,c) the pristine graphene slit pores;(b,d) the hydroxyl modified graphene slit pores

對于受限水的運動性質,本文進一步分析孔道內水分子在z方向的位置隨時間的變化.其圖形54,55主要根據輸出的分子運動軌跡數據,通過隨機選取受限孔道內一定水分子,對其z方向的位置隨時間的變化進行數據處理,根據該圖形(圖6(c,d)),我們可以直接觀察水分子沿垂直方向的運動行為隨著孔寬的變化.從圖形可知,水在該方向的運動軌跡與前面的密度分布與孔寬的變化行為一致,表明其呈現明顯的層狀結構,并且隨著狹縫孔寬的增大,水分子的運動變劇烈.

3.2.2 氫 鍵

分子在受限體系下的動力學性質不僅與液-固之間的相互作用有關,液-液之間的氫鍵(HB)作用也占主導地位,圖7(a,b)為水分子沿著表面法向(z軸)的氫鍵分布.氫鍵的定義是依據幾何標準56提出的:O···O ―H的夾角不超過30°,O ···O的距離不超過0.35 nm.從圖7中可以看出氫鍵分布基本上關于狹縫中心對稱.在原態體系中,平均氫鍵數目介于3.13 –3.47,其中H=0.7 nm形成的平均氫鍵數目最多,而H=0.8 nm形成的平均氫鍵數目最少.對于H=0.7 nm,從構型圖可知形成四元環形結構從而容易形成氫鍵,而這種氫鍵網絡結構不利于單個水分子在納米通道中熱運動.34,57在0.8 nm的納米通道中水分子間相對少的氫鍵數目對應其較大的擴散系數;對于羥基化石墨烯孔道,在固-液界面層,水分子之間形成的氫鍵數目明顯低于原態石墨烯孔道,特別是在小于0.9 nm的羥基化石墨烯孔道.這主要是由于表面羥基與水分子更容易形成氫鍵,從而削弱了受限水分子之間形成氫鍵的能力.同時可以看出受限于羥基修飾的石墨烯狹縫中水分子平均氫鍵數目明顯比主體相的平均氫鍵數目低,隨著狹縫寬度的增大,平均的氫鍵數目也相應增大.

3.2.3 氫鍵相關函數

我們進一步分析了受限體系氫鍵穩定性的動力學行為,這里采用間歇性相關函數進行分析,其定義如下方程表示.58,59

在式(2)中如果某兩個水分子在t時刻形成氫鍵則h(t)=1,否則,h(t)=0.CHB(t)能夠反映氫鍵的結構松弛動力學性質.由圖得知受限水的CHB(t)衰減比主體相水慢.同時從圖7(c)中發現狹縫寬度為0.7 nm時,CHB(t)衰減非常緩慢,可認為主要是由于0.7 nm狹縫中水分子之間形成的四元環結構使得氫鍵十分牢靠.同樣,在羥基修飾石墨烯狹縫孔道中,寬度為0.6、0.7、0.8 nm時,CHB(t)表現出類似的性質,同時衰減更為緩慢.總體看,對于兩種石墨烯表面都表現為隨著狹縫寬度的增大,受限水的CHB(t)衰減變得越來越快,這表明CHB(t)的衰減很大程度上取決于孔道寬度.CHB(t)衰減的越慢,形成的氫鍵越穩定,越不利于水分子在狹縫中的運動.

圖7 (a,b)狹縫孔道內水分子平均氫鍵數沿z軸的分布;(c,d)受限水及體態水氫鍵的間歇性時間相關函數(CHB(t))Fig.7 (a,b) Average number of hydrogen bond profiles along the z axis for water inside slit pores;(c,d) the intermittent time correlation function(CHB(t)) for the H2O-H2O hydrogen bond numbers inside the nanochannnel as well as bulk water (a,c) the pristine graphene slit pores;(b,d) the hydroxyl modified graphene slit pores

4 結 論

利用分子動力學模擬的方法,研究了孔寬為0.6–1.5 nm范圍內,原態石墨烯和羥基化石墨烯狹縫孔道中水分子結構性質及動力學行為.通過研究發現水分子不能自發進入0.6 nm原態石墨烯狹縫孔內,而在同等寬度的羥基化狹縫中,由于親水性毛細管推動力的作用,水分子能夠自發滲入.通過密度及取向分析,狹縫孔內水分子形成明顯的層狀結構,并且在原態石墨烯狹縫孔道內水分子以平行于碳板表面的形式存在,而在0.6–0.7 nm改性的石墨烯狹縫孔道內,部分水分子氫鍵指向碳板平面,與表面的羥基官能團形成氫鍵作用.在超微(0.6–0.8 nm)狹縫孔內,水分子之間可以形成四元環、六元環等類似固態一樣的有序結構,這種特殊構型促進水分子之間形成氫鍵,從而不利于水分子的擴散運動.其次,在改性的石墨烯狹縫孔道內,由于羥基的親水性及空間效應,水分子的擴散能力比同等原態石墨烯孔道要低1–2個數量級,而在0.8 nm原態石墨烯狹縫孔道內,由于水分子之間形成的平均氫鍵數目最少,水分子表現出相對較快的擴散能力.相比于主體相,在受限體系下,水分子之間形成的氫鍵更穩定,并且隨著狹縫寬度的減小,氫鍵的CHB(t)衰減越慢.相比較,羥基化石墨烯受限孔道中水分子形成氫鍵結構穩定性更強.綜上所述,石墨烯狹縫孔道的尺寸效應及表面性質對受限水分子的結構性質及動力學行為有著直接的影響.

(1)Pan,Y.S.;Birkedal,H.;Pattison,P.;Brown,D.;Chapuis,G.J.Phys.Chem.B 2004,108(20),6458.doi:10.1021/jp037219v

(2)Newsome,D.A.;Sholl,D.S.J.Phys.Chem.B 2005,109(15),7239.doi:10.1021/jp044247k

(3)Milischuk,A.A.;Ladanyi,B.M.J.Chem.Phys.2011,135(17),174709.doi:10.1063/1.3657408

(4)Qiao,Y.;Xu,X.;Li,H.Appl.Phys.Lett.2013,103(23),233106.doi:10.1063/1.4839255

(5)Han,S.;Choi,M.Y.;Kumar,P.;Stanley,H.E.Nat.Phys.2010,6(9),685.doi:10.1038/nphys1708

(6)Du,F.;Qu,L.T.;Xia,Z.H.;Feng,L.F.;Dai,L.M.Langmuir 2011,27(13),8437.doi:10.1021/la200995r

(7)Strauss,I.;Chan,H.;Král,P.J.Am.Chem.Soc.2014,136(4),1170.doi:10.1021/ja4087962

(8)Cicero,G.;Grossman,J.C.;Schwegler,E.;Gygi,F.;Galli,G.J.Am.Chem.Soc.2008,130(6),1871.doi:10.1021/ja074418+

(9)Thomas,J.A.;McGaughey,A.J.H.Nano Lett.2008,8(9),2788.doi:10.1021/nl8013617

(10)Mashl,R.J.;Joseph,S.;Aluru,N.R.;Jakobsson,E.Nano Lett.2003,3(5),589.doi:10.1021/nl0340226

(11)Liu,Y.C.;Wang,Q.;Lü,L.H.;Zhang,L.Z.Acta Phys.-Chim.Sin.2005,21(1),63.[劉迎春,王 琦,呂玲紅,章連眾.物理化學學報,2005,21(1),63.] doi:10.3866/ PKU.WHXB 20050113

(12)Iiyama,T.;Nishikawa,K.;Otowa,T.;Kaneko,K.J.Phys.Chem.1995,99(25),10075.doi:10.1021/j100025a004

(13)Koga,K.;Gao,G.T.;Tanaka,H.;Zeng,X.C.Nature 2001,412(6849),802.doi:10.1038/35090532

(14)Stoller,M.D.;Park,S.;Zhu,Y.W.;An,J.H.;Ruoff,R.S.Nano Lett.2008,8(10),3498.doi:10.1021/nl802558y

(15)Chandra,V.;Park,J.;Chun,Y.;Lee,J.W.;Hwang,I.C.;Kim,K.S.ACS Nano 2010,4(7),3979.doi:10.1021/ nn1008897

(16)Zhang,H.;Lv,X.J.;Li,Y.M.;Wang,Y.;Li,J.H.ACS Nano 2010,4(1),380.doi:10.1021/nn901221k

(17)Cohen-Tanugi,D.;Grossman,J.C.Nano Lett.2012,12(7),3602.doi:10.1021/nl3012853

(18)Hu,Y.J.;Jin,J.;Zhang,H.;Wu,P.;Cai,C.X.Acta Phys.-Chim.Sin.2010,26(8),2073.[胡耀娟,金 娟,張 卉,吳 萍,蔡稱心.物理化學學報,2010,26(8),2073.] doi:10.3866/PKU.WHXB20100812

(19)Chen,H.Q.;Müeller,M.B.;Gilmore,K.J.;Wallace,G.G.;Li,D.Adv.Mater.2008,20(18),3557.doi:10.1002/adma.200800757

(20)Li,D.;Mueller,M.B.;Gilje,S.;Kaner,R.B.;Wallace,G.G.Nat.Nanotechnol.2008,3(2),101.doi:10.1038/nnano.2007.451

(21)Han,Y.;Xu,Z.;Gao,C.Adv.Funct.Mater.2013,23(29),3693.doi:10.1002/adfm.v23.29

(22)Mi,B.X.Science 2014,343(6172),740.doi:10.1126/ science.1250247

(23)Joshi,R.K.;Carbone,P.;Wang,F.C.;Kravets,V.G.;Su,Y.;Grigorieva,I.V.;Wu,H.A.;Geim,A.K.;Nair,R.R.Science 2014,343(6172),752.doi:10.1126/science.1245711

(24)Nair,R.R.;Wu,H.A.;Jayaram,P.N.;Grigorieva,I.V.;Geim,A.K.Science 2012,335(6067),442.doi:10.1126/ science.1211694

(25)Sun,P.Z.;Zhu,M.;Wang,K.L.;Zhong,M.L.;Wei,J.Q.;Wu,D.H.;Xu,Z.P.;Zhu,H.W.ACS Nano 2013,7(1),428.doi:10.1021/nn304471w

(26)Sun,P.Z.;Zheng,F.;Zhu,M.;Song,Z.G.;Wang,K.L.;Zhong,M.L.;Wu,D.H.;Little,R.B.;Xu,Z.P.;Zhu,H.W.ACS Nano 2014,8(1),850.doi:10.1021/nn4055682

(27)Hu,M.;Mi,B.X.Environ.Sci.Technol.2013,47(8),3715.doi:10.1021/es400571g

(28)Xu,L.;Hu,Y.Z.;Ma,T.B.;Wang,H.Nanotechnology 2013,24(50),505504.doi:10.1088/0957-4484/24/50/505504

(29)Kolesnikov,A.I.;Zanotti,J.M.;Loong,C.K.;Thiyagarajan,P.;Moravsky,A.P.;Loutfy,R.O.;Burnham,C.J.Phys.Rev.Lett.2004,93(3),035503.doi:10.1103/PhysRevLett.93.035503

(30)Fernández-Serra,M.V.;Artacho,E.Phys.Rev.Lett.2006,96(1),016404.doi:10.1103/PhysRevLett.96.016404

(31)Gao,W.X.;Wang,H.L.;Li,S.M.Acta Phys.-Chim.Sin.2014,30(9),1625.[高文秀,王洪磊,李慎敏.物理化學學報,2014,30(9),1625.] doi:10.3866/PKU.WHXB201407031

(32)Xiong,W.;Liu,J.Z.;Ma,M.;Xu,Z.P.;Sheridan,J.;Zheng,Q.S.Phys.Rev.E 2011,84(5),056329.doi:10.1103/ PhysRevE.84.056329

(33)Falk,K.;Sedlmeier,F.;Joly,L.;Netz,R.R.;Bocquet,L.Nano Lett.2010,10(10),4067.doi:10.1021/nl1021046

(34)Mosaddeghi,H.;Alavi,S.;Kowsari,M.H.;Najafi,B.J.Chem.Phys.2012,137(18),184703.doi:10.1063/1.4763984

(35)Kumar,P.;Buldyrev,S.V.;Starr,F.W.;Giovambattista,N.;Stanley,H.E.Phys.Rev.E 2005,72(5),051503.doi:10.1103/ PhysRevE.72.051503

(36)Hirunsit,P.;Balbuena,P.B.J.Phys.Chem.C 2007,111(4),1709.doi:10.1021/jp063718v

(37)Warner,J.H.;Mukai,M.;Kirkland,A.I.ACS Nano 2012,6(6),5680.doi:10.1021/nn3017926

(38)Argyris,D.;Tummala,N.R.;Striolo,A.;Cole,D.R.J.Phys.Chem.C 2008,112(35),13587.doi:10.1021/jp803234a

(39)Liu,L.;Zhang,L.;Sun,Z.G.;Xi,G.Nanoscale 2012,4(20),6279.doi:10.1039/c2nr31847d

(40)Mark,P.;Nilsson,L.J.Phys.Chem.A 2001,105(43),9954.doi:10.1021/jp003020w

(41)Cheng,A.;Steele,W.A.J.Chem.Phys.1990,92(6),3858.doi:10.1063/1.458562

(42)Wei,N.;Lv,C.J.;Xu,Z.P.Langmuir 2014,30(12),3572.doi:10.1021/la500513x

(43)Jane?ek,J.;Netz,R.R.Langmuir 2007,23(16),8417.doi:10.1021/la700561q

(44)Cornell,W.D.;Cieplak,P.;Bayly,C.I.;Gould,I.R.;Merz,K.M.;Ferguson,D.M.;Spellmeyer,D.C.;Fox,T.;Caldwell,J.W.;Kollman,P.A.J.Am.Chem.Soc.1995,117(19),5179.doi:10.1021/ja00124a002

(45)Gotzias,A.;Tylianakis,E.;Froudakis,G.;Steriotis,T.Microporous Mesoporous Mat.2012,154,38.doi:10.1016/ j.micromeso.2011.10.011

(46)Zhu,Y.D.;Guo,X.J.;Shao,Q.;Wei,M.J.;Wu,X.M.;Lu,L.H.;Lu,X.H.Fluid Phase Equilibr.2010,297(2),215.doi:10.1016/j.fluid.2010.05.005

(47)Eun,C.S.;Berkowitz,M.L.J.Phys.Chem.B 2010,114(42),13410.doi:10.1021/jp1072654

(48)Lum,K.;Chandler,D.;Weeks,J.D.J.Phys.Chem.B 1999,103(22),4570.doi:10.1021/jp984327m

(49)Ren,X.P.;Zhou,B.;Wang,C.L.J.Chem.Phys.2012,137(2),024703.doi:10.1063/1.4733719

(50)Boukhvalov,D.W.;Katsnelson,M.I.;Son,Y.W.Nano Lett.2013,13(8),3930.doi:10.1021/nl4020292

(51)Deshmukh,S.A.;Kamath,G.;Baker,G.A.;Sumant,A.V.;Sankaranarayanan,S.K.R.S.Surf.Sci.2013,609,129.doi:10.1016/j.susc.2012.11.017

(52)Wei,N.;Peng,X.S.;Xu,Z.P.ACS Appl.Mater.Inter.2014,6(8),5877.doi:10.1021/am500777b

(53)Pertsin,A.;Grunze,M.J.Phys.Chem.B 2004,108(4),1357.doi:10.1021/jp0356968

(54)Hub,J.S.;Winkler,F.K.;Merrick,M.;de Groot,B.L.D.J.Am.Chem.Soc.2010,132(38),13251.doi:10.1021/ jp0356968

(55)Zang,J.;Konduri,S.;Nair,S.;Sholl,D.S.ACS Nano 2009,3(6),1548.doi:10.1021/nn9001837

(56)Luzar,A.;Chandler,D.Nature 1996,379(6560),55.doi:10.1021/jp044247k

(57)Striolo,A.Nano Lett.2006,6(4),633.doi:10.1038/379055a0

(58)Martí,J.;Sala,J.;Guàrdia,E.J.Mol.Liq.2010,153(1),72.doi:10.1016/j.molliq.2009.09.015

(59)Jorgensen,W.L.;Chandrasekhar,J.;Madura,J.D.;Impey,R.W.;Klein,M.L.J.Chem.Phys.1983,79(2),926.doi:10.1063/1.445869

Molecular Dynamics Simulation of Water Molecules in Confined Slit Pores of Graphene

ZHAO Meng-Yao YANG Xue-Ping YANG Xiao-Ning*
(State Key Laboratory of Materials-Oriented Chemical Engineering,College of Chemistry and Chemical Engineering,Nanjing Tech University,Nanjing 210009,P.R.China)

Graphene has potential applications in many fields.In particular,two-dimensional graphene nanochannels assembled from graphene sheets can be used for filtration and separation.In this work,molecular dynamics simulations were performed to investigate the microscopic structural and dynamical properties of water molecules confined in pristine and hydroxyl-modified graphene slit pores with widths of 0.6–1.5 nm.The simulation results indicate that water molecules have layered structure distributions within the graphene nanoscale channels.The special ordered ring structure can be formed for water confined in the subnanometer pores(0.6–0.8 nm).Graphene surfaces are able to induce distinctive molecular interfacial orientations of water molecules.In the graphene slits,the diffusion of water molecules was slower than that in bulk water,and the hydroxyl-modified graphene pores could lead to more reduced water diffusion ability.For the hydroxyl-modified graphene pores,water molecules spontaneously permeated into the 0.6 nm slit pore.According to the simulation results,the dynamic behavior of confined water is associated with the ordered water structures confined within the graphene-based nanochannels.These simulation results will be helpful in understanding the penetration mechanism of water molecules through graphene nanochannels,and will provide a guide for designing graphene-based membrane structures.

November 26,2014;Revised:May 28,2015;Published on Web:June 1,2015.

O641

icle]

10.3866/PKU.WHXB201506011 www.whxb.pku.edu.cn

*Corresponding author.Email:Yangxia@njut.edu.cn.

The project was supported by the National Natural Science Foundation of China(21376116).

國家自然科學基金(21376116)資助項目? Editorial office of Acta Physico-Chimica Sinica