物理截?cái)嗯c電子局域函數(shù)結(jié)合法研究Ge2Sb2Te5非晶態(tài)原子成鍵

王鑫洋+陳念科+王雪鵬+張斌+陳志紅+李賢斌+劉顯強(qiáng)

【摘 要】采用物理截?cái)嗯c電子局域函數(shù)相結(jié)合的方法，可以提高非晶態(tài)結(jié)構(gòu)分析的準(zhǔn)確性，并已在GeTe二元非晶結(jié)構(gòu)分析中得到應(yīng)用。本文將該方法應(yīng)用于三元系Ge2Sb2Te5合金非晶結(jié)構(gòu)的分析中，確定了分析Ge2Sb2Te5非晶結(jié)構(gòu)合理的物理截?cái)嚅L度為3.20，電子局域函數(shù)閾值為0.63。在此基礎(chǔ)上對(duì)Ge和Sb原子的配位數(shù)和鍵角分布等結(jié)構(gòu)信息進(jìn)行分析。結(jié)果顯示：Ge原子以3、4配位為主，對(duì)應(yīng)的局域結(jié)構(gòu)分別為缺陷八面體和四面體；Sb原子同樣以3、4配位為主，對(duì)應(yīng)的局域結(jié)構(gòu)主要為缺陷八面體。

【關(guān)鍵詞】物理截?cái)?電子局域函數(shù) 配位數(shù) 鍵角分布

【Abstract】The method of combining cutoff distance with electron localization function （ELF） can improve the accuracy of the analysis for the amorphous structure， and has been applied to analyse the amorphous structure of the phase change material Ge2Sb2Te5， determining the reasonable value of cutoff distance is 3.20 and the value of ELF Threshold is 0.63. On this basis， the Ge and Sb atomic coordination number and bond angle distribution of structural information are discussed. The results show that Ge atoms are in 3， 4 coordination， which the local structure correspond to the defect octahedral forms and tetrahedron forms； the Sb atoms are also in 3， 4 coordinaion， and the corresponding local structures are mainly the defect octahedral forms.

【Key words】cutoff distance； electron localization function； coordinate numbers；bond angle distribution

1 引言

非晶結(jié)構(gòu)具有短程有序，長程無序的特征，其結(jié)構(gòu)解析一直是困擾科學(xué)界的難題之一[1-3]。在研究非晶模型時(shí)，首先要確定模型中各原子的成鍵信息。目前在確定非晶態(tài)結(jié)構(gòu)中各原子成鍵情況時(shí)，最常用的是物理截?cái)喾椒?。該方法是以某一個(gè)原子為中心選取一個(gè)經(jīng)驗(yàn)長度，當(dāng)兩個(gè)原子之間的距離不超過該長度時(shí)就認(rèn)為這兩個(gè)原子成鍵[4，5]；反之，則認(rèn)為不成鍵。該方法由于受到人為因素影響，致使其確定非晶結(jié)構(gòu)中的原子成鍵信息與實(shí)際情況存在一定差異[6]。電子局域函數(shù)（Electron Localization Function，ELF）常用來分析共價(jià)系統(tǒng)中原子之間的成鍵情況，可以給出結(jié)構(gòu)模型中電子的局域化程度，并反映兩原子之間所成化學(xué)鍵的相對(duì)強(qiáng)度[7，8]。由于ELF方法是通過原子周圍電子的局域化程度來判斷兩原子是否成鍵，因此采用該方法分析非晶模型中各原子的成鍵情況要比物理截?cái)喾椒ǜ泳_。

相變存儲(chǔ)器被認(rèn)為是下一代非易失存儲(chǔ)技術(shù)的最優(yōu)解決途徑之一，具有廣闊的應(yīng)用和發(fā)展前景[9]。在相變材料中，GeTe-Sb2Te3系列合金是迄今研究最多、最成熟的一類[10，11]，該材料已經(jīng)廣泛的應(yīng)用于硬盤、DVD和藍(lán)光光盤中[12]。其工作原理在于利用該類材料在非晶態(tài)和晶態(tài)之間的快速相變，導(dǎo)致兩種狀態(tài)的光學(xué)和電學(xué)性質(zhì)發(fā)生巨大變化，以此來實(shí)現(xiàn)數(shù)據(jù)的快速擦寫與存儲(chǔ)[13]。

