史若宇 王林鋒 高磊 宋愛生 劉艷敏 胡元中 馬天寶?

1)(清華大學,摩擦學國家重點實驗室,北京 100084)

2)(南京航空航天大學航天學院仿生結(jié)構(gòu)與材料防護研究所,南京 210016)

3)(北京科技大學腐蝕與防護中心,北京 100083)

基于滑動勢能面的二維材料原子尺度摩擦行為的量化計算?

史若宇1)王林鋒2)高磊3)宋愛生1)劉艷敏1)胡元中1)馬天寶1)?

1)(清華大學,摩擦學國家重點實驗室,北京 100084)

2)(南京航空航天大學航天學院仿生結(jié)構(gòu)與材料防護研究所,南京 210016)

3)(北京科技大學腐蝕與防護中心,北京 100083)

近年來,二維材料優(yōu)異的摩擦特性成為人們關(guān)注的焦點,然而目前缺乏理論上對其摩擦力進行快速、有效、精確的計算預測方法.本文提出采用密度泛函理論計算真實體系的滑動勢能面,利用得到的“數(shù)值型勢能面”替代傳統(tǒng)的解析勢函數(shù),并結(jié)合Prandtl-Tomlinson模型,量化求解具有復雜形狀勢能面的真實二維材料體系的摩擦行為.基于該方法,揭示了原子力顯微鏡實驗中觀察到的石墨烯Moiré紋超晶格結(jié)構(gòu)的雙周期“黏-滑”摩擦現(xiàn)象;理論預測了二維材料異質(zhì)結(jié)構(gòu)的層間超低摩擦現(xiàn)象,相對于同質(zhì)材料,其靜摩擦力和滑動摩擦力均成數(shù)量級降低,發(fā)現(xiàn)勢能面起伏和驅(qū)動彈簧剛度均會影響層間相對滑動路徑,進而對層間的摩擦行為產(chǎn)生影響.該方法同樣可拓展到其他van der Waals作用主導的界面摩擦體系.

原子級摩擦,Prandtl-Tomlinson模型,二維材料,層間滑動

1 引 言

隨著科學研究尺度的微觀化和技術(shù)應用的納米化,原子尺度的摩擦行為一直受到科研人員的廣泛關(guān)注.原子力顯微鏡(atomic force microscopy,AFM)是研究原子級摩擦行為的重要手段,探針基座通過懸臂梁帶動探針針尖在樣品表面滑動,通過測量和記錄懸臂梁的扭轉(zhuǎn)變形得到針尖和樣品之間的摩擦力.Prandtl-Tomlinson模型［1,2］(PT模型)可以在理論上很好地描述AFM的實驗過程,將探針基座、懸臂梁、針尖分別簡化成驅(qū)動塊、彈簧以及振子,將樣品表面原子的勢能起伏簡化成三角函數(shù)的形式,那么這一過程就可以被抽象成一個可以解析的動力學過程即PT模型:一個勻速運動的驅(qū)動塊通過彈簧帶動一個振子在三角函數(shù)的勢能面上滑動,振子感受到的勢能包括兩個部分,一是振子受到基底作用的勢能,二是彈簧中的彈性勢能.對PT模型建立動力學方程并求解,可以分析體系的摩擦行為,即在給定的勢能面、彈簧剛度、溫度等條件下,能夠獲得每一時刻振子的位置、驅(qū)動塊所受到的力以及系統(tǒng)各部分的能量等詳細信息,從而為人們研究原子級摩擦行為和能量耗散機制提供指導［3,4］.

