賈汝娟 王蒼龍 楊陽(yáng) 茍學(xué)強(qiáng) 陳建敏 段文山3)?

2)(中國(guó)科學(xué)院蘭州化學(xué)物理研究所,固體潤(rùn)滑國(guó)家重點(diǎn)實(shí)驗(yàn)室,蘭州 730000)

3)(甘肅省原子分子物理與功能材料重點(diǎn)實(shí)驗(yàn)室,蘭州 730070)

(2012年10月7日收到;2012年11月5日收到修改稿)

1 引言

凝聚態(tài)物理中有許多非線性現(xiàn)象都可以用一個(gè)處于周期外勢(shì)的原子鏈來(lái)描述,這就是著名的Frenkel-Kontorova(FK)模型[1].FK模型在研究非平衡性質(zhì)及其他物理領(lǐng)域,特別是在凝聚態(tài)物理中被廣泛應(yīng)用,如在超導(dǎo)體中的旋渦晶格[2,3]、電荷密度波(CDW)[4]、膠體[5]、熱傳導(dǎo)[6,7]、公度-不公度(CI)相變[8-10]等等,尤其是對(duì)固體摩擦現(xiàn)象[11-13].FK模型越來(lái)越多地受到研究者們的關(guān)注,因?yàn)樗鳛樯钊胙芯考{米摩擦學(xué)領(lǐng)域復(fù)雜問(wèn)題的一種理論工具[14],可以使人們更容易理解納米摩擦學(xué)的機(jī)理[15,16].雖然目前有許多一維FK模型的理論研究[17-21],但真實(shí)的物理系統(tǒng)并沒(méi)有這么簡(jiǎn)單,所以將一維FK模型推廣到二維FK模型很有必要.

近幾年,在納米摩擦學(xué)的研究中,從相互接觸的兩層原子間摩擦力的觀測(cè)中可知,錯(cuò)合角對(duì)摩擦力有著明顯的影響,對(duì)于一定大小的錯(cuò)合角,可以產(chǎn)生超潤(rùn)滑[9,22,23].這些結(jié)果和超導(dǎo)電性[24-27]中發(fā)現(xiàn)的結(jié)果相似,這就促使了進(jìn)一步的理論研究,對(duì)摩擦力產(chǎn)生的微觀機(jī)理的認(rèn)識(shí),對(duì)生物醫(yī)學(xué)和工程材料的巧妙設(shè)計(jì)具有理論指導(dǎo)作用[21,28,29].

基于以上提到的理論[30]和實(shí)驗(yàn)[31,32]研究,我們討論了二維FK模型,基底勢(shì)采用周期性六角對(duì)稱結(jié)構(gòu),系統(tǒng)從locked到sliding態(tài)的相變.運(yùn)用分子動(dòng)力學(xué)模擬方法[33,34],可以得到,隨著外驅(qū)動(dòng)力的增加,在某一臨界力處,系統(tǒng)開(kāi)始運(yùn)動(dòng),原子的運(yùn)動(dòng)方向和外驅(qū)動(dòng)力的方向不同.隨著外驅(qū)動(dòng)力的進(jìn)一步增加,在另一臨界力處,原子開(kāi)始在外驅(qū)動(dòng)力的方向上運(yùn)動(dòng).我們將這兩個(gè)臨界力分別定義為最大靜摩擦力和動(dòng)摩擦力,為了更清楚地研究這兩個(gè)力的特征,我們討論了系統(tǒng)中外驅(qū)動(dòng)力的大小和方向、勢(shì)壘高度的大小、上層原子間的耦合系數(shù),尤其是錯(cuò)合角θ對(duì)這兩個(gè)臨界力的影響.

2 模型

二維FK模型由上下兩層原子層組成,上層原子采用六角對(duì)稱結(jié)構(gòu),在外驅(qū)動(dòng)力作用下,任意一個(gè)原子不僅受到最近鄰6個(gè)粒子的作用,而且還受到下層二維周期性墊底勢(shì)的作用.模擬中,基底勢(shì)函數(shù)采用近似的六角對(duì)稱結(jié)構(gòu)：

其中,f為勢(shì)函數(shù)的勢(shì)壘高度,b為基底勢(shì)x方向的周期.上下兩層間的錯(cuò)合角為θ,因此墊底勢(shì)可改寫為

上層原子中的任意一個(gè)原子(n,m),它受到最近鄰6個(gè)原子的作用,原子之間的相互作用勢(shì)采用簡(jiǎn)單的簡(jiǎn)諧形式：(yi,j-yn,m-a)2],其中,K為耦合系數(shù),a為上層原子間的平衡距離,任意原子(n,m)的位置可表示為(xn,m,yn,m),它的位移矢量rn,m(xn,m,yn,m)滿足運(yùn)動(dòng)學(xué)方程：