在我們先前的研究中，為更加準(zhǔn)確的研究相變材料非晶態(tài)中原子成鍵與結(jié)構(gòu)信息，采用物理截?cái)嗪虴LF相結(jié)合的方法對(duì)GeTe非晶態(tài)結(jié)構(gòu)進(jìn)行了研究。結(jié)果發(fā)現(xiàn)，總成鍵數(shù)曲線會(huì)隨著ELF閾值的增加而單調(diào)的減少，在減少過程中曲線出現(xiàn)了明顯的拐點(diǎn)，該拐點(diǎn)所對(duì)應(yīng)的ELF值就是分析GeTe非晶態(tài)時(shí)所需要的ELF閾值。并且通過ELF方法還可以確定最佳的物理截?cái)嚅L度。在本文中，我們將該方法應(yīng)用于三元系Ge2Sb2Te5合金非晶結(jié)構(gòu)的分析，給出了用于Ge2Sb2Te5非晶結(jié)構(gòu)分析的合理物理截?cái)嚅L度和ELF閾值，并進(jìn)一步對(duì)Ge和Sb原子的配位及局域結(jié)構(gòu)等信息進(jìn)行了分析。

2 實(shí)驗(yàn)方法

本文第一性原理計(jì)算采用基于密度泛函理論結(jié)合投影綴加平面波方法的Vienna ab-initio simulation package（VASP）軟件包[14]。其中投影綴加平面波用來描述離子實(shí)與價(jià)電子之間的相互作用，而電子與電子之間的交換關(guān)聯(lián)勢則采用廣義梯度近似（general gradient approximate，GGA）泛函進(jìn)行處理[15]。在模擬的過程中使用正則系綜（NVT）方法[6]，采用Nose恒溫法控制溫度[16]，并且僅在該體系的Γ點(diǎn)處進(jìn)行模擬。首先搭建Ge2Sb2Te5的晶體結(jié)構(gòu)模型[17]，該模型結(jié)構(gòu)為立方相，初始模型為333超胞，其中共包含194個(gè)原子（Ge：42個(gè)，Sb：44個(gè)，Te：108個(gè)）。先對(duì)其進(jìn)行結(jié)構(gòu)優(yōu)化，然后進(jìn)行分子動(dòng)力學(xué)模擬，選取截?cái)嗄転?20eV，每一個(gè)離子步的步長為3fs。將優(yōu)化后的Ge2Sb2Te5晶體瞬間升溫至3000K，并保溫9ps；隨后以20K/ps的降溫速率冷卻到300K，并在300K下保溫9ps使其充分弛豫，最終得到Ge2Sb2Te5非晶結(jié)構(gòu)模型[6]。

電子局域函數(shù)反映了電子在結(jié)構(gòu)中的局域化程度和相對(duì)化學(xué)鍵強(qiáng)度，其數(shù)學(xué)表述為[7，8]：

其中用來衡量Pauli排斥力，它表示找到其他相同自旋電子的概率密度。則代表在均勻的電子氣中具有相同局域自旋密度的值，并且也是的歸一化系數(shù)。ELF值是介于0到1之間無量綱的常數(shù)。當(dāng)ELF=0.5時(shí)代表均勻電子氣；當(dāng)該值大于0.5時(shí)，說明電子更加局域在某一個(gè)位置。通過VASP軟件計(jì)算，可以直接得到非晶結(jié)構(gòu)中每個(gè)格點(diǎn)的ELF值。為了計(jì)算兩相鄰原子之間的ELF值，在它們的連線上取21個(gè)點(diǎn)，以每個(gè)點(diǎn)為中心做一定半徑的圓柱形管，該點(diǎn)的ELF值為管內(nèi)所包含格點(diǎn)的ELF值的平均值。在21個(gè)點(diǎn)的數(shù)值中存在一個(gè)極小值，該數(shù)值可以反映出兩原子之間所成化學(xué)鍵的相對(duì)強(qiáng)度[6]。

3 實(shí)驗(yàn)結(jié)果與討論

為了確定分析Ge2Sb2Te5非晶態(tài)結(jié)構(gòu)時(shí)所需要的物理截?cái)嚅L度和ELF閾值，我們分別選取不同的物理截?cái)嚅L度和ELF閾值，對(duì)該結(jié)構(gòu)的成鍵情況進(jìn)行了分析。如圖1（a）所示，總成鍵數(shù)由物理截?cái)嚅L度和ELF閾值共同決定。隨著ELF閾值增大，總成鍵數(shù)單調(diào)下降，曲線在下降的過程中分為兩個(gè)階段。第一階段，當(dāng)ELF閾值小于0.63時(shí)，總成鍵數(shù)減小比較緩慢。第二階段，當(dāng)ELF閾值大于0.63時(shí)，成鍵數(shù)量快速減少。這是因?yàn)椋?dāng)ELF閾值小于0.63時(shí)，原子間強(qiáng)度較弱的化學(xué)鍵將不計(jì)入統(tǒng)計(jì)；而當(dāng)該值大于0.63時(shí)，會(huì)將一部分強(qiáng)度較強(qiáng)的真實(shí)存在的化學(xué)鍵濾除。利用我們?cè)诜治鯣eTe非晶態(tài)時(shí)所得到的結(jié)論：隨ELF數(shù)值增加而出現(xiàn)成鍵數(shù)快速下降的拐點(diǎn)就是分析材料結(jié)構(gòu)時(shí)所需要的ELF閾值。因此，0.63應(yīng)是分析Ge2Sb2Te5體系時(shí)較為合理的ELF閾值。

