第一性原理計算Au-Ag共格界面能

張慧玲， 房 勇

（1.太原學院 材料與化學工程系，山西 太原030032； 2.四川輕化工大學 材料科學與工程學院，四川 自貢643000）

納米金屬多層膜是兩種或兩種以上金屬交替沉積，最終形成層間界明顯、組分周期性變化的薄膜材料［1］.目前，人們已制備出的納米金屬多層膜主要有Cu／Ag、Cu／Ni、Co／Au、Au／Ni等［2-4］.研究表明，納米金屬多層膜在機械性能、光學性能、電學和電化學性能等許多方面有著獨特的優勢，在某些性能方面明顯優于其單一組分金屬薄膜.界面是晶體中的面缺陷，對晶體材料的性質和發生轉變的過程有重要影響［1］.近年來，新興的納米多層膜研究更是將界面的研究推向了一個新的階段［2-4］.界面處原子排列混亂而使系統能量升高，額外升高的自由能部分稱為界面能，界面能是描述金屬界面基本屬性的重要參數，它在解釋各種物理化學過程中起著決定性作用.實驗上測定界面能非常復雜，通常采用理論方法進行估算，其中分子動力學方法在研究同質界面上起著重要作用［5-6］，但由于方法中對勢函數的依賴性，對異質界面研究很少.第一性原理模擬方法能得到準確的表面能數值［7-8］，但由于該方法對計算模型大小的限制，很難對異質界面展開研究.

本文運用密度泛函理論的第一性原理理論，首先采用廣義梯度近似方法計算了FCT結構Au、Ag的晶格常數，隨后建立Au-Ag以（100）面結合的異質界面模型，定量計算了該界面能.

1 計算模型及計算參數

Au、Ag為面心立方結構金屬，晶格常數分別為0.407 83、0.408 57 nm.為了獲得合理的Au-Ag共格界面模型，建立FCT結構Au、Ag晶體（a＝b≠c，α＝β＝γ＝90°），模型如圖1所示，其中a＝b＝0.408 2 nm，以c為自變量.隨后采用層晶模型［9］（Slab模型）建立Au、Ag（100）自由表面模型，及Au-Ag以（100）面結合的共格異質界面模型.共格異質界面模型如圖2所示，模型中Au、Ag原子層數分別為8、7層，中間接觸距離為0.437 5 nm，體系共有60個原子.

圖1 FCT晶體結構Fig.1 The crystal structure of FCT

圖2 Au-Ag以（100）晶面為接觸的共格界面模型，深色為金原子，淺色為銀原子Fig.2 The model of Au-Ag coherent interface matched along（100）crystallographic plane，the ball with dark color is Au，with light color is Ag

所有的計算工作都是在Material studio 4.0軟件包中的Castep計算模塊下完成，計算采用周期性邊界條件，用廣義梯度近似（GGA）的PBE處理電子間的交換關聯能，自洽收斂精度設置為2×10－5eV／atom，原子間相互作用力收斂標準為0.05 eV／A，Brillouin區積分K網格點采用4×4×10.測試計算發現增加層晶厚度、真空層厚度或增加表面原子數目對計算結果影響很小.

2 結果及分析

為了驗證數據的有效性，采用公式

計算FCC結構Ag的低密勒指數面表面能，其中，Esurf為表面能，Eslab為層晶總能量，Ebulk為FCC結構Ag能量，n為層晶中原子數目.如表1所示，比較發現，計算結果與文獻［9］結果相近，說明我們計算的結果是可信的.

表1 FCC結構Ag低密勒指數表面表面能Tab.1 Surface Energy of FCC structure Ag low miller index planes

為了獲得合理的共格界面模型，首先建立如圖1所示FCT Ag晶體結構（a＝b＝0.408 2 nm），計算體系總能隨c值的變化關系，結果如圖3所示，可知該結構最低能量所對應的c值為0.428 nm，同樣的方法可得Au所對應的c值為0.446 nm.

圖3 a＝0.408 2 nm FCT結構Ag晶體總能隨c變換關系圖Fig.3 Dependence of the total energy of FCT Ag on the numbers of the lattice constants，a＝0.408 2 nm

隨后采用層晶模型，建立Au-Ag自由（100）表面及Au-Ag以（100）面結合的共格異質界面.對FCT結構的Au、Ag（100）表面進行幾何結構優化并計算表面能，得出Au、Ag（100）自由表面能分別為1.132 5、0.790 4 J／m2.共格異質界面模型如圖2所示，模型中將Au-Ag（100）相互接觸面間距設為自變量，計算slab總能隨面間距的變化關系，得到最低能量所對應間距（0.437 5 nm），隨后進行幾何結構優化并計算層晶總能.界面能采用公式

其中，Einter為界面能，n1、n2分別為模型中FCT結構Au、Ag的原子數目，Esurf為FCT結構表面能，最終得到的 Au-Ag（100）面結合界面能為0.074 1 J／m2.

由計算結果可知，Au-Ag共格異質界面的界面能非常低，一個原因是共格界面并不像自由表面一樣需要斷鍵，另一個原因是因為FCT結構Au、Ag的晶格常數相近，結合成共格界面后彈性應變能較小.如果將FCC到FCT結構轉變所需能量定義為彈性應變能，則需建立復雜模型進行進一步計算.

3 結論

采用近似方法建立共格Au-Ag（100）異質界面模型，并定量計算了該界面能數值為0.074 1 J／m2，為后期建立復雜界面模型打下了基礎.