直流電場(chǎng)下水中石墨烯定向行為研究?

董若宇 曹鵬 曹桂興 胡幗杰 曹炳陽(yáng)?

1)(清華大學(xué)航天航空學(xué)院,熱科學(xué)與動(dòng)力工程教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)2)(中國(guó)空間技術(shù)研究院通信衛(wèi)星事業(yè)部,北京 100094)(2016年7月26日收到;2016年10月13日收到修改稿)

直流電場(chǎng)下水中石墨烯定向行為研究?

董若宇1)曹鵬2)曹桂興2)胡幗杰2)曹炳陽(yáng)1)?

1)(清華大學(xué)航天航空學(xué)院,熱科學(xué)與動(dòng)力工程教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)2)(中國(guó)空間技術(shù)研究院通信衛(wèi)星事業(yè)部,北京 100094)(2016年7月26日收到;2016年10月13日收到修改稿)

納米顆粒在液體環(huán)境中的定向控制與系統(tǒng)物理性能的調(diào)控及新型納米器件的制備等應(yīng)用領(lǐng)域密切相關(guān).本文使用分子動(dòng)力學(xué)模擬方法,研究了水中單片不帶電矩形石墨烯在直流電場(chǎng)下的定向行為.結(jié)果發(fā)現(xiàn)石墨烯平面趨向平行于電場(chǎng)方向且隨著電場(chǎng)強(qiáng)度增大定向性增強(qiáng),其主要原因在于極性水分子在電場(chǎng)下的響應(yīng)以及水合作用;減小石墨烯長(zhǎng)寬比,石墨烯法向矢量和長(zhǎng)邊矢量的定向性減弱,定量結(jié)果表明法向和長(zhǎng)邊定向度分別與繞長(zhǎng)邊和法向的轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)存在負(fù)相關(guān)關(guān)系.

石墨烯,定向,直流電場(chǎng),分子動(dòng)力學(xué)模擬

1引 言

石墨烯在21世紀(jì)之初被發(fā)現(xiàn),并一舉成為21世紀(jì)新興材料中的翹楚.作為典型的二維材料,石墨烯擁有極為優(yōu)異的熱學(xué)、力學(xué)、電學(xué)等物理性能,例如其沿面向的熱導(dǎo)率可高達(dá)約103W/mK.更為重要的是,這些特殊的性能還能移植到以傳統(tǒng)有機(jī)物或無(wú)機(jī)物作為介質(zhì)的復(fù)合系統(tǒng)中,從而在微納材料領(lǐng)域擁有廣泛的應(yīng)用前景.

在外場(chǎng)的作用下,液體環(huán)境中的非球形納米顆粒會(huì)受到外場(chǎng)和隨機(jī)布朗運(yùn)動(dòng)兩者的共同影響,呈現(xiàn)某些確定性的運(yùn)動(dòng)形式,其中之一是沿某方向的定向行為.對(duì)顆粒的定向控制一方面可應(yīng)用于改善復(fù)合系統(tǒng)的物理性能[1?3],例如,納米流體[4]中添加的碳納米管或石墨烯顆粒分布雜亂無(wú)章,實(shí)現(xiàn)沿某方向的定向?qū)R分布可有效提高沿定向方向的系統(tǒng)熱導(dǎo)率;另一方面可基于定向獲得納米自組裝結(jié)構(gòu)并實(shí)現(xiàn)各種復(fù)雜的功能[5],例如制備獲得納米導(dǎo)線[6]、光開關(guān)[7]、自旋電子元件等[8],進(jìn)一步作為芯片實(shí)驗(yàn)室和微納機(jī)電系統(tǒng)中的關(guān)鍵組成模塊.

