岳 莉

（凱里學(xué)院，貴州凱里 556011）

合金團(tuán)簇與純金屬團(tuán)簇相比，電子結(jié)構(gòu)更加復(fù)雜，其性質(zhì)與團(tuán)簇的大小和化學(xué)成分有關(guān)［1］.鎳，是一種硬而有延展性并具有鐵磁性的金屬，它能夠高度磨光和抗腐蝕，在軍工、航空航天、機(jī)械、不銹鋼等領(lǐng)域有著重要應(yīng)用，它的化學(xué)性質(zhì)較活潑，主要用于合金及催化劑，對(duì)Ni 團(tuán)簇和以Ni 為基體的合金團(tuán)簇研究已有許多［2?4］.由于過(guò)渡金屬原子中含有d 電子，形成的合金團(tuán)簇就具有更為特殊的電子結(jié)構(gòu)和物理化學(xué)性質(zhì)［5］.鉑是過(guò)渡金屬，有良好的延展性、導(dǎo)熱性和導(dǎo)電性，化學(xué)性質(zhì)不活潑，有很高的化學(xué)穩(wěn)定性和催化活性，對(duì)以Pt 為基體的合金團(tuán)簇研究也有不少［5?6］，而合金團(tuán)簇的物理化學(xué)性質(zhì)是可以通過(guò)團(tuán)簇的組成、尺寸和化學(xué)順序進(jìn)行調(diào)節(jié)的，目前，對(duì)NixPt團(tuán)簇結(jié)構(gòu)穩(wěn)定性與磁性的研究還未見(jiàn)報(bào)道.為此，本文采用基于密度泛函理論的第一性原理方法，對(duì)NixPt（x=1?9）團(tuán)簇結(jié)構(gòu)穩(wěn)定性與磁性進(jìn)行研究.

1 計(jì)算方法

NixPt（x=1?9）合金團(tuán)簇構(gòu)型設(shè)計(jì)時(shí)，對(duì)于總原子數(shù)N=x+1 ≤8 的團(tuán)簇，幾乎考慮到每個(gè)尺寸的所有構(gòu)型，同時(shí)參考Nix+1團(tuán)簇構(gòu)型.由于團(tuán)簇結(jié)構(gòu)的簡(jiǎn)并電子態(tài)會(huì)對(duì)其穩(wěn)定性造成影響，因此，對(duì)經(jīng)篩選后的NixPt穩(wěn)定結(jié)構(gòu)設(shè)置不同自旋多重度進(jìn)行優(yōu)化，再比較NixPt團(tuán)簇總能量，能量最低的構(gòu)型即為NixPt 基態(tài)結(jié)構(gòu).設(shè)計(jì)NixPt（x=1?9）團(tuán)簇結(jié)構(gòu)的具體方法如下：首先在穩(wěn)定、次穩(wěn)定、亞穩(wěn)定Nix+1和同族同尺寸團(tuán)簇非對(duì)稱位置上用一個(gè)鉑原子替代一個(gè)鎳原子；接著在穩(wěn)定、次穩(wěn)定、亞穩(wěn)定Nix團(tuán)簇的不同位置上添加一個(gè)鉑原子；最后在Nix-1Pt 團(tuán)簇非對(duì)稱位置上添加一個(gè)鎳原子.同時(shí)，為了排除部分局域最小值而更加接近全局最小值，在對(duì)每個(gè)初始結(jié)構(gòu)優(yōu)化的過(guò)程中均沒(méi)有考慮對(duì)稱性.

為了能更精確地描述積分和電子自旋的貢獻(xiàn)，計(jì)算過(guò)程中考慮自旋非限制，角動(dòng)量多級(jí)擴(kuò)展使用十六偶極機(jī)制和雙數(shù)值極化基組.本文采用廣義梯度（GGA）近似下的BLPY 關(guān)聯(lián)函數(shù)，在總能、原子位移和能量梯度的收斂精度設(shè)為1×10?6eV、5×10?4? 和1×10?5eV/?下進(jìn)行計(jì)算.

2 結(jié)果與分析

2.1 團(tuán)簇的最低能量結(jié)構(gòu)

圖1和2分別顯示了Nix+1（x=1?9）團(tuán)簇和NixPt（x=1?9）合金團(tuán)簇的最低能量結(jié)構(gòu)，其中黑色小球代表Pt 原子，表1、表2分別列出了Nix+1（x=1?9）團(tuán)簇和NixPt（x=1?9）合金團(tuán)簇的計(jì)算結(jié)果.