近年來,二維材料呈現(xiàn)出許多新奇而優(yōu)異的摩擦特性,例如重構(gòu)表面或二維材料異質(zhì)結(jié)構(gòu)會出現(xiàn)摩擦力的調(diào)制現(xiàn)象,顯著影響體系的摩擦特性和能量耗散［5?7］,特別是二維材料層間的超滑狀態(tài)［8,9］,使其成為納米摩擦學研究中的熱點.二維異質(zhì)材料或二維材料層間的滑動勢能面往往難以用簡單的解析型函數(shù)(如三角函數(shù))來描述,因而限制了人們采用經(jīng)典PT模型對其摩擦特性進行理論計算.另一方面,分子動力學方法可以直接模擬原子尺度的材料摩擦行為,特別是在石墨烯體系中得到了廣泛應用.然而分子動力學方法依賴于層間作用勢的合理選取,對于實際中大量二維材料體系尚缺乏有效的經(jīng)驗勢函數(shù)描述其層間相互作用.本文采用由密度泛函理論(DFT-D)計算的“數(shù)值型勢能面”替代簡單的解析勢函數(shù),實現(xiàn)對二維異質(zhì)材料表面和層間摩擦行為的求解.利用DFT-D計算可以精確得到任意體系滑動勢能面離散化的數(shù)據(jù),將其處理為可用于PT模型計算的“數(shù)值型勢能面”,從而可以計算其摩擦行為.我們采用這一方法,研究了單原子在石墨烯和金屬釕構(gòu)成的異質(zhì)結(jié)構(gòu)表面的摩擦行為,發(fā)現(xiàn)了與AFM實驗中相近的雙周期黏滑現(xiàn)象.利用第一性原理得到的二維材料層間滑動勢能面,研究了二維材料層間的摩擦行為,發(fā)現(xiàn)二維異質(zhì)材料的層間靜摩擦力和動摩擦力遠小于同質(zhì)材料之間的摩擦力,詳細分析了石墨烯層間滑動時的滑移路徑以及剛度對摩擦力的影響規(guī)律.這一基于數(shù)值型勢能面的差分法求解PT模型的研究手段,可以適用于各種二維材料,以及層間以van der Waals作用為主的界面摩擦問題,同時有助于深層次揭示原子尺度摩擦和能量耗散的物理本質(zhì).

2 方 法

2.1 PT模型計算

PT模型中,摩擦副的上表面簡化為速度為v的驅(qū)動塊,通過剛度為k的彈簧帶動質(zhì)量為m的振子,下表面一般簡化為三角函數(shù)形式的勢能面于是滑動的過程中體系的總勢能為

其中,x為t時刻振子的位置.于是,振子的動力學過程可以由朗之萬(Langevin)方程表達:

式中m為振子質(zhì)量,μ為系統(tǒng)的阻尼系數(shù),ξ(t)為溫度項.基于MATLAB采用四階的龍格-庫塔方法對方程進行求解,由于本文不討論溫度的影響,求解過程中將溫度設為0 K,ξ(t)=0.對方程(2),求解每一時刻的振子位置x,進一步可以得到每一時刻的摩擦力F,即驅(qū)動塊受到的彈簧力:

但在很多情況下,實際的勢能面不是簡單的三角函數(shù)形式,甚至難以準確寫出其數(shù)學表達式,也就不能進行計算.為了突破這一方法的局限性,等間隔地建立勢能-位置的“數(shù)值型勢能面”:U(x)=(U1,U2,···Un),這一勢能數(shù)據(jù)庫精度足夠高(間隔?x足夠小)時,差分可以得到勢能導數(shù)-位置的數(shù)據(jù)庫:

每一步求解過程中,利用上一步得到的振子位置在數(shù)據(jù)庫中讀取勢能導數(shù)的數(shù)值,代入本步進行求解.本文將這種解法稱為差分法,而將以往所用的方法稱為微分法.

圖1 (網(wǎng)刊彩色)(a)勢能面、勢能庫和勢能導數(shù)庫示意圖;(b)差分法得到的黏滑曲線;(c)微分法得到的黏滑曲線Fig.1.(color online)(a)Schematics of potential surface,potential data base and potential derivative data base;(b)friction obtained via fi nite di ff erence method;(c)friction obtained via di ff erential method.

經(jīng)驗證,采用差分法和微分法在相同的時間步長下得到的結(jié)果基本相同.采用一個正弦函數(shù)形式的勢能面進行測試,勢能面的周期為3 ?,幅值為5 eV,分別用差分法和微分法在同樣的條件下(彈簧剛度k=30 N/m,溫度T=0 K,阻尼系數(shù)0.6,驅(qū)動塊速度v=10?5m/s,振子質(zhì)量m=10?12kg,對于差分法勢能庫格子取0.01 ?)進行計算,得到了如圖1(b)和圖1(c)所示的黏滑曲線,對比發(fā)現(xiàn)兩者平均差距僅0.12%.

差分法也可以應用于二維的勢能面,二維的動力學方程為:

此時勢能導數(shù)也要擴展至二維:

振子受到兩個方向上的力:

其中kx,ky代表彈簧在x,y兩個方向上的剛度,在AFM實驗中分別是垂直探針懸臂長軸的剛度和沿著長軸方向的剛度.AFM的摩擦實驗中針尖一般沿著垂直探針懸臂長軸的方向滑動,因此將探針懸臂受到的扭轉(zhuǎn)力作為摩擦力,本文二維勢能面的求解中,也采用驅(qū)動塊在x方向上的彈簧力Fx作為摩擦力.

2.2 DFT計算

PT模型中所需的數(shù)值型勢能面由DFT-D計算得到.所用方法為具有長程色散作用校正的密度泛函理論計算方法［10?12］,交換關(guān)聯(lián)泛函采用廣義梯度近似中的Perdew-Burke-Ernzerhof泛函［13］.價電子和離子實之間的作用采用投影綴加波方法［14,15］描述,截斷能為400 eV.總能的收斂精度為10?4eV.計算優(yōu)化并記錄兩層材料處于不同的相對位置時的體系勢能,得到體系的滑動勢能面.對勢能面插值擴展之后可用于PT模型差分求解摩擦行為.

3 結(jié)果與討論

3.1 針尖在異質(zhì)結(jié)構(gòu)表面的摩擦行為

兩種不同的二維材料之間若具有晶格失配或者相對轉(zhuǎn)角,疊在一起會形成一種稱為Moiré紋(Moiré pattern)的周期結(jié)構(gòu)［16］. 事實上,Moiré紋是一種廣泛存在的超晶格結(jié)構(gòu),不僅限于二維材料之間.最近Zhu等［17］模擬發(fā)現(xiàn),表面電荷可以改變水分子的偶極取向,從而調(diào)控水的結(jié)構(gòu)相變.研究表明,水分子在冰相下呈現(xiàn)六方結(jié)構(gòu),并可與石墨烯形成公度或者不公度結(jié)構(gòu),并形成Moiré紋結(jié)構(gòu).這樣,在異質(zhì)結(jié)構(gòu)表面就形成了兩種周期性勢能面的疊加,一種是原子晶格周期的勢能面,另一種是Moiré紋周期的勢能面.由于勢能面比較復雜,這種結(jié)構(gòu)會呈現(xiàn)出一些新奇的摩擦行為［6,7］.石墨烯Gr晶格周期為aa=2.46 ?,金屬Ru晶格周期為2.71 ?,Gr和Ru構(gòu)成異質(zhì)體系便會出現(xiàn)Moiré紋［18,19］,并且AFM實驗中發(fā)現(xiàn)了雙周期“黏-滑”摩擦現(xiàn)象,如圖2(a)所示［20］.Gr蜂窩狀的原子結(jié)構(gòu)使得其勢能面分布也是蜂窩狀,再疊加上一個Moiré紋周期的勢能面,這種二維的勢能面極為復雜,很難準確寫出表面勢能的數(shù)學表達式,采用微分法計算比較困難,于是我們用差分法對這一體系的摩擦行為進行研究.

首先通過第一性原理軟件(VASP)建立并優(yōu)化Gr和Ru的異質(zhì)結(jié)構(gòu),二者構(gòu)成的Moiré紋周期為am=3 nm［21］.導出表層Gr的原子位置和Moiré紋中心的位置.然后以每個碳原子位置為中心生成一個高度為Ua0、半徑為aa的原子級勢壘:

以每個Moiré紋中心位置生成一個高度為Um0半徑為am的勢壘:

其中Ua0,Um0根據(jù)實驗中Fmax確定,這樣就構(gòu)成了整個異質(zhì)結(jié)構(gòu)體系的勢能庫,如圖2(d)所示.

這里勢能庫x和y方向的格子均為0.05 ?,Ua0=16 eV,Um0=120 eV,aa=0.246 nm,am=3 nm,x和y方向的彈簧剛度分別為kx=35 N/m和ky=150 N/m,溫度為0 K.得到的摩擦力圖如圖2(e)所示,從中可以看出原子晶格的周期性特征和Moiré紋大小的周期性特征,從圖2(g)的側(cè)向力曲線中也可以看到原子晶格周期和Moiré紋周期的兩種黏滑特征.Gr和Ru之間的作用力比較強,二者所形成的異質(zhì)結(jié)構(gòu)表面Moiré紋級的勢壘較大,會發(fā)生雙周期的黏滑.除此之外,Gr和鍺Ge［6］、鉑Pt形成的異質(zhì)結(jié)構(gòu)［22］,KBr和NaCl形成的異質(zhì)結(jié)構(gòu)［7］,以及表面重構(gòu)的金Au［5］,表面都會出現(xiàn)或強或弱的摩擦力調(diào)制現(xiàn)象,均可用差分法PT模型進行摩擦行為的分析.