外加電場(chǎng)是控制液體中納米顆粒定向的有效方法:在電場(chǎng)下,懸浮液中介電質(zhì)的碳納米管及石墨烯可發(fā)生極化現(xiàn)象,其結(jié)構(gòu)的各向異性導(dǎo)致了極化率的各向異性,在非均勻的電場(chǎng)強(qiáng)度作用下可形成沿電場(chǎng)方向的定向分布及連接正負(fù)電極的網(wǎng)狀結(jié)構(gòu)[9?11].文獻(xiàn)[12]發(fā)現(xiàn)碳納米管顆粒在交流電場(chǎng)下的定向程度優(yōu)于直流電場(chǎng),且利用交流電場(chǎng)可避免由于顆粒表面所帶電荷而引起的向一端電極遷移的電泳現(xiàn)象.盡管實(shí)驗(yàn)顯示顆粒極化后受到介電泳力和介電泳轉(zhuǎn)矩是形成定向分布的主要原因[13],但分子動(dòng)力學(xué)模擬同時(shí)指出在水中表面不帶有凈電荷且未考慮極化效應(yīng)的納米顆粒也可在直流電場(chǎng)下定向,其原因可歸結(jié)為極性水分子的水合作用[5,14].

然而,現(xiàn)有文獻(xiàn)集中在對(duì)棒狀顆粒如碳納米管在電場(chǎng)下的定向行為研究[10?12],針對(duì)片狀顆粒如石墨烯的研究還很少.實(shí)際上,棒狀顆粒只需選取惟一方向矢量(軸線方向)研究其定向性規(guī)律,而不同形狀的片狀顆粒除需要考慮垂直于其平面的方向矢量,還需考慮處于平面內(nèi)的方向矢量的定向情況,因而增加了描述問題的復(fù)雜性.在我們之前的研究中,利用剪切流場(chǎng)實(shí)現(xiàn)了對(duì)碳納米管的定向控制,特別是利用轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)[15]定量表征了布朗轉(zhuǎn)動(dòng)的強(qiáng)弱并與定向性建立關(guān)系[16].本文將采用分子動(dòng)力學(xué)模擬研究在直流電場(chǎng)下單片矩形石墨烯顆粒的三維定向行為,建立垂直于石墨烯平面和平行于石墨烯長(zhǎng)邊方向的方向矢量,分別研究沿這兩個(gè)矢量方向的定向特征,通過(guò)將定向性與相應(yīng)的轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)建立對(duì)應(yīng)關(guān)系,明確石墨烯的形狀效應(yīng)對(duì)定向性的影響.

2模擬細(xì)節(jié)

利用分子動(dòng)力學(xué)模擬研究單片石墨烯在水中的轉(zhuǎn)動(dòng)和定向行為.水的勢(shì)能模型采用SPC/E[17],

不同水分子的氧原子與氧原子間的相互作用包括靜電相互作用和Lennard-Jones(LJ)勢(shì)能部分:

不同水分子間的氧原子與氫原子、氫原子與氫原子間只考慮靜電相互作用,

其中σO-O=3.166 ?,εO-O=0.1553kcal/mol,氧原子所帶電荷為qO=?0.8476e,氫原子所帶電荷為qH=?0.4238e.利用SHAKE算法[18]保持氧原子與氫原子間的鍵長(zhǎng)為1.0 ?,H—O—H的鍵角為109.47°.

石墨烯中的碳原子不帶電且不考慮極化效應(yīng),整體作為剛體處理.氧原子和碳原子之間的相互作用采用L-J勢(shì)能模型描述:

其中σO-C=3.275 ?,εO-C=0.114kcal/mol,不考慮碳原子與氫原子之間的相互作用[5,19?21].(1)和(4)式中LJ勢(shì)能截?cái)喟霃饺?.2nm.

模擬系統(tǒng)建立在直角坐標(biāo)系x,y,z上,三個(gè)方向使用周期性邊界條件.圖1是矩形石墨烯在水中初始位置示意圖,長(zhǎng)邊方向與y軸平行,短邊方向與x軸平行,石墨烯可自由地在水中平動(dòng)和轉(zhuǎn)動(dòng).分子動(dòng)力學(xué)計(jì)算利用LAMMPS軟件包完成[22],模擬過(guò)程采用NPT系綜,即系統(tǒng)恒粒子數(shù)、恒壓力和恒溫度,使用Nose-Hoover熱浴控溫,系統(tǒng)溫度設(shè)為T=300 K;Nose-Hoover壓浴控壓,系統(tǒng)壓強(qiáng)設(shè)為P=1atm,水分子總數(shù)N=4180.使用particle-mesh Ewald方法處理水分子間的長(zhǎng)程靜電相互作用[23],時(shí)間步長(zhǎng)2fs,模擬總時(shí)長(zhǎng)6.2ns,使用0.2ns使系統(tǒng)達(dá)到穩(wěn)定,其余6ns進(jìn)行數(shù)據(jù)統(tǒng)計(jì).另外,給模擬系統(tǒng)中施加一直流電場(chǎng),電場(chǎng)方向沿x方向,相當(dāng)于在氫原子和氧原子上施加一外力作用:

