馮翠菊 劉 英

(1.華北科技學院基礎部，北京 東燕郊 101601;2.河北師范大學物理與信息工程學院，河北 石家莊 050800)

基于第一性原理計算CuxNi19-x(x＜19)混合團簇的結構研究①

馮翠菊1②劉 英2

(1.華北科技學院基礎部，北京 東燕郊 101601;2.河北師范大學物理與信息工程學院，河北 石家莊 050800)

采用基于第一性原理計算的Chen-M?bius三維晶格反演基礎上獲得的原子間相互作用勢參數，運用最陡下降法和共軛梯度法結合起來的最優(yōu)化計算Smart方法得到了CuxNi19-x(x＜19)二元過渡金屬混合團簇的穩(wěn)態(tài)構型，發(fā)現Cu-Ni混合團簇中有明顯Cu元素偏析現象，即Cu和Ni并沒有互相混合形成有序結構，而是分別聚集在一起，所形成的體系對稱性很低，呈割據狀態(tài)。

晶格反演，Cu-Ni混合團簇，穩(wěn)態(tài)結構

0 引言

團簇結構和性質的研究對于理解物質從微觀到宏觀的過渡具有重要作用，團簇結構研究的基本問題是弄清團簇如何由原子、分子一步步發(fā)展而成，以及隨著這種發(fā)展，團簇的結構和性質如何變化。混合團簇可以極大的擴展單質團簇的性質，而且通過混合可以得到和加強團簇的某些特殊性質。由于組成、結構及性質的豐富多樣性，混合團簇在電子、工程、催化等領域具有廣泛的應用前景。

團簇在理論計算中存在著多種計算方法。最精確的從頭計算方法，由于其計算量與基函數數目的五次方或六次方成正比，限制了其在電子數目較多的體系(即較大尺寸團簇中)的應用。局域密度泛函方法在小團簇計算中取得了巨大的成功，但當團簇尺寸增大到20個原子以上時，該方法同樣受到了計算量的限制，對大中型團簇顯得無能為力。目前對于金屬團簇的研究主要是基于Gupta 、Sutton-Chen、EAM 等經驗勢，Gupta勢已成功研究了許多團簇的結構，現在已有關于Cu、Ni、Ag、Au 等團簇［1］及 Ni-Al、Cu-Co、Cu-Au 團簇［2，3］的報道。

基于第一性原理計算的材料的結合能，運用Chen-M?bius反演［4-6］，可以嚴格導出的原子間相互作用對勢，該方法廣泛地應用于研究材料的靜力學和動力學特性［7-10］。本文將Chen-M?bius晶格反演應用推廣到過渡金屬二元團簇的研究，運用最優(yōu)化方法計算得到了CuxNi19-x(x＜19)二元過渡金屬團簇的穩(wěn)態(tài)結構，由此開辟了一條基于第一性原理計算的、較高效率地研究二元金屬團簇的結構及其特性的途徑。

1 勢函數的獲得與計算方法

1.1 勢函數的獲得

基于第一性原理計算的Chen-M?bius晶格反演獲得的原子間相互作用對勢，不依賴于經驗的形式和方法，直接由金屬和金屬間化合物的第一性原理自洽能帶計算和嚴格的晶格反演過程得出，其準確性僅僅取決于第一性原理計算的準確度和對勢近似的可靠性，因此它在材料的物理性質的理論計算中得到了廣泛的應用。Chen-M?bius反演的理論要點如下：

對于同種原子構成的晶體。設從第一性原理算出的晶體結合能函數可表示為：

式中x為最近鄰原子間距，r0(n)為第n階近鄰的配位數，b0(n)是以最近鄰距離為單位的、第n級近鄰到參考原子的距離，Φ(x)為對勢函數。為了理論上處理的方便，通過{b0(n)}的自乘得到{b(n)}，{b(n)}可認為構成乘法半群，這里多出的虛格點相應的虛配位數均為零。這時有

這里

由此可反演出的原子間相互作用勢普遍公式

式中I(n)可通過下式求出

這里求和號下b(d)|b(n)表示對所有這樣的b(d)求和：即b(d)∈{b(n)}，且可以在{b(n)}中找到一個{b(k)}，使得b(d)b(k)=b(n)。

可以證明，(2)式是(1)式的解：