圖2 (網(wǎng)刊彩色)(a)AFM實驗中的黏滑曲線;(b)只有原子晶格勢壘的勢能面;(c)只有Moiré紋勢壘的勢能面;(d)異質(zhì)結(jié)構(gòu)的勢能面,由(b)和(c)疊加而成;(e)由異質(zhì)結(jié)構(gòu)的勢能面得到的摩擦力圖;(f)圖(d)中黑色虛線上的勢能分布;(g)圖(e)中黑色虛線上的黏滑曲線［20］Fig.2.(color online)(a)Experimental Friction loop;(b)potential surface with atomic-level barrier;(c)potential surface with superlattice-level barrier;(d)potential surface with both atomic-level and superlattice-level barriers;(e)friction map corresponding to the potential surface in(d);(f)potential energy pro fi le;(g)calculated friction loop corresponding to(f)［20］.

3.2 二維材料層間的摩擦力

兩層二維材料在不公度接觸的情況下,每個原子受到來自不同方向不同大小的力相互抵消,層間滑動勢壘很小,有可能出現(xiàn)超滑行為［23］.這里,我們利用第一性原理計算了二維材料的層間滑動勢能面［9］,并輸入差分法PT模型,計算各個體系的摩擦力,研究二維材料層間的摩擦行為.

3.2.1 不同二維材料之間的摩擦行為

圖3(a)—(e)是第一性原理計算得到的各個材料層間滑動勢能面,相對于二硫化鉬(MoS2)、石墨烯(Gr)、氟化石墨烯(FG)和各自同質(zhì)的材料層間滑動勢能面,FG和MoS2以及Gr和MoS2之間因不公度性的存在,勢壘很小.

基于上述的滑動勢能面,我們通過差分法計算了在這幾種的體系下各自的摩擦力圖,如圖3(f)—(j)所示,類似AFM實驗中摩擦力圖的做法,將勢能面等分為若干條水平線,沿著水平線從左向右滑動,得到這些水平線上的摩擦力,進而將各條線組合成一幅摩擦力圖.從圖3中可以看出,在同質(zhì)結(jié)構(gòu)的體系里,MoS2/MoS2體系的勢能面起伏最大,相應的摩擦力也較大,而FG/FG體系的勢壘和摩擦力相對較小.而對比同質(zhì)的體系,異質(zhì)材料之間摩擦力則小很多.圖4中定量地對比了各個體系的勢壘高度和摩擦力大小,FG/MoS2和Gr/MoS2的摩擦力量級在10?4nN,與同質(zhì)體系相比可以忽略不計.由于這兩種異質(zhì)體系能壘非常小,二者進入了一種稱為連續(xù)滑動的狀態(tài)［5］,這時兩者之間不發(fā)生黏滑失穩(wěn),去程和回程的摩擦力曲線基本重合,因此摩擦力和能量耗散均極小,發(fā)生超滑.而另一個值得注意的是,盡管MoS2/MoS2體系勢壘比Gr/Gr體系大,但MoS2/MoS2體系周期較大,使得兩種勢能面梯度相近,在PT模型中,最大摩擦力和勢能面梯度直接相關(guān):Fmax=πU/a,因此二者摩擦力基本相同.

圖4 (網(wǎng)刊彩色)MoS2/MoS2,Gr/Gr,FG/FG,FG/MoS2和Gr/MoS2體系的勢壘和摩擦力對比Fig.4.(color online)Potential barrier and friction of MoS2/MoS2,Gr/Gr,FG/FG,FG/MoS2and Gr/MoS2.

3.2.2 摩擦過程中振子的運動路徑

與一維不同的是,本文涉及x,y兩個方向上的運動自由度,因此PT模型描述的驅(qū)動塊在沿著x方向運動時,振子不僅在x方向上運動,在y方向上也會和驅(qū)動塊的掃描路徑發(fā)生偏移,因此振子的運動路徑將是一條在平面上尋求最低總能量的曲線.對于層間滑動過程而言,二維PT模型中的振子軌跡代表了層間滑動時實際運動的相對路徑,而驅(qū)動塊路徑是指外部驅(qū)動力的方向.為了對此做進一步的探討,我們以Gr/Gr體系為例,計算了不同掃描路徑下振子的運動軌跡.