其中,E為電場(chǎng)強(qiáng)度,q為氫原子或氧原子所帶電荷量.

圖1 (網(wǎng)刊彩色)模擬系統(tǒng)初始狀態(tài)示意圖Fig.1.(color online)Schematic diagram of the initial con fi guration of the simulation system.

3結(jié)果與討論

首先,如圖2定義方向矢量L和N,其中L平行于石墨烯長(zhǎng)邊,N垂直于石墨烯平面.在圖3中對(duì)比了平衡態(tài)未施加電場(chǎng)和電場(chǎng)強(qiáng)度E=1.0 V/nm時(shí),方向矢量L和N與電場(chǎng)方向夾角φN-E,φL-E隨時(shí)間的變化關(guān)系,此時(shí)石墨烯尺寸取為長(zhǎng)L=3.26nm,寬W=0.74nm,長(zhǎng)寬比L/W=4.4.可以發(fā)現(xiàn)在平衡態(tài)下,兩個(gè)角度隨時(shí)間變化呈現(xiàn)無(wú)規(guī)則的隨機(jī)特征,而在強(qiáng)度為1.0 V/nm的電場(chǎng)作用下,石墨烯的轉(zhuǎn)動(dòng)行為體現(xiàn)出良好的定向性:其中,法向矢量N與電場(chǎng)方向夾角φN-E分布在90°附近,表明石墨烯平面平行于電場(chǎng)方向;而長(zhǎng)邊矢量L與電場(chǎng)方向夾角φL-E大致分布在150°—180°之間,考慮到0—90°和180°—90°的等價(jià)性,表明此時(shí)石墨烯長(zhǎng)邊與電場(chǎng)方向夾角在0°—3°之間.

圖2 石墨烯單位方向矢量定義Fig.2.The de fi nition of the unit vectors of graphene.

進(jìn)一步改變電場(chǎng)強(qiáng)度,如圖4所示為φN-E和φL-E的角度概率密度分布.可以發(fā)現(xiàn)φN-E概率密度峰值分布在90°附近(圖4(a)),且隨著電場(chǎng)強(qiáng)度增大,峰值增大,定向性增強(qiáng);而φL-E概率分布的規(guī)律性卻并不顯著 (圖4(b):由于φL-E的取值會(huì)出現(xiàn)跨越0°和180°的情況,為了不使概率分布出現(xiàn)兩個(gè)等價(jià)的峰值,將角度范圍取為0°—90°),在E=0.01 V/nm時(shí),φL-E分布隨機(jī),而在更大的電場(chǎng)強(qiáng)度下,φL-E出現(xiàn)了沿某角度的定向現(xiàn)象,但該角度卻并不惟一,取值范圍大致處于5°—50°的區(qū)間.

不帶電且不進(jìn)行極化處理的石墨烯在電場(chǎng)下呈現(xiàn)定向行為,其主要機(jī)理應(yīng)為帶電水分子在電場(chǎng)下的響應(yīng),而石墨烯通過(guò)和水分子之間的相互作用來(lái)“感知”外加電場(chǎng),并產(chǎn)生運(yùn)動(dòng)行為的改變.實(shí)際上,對(duì)石墨烯顆粒作用最強(qiáng)的水分子分布在第一水合層中,如圖5所示(未施加電場(chǎng),石墨烯長(zhǎng)寬比為4.4),氧原子在石墨烯法向分層分布,處在第一個(gè)概率密度波谷以內(nèi)的部分稱為第一水合層[5],而下面將研究第一水合層中水分子在電場(chǎng)下的分布規(guī)律.