式中,Fext=(Fextcosα,Fextsinα)是外驅(qū)動(dòng)力,α為外驅(qū)動(dòng)力Fext與x軸之間的夾角,γ為黏性阻尼系數(shù).模擬中,為了計(jì)算方便,對(duì)運(yùn)動(dòng)方程(3)進(jìn)行無(wú)量綱化處理,并設(shè)每個(gè)原子的質(zhì)量Mn,m=1,b=1,f=0.01,K=1,γ=0.1.

在數(shù)值模擬中,我們采用四階龍格-庫(kù)塔法求解運(yùn)動(dòng)方程(3),初始時(shí)刻所有原子處在能量的最低點(diǎn),外力Fext在絕熱條件下逐漸增大,且對(duì)于每個(gè)Fext值都要給一個(gè)足夠長(zhǎng)的弛豫時(shí)間使系統(tǒng)達(dá)到一個(gè)穩(wěn)定的狀態(tài).

3 數(shù)值結(jié)果及其分析

3.1 摩擦力(Fs和Fc)與外驅(qū)動(dòng)力Fext的關(guān)系

圖1描述了上層原子的平均速度ˉυ隨外驅(qū)動(dòng)力Fext的變化曲線,初始時(shí)刻上層的每個(gè)原子均勻地分布在下層周期勢(shì)阱底,在外驅(qū)動(dòng)力的作用下,上層原子的平均速度逐漸地增大,系統(tǒng)從locked態(tài)到sliding態(tài)的過(guò)程中會(huì)產(chǎn)生一個(gè)臨界力,當(dāng)外驅(qū)動(dòng)力小于這個(gè)臨界力時(shí),上層原子的平均速度為零,此時(shí)系統(tǒng)處于靜止?fàn)顟B(tài).反之,當(dāng)外驅(qū)動(dòng)力大于該臨界力時(shí),系統(tǒng)會(huì)產(chǎn)生相對(duì)運(yùn)動(dòng).我們稱該臨界力為系統(tǒng)的最大靜摩擦力Fs(簡(jiǎn)稱靜摩擦力),由于外驅(qū)動(dòng)力沿著不同的方向α作用到原子上時(shí),系統(tǒng)受到的靜摩擦力Fs的大小也不同,所以系統(tǒng)的靜摩擦力Fs會(huì)受到外驅(qū)動(dòng)力方向的影響.β為原子的平均速度與x軸之間的夾角,把上層原子的平均速度方向與外驅(qū)動(dòng)力的方向相同時(shí)對(duì)應(yīng)的外驅(qū)動(dòng)力的值定義為動(dòng)摩擦力Fc,即α的變化從α/=β到α=β時(shí)的外驅(qū)動(dòng)力的值.其中,Fc也受到外驅(qū)動(dòng)力的方向α的影響.

圖1 上層原子的平均速度隨外驅(qū)動(dòng)力的變化曲線 模型參數(shù)：f=0.01,K=1,θ=0°

圖2 所示為Fs和Fc受外驅(qū)動(dòng)力方向α的影響,當(dāng)θ=0°時(shí),曲線Fs和Fc關(guān)于α=60°對(duì)稱,這與我們選擇的基底勢(shì)函數(shù)的形狀有關(guān).根據(jù)圖2所示,我們把參數(shù)空間分為三個(gè)不同的區(qū)域：

1)區(qū)域AA(arbitrary angle)Fext＜Fs,系統(tǒng)的平均速度ˉυ=0;

2)區(qū)域 CA(constant angle)Fs＜Fext＜Fc,β/=α;

3)區(qū)域SA(same angle)Fext＞Fc,β=α.

在區(qū)域AA(ˉυ=0)中,上層系統(tǒng)處于靜止?fàn)顟B(tài).在區(qū)域CA中,ˉυ的值不為零,系統(tǒng)開(kāi)始運(yùn)動(dòng),但運(yùn)動(dòng)方向和外驅(qū)動(dòng)力的方向不相同.同時(shí),靜摩擦力 Fs關(guān)于α=60°對(duì)稱,并且在α=0°,α=60°,α=120°處達(dá)到了最大值,然而在α=30°和α=90°處Fs有最小值.