圖5(a)展示了驅(qū)動塊沿著三條不同的路徑運動時,振子每一時刻在勢能面上的位置.路徑1○是沿著晶格的密排方向跨過勢能面最高點的直線,路徑2○是沿著勢能面較低位置的直線,路徑3○是介于1○,2○之間的直線,經(jīng)過勢壘的半山腰.圖中虛直線代表驅(qū)動塊的運動軌跡.每隔固定的時間間隔,用圓點記錄下振子的位置,圓點連線是振子的運動軌跡.由于這里的采樣密度比較大,很多圓點出現(xiàn)了重合,因此我們將每條路徑上的振子軌跡進行了不同程度的放大,分別列在圖5(c)—(e)中,路徑1○,2○,3○上的摩擦力如圖5(f)—(h)所示.

從記錄的振子路徑數(shù)據(jù)上可以看出,振子大多時間在勢能低點.在路徑1○,振子路徑雖然和驅(qū)動塊在y方向上有偏移,但基本是一條直線,由于驅(qū)動塊路徑距離勢能最低點較遠,在彈簧力的作用下,振子不能處于整個勢能面的最低點.此時,每一個黏滑周期中振子跳過一個晶格大小.而在路徑2○,由于整個勢能面最低點距離驅(qū)動塊路線較近,彈簧對振子作用較小,因此振子可以從一個勢阱跳躍到另一個勢阱,并且在一個黏滑周期中,振子跳過半個晶格大小.由于路徑2○上的勢能起伏相對較弱,此時的摩擦力也較小.對于路徑3○,振子軌跡也基本是一條直線,摩擦力曲線和路徑1○比較接近.

3.2.3 y方向剛度對振子路徑及摩擦行為的影響

振子的軌跡是勢能面和彈簧的彈性勢能共同作用的結(jié)果.這里涉及勢能面的形狀和彈簧的剛度等參數(shù),前文已經(jīng)討論了勢能面的形狀對振子位置和摩擦力的影響,而彈簧剛度可以分為x,y兩個方向,在AFM實驗中分別是垂直探針懸臂長軸的剛度和沿著長軸方向的剛度,不同的AFM探針有不同的x,y方向剛度,以往的文獻對x方向的彈簧剛度kx研究較多,這里借助二維PT模型,對y方向彈簧剛度ky和摩擦行為之間的關(guān)系進行進一步探討.

圖5 (網(wǎng)刊彩色)(a)Gr/Gr體系的滑動勢能面,箭頭和虛直線標記了驅(qū)動塊的運動路徑,圓點標記了振子的實際軌跡,相鄰圓點間等時間間隔;(b)圖(a)中勢能面對應的摩擦力圖;(c)—(e)圖(a)中三個振子路徑的放大,橫坐標為x方向的位置,y坐標為偏離各自驅(qū)動塊路徑的距離;(f)—(h)三條路徑各自對應的摩擦力曲線Fig.5.(color online)(a)Sliding potential surface of Gr/Gr,arrows and dashed lines indicate the direction of the driver,black dots indicate the trajectory of the oscillator,time intervals are constant between adjacent dots;(b)corresponding friction map;(c)–(e)trajectory after zooming in on(a),axis X indicates the driver’s position in x direction and axis Y indicates the oscillator’s shift in y direction away from the trace of the driver;(f)–(h)friction loop on di ff erent trajectories.

圖6 (網(wǎng)刊彩色)(a)—(d)在Gr/Gr體系滑動勢能面上分別標注了y方向剛度為200 N/m,50 N/m,10 N/m,1 N/m時驅(qū)動塊不同路徑下振子的軌跡;(e)—(h)y方向剛度分別為200 N/m,50 N/m,10 N/m,1 N/m時對應的摩擦力圖;(i)—(l)分別對應(a)—(d)路徑3○上的摩擦力曲線Fig.6.(color online)(a)–(d)Trajectory of the oscillator when spring constant on y direction is 200 N/m,50 N/m,10 N/m,1 N/m;(e)–(h)friction map corresponding to(a)–(d);(i)–(l)friction loop on trace3○.