圖3 (網(wǎng)刊彩色)長(zhǎng)寬比為L(zhǎng)/W=4.4的石墨烯方向矢量N,L與電場(chǎng)方向(未施加電場(chǎng)時(shí)為x方向)的夾角(a)φN-E,(b)φL-E隨時(shí)間的變化Fig.3.(color online)The time-varying angles(a)φN-E,(b)φL-Ebetween the vectors N,L and the electric direction(x direction when the electric fi eld is o ff)for the graphene with a length-width ratio of L/W=4.4.

圖4 (網(wǎng)刊彩色)改變電場(chǎng)強(qiáng)度時(shí)的角度概率分布 (a)φN-E;(b)φL-EFig.4.(color online)The angle probability distributions at di ff erent electric strengths:(a)φN-E;(b)φL-E.

圖5 氧原子沿石墨烯法向的概率密度分布Fig.5.The probability density distribution of the oxygen atoms along the normal direction of graphene.

在第一水合層中水分子的偶極矩矢量與電場(chǎng)方向及石墨烯法向的夾角φD-E和φD-N的概率密度分布如圖6所示,此處的概率分布為既對(duì)空間也對(duì)時(shí)間進(jìn)行統(tǒng)計(jì)平均的結(jié)果.其中,偶極矩矢量方向定義為D,它與電場(chǎng)方向的夾角φD-E在0°附近的概率密度取為最大值(圖6(a)),這是因?yàn)樗肿釉陔妶?chǎng)下產(chǎn)生規(guī)則性排布,即偶極矩方向與電場(chǎng)方向平行,且隨著電場(chǎng)強(qiáng)度的增大,概率分布在0°附近的峰值增大,表明水分子在強(qiáng)電場(chǎng)下的響應(yīng)增強(qiáng).而在圖6(b)中,偶極矩矢量D與石墨烯法向的夾角φD-N在90°時(shí),概率密度取為最大值,表明絕大部分水分子與石墨烯表面平行.這一現(xiàn)象在未施加電場(chǎng)(E=0 V/nm)時(shí)已經(jīng)存在,且隨著電場(chǎng)強(qiáng)度的增大,φD-N在90°附近的概率密度峰值隨之增大.其主要物理機(jī)理在于水分子間存在氫鍵相互作用,當(dāng)在水分子中加入石墨烯納米顆粒時(shí),大量氫鍵遭到破壞,造成系統(tǒng)焓的增加[24].因而為了保持較大的氫鍵數(shù)目,水分子應(yīng)和納米顆粒表面平行,而形成覆蓋納米顆粒表面的水分子層,這即為水合作用.因而,在以上水合作用和直流電場(chǎng)的綜合影響下,水分子在石墨烯顆粒表面排布如圖7所示,石墨烯便會(huì)趨向于和電場(chǎng)方向平行,且如圖4(a)所示定向性隨電場(chǎng)強(qiáng)度增大而增強(qiáng).

保持電場(chǎng)強(qiáng)度E=1.0 V/nm不變,改變石墨烯長(zhǎng)寬比L/W,如圖8所示.與改變電場(chǎng)強(qiáng)度時(shí)的規(guī)律相似,φN-E在90°附近的定向分布仍然非常明顯(圖8(a)),且隨著L/W減小,峰值減小,定向性減弱;而φL-E隨著L/W減小(圖8(b)),峰值減小,角度分布的幅度變寬,且峰值處所對(duì)應(yīng)的角度也不惟一,當(dāng)L/W=4.4,2.1時(shí),峰值出現(xiàn)在5°—20°附近,而當(dāng)L/W=1.,峰值卻出現(xiàn)在45°附近.實(shí)際上,當(dāng)L/W=1.0時(shí),石墨烯為正方形,此時(shí)“長(zhǎng)邊”與電場(chǎng)強(qiáng)度夾角φL-E的概率分布出現(xiàn)在45°附近的對(duì)稱分布是較為合理的.