從圖1（b）可以看出，ELF取0.63時(shí)，雖然對(duì)于不同的物理截?cái)嚅L度會(huì)得到不同的總成鍵數(shù)目，但是當(dāng)物理截?cái)嚅L度處于3.20～3.30時(shí)，總成鍵數(shù)基本相同，而選取3.0～3.15時(shí)總成鍵數(shù)相對(duì)偏小。這說明，當(dāng)物理截?cái)嚅L度為3.0～3.15時(shí)，由于截?cái)嚅L度選擇過短會(huì)將一些真實(shí)的成鍵截?cái)?。而從圖1（a）可以看出，當(dāng)截?cái)嚅L度處于3.20～3.30時(shí)，總成鍵數(shù)幾乎只取決于ELF閾值，而與截?cái)嚅L度的選取無關(guān)。圖1（a）所示，相較于物理截?cái)?.20～3.30范圍內(nèi)的成鍵總數(shù)曲線，當(dāng)取值為3.20時(shí)，ELF閾值在0.58～0.63之間曲線的斜率最小，說明在該范圍內(nèi)隨ELF值增加而被截?cái)嗟脑渔I是相對(duì)最少的。換言之，此時(shí)采用物理截?cái)喾椒ㄅc采用ELF方法得到的成鍵分析結(jié)果是最為接近。因此，對(duì)于不含ELF數(shù)據(jù)的Ge2Sb2Te5非晶結(jié)構(gòu)模型，取物理截?cái)嚅L度為3.20時(shí)得到的結(jié)果將最接近真實(shí)的成鍵數(shù)。至此，通過采用物理截?cái)嗯cELF相結(jié)合的方法，精確的確定出了分析 Ge2Sb2Te5非晶態(tài)結(jié)構(gòu)時(shí)所需要的物理截?cái)嚅L度為3.20，ELF閾值為0.63。

通過比較Ge2Sb2Te5與GeTe，其ELF閾值相同，而Ge2Sb2Te5物理截?cái)嚅L度卻略大（≈0.15），其原因是Ge2Sb2Te5的亞穩(wěn)相結(jié)構(gòu)為NaCl結(jié)構(gòu)，其中Te原子占據(jù)一套晶格，Ge、Sb和空位隨機(jī)分布在另一套晶格上[18]。正是由于Sb原子的加入，以及空位含量的增加，導(dǎo)致了在分析Ge2Sb2Te5的非晶態(tài)時(shí)需要的物理截?cái)嚅L度要大于分析GeTe非晶態(tài)時(shí)所需要的物理截?cái)嚅L度。

在物理截?cái)嚅L度取3.20的情況下，我們進(jìn)一步對(duì)比分析了ELF閾值分別取0.58和0.63時(shí)Ge2Sb2Te5非晶結(jié)構(gòu)模型中Ge和Sb原子的配位情況（配位數(shù)和鍵角分布）。圖2給出了物理截?cái)嚅L度為 3.20，ELF閾值分別為0.58和0.63時(shí)Ge和Sb原子的配位數(shù)。從圖2（a）中可看出，當(dāng)ELF閾值取0.58時(shí)Ge原子的配位數(shù)主要集中在3、4配位，Sb原子的配位數(shù)集中在3、4以及5配位。從圖2（b）中可以看出，當(dāng)ELF閾值取0.63時(shí)，Ge原子的配位數(shù)仍以3、4配位為主，但增加了少量2配位，這些2配位的出現(xiàn)是由于截?cái)嗔?配位和4配位中部分強(qiáng)度較弱的化學(xué)鍵。一般的，Sb原子在化合物中的常見價(jià)態(tài)為+3和+5價(jià)。如圖2所示，當(dāng)ELF閾值增大時(shí)，3配位的Sb原子所占的比例明顯增加，4配位的Sb原子所占的比例減少，而5配位的Sb原子基本保持不變。同樣，2配位的Sb原子的出現(xiàn)是由于截?cái)嗔?配位和4配位中強(qiáng)度較弱的化學(xué)鍵。綜上所述，當(dāng)物理截?cái)噙x取3.20，ELF閾值取0.63時(shí)所得到的結(jié)果更加趨向于真實(shí)的原子成鍵情況。