關于異種原子間相互作用勢問題可以通過以下的方案解決［6］，即認為通過純元素晶體反演出的原子間相互作用勢在二元合金中仍然有效，如此，可以從合金的總結合能中扣除同種原子之間相互作用對結合能的貢獻而得出異種原子的“偏結合能”曲線，從偏結合能曲線反演出異種原子之間相互作用勢。

表一中給出了Chen-M?bius晶格反演得到的原子間相互作用對勢，然后利用Morse勢形式擬合得到的參數，其中Rcut表示計算原子間相互作用對勢的截止距離。

1.2 計算方法

基于原子間相互作用勢，尋找團簇的穩(wěn)態(tài)結構可以通過團簇的能量最小化過程實現。團簇的能量最小化過程實質上是一個多元函數的極值問題，我們采用最陡下降法(Steepest descents)和共軛梯度方法(Conjugate gradient)有機結合起來的Smart方法［11］。

由于在極值點附近能量的梯度趨于零，最陡下降法在系統(tǒng)趨近極值點時收斂變慢，盡管如此，這種方法由于運算穩(wěn)定，仍被廣泛采用或作結構能量最小化的預處理方法。共軛梯度方法以及各種Newton-Raphson方法對于偏離局域極值較大的系統(tǒng)構型(即勢能曲面偏離二次曲面較大)是不穩(wěn)定的。另外，Newton-Raphson方法需要的存儲空間較大，不適于大體系計算。

隨著團簇的尺度增加，團簇可能的構型數目增加非?？?，顯然我們不能用窮舉的方法得到團簇的基態(tài)。實際計算中，判斷系統(tǒng)最小化后的構型是局域極小還是全域極小，是一個難以回答的問題。在我們的計算中，采用不同的系統(tǒng)構型進行多次能量最小化，由各極值點的能量值判斷能量最優(yōu)化結構。對每一個固定原子數的團簇，我們隨機地選取系統(tǒng)1000～2000個不同的初始構型進行能量最小化。本文首先采用最陡下降法對團簇的結構進行預處理，然后采用共軛梯度方法計算團簇的穩(wěn)態(tài)結構，從而找到二元金屬團簇的最小能量結構。當然所找到的每一固定原子數的穩(wěn)態(tài)結構不止一種，每種原子數的團簇都分別對應著多種穩(wěn)態(tài)結構，但其中有些是亞穩(wěn)態(tài)，通過對它們的能量比較，我們可以找到(近似的)最小能量結構，也就是要找的最穩(wěn)態(tài)結構

表1 Chen-M?bius晶格反演得到的原子間相互作用對勢參數

2 結果與討論

圖1列出了CuxNi19-x(x＜19)各種成分的穩(wěn)態(tài)構型，從中可以看到隨著Cu原子數的逐漸增大，混合團簇CuxNi19-x的穩(wěn)態(tài)構型都是三個五邊形金字塔結構互相交叉構成的對稱結構，即雙二十面體套接結構，但由于Cu原子和Ni原子的半徑不同而發(fā)生了輕微的形變。在CuNi18中，唯一的Cu原子占據了雙二十面體的一個頂點，而不是占據中心位置，這是由于兩種元素具有不同的表面能造成的，Cu的表面能比Ni的表面能低，在Cu原子和Ni原子形成混合團簇時，Cu原子會自動占據團簇的表面位置，而Ni卻占據了內部位置，形成了偏析現象。同樣的現象也發(fā)生在x=2、3、4，以及 x=12～18的情況下，即在Cu、Ni原子數混合比例相差較大時，有明顯的Cu偏析現象。當兩者的混合比例相差不大時，我們同樣觀察到偏析現象，Cu和Ni并沒有互相混合形成有序結構，而是分別聚集在一起，所形成的體系對稱性很低，呈割據狀態(tài)。這種割據狀態(tài)的形成，可以看成Cu、Ni合金微觀不均勻性的具體體現，這導致了宏觀上Cu-Ni固溶體的體現。

圖1 CuXNi19-x的穩(wěn)態(tài)構型 深色球代表Ni，淺色代表Cu

同樣我們也計算了CuxNi55-x的穩(wěn)態(tài)構型，圖2我們給出了Cu8Ni47、Cu27Ni28的穩(wěn)態(tài)構型，這些穩(wěn)態(tài)構型都是二十面體結構，一個原子占據了二十面體的中心，外殼層包含了兩層原子，最外面一層包含了42個原子，次外層包含了12個原子。由圖可以明顯看出，Cu、Ni原子混合比例相差較大時，Cu原子明顯的偏析，混合比例相近時，Cu也偏析呈明顯的割據狀態(tài)，這一計算結果與J.M.Montejano-Carrizales等人［12］計算結果相一致，這說明了我們計算的可靠性。