圖6 (網(wǎng)刊彩色)改變電場(chǎng)強(qiáng)度時(shí),偶極矩矢量與以下方向夾角的角度概率分布 (a)電場(chǎng)方向;(b)石墨烯法向Fig.6.(color online)The probability distributions of the angles between the water dipole moment and the following directions at di ff erent electric strengths:(a)Electric direction;(b)normal direction of graphene.

圖7 (網(wǎng)刊彩色)水合作用和直流電場(chǎng)綜合作用下水分子在石墨烯表面排布示意圖Fig.7.(color online)The schematic diagram of the water orientation distribution on the graphene surface under the combined e ff ect of the hydration and electric fi eld.

圖8 (網(wǎng)刊彩色)改變長(zhǎng)寬比時(shí)的角度概率分布 (a)φN-E;(b)φL-EFig.8.(color online)The angle probability distributions at di ff erent length-width ratios of graphene:(a)φN-E;(b)φL-E.

圖9 (網(wǎng)刊彩色)改變石墨烯長(zhǎng)寬比時(shí),偶極矩矢量與以下方向夾角的角度概率分布 (a)電場(chǎng)方向;(b)石墨烯法向Fig.9.(color online)The probability distributions of the angles between the water dipole moment and the following directions at di ff erent length-width ratios of graphene:(a)Electric direction;(b)normal direction of graphene.

圖9顯示當(dāng)石墨烯長(zhǎng)寬比變化時(shí),第一水合層中水分子的偶極矩矢量與電場(chǎng)方向及石墨烯法向的夾角φD-E和φD-N的概率密度分布.可以發(fā)現(xiàn)圖中不同長(zhǎng)寬比下的概率分布曲線基本重合,φD-E和φD-N的概率分布都基本不隨長(zhǎng)寬比的變化而變化,這表明以上水合作用以及水分子在電場(chǎng)下的規(guī)則性排布與顆粒的形狀尺寸無(wú)明顯依賴關(guān)系.

在我們之前的工作中[15,25],基于統(tǒng)計(jì)理論推導(dǎo)獲得了計(jì)算非球形顆粒轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)Dr的Einstein關(guān)系式,可表示為散系數(shù)Dr_L及Dr_N,它們定量表征了顆粒繞上述轉(zhuǎn)動(dòng)軸布朗轉(zhuǎn)動(dòng)的強(qiáng)弱.石墨烯在電場(chǎng)下的轉(zhuǎn)動(dòng)行為,可由方向矢量L和N的取向所決定,當(dāng)選定電場(chǎng)強(qiáng)度不變時(shí),L和N的定向性應(yīng)該與布朗轉(zhuǎn)動(dòng)的強(qiáng)弱有一定的定量關(guān)系.

為了確定此定量關(guān)系,下面利用定向度張量(7)式獲得描述石墨烯在電場(chǎng)下定向性強(qiáng)弱的定量參數(shù)定向度[26]:

其中,N是所有瞬時(shí)的石墨烯數(shù)量,l(n)是單位軸向矢量,下標(biāo)i和j在x,y和z三個(gè)方向之間選取,因而T是一個(gè)3×3的對(duì)稱矩陣.該張量的最大特征值即為定向度λ,當(dāng)λ=1.時(shí),此時(shí)定向性最強(qiáng),所有瞬時(shí)的N或L皆沿某一方向定向,此方向即為該矩陣的特征向量,對(duì)應(yīng)某個(gè)三維定向方向.定向度λ越小,表明定向性越弱.

其中,〈···〉代表系綜平均,Δ?為以方向矢量L或N 作為轉(zhuǎn)動(dòng)軸而進(jìn)行轉(zhuǎn)動(dòng)的角位移,t為時(shí)間.在平衡態(tài)下,統(tǒng)計(jì)石墨烯顆粒分別繞L和N方向轉(zhuǎn)動(dòng)的均方角位移(MSD)〈Δ?2〉隨時(shí)間的變化,如圖10所示.通過(guò)(6)式可分別計(jì)算獲得三維轉(zhuǎn)動(dòng)擴(kuò)