圖6展示了ky不同時振子路徑的變化,圖6(a)—(d)分別是ky為200,50,10,1 N/m時振子的軌跡,而kx一直保持10 N/m.當ky很大時,隨著振子在y方向上的偏移,彈簧的彈性勢能積攢很快,總能量迅速增加,因此振子在y方向上的偏移很小,振子的軌跡近似為沿著驅(qū)動塊路徑的直線,并且不論驅(qū)動塊沿著哪條路徑運動,振子軌跡都是如此.ky的減小給振子y方向的偏移提供了條件,ky=50 N/m時,路徑2○對應的振子軌跡開始發(fā)生偏移,振子的停留位置大多在勢能最低點,而路徑1○,3○對應的振子軌跡仍然保持直線.隨著ky繼續(xù)減小,振子的偏移也逐漸增加,當ky非常小時,振子受到彈簧y方向的約束很弱,此時,彈簧積累的彈性勢能很小,下表面勢能面的影響占據(jù)主導作用,從而振子被迫尋找勢能面的最低點,因此即使在路徑1○和3○上,振子也會沿著勢能面上勢能最低的路徑運動,如圖6(d)所示.而由于y方向剛度影響了振子的運動軌跡,投影在x方向上的振子軌跡也會發(fā)生改變,從而影響了體系的摩擦力.以路徑3○對應的摩擦力曲線為例(圖6(i)—(l),在ky較大時,黏滑周期和勢能面的晶格周期相同,是比較規(guī)則的黏滑曲線.而ky=10 N/m時,由于在一個黏滑周期中,振子在更遠處的一個勢能低點有短暫的停留,因此黏滑曲線的突滑階段并沒有直接掉落至最低點,而是有一個微小的上升,之后才繼續(xù)降低,最終進入另一個黏著階段.當ky非常小時,無論驅(qū)動塊沿著哪個路徑運動,振子的軌跡基本相同,都是從一個勢阱跳躍至下一個臨近的勢阱,因此每條線上的摩擦力Fx差別不大,都是在同樣的位置發(fā)生黏著和突滑,也就出現(xiàn)了摩擦力圖條紋狀的特征.

4 結(jié) 論

利用DFT-D得到的“數(shù)值型勢能面”替代簡化的解析勢函數(shù),將PT模型計算摩擦行為的適用范圍拓展至任意勢能面,而非局限于規(guī)則的周期性簡單勢能函數(shù),可量化求解具有復雜形狀勢能面的真實二維材料體系的摩擦行為,給原子級摩擦的研究提供了一個快捷可靠的計算分析手段.以此為基礎(chǔ),探討了針尖和異質(zhì)材料表面之間的摩擦行為,揭示了實驗中的雙周期黏滑現(xiàn)象的機理.同時,詳細計算分析了不同二維材料層間的摩擦行為,發(fā)現(xiàn)異質(zhì)材料之間摩擦力可以成數(shù)量級降低,討論了層間滑動路徑的變化規(guī)律,指出y方向剛度和摩擦行為之間有密切的聯(lián)系,為二維材料間摩擦行為的預測提供了理論指導.上述通過差分法PT模型的計算和分析,為摩擦行為調(diào)控、超滑系統(tǒng)設計等提供了理論依據(jù).

本文計算部分在清華信息科學與技術(shù)國家實驗室“探索100”上進行.

［1］Tomlinson G A 1929Phil.Mag.7 905

［2］Zhao Y P 2012Surface and Interface Mechanics(Beijing:Science Press)p301(in Chinese)［趙亞溥2012表面與界面物理力學(北京:科學出版社)第301頁］

［3］Dong Y,Vadakkepatt A,Martini A 2011Tribol.Lett.44 367

［4］Wang Z,Ma T,Hu Y,Xu L,Wang H 2015Friction3 170

［5］Li Q,Dong Y,Martini A,Carpick R W 2011Tribol.Lett.43 369

［6］Zheng X,Gao L,Yao Q,Li Q,Zhang M,Xie X,Qiao S,Wang G,Ma T,Di Z,Luo J,Wang X 2016Nat.Commun.7 13204

［7］Maier S,Gnecco E,Barato ffA,Bennewitz R,Meyer E 2008Phys.Rev.B78 045432

［8］Liu S,Wang H,Xu Q,Ma T,Yu G,Zhang C,Geng D,Yu Z,Zhang S,Wang W,Hu Y,Wang H,Luo J 2017Nat.Commun.8 14029

［9］Wang L F,Ma T B,Hu Y Z,Zheng Q,Wang H,Luo J 2014Nanotechnology25 385701

［10］Grimme S 2006J.Comput.Chem.27 1787

［11］Kresse G,Hafner J 1994Phys.Rev.B49 14251

［12］Kresse G,Furthmüller J 1996Comput.Mater.Sci.6 15

［13］Perdew J P,Burke K,Ernzerhof M 1996Phys.Rev.Lett.77 3865

［14］Bl?chl P E 1994Phys.Rev.B50 17953

［15］Kresse G,Joubert D 1999Phys.Rev.B59 1758

［16］Geim A K,Grigorieva I V 2013Nature499 419

［17］Zhu X,Yuan Q,Zhao Y P 2014Nanoscale6 5432

［18］Pan Y,Zhang H,Shi D,Sun J,Du S,Liu F,Gao H J 2009Adv.Mater.21 2777

［19］Pan Y,Shi D,Gao H J 2007Chin.Phys.16 3151

［20］Shi R,Gao L,Lu H,Li Q,Ma T,Guo H,Du S,Feng X,Zhang S,Liu Y,Cheng P,Hu Y,Gao H,Luo J 20172D Mater.4 025079

［21］Gao L,Liu Y,Ma T,Shi R,Hu Y,Luo J 2016Appl.Phys.Lett.108 261601

［22］Filleter T,Bennewitz R 2010Phys.Rev.B81 155412

［23］Hirano M,Shinjo K 1990Phys.Rev.B41 11837

Quantitative calculation of atomic-scale frictional behavior of two-dimensional material based on sliding potential energy surface?

Shi Ruo-Yu1)Wang Lin-Feng2)Gao Lei3)Song Ai-Sheng1)Liu Yan-Min1)Hu Yuan-Zhong1)Ma Tian-Bao1)?
1)(State Key Laboratory of Tribology,Tsinghua University,Beijing 100084,China)
2)(Institute of Bio-inspired Structure and Surface Engineering,College of Astronautics,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China)
3)(Corrosion and Protection Center,Key Laboratory for Environmental Fracture(MOE,University of Science and Technology Beijing),Beijing 100083,China)