當(dāng)θ=0°時(shí),Fc隨α的變化分別在α=25°,α=35°,α=85°和α=95°處有四個(gè)相同的峰值,在α=30°和α=90°處,Fc有最小值,Fc也關(guān)于α=60°對(duì)稱.所以在區(qū)域CA中,Fs和Fc受α的影響呈現(xiàn)出周期性的變化,數(shù)值模擬結(jié)果表明,靜摩擦力周期性的變化與模型中采用的基底勢(shì)函數(shù)的結(jié)構(gòu)有關(guān),即六角對(duì)稱結(jié)構(gòu).在區(qū)域SA,β=α,原子的運(yùn)動(dòng)方向和外驅(qū)動(dòng)力的方向相同.

圖2 Fs和Fc依賴于外驅(qū)動(dòng)力的大小和方向 在 f=0.01,K=1,θ=0°;在區(qū)域AA,系統(tǒng)的平均速度為零;在區(qū)域CA,β/=α;在區(qū)域SA,β=α

圖3 Fs和Fc依賴于外驅(qū)動(dòng)力的大小和方向 在 f=0.01,K=1,θ=30°;在區(qū)域AA,系統(tǒng)的平均速度為零;在區(qū)域CA,β/=α;在區(qū)域SA,β=α

對(duì)于θ=30°的特殊情況,如圖3所示.與圖2相比,也有三個(gè)區(qū)域：AA,CA和SA,分別與Fext＜Fs,Fs＜Fext＜Fc和Fext＞Fc相對(duì)應(yīng).隨著α的變化,在 α=55°,α=65°和 α=115°處,Fc有三個(gè)相同的峰值;在 α=0°,α=60°,α=120°處,Fc有最小值,但不再關(guān)于α=60°對(duì)稱.而Fs在θ/=0°比θ=0°時(shí)的值小,幾乎為零,不隨α的變化而變化.因此,這兩種臨界力也依賴于錯(cuò)合角θ.

3.2 摩擦力(Fs和Fc)與錯(cuò)合角θ的關(guān)系

臨界力Fs和Fc隨著外驅(qū)動(dòng)力的大小和錯(cuò)合角θ的變化見(jiàn)圖4.根據(jù)圖2我們可以劃分為三個(gè)區(qū)域AA,CA和SA,分別對(duì)應(yīng)于Fext＜Fs,Fs＜Fext＜Fc和Fext＞Fc.

圖4 Fs和Fc依賴于錯(cuò)合角θ,f=0.01,K=1,α=0°

在區(qū)域AA,原子的平均速度為零.在區(qū)域CA中,原子的平均速度為有限的值,但是原子的運(yùn)動(dòng)方向和外驅(qū)動(dòng)力的方向不同.在θ=0°,60°,120°時(shí)處,Fs有最大值,且關(guān)于θ=60°對(duì)稱.而在θ=30°,90°時(shí),Fc有最小值,也關(guān)于θ=60°對(duì)稱,同時(shí),在θ=25°,35°,85°,95°時(shí),Fc有四個(gè)幾乎相等的峰值.在區(qū)域SA中,原子在外驅(qū)動(dòng)力的方向上運(yùn)動(dòng).從圖4中可以看到,這兩種臨界力Fs和Fc明顯地依賴于錯(cuò)合角θ.

圖5 Fs和Fc依賴于錯(cuò)合角θ,f=0.01,K=1,α=30°

圖5 中,對(duì)于α=30°的情況,與圖4相比較,Fs和Fc隨著θ的變化曲線不再關(guān)于θ=60°對(duì)稱.隨著θ值的變化,Fs和Fc的曲線呈現(xiàn)出周期性變化.顯然α也對(duì)Fs和Fc有影響.

從以上的比較分析中,可以得到,Fs和Fc依賴于錯(cuò)合角θ的同時(shí),α對(duì)Fc的影響也很大,進(jìn)而Fs和Fc隨著θ和α的變化而變化的情況比較復(fù)雜,這里我們只分析了α=0°,θ/=0°或α/=0°,θ=0°的情況.

3.3 摩擦力(Fs和Fc)與原子間的耦合系數(shù)K的關(guān)系

對(duì)于給定不同的θ和α的值,Fs和Fc依賴于原子的耦合系數(shù)K,見(jiàn)圖6.