圖10 (網(wǎng)刊彩色)平衡態(tài)下,L/W=4.4,以L和N為轉(zhuǎn)動(dòng)軸轉(zhuǎn)動(dòng)的石墨烯MSD〈Δ?2〉隨時(shí)間的變化Fig.10. (color online)At equilibrium,and when L/W=4.4,the time-varying mean squared displacements〈Δ?2〉of graphene around the rotation axis L,N.

對(duì)于圖8中的情況,電場(chǎng)強(qiáng)度保持為E=1.0 V/nm,石墨烯平面接近平行于電場(chǎng)方向,且長(zhǎng)邊方向與電場(chǎng)方向夾角較小,集中分布在5°—45°之間.此時(shí)法向方向N取向的改變可近似看成由繞石墨烯長(zhǎng)邊L的回轉(zhuǎn)運(yùn)動(dòng)所導(dǎo)致(圖11(a)),因而討論N定向度λN和Dr_L間的關(guān)系;而長(zhǎng)邊方向N 取向的改變可近似看成由繞石墨烯法向N的平面轉(zhuǎn)動(dòng)所導(dǎo)致(圖11(b)),因而討論L定向度λL和Dr_N間的關(guān)系.

圖11 (a)石墨烯繞長(zhǎng)邊方向L回轉(zhuǎn)運(yùn)動(dòng)示意圖;(b)石墨烯繞法向方向N平面轉(zhuǎn)動(dòng)示意圖Fig.11.(a)Schematic diagram of the graphene rotation around the long-side vector L;(b)schematic diagram of the graphene rotation around the normal vector N.

圖12所示是法向矢量N定向度λN隨Dr_L以及長(zhǎng)邊矢量L定向度λL隨Dr_N的變化關(guān)系.以上兩組定向度和轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)λN-Dr_L,λLDr_N之間的定量關(guān)系皆為負(fù)相關(guān)關(guān)系,即轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)越大,相對(duì)應(yīng)的定向度越小,這與我們之前發(fā)現(xiàn)的在線性剪切流場(chǎng)作用下棒狀顆粒碳納米管的定向規(guī)律相一致[16](不改變碳納米管橫截面積時(shí)).在此,我們通過(guò)將轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)和定向度建立定量關(guān)系,對(duì)外場(chǎng)不變、顆粒形狀變化時(shí)的二維片狀石墨烯的定向規(guī)律成因給出了初步分析.

同時(shí),以上結(jié)果也對(duì)實(shí)際應(yīng)用中利用直流電場(chǎng)調(diào)控顆粒的定向具有指導(dǎo)意義:當(dāng)前模擬中石墨烯尺寸較小,在納米量級(jí);而實(shí)際中石墨烯、碳納米管等顆粒尺寸可達(dá)微米量級(jí).隨著顆粒尺寸增大,相應(yīng)的轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)減小,由于前述的負(fù)相關(guān)關(guān)系,定向度增大.因而,使用更小的直流電場(chǎng)強(qiáng)度就可以觀察到和模擬中相似的定向結(jié)果.例如,對(duì)于管長(zhǎng)在102μm量級(jí)的碳納米管,可在104V/m(遠(yuǎn)小于模擬中的109V/m)的直流電場(chǎng)強(qiáng)度下觀察到顯著的定向行為[12].

圖12 (網(wǎng)刊彩色)法向矢量N定向度λN隨Dr_L的變化,長(zhǎng)邊矢量L定向度λL隨Dr_N的變化Fig.12.(color online)The variation of the orientation order of the normal vector N with Drmr_Land the variation of the orientation order of the long-side vector L with Dr_N.

4結(jié) 論