表1 Nix+1（x=1-9）團(tuán)簇基態(tài)結(jié)構(gòu)的對(duì)稱性、平均束縛能、HOMO-LUMO能隙、自旋總磁矩

表2 NixPt（x=1-9）團(tuán)簇基態(tài)結(jié)構(gòu)的對(duì)稱性、平均束縛能、能隙、鉑原子電荷、自旋總磁矩

圖1 Nix+1（x=1?9）團(tuán)簇最低能量結(jié)構(gòu)

從圖1 的Nix+1（x=1?9）和圖2 的NixPt（x=1?9）團(tuán)簇最低能量結(jié)構(gòu)可以看出，Ni6Pt、Ni7Pt 分別與Ni7、Ni8不僅形狀沒(méi)變而且對(duì)稱性也沒(méi)變，說(shuō)明對(duì)于這些尺寸，Pt原子的摻雜并沒(méi)有改變純Ni團(tuán)簇幾何構(gòu)型.另外，NiPt、Ni2Pt、Ni3Pt 分別與Ni2、Ni3、Ni4雖然有相似的最低能量結(jié)構(gòu)形狀，但有的對(duì)稱性卻發(fā)生了變化，如純Ni3金屬團(tuán)簇的最低能量結(jié)構(gòu)是等邊三角形而合金Ni2Pt 團(tuán)簇的最低能量結(jié)構(gòu)卻是等腰三角形，純Ni4金屬團(tuán)簇的最低能量結(jié)構(gòu)是各邊均相等的三棱錐，合金Ni3Pt團(tuán)簇的最低能量結(jié)構(gòu)卻是底為等邊三角形的三棱錐，合金Ni4Pt、Ni5Pt、Ni9Pt團(tuán)簇的最低能量結(jié)構(gòu)與純Ni5、Ni6、Ni10團(tuán)簇的最低能量結(jié)構(gòu)形狀和對(duì)稱性都不一樣，如純Ni5團(tuán)簇是底為四邊形的三棱錐而合金Ni4Pt團(tuán)簇卻是共等邊三角形底的兩個(gè)三棱錐.

圖2 NixPt（x=1?9）團(tuán)簇最低能量結(jié)構(gòu)

2.2 團(tuán)簇的能級(jí)與穩(wěn)定性

團(tuán)簇的平均束縛能和HOMO?LUMO 能隙是反映其穩(wěn)定性的主要依據(jù)，平均束縛能越大結(jié)構(gòu)越穩(wěn)定，HOMO?LUMO 能隙越大，團(tuán)簇化學(xué)性子越穩(wěn)定.為了描述穩(wěn)定性，作者計(jì)算了團(tuán)簇的平均束縛能和HOMO?LUMO能隙，并將其繪制于圖3和圖4中.

圖3 團(tuán)簇平均結(jié)合能隨尺寸增加的變化規(guī)律

圖4 團(tuán)簇能隙隨尺寸增加的變化規(guī)律

從圖3 的團(tuán)簇平均結(jié)合能隨尺寸增加的變化規(guī)律可以看出，NixPt（x=1?9）團(tuán)簇的平均束縛能隨著原子總數(shù)N=x+1的增多而增大，說(shuō)明NixPt（x=1?9）團(tuán)簇的穩(wěn)定性隨尺寸增加而加強(qiáng).與純Nix+1（x=1?9）團(tuán)簇的平均束縛能相比較看出，合金NixPt團(tuán)簇的平均束縛能與純Nix+1團(tuán)簇的走勢(shì)大致相同，都隨著原子總數(shù)的增多而增大，但合金NixPt 團(tuán)簇的平均束縛能都比相同原子總數(shù)的純Nix+1團(tuán)簇的高，說(shuō)明Pt原子的摻雜增加了團(tuán)簇的穩(wěn)定性，其原因是原子的配位數(shù)隨原子總數(shù)增加而增多，合金NixPt團(tuán)簇大于純Nix+1團(tuán)簇，使得合金NixPt團(tuán)簇結(jié)構(gòu)總能量降低較多，釋放的能量較多，合金團(tuán)簇的穩(wěn)定性稍高于同尺寸的純團(tuán)簇.