接下來，我們將對(duì)該模型中Ge和Sb原子的鍵角分布進(jìn)行分析。如圖3（a）所示，當(dāng)物理截?cái)嚅L度取3.20，ELF閾值取0.63時(shí)，3配位的Ge原子的鍵角分布主要集中在90°附近，說明其局域結(jié)構(gòu)為缺陷八面體結(jié)構(gòu)；而4配位的Ge原子的鍵角分布主要集中在109°附近，說明其局域結(jié)構(gòu)為四面體[18]。對(duì)于Sb原子，3、4配位Sb原子的鍵角分布都集中在90°附近，說明3、4配位的Sb原子的局域結(jié)構(gòu)均為缺陷八面體。如圖3（b）所示，當(dāng)ELF閾值從0.58增大至0.63時(shí)，雖然截?cái)嗔艘恍┗瘜W(xué)鍵，但是總體上Ge和Sb原子的鍵角分布并沒有發(fā)生明顯改變。這也從側(cè)面證明了通過增加ELF閾值所截?cái)嗟哪切┗瘜W(xué)鍵是一些強(qiáng)度相對(duì)較弱的化學(xué)鍵。正是由于這些被截掉的化學(xué)鍵之間的兩個(gè)原子的相互作用力很弱，才沒有對(duì)Ge和Sb原子的局域結(jié)構(gòu)產(chǎn)生明顯的影響。

4 結(jié)語

我們將物理截?cái)嗯cELF相結(jié)合的方法應(yīng)用于三元系Ge2Sb2Te5非晶態(tài)結(jié)構(gòu)的分析中，確定了在Ge2Sb2Te5非晶結(jié)構(gòu)中判斷成鍵的合理?xiàng)l件為：物理截?cái)嚅L度取3.20，ELF閾值取0.63。對(duì)Ge和Sb原子配位數(shù)和其鍵角分布等結(jié)構(gòu)信息分析，發(fā)現(xiàn)Ge原子以3、4配位為主，其中3配位的Ge原子主要以缺陷八面體形式存在，而4配位則主要以四面體的形式存在；Sb原子同樣以3、4配位為主，而其局域結(jié)構(gòu)均為缺陷八面體。采用物理截?cái)嗪虴LF相結(jié)合的方法對(duì)Ge2Sb2Te5非晶結(jié)構(gòu)進(jìn)行分析，不但驗(yàn)證了采用該方法分析非晶態(tài)結(jié)構(gòu)所得結(jié)果的合理性，同時(shí)也證明了該方法可以廣泛的適用于其他具有共價(jià)性質(zhì)的非晶態(tài)材料的研究中。

參考文獻(xiàn)：

[1]Zallen R 1983 The Physics of Amorphous Solids （New York： Wiley） pp10-16.

[2]Ziman J M 1979 Models of Disorder： The Theoretical Physics of Homogeneously Disordered Systems （Cambridge ：Cambridge University Press） pp51-56.

[3]Yonezawa F， Ninomiya T 1983 Topological Disorder in Condensed Matter （Berlin：Springer） pp30-39.

[4]Akola J， Jones R O 2007 Phys. Rev. B 76 235201.

[5]Xu M， Cheng Y Q， Wang L， Sheng H W， Meng Y， Yang W G， Han X D， Ma E 2012 Proc. Natl. Acad. Sci. U. S. A. 109 E1055-62.

[6]Xu M， Cheng Y Q， Sheng H W， Ma E 2009 Phys. Rev. Lett. 103 195502.

[7]Hughbanks T， Hoffmann R 1983 J. Am. Chem. Soc. 105 3528-37.

[8]Silvi B， Savin A 1994 Nature 371 683-6.

[9]Ovshinsky S R 1968 Phys. Rev. Lett. 21 1450-3.

[10]Yoon S M， Choi K J， Lee N Y， Jung S W， Lee S Y， Park Y S， Yu B G， Lee S J， Yoon S G 2008 J. Electrochem. Soc. 155 H421-H5.

[11]Wang K， Steitner C， Warnwangi D， Ziegler S， Wuttig M， Tomforde J， Bensch W 2007 Microsyst. Technol. 13 203-6.

[12]Welnic W， Wuttig M 2008 Mater. Today 11 20-7.

[13]Wuttig M， Yamada N 2007 Nature Mater. 6 824-32.

[14]Kresse G， Furthmuller J 1996 Comput. Mater. Sci. 6 15-50.

[15]Ernzerhof M， Scuseria G E 1999 J. Chem. Phys. 110 5029-36.

[16]Tuckerman M， Berne B J， Martyna G J 1992 J. Chem. Phys. 97 1990-2001.

[17]Yamada N， Matsunaga T 2000 J. Appl. Phys. 88 7020-7028.

[18]Caravati S， Bernasconi M， K hne T D， Krack M， Parrinello M 2007 Appl. Phys. Lett. 91 171906.