利用分子動(dòng)力學(xué)模擬發(fā)現(xiàn)不帶電且未考慮極化效應(yīng)的矩形石墨烯在水中受到直流電場(chǎng)調(diào)控可出現(xiàn)定向行為,其中石墨烯平面趨向平行于電場(chǎng)方向,而長(zhǎng)邊方向與電場(chǎng)方向的定向夾角隨電場(chǎng)強(qiáng)度及長(zhǎng)寬比的變化取值分布在5°—50°范圍內(nèi).石墨烯平面平行于電場(chǎng)方向且隨著電場(chǎng)強(qiáng)度增大而定向性增強(qiáng)的主要原因在于第一水合層中水分子偶極矩在電場(chǎng)下出現(xiàn)響應(yīng),并趨向平行于電場(chǎng)方向,且由于水合作用偶極矩趨向平行于石墨烯平面,以上偶極矩的取向性會(huì)隨著電場(chǎng)強(qiáng)度的增大而增強(qiáng),水分子通過(guò)和石墨烯的相互作用而影響石墨烯的定向行為.減小石墨烯長(zhǎng)寬比時(shí),石墨烯法向和長(zhǎng)邊定向性減弱,利用定向度張量得到定量描述定向性的參數(shù)定向度,并使用Einstein關(guān)系式計(jì)算獲得平衡態(tài)下石墨烯繞長(zhǎng)邊方向和法向方向轉(zhuǎn)動(dòng)的三維轉(zhuǎn)動(dòng)擴(kuò)散系數(shù),發(fā)現(xiàn)石墨烯法向、長(zhǎng)邊定向度分別與繞長(zhǎng)邊和法向的轉(zhuǎn)動(dòng)擴(kuò)散系數(shù)存在負(fù)相關(guān)關(guān)系.片狀顆粒石墨烯在直流電場(chǎng)下的定向行為本質(zhì)上來(lái)源于外場(chǎng)和布朗轉(zhuǎn)動(dòng)的共同作用,增大外場(chǎng)強(qiáng)度或減小轉(zhuǎn)動(dòng)擴(kuò)散系數(shù),皆可有效增強(qiáng)定向性.

[1]Huang H,Liu C H,Wu Y,Fan S S 2005 Adv.Mater.17 1652

[2]Liang Q,Yao X,Wang W,Liu Y,Wong C P 2011 ACS Nano 5 2392

[3]Behabtu N,Young C C,Tsentalovich D E,Kleinerman O,Wang X,Ma A W K,Bengio E A,ter Waarbeek R F,de Jong J J,Hoogerwerf R E,Fairchild F B,Ferguson J B,Maruyama B,Kono J,Talmon Y,Cohen Y,Otto M J,Pasquali M 2013 Science 339 182

[4]Liu M S,Lin M C C,Huang I T,Wang C C 2005 Int.Commun.Heat Mass Trans.32 1202

[5]Guo X,Su J,Guo H 2012 Soft Matter 8 1010

[6]Hermanson K D,Lumsdon S O,Williams J P,Kaler E W,Velev O D 2001 Science 294 1082

[7]Hsu H Y,Sharma N,Ruo ffR S,Patankar N A 2005 Nanotechnology 16 312

[8]Alexandrou I,Ang D K H,Mathur N D,Haq S,Amaratunga G A J 2004 Nano Lett.4 2299

[9]van der Zande B M I,Koper G J M,Lekkerkerker H N W 1999 J.Phys.Chem.B 103 5754

[10]Ma C,Zhang W,Zhu Y,Ji L,Zhang R,Koratkar N,Liang J 2008 Carbon 46 706

[11]Li J,Zhang Q,Peng N,Zhu Q 2005 Appl.Phys.Lett.86 153116

[12]Martin C A,Sandler J K W,Winder A H,Schwarz M K,Bauhofer W,Schulte K,Sha ff er M S P 2005 Polymer 46 877

[13]Oliveira L,Saini D,Gaillard J B,Podila R,Rao A M,Serkiz S M 2015 Carbon 93 32

[14]Daub C D,Bratko D,Ali T,Luzar A 2009 Phys.Rev.Lett.103 207801

[15]Cao B Y,Dong R Y 2014 J.Chem.Phys.140 34703

[16]Dong R Y,Cao B Y 2014 Sci.Rep.4 6120

[17]Song Y,Dai L L 2010 Mol.Simulat.36 560

[18]Ryckaert J P,Cicotti G,Berendsen H J C 1977 J.Comput.Phys.23 327

[19]Won C Y,Joseph S,Aluru N R 2006 J.Chem.Phys.125 114701

[20]Werder T,Walther J H,Ja ff e R L,Halicioglu T,Noca F,Koumoutsakos P 2001 Nano Lett.1 697