The excellent tribological characteristics of two-dimensional(2D)materials have

great attention,however,how to e ff ectively predict their frictions is still lacking.Here,we propose to obtain the sliding potential energy surface by density functional theory calculations,instead of simpli fi ed potential energy function.Thus it is able to solve the frictional behaviors of 2D materials with irregular complex potential energy surfaces.Firstly,we reveal the mechanism of dual-scale stick-slip behavior between a tip and a graphene/Ru(0001)heterostructure.With a dual-wavelength potential energy surface,we observe a similar frictional behavior to those captured in atomic force microscopy experiments,in which a signi fi cant long-range stick-slip sawtooth modulation emerges with a period coinciding with the Moiré superlattice structure.Secondly,we discuss the interlayer frictions of 2D materials,including graphene/graphene, fl uorinated graphene/ fl uorinated graphene,MoS2/MoS2,graphene/MoS2and fl uorinated graphene/MoS2.With sliding potential energy surface obtained by density functional theory calculations,the interlayer friction is estimated according to the Prandtl-Tomlinson model calculation method.Compared with the friction between homostructures,the friction between heterostructures is lowered by orders of magnitude,which could be attributed to its ultralow sliding potential barrier.The stick-slip instability could be observed in homostructure,while heterostructure exihibits smooth friction loops.The 2D sliding path between the layers is recorded in the sliding process,showing its dependence on both the potential energy barrier and the spring constant.The sliding path shift increases with the increase of potential energy barrier and the decrease of spring constant in the y direction.This method is also applicable to tribological systems with dominated interfacial van der Waals interaction.

atomic scale friction,Prandtl-Tomlinson model,two-dimensional materials,interlayer sliding

4 June 2017;revised manuscript received 11 July 2017)

(2017年6月4日收到;2017年7月11日收到修改稿)

10.7498/aps.66.196802

?國家自然科學基金(批準號:51422504,51505217,51527901)資助的課題.

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

?2017中國物理學會Chinese Physical Society

PACS:68.35.Af,81.05.ue

10.7498/aps.66.196802

*Project supported by the National Natural Science Foundation of China(Grant Nos.51422504,51505217,51527901).

?Corresponding author.E-mail:mtb@mail.tsinghua.edu.cn