圖 6(a)中,α=0°,θ=0°,Fs不受 K 的影響,但Fc明顯地受到參數(shù)K的影響.此時(shí)有三個(gè)區(qū)域：AA,CA和SA.在區(qū)域AA內(nèi),原子幾乎不動(dòng).在區(qū)域CA內(nèi),原子的運(yùn)動(dòng)方向和外驅(qū)動(dòng)力的方向不同.在區(qū)域SA中,原子在外驅(qū)動(dòng)力的方向上移動(dòng).圖 6(b)中,α=10°,θ=0°時(shí),Fs和 Fc的值也不相同,只有Fc依賴于耦合系數(shù)K,也有三個(gè)區(qū)域AA,CA和SA,分別與Fext＜Fs,Fs＜Fext＜Fc和Fext＞Fc相對(duì)應(yīng).在θ=0°時(shí)對(duì)于不同的α值,Fs的值幾乎沒(méi)有變化且不為零.對(duì)于θ/=0°,在圖6(c)和圖 6(d)中,分別為 α=0°,θ=10°和 α=10°,θ=10°.隨著耦合系數(shù)K的增加,Fc的值的變化也很大.但是Fs的值變化很小,幾乎為零,所以Fs與耦合系數(shù)K是相互獨(dú)立的,即Fs不依賴于耦合系數(shù)K,在這種情況下,超潤(rùn)滑可能會(huì)產(chǎn)生.

圖 6 Fs和 Fc隨著耦合系數(shù) K 的變化曲線 f=0.01;α,θ值：(a)α=0°,θ=0°;(b)α=10°,θ=0°;(c)α=0°,θ=10°;(d)α=10°,θ =10°

3.4 摩擦力(Fs和Fc)與勢(shì)壘高度 f的關(guān)系

Fs和Fc隨著勢(shì)壘高度的變化曲線,如圖7所示：(a)α=0°,θ=0°;(b)α=10°,θ=0°;(c)α=0°,θ=10°;(d)α=10°,θ=10°.在圖 7(a)中,當(dāng)外驅(qū)動(dòng)力足夠大時(shí),可以看成有兩個(gè)區(qū)域：AA和SA,此時(shí),Fs=Fc.當(dāng) f繼續(xù)增加時(shí),Fs和Fc也隨著增加.對(duì)于圖7(b),(c),(d),我們可以看到有三個(gè)區(qū)域：AA,CA和SA且分別與Fext＜Fs,Fs＜Fext＜Fc,Fext＞Fc相對(duì)應(yīng).Fs和Fc也隨著 f增加而增加.其中在圖7(c)和圖7(d)中,當(dāng)θ/=0°時(shí),靜摩擦力很小,當(dāng) f接近于零時(shí),靜摩擦力也接近于零.

圖 7 Fs和 Fc隨著勢(shì)壘高度 f的變化曲線 K=1;α,θ 值：(a)α=0°,θ=0°;(b)α=10°,θ=0°;(c)α=0°,θ=10°;(d)α=10°,θ =10°

根據(jù)以上的分析,我們要獲得超潤(rùn)滑,必須采用勢(shì)壘高度 f小的材料.同時(shí),為了得到很小的摩擦力,必須選擇合適的錯(cuò)合角θ.

4 結(jié)論

在具有六角晶格對(duì)稱結(jié)構(gòu)的二維FK模型中,我們研究了系統(tǒng)從locked態(tài)到sliding態(tài)的相變,數(shù)值模擬了具有六角晶格對(duì)稱結(jié)構(gòu)的摩擦行為.隨著外驅(qū)動(dòng)力的增加,在某一個(gè)臨界值處,原子開(kāi)始運(yùn)動(dòng).隨著外力的進(jìn)一步增大,在另一臨界力處,原子在外驅(qū)動(dòng)力的方向上運(yùn)動(dòng).把這兩個(gè)臨界力定義為兩種不同的摩擦力,靜摩擦力和動(dòng)摩擦力,它們依賴于外驅(qū)動(dòng)力的大小和方向,勢(shì)壘高度,耦合系數(shù)和錯(cuò)合角.對(duì)于一定大小的錯(cuò)合角,摩擦力的值很小.特別是當(dāng)錯(cuò)合角θ/=0°的情況,系統(tǒng)容易出現(xiàn)超潤(rùn)滑現(xiàn)象,我們將選擇上層原子耦合系數(shù)較大和勢(shì)壘高度較小的材料.同時(shí),為了得到較小的摩擦力必須選擇合適的錯(cuò)合角θ.

[1]Kontorova T A,Frenkel Y I 1938 Zh.Eksp.Teor.Fiz.8 1340

[2]Blatter G,Feigel’man M V,Geshkenbein V B,Larkin A I,Vinokur V M 1994 Rev.Mod.Phys.66 1125