[21]Shiomi J,Maruyama S 2009 Nanotechnology 20 055708

[22]Plimpton S 1995 J.Comput.Phys.7 1

[23]Hockney R W,Eastwood J W 1988 Computer Simulation Using Particles(New York:Taylor&Francis Group)pp267–304

[24]Djikaev Y S,Ruckenstein E 2012 J.Phys.Chem.B 116 2820

[25]Dong R Y,Cao B Y 2015 J.Nanosci.Nanotechnol.15 2984

[26]B?rzs?nyi T,Szabó B,T?r?s G,Wegner S,T?r?k J,Somfai E,Bien T,Stannarius R 2012 Phys.Rev.Lett.108 228302

PACS:47.11.Mn,81.05.ue,66.10.cg,87.16.dpDOI:10.7498/aps.66.014702

*Project supported by the National Natural Science Foundation of China(Grant No.51322603)and Institute of Telecommunication Satellite Innovation Fund,China.
?Corresponding author.E-mail:caoby@tsinghua.edu.cn

DC electric fi eld induced orientation of a graphene in water?

Dong Ruo-Yu1)Cao Peng2)Cao Gui-Xing2)Hu Guo-Jie2)Cao Bing-Yang1)?

1)(Key Laboratory for Thermal Science and Power Engineering of Ministry of Education Department of Engineering Mechanics Tsinghua University,Beijing 100084,China)2)(Institute of Telecommunication Satellite China Academy of Space Technology,Beijing 100094,China)(Received 26 July 2016;revised manuscript received 13 October 2016)

Graphene,as a classical two-dimensional material,has various excellent physical properties,which can be further transferred into its nanocomposite.Under external fi elds,the nonspherical nanoparticles in liquid environment will exhibit various deterministic movements,among them is the orientation behavior.By realizing the orientation control of nanoparticles,we can,on one hand,increase the thermal conductivity of the system along the oriented direction,and on the other hand,fabricate novel nano-devices based on the nanoscale self-assembly,which may become the key components in NEMS and Lab-on-a-chip architectures.However,current studies mainly focus on the orientations of one-dimensional rod-shaped particles,like carbon nanotubes.For a two-dimensional nanoparticle,like graphene,the situation is more complex than the one-dimensional one,because two unit vectors should be de fi ned to monitor the orientation behaviors.As far as we know,this part of research has not been extensively carried out.Thus,in this paper,the molecular dynamics method is used to study the orientation of a single uncharged rectangular graphene in water,induced by DC electric fi elds.We track the orientations of the normal and long-side vectors of graphene.The results show that at a relatively high electric strength of 1.0 V/nm,the graphene is preferred to orient its normal vector perpendicular and its long-side vector with a small angle(located between 0?and 30?)with respect to the electric direction,respectively.With the increase of the electric fi eld strength,the orientation preference of the normal vector along the electric direction is increased.To explain this phenomenon,we calculate the orientation distribution of water molecules in the fi rst hydration shell.The dipoles tend to be parallel to the electric direction,and the surfaces of water molecules tend to be parallel to the surface of graphene.These two combined e ff ects result in the above orientation behavior of the normal vector.Another interesting phenomenon is that the decrease of the length to width ratio of graphene will cause both the orientation preferences of the normal vector and the long-side vector to decrease.By utilizing the Einstein relation,we can obtain the rotational di ff usion coefficients of graphene around the normal vector and long-side vector.The qualitative results show that the orientation orders of the normal vector and long-side vector respectively have negative correlations with the rotational di ff usion coefficients of the rotation around the long-side vector and the normal vector.The orientation behavior of the platelike graphene actually comes from the competing e ff ects between its rotational Brownian motion and the external fi eld.Increasing the strength of the external fi eld or reducing the rotational di ff usivity will both lead to an increased orientation order of the nonspherical nanoparticle.

graphene,orientation,DC electric fi eld,molecular dynamics

10.7498/aps.66.014702

?國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào):51322603)和通信衛(wèi)星系統(tǒng)創(chuàng)新基金資助的課題.

?通信作者.E-mail:caoby@tsinghua.edu.cn