圖2 Cu8Ni47Cu27Ni28的穩(wěn)態(tài)構型深色球代表Ni，淺色代表Cu

3 結論

本文采用基于第一性原理計算的Chen-M?bius三維晶格反演基礎上獲得的原子間相互作用勢參數，運用最陡下降法和共軛梯度法結合起來的最優(yōu)化計算Smart方法得到了Cu-Ni二元過渡金屬混合團簇的穩(wěn)態(tài)構型，發(fā)現CuNi混合團簇的穩(wěn)態(tài)構型與純金屬Ni團簇的基態(tài)構型完全一樣，只是由于Cu、Ni原子半徑的不同而略有畸變;混合團簇中Cu原子趨于占據團簇的表面位置，而Ni原子趨于占據中心位置，即Cu-Ni混合團簇中有明顯Cu元素偏析現象，Cu和Ni并沒有互相混合形成有序結構，而是分別聚集在一起，所形成的體系對稱性很低，呈割據狀態(tài)。

［1］ J.Oviedo，R.E.Palmer，J.Chem.Phys.，117(2002)，9548.

［2］ J.Jellinek，E.B.Krissinel，Chem Phys，lett.258(1996)，283; ib.id， Chem， Phys.lett.272(1997)，30

［3］ J.Wang，F.Ding，el.al.，Solid State Commun，119(2001)，13.

［4］ M.Li and N.X.Chen，Phys.Rev.B52(2)(1995)997-1003.

［5］ M.Li and N.X.Chen，Z.Phys.，B100(2)(1996)169 - 172; M.Li， S.J.Liu， N.X.Chen，Y.G.Teng， W.X.Zhang， and S.Z.Zhang，Phys.Lett.，A169(5)(1992)364-370.

［6］ M.Li，S.J.Liu，and N.X.Chen，Phys.Lett.，A177(2)(1993)134-138.

［7］ N.X.Chen，Z.D.Chen，Y.N.Shen，S.J.Liu，M.Li，Phys.Lett.，A184(4-5)(1994)347-351.

［8］ M.Li and N.X.Chen，Phys.Rev.B52(2)(1995)997-1003.

［9］ M.Li and N.X.Chen，Z.Phys.，B100(2)(1996) 169 - 172; M.Li， S.J.Liu，N.X.Chen，Y.G.Teng，W.X.Zhang，and S.Z.Zhang， Phys.Lett.， A169 (5)(1992)364-370.

［10］ M.Li，S.J.Liu，and N.X.Chen，Phys.Lett.，A177(2)(1993)134-138.

［11］ J.P.Perdew，J.A.Chevary，S.H.Vosko，K.A.Jackson，M.R.Pederson，D.J.Singh，C.Fiolhais，Phys.Rev.，B46(11)(1992)6671-6687;J.P.Perdew，Physics Letters，A165(1)(1992)79-82;J.P.Perdew， Y.Wang，Phys.Rev.B45(23)(1992)13244-13249.

［12］ J.M.Montejano-CarrizalesM.P.Iniguer and.A.Alonso，Phys.Rev.B49(1995)16649-16658.

First-principles study of the stable structures of CuxNi19-x(x＜19)bimetallic clusters

FENG CuiJu

(Department of basic courses，North China institute of Science And Technology，Yanjiao Beijing-East 101601)

Using first principles total-energy calculations and the Mobius inversion method in metal Cu，Ni and Cu-Ni，We obtained the pair-potential from cohesive energy.Structural of Cu-Ni bimetallic clusters derived from the smart method are presented and discussed.In Cu-Ni clusters，the Cu atoms preferentially occupy the surface sites，whereas the Ni atoms occupy the interior of the cluster.Cu segregates to the surface in the Cu-Ni clusters.

Chen’s 3D lattice inverse method;bimetallic clusters;stable structure

O469

A

1672-7169(2011)01-0071-05

2011-01-04?；痦椖浚喝A北科技學院科研基金(B09017)。

馮翠菊(1970-)，女，河北藁城人，碩士，華北科技學院基礎部講師，研究方向：材料的模擬與設計。