[3]Marley A C,Higgins M J,Bhattacharya S 1995 Phys.Rev.Lett.74 3029

[4]Gr¨uner G 1988 Rev.Mod.Phys.60 1129

[5]Reichhardt C,Olson Reichhardt C J 2004 Phys.Rev.Lett.92 108301

[6]Hu B,Yang L 2005 Chaos 15 015119

[7]Shao Z G,Yang L,Zhong W R,He D H,Hu B 2008 Phys.Rev.E 78 061130

[8]Bak P 1982 Rep.Prog.Phys.45 587

[9]Yang Y,Wang C L,Duan W S,Shi Y R,Chen J M 2012 Acta Phys.Sin.61 130501(in Chinese)[楊陽(yáng),王蒼龍,段文山,石玉仁,陳建敏2012物理學(xué)報(bào)61 130501]

[10]Li R T,Duan W S,Yang Y,Wang C L,Chen J M 2011 Euro.Phys.Lett.94 56003

[11]Xu Z M,Huang P 2006 Acta Phys.Sin.55 2427(in Chinese)[許中明,黃平2006物理學(xué)報(bào) 55 2427]

[12]Li X L,Liu F,Lin M M,Chen J M,Duan W S 2010 Acta Phys.Sin.59 2589(in Chinese)[李曉禮,劉鋒,林麥麥,陳建敏,段文山2010物理學(xué)報(bào) 59 2589]

[13]Yang Y,Duan W S,Chen J M,Yang L,Teki′c J,Shao Z G,Wang C L 2010 Phys.Rev.E 82 051119

[14]Persson B N J 1999 Surf.Sci.Rep.33 83

[15]Qian L M,Luo J B,Wen S Z,Xiao X D 2000 Acta Phys.Sin.49 2240(in Chinese)[錢林茂,雒建斌,溫詩(shī)鑄,蕭旭東2000物理學(xué)報(bào)49 2240]

[16]Wen S Z 1998 Nano Tnbology(Beijing：Tisinghua University Press)(in Chinese)[溫詩(shī)鑄1998納米摩擦學(xué)(北京：清華大學(xué)出版社)]

[17]Zheng Z G 2001 Commun.Theor.Phys.36 37

[18]Hirano M 2003 Wear 254 932

[19]Xu H B,Wang G R,Chen S G 2000 Chin.Phys.9 0611

[20]Braun O M,Zhang H,Hu B,Teki′c J 2003 Phys.Rev.E 67 066602

[21]Li H X,Xu T,Wang C B,Chen J M,Zhou H D,Liu H W 2007 Tribol.Int.40 132

[22]Lin M M,Duan W S,Chen J M 2010 Chin.Phys.B 19 026201

[23]Wang C L,Duan W S,Hong X R,Chen J M 2008 Appl.Phys.Lett.93 153116

[24]Reichhardt C,Gr?nbech-Jensen N 2001 Phys.Rev.B 63 054510

[25]Reichhardt C,Zim′anyi G T,Gr?nbech-Jensen N 2001 Phys.Rev.B 64 014501

[26]Reichhardt C,Zim′anyi G T,Scalettar R T,Hoffmann A,Schuller I K 2001 Phys.Rev.B 64 052503

[27]Reichhardt C,Olson C J,Hastings M B 2002 Phys.Rev.Lett.89 024101

[28]Persson B N J,Albolu O,Tartaglino U,Volokitin A I,Tosatti E 2005 J.Phys.：Condens.Matter 17 R01

[29]Yang Y,Duan W S,Yang L,Chen J M,Lin M M 2011 Euro.Phys.Lett.93 16001

[30]Teki′c J,Braun O M,Hu B 2005 Phys.Rev.E 71 026104

[31]Yaron U,Gammel P L,Huse D A,Kleiman R N,Oglesby C S 1995 Nature 376 753

[32]Pardo F,de la Cruz F,Gammel P L,Bucher E,Bishop D J 1998 Nature 396 348

[33]Zhang Z H,Han K,Li H P,Tang G,Wu Y X,Wang H T,Bai L 2008 Acta Phys.Sin.57 3160(in Chinese)[張兆慧,韓奎,李海鵬,唐鋼,吳玉喜,王洪濤,白磊2008物理學(xué)報(bào) 57 3160]

[34]Gong B Z,Zhang B J 2009 Acta Phys.Sin.58 1504(in Chinese)[龔博致,張秉堅(jiān)2009物理學(xué)報(bào)58 1504]