從圖4 的團(tuán)簇能隙隨尺寸增加的變化規(guī)律可以看出，NixPt（x=1?9）團(tuán)簇的HOMO?LUMO能隙隨著原子總數(shù)N=x+1的增多在振蕩過(guò)程有增大有減小.HOMO?LUMO能隙直接影響著團(tuán)簇的化學(xué)活性和穩(wěn)定性，它的大小反映了電子從占據(jù)軌道向空軌道躍遷的能力，代表了參與化學(xué)反應(yīng)的能力，HOMO?LUMO 能隙值較大時(shí)，電子不容易發(fā)生躍遷，則被認(rèn)為具有較強(qiáng)的化學(xué)穩(wěn)定性，HOMO?LUMO 能隙值較小時(shí)，電子躍遷能力強(qiáng)，則被認(rèn)為具有較強(qiáng)的化學(xué)活性.與純Nix+1（x=1?9）團(tuán)簇的HOMO?LUMO 能隙相比較看出，合金NixPt團(tuán)簇的HOMO?LUMO 能隙與純Nix團(tuán)簇的走勢(shì)不一致，但在N=4 時(shí)都有最大值，說(shuō)明合金Ni3Pt 團(tuán)簇與純Ni4團(tuán)簇一樣有較強(qiáng)的化學(xué)穩(wěn)定性，合金Ni9Pt 團(tuán)簇比其周圍合金團(tuán)簇的化學(xué)穩(wěn)定性強(qiáng)，這與具有相同原子總數(shù)的純Ni10團(tuán)簇相反.

2.3 團(tuán)簇的磁性

團(tuán)簇的磁矩是宏觀之下磁力的來(lái)源，為此作者計(jì)算了團(tuán)簇的自旋總磁矩、鉑原子電荷和鉑原子自旋磁矩.

從圖5 的團(tuán)簇自旋總磁矩隨尺寸增加的變化規(guī)律可以看出，NixPt（x=1?9）團(tuán)簇的自旋總磁矩隨著原子總數(shù)N=x+1 的增多迅速增大，最后趨近飽和值7.792μB.與純NixPt（x=1?9）團(tuán)簇的自旋總磁矩相比較看出，合金NixPt 團(tuán)簇的自旋總磁矩與純Nix+1團(tuán)簇的走勢(shì)相似，合金NixPt 團(tuán)簇的自旋總磁矩總體比純Nix+1團(tuán)簇的小，但隨著原子總數(shù)N的增多最后趨于相等，且純Nix+1團(tuán)簇的自旋總磁矩在N=6時(shí)達(dá)到飽和，而合金NixPt團(tuán)簇是在N=7才達(dá)到飽和.

圖5 團(tuán)簇自旋總磁矩隨尺寸增加的變化規(guī)律

從圖6 的Pt 原子電荷與自旋磁矩隨尺寸增加的變化規(guī)律可以看出，隨著NixPt（x=1?9）團(tuán)簇原子總數(shù)N=x+1 的增多，Pt 原子獲得電子逐漸增加，Pt 原子自旋磁矩在振蕩過(guò)程中急劇減小，說(shuō)明隨著原子總數(shù)N的增多，NixPt團(tuán)簇中用于成鍵的電子逐漸增多，當(dāng)N=5、10 時(shí)，Pt原子獲得電子達(dá)到較大值、Pt 原子自旋磁矩有較小值，表明Ni4Pt、Ni9Pt 團(tuán)簇的穩(wěn)定性較好，這與平均束縛能部分分析的結(jié)果吻合的很好.

圖6 Pt原子電荷與自旋磁矩隨尺寸增加的變化規(guī)律

3 結(jié)論

利用廣義梯度（GGA）近似下的BLYP關(guān)聯(lián)函數(shù)對(duì)NixPt（x=1?9）團(tuán)簇及純NixPt（x=1?9）團(tuán)簇的結(jié)構(gòu)穩(wěn)定性與磁性進(jìn)行系統(tǒng)的研究，得到以下結(jié)論：隨著團(tuán)簇原子總數(shù)N的增加Pt原子的摻雜對(duì)純Ni團(tuán)簇的幾何結(jié)構(gòu)改變不大，但是，平均束縛能卻改變?cè)S多，相對(duì)于純團(tuán)簇而言，Pt原子的摻雜使所有研究尺寸團(tuán)簇的穩(wěn)定性都得到增強(qiáng)；NixPt（x=1?9）團(tuán)簇的HOMO?LUMO 能隙隨著原子總數(shù)N=x+1的增多出現(xiàn)振蕩行為；隨著團(tuán)簇總原子數(shù)的增加，合金NixPt團(tuán)簇的自旋總磁矩總體比純Nix+1團(tuán)簇的小，但最后都趨于相近，而純Nix+1團(tuán)簇的自旋總磁矩在N=6時(shí)達(dá)到飽和，而合金NixPt團(tuán)簇是在N=7才達(dá)到飽和，整個(gè)過(guò)程中，Pt原子得到的電子數(shù)逐漸增大.