趙 毅 王 棟 徐晨皓 吳江桅 王 洋 謝耀平

(1. 中國核動力研究設(shè)計院 反應(yīng)堆燃料及材料重點實驗室，四川 成都 610213； 2. 上海大學(xué)材料研究所，上海 200072；3. 上海大學(xué)計算機學(xué)院高性能計算中心，上海 200444)

鋯合金是重要的核燃料包殼材料，耐腐蝕性是其關(guān)鍵性能指標(biāo)，而鋯合金氧化膜的性質(zhì)決定著鋯合金耐腐蝕性的好壞，因此研究鋯合金氧化膜的性質(zhì)具有重要意義。許多商用鋯合金屬于Zr- Sn- Nb系。Zr- Sn- Nb系合金氧化膜主要由單斜二氧化鋯(monoclinic ZrO2, m- ZrO2)組成，還存在少量四方二氧化鋯(tetragonal ZrO2, t- ZrO2)以及Sn和Nb的氧化物，如SnO、SnO2、NbO2和Nb2O5等。Preuss等[1- 2]的研究表明：鋯合金氧化膜中不同位置t- ZrO2含量不同，距離金屬/氧化膜(oxide/metal, O/M)界面越遠，t- ZrO2含量越少；還通過電子背散射衍射研究了2種不同Sn含量鋯合金(ZIRLO和Zr- 1.0Nb- 0.1Fe)氧化膜的相組成，發(fā)現(xiàn)氧化膜大部分物相為m- ZrO2，ZIRLO合金氧化膜中t- ZrO2含量比Zr- 1.0Nb- 0.1Fe合金的多。Wei等[3]研究了Sn含量對鋯合金氧化膜相組成的影響，發(fā)現(xiàn)，在腐蝕過程中高Sn含量合金O/M界面形成了較多的t- ZrO2，隨著腐蝕過程的進行，ZrO2向外擴展，形成了較多的m- ZrO2，而低Sn含量合金形成了較少的t- ZrO2。Frankel等[4]研究發(fā)現(xiàn)，在腐蝕過程中，ZIRLO合金、A- 0.6Sn合金以及A- 0.0Sn合金氧化膜中t- ZrO2的質(zhì)量分?jǐn)?shù)分別為4%～8%，4%～7%，2%～5%。

鋯合金的耐腐蝕性與Sn含量密切相關(guān)，Sn含量降低，合金的耐腐蝕性提高[2]。Sn在鋯合金中一般不以第二相的形式析出，通常偏聚在合金表面，不同初始取向鋯合金表面Sn的偏聚程度也不同。Sn在氧化膜中有單質(zhì)Sn、SnO2、SnO等存在形式[5- 6]。Takeda等[7]通過透射電子顯微鏡觀察含Sn鋯合金的氧化膜，發(fā)現(xiàn)m- ZrO2附近晶界存在SnO。

Nb在鋯合金氧化膜中主要有2種存在形式：一是固溶于α- Zr或β- Zr；二是以第二相β- Nb存在。β- Nb是體心立方(body- centered cubic，BCC)結(jié)構(gòu)，其氧化速率小于基體α- Zr[8]。β- Nb會氧化為非晶或四方相的二氧化鈮(tetragonal NbO2, t- NbO2)[9- 10]，或氧化為單斜相的五氧化二鈮(monoclinic Nb2O5, m- Nb2O5)[11]。添加質(zhì)量分?jǐn)?shù)為0.1%的Fe的鋯合金氧化膜中會出現(xiàn)Laves結(jié)構(gòu)的第二相Zr(Nb,Fe)2，在距O/M界面0.8 μm處Laves相開始氧化為非晶態(tài)，距O/M界面1.5 μm處β- Nb開始氧化，這表明β- Nb具有更好的抗氧化性；Nb在靠近O/M界面為+2價，而靠近氧化物/水界面為+5價[12]。

目前有關(guān)Zr- Sn- Nb合金氧化膜性質(zhì)的研究報道很多，但對氧化膜中Zr、Sn、Nb氧化物的形成與演化機制尚未完全解釋清楚。因此，本文首先基于第一性原理計算了幾種氧化物的形成焓，研究不同O化學(xué)勢下氧化物的相對穩(wěn)定性；然后結(jié)合熱力學(xué)計算獲得不同溫度和壓力下O化學(xué)勢，得到氧化物的穩(wěn)定性與壓力和溫度之間的關(guān)系；最后系統(tǒng)地確定了Zr、Sn、Nb氧化物形成的熱力學(xué)驅(qū)動力隨溫度和壓力的變化關(guān)系。

1 計算方法與模型構(gòu)建

1.1 第一性原理計算方法

利用基于密度泛函理論(density functional theory, DFT)的第一性原理計算方法(first- principles method)，對Zr、Sn、Nb氧化物的相結(jié)構(gòu)總能進行計算，所有計算均采用vienna ab- initio simulation package(VASP)軟件包實現(xiàn)[13]。采用投影綴加波(projector- augmented wave, PAW)描述[14- 15]離子與價電子之間的相互作用，采用廣義梯度近似(generalized gradient approximation, GGA)下的Perdew- Burke- Ernzerhof(PBE)泛函處理電子間交換關(guān)聯(lián)。k點網(wǎng)絡(luò)采用Monkhorst- Pack方法產(chǎn)生[16]。計算過程中，對超胞的矢量和所有原子都進行弛豫，能量收斂標(biāo)準(zhǔn)為0.1 eV/nm。

1.2 參數(shù)測試

對平面波展開的能量截斷進行計算。對于單質(zhì)Zr，以密排六方(hexagonal close-packed, HCP)- Zr為對象，計算體系總能，其能量截斷取100～450 eV，計算結(jié)果如圖1(a)所示。當(dāng)能量截斷大于150 eV時，繼續(xù)增加能量截斷使總能的變化小于0.01 eV/atom。進一步測試能量截斷對Zr- O體系形成焓的影響。根據(jù)O2、HCP- Zr和ZrO2在不同能量截斷下的總能，計算m- ZrO2的形成焓，結(jié)果列于圖1(b)。由圖1(b)可以看出，當(dāng)能量截斷達到500 eV時，繼續(xù)增加能量截斷，形成焓的變化小于0.01 eV/atom。

圖1 HCP- Zr總能與m- ZrO2形成焓隨能量截斷的變化

采用Monhkorst- Packk點取樣方法對布里淵區(qū)進行積分，根據(jù)HCP- Zr(a=0.323 nm,c=0.515 nm)晶胞形狀，分別取(2×2×1)、(4×4×2)、(6×6×4)、(8×8×6)、(10×10×6)、(12×12×8)、(14×14×8)等k網(wǎng)格點，測得不同k點網(wǎng)格密度下的總能如圖2所示。由圖2可以看出，當(dāng)k點網(wǎng)格密度達到(6×6×4)時，繼續(xù)增加網(wǎng)格密度，總能變化小于0.01 eV/atom。因此，參照HCP- Zr體系的結(jié)果，本文所有k點取樣均按照密度相等的原則進行設(shè)置。

圖2 HCP- Zr總能隨k點網(wǎng)格密度的變化

2 計算結(jié)果與討論

2.1 典型氧化物的晶體結(jié)構(gòu)與形成焓

Zr- Sn- Nb系合金氧化膜中典型氧化物有m- ZrO2、t- ZrO2、t- SnO、t- SnO2、NbO2和m- Nb2O5，氧化物原胞如圖3所示。計算選取截斷能為500 eV，k點網(wǎng)格分別為7×7×7、9×9×7、2×2×1、1×1×1、1×1×2、1×1×1。通過晶格優(yōu)化獲得的氧化物晶格常數(shù)如表1所示，可知m- ZrO2和t- ZrO2的晶格常數(shù)計算值與試驗值吻合較好。Sn和Nb的氧化物的晶格常數(shù)與m- ZrO2相差較大，因此氧化物界面為非共格關(guān)系，界面能對微結(jié)構(gòu)的影響較大。形成焓(Hf)可用來反映氧化物的穩(wěn)定性，其計算公式[17]為：

表1 Zr- Sn- Nb系合金氧化膜中氧化物的晶格常數(shù)a, b, c及晶胞矢量夾角α, β, γ

圖3 Zr- Sn- Nb系合金氧化膜中氧化物原胞(紅色小球代表O原子，其余小球代表金屬原子)

(1)

式中：Etot為晶胞總能；EM(M=Zr,Nb,Sn)和EO分別為金屬元素M和O基態(tài)時的平均單個原子能量，基態(tài)下Zr為HCP結(jié)構(gòu)，Nb為BCC結(jié)構(gòu)，Sn為面心立方(face- centered cubic, FCC)結(jié)構(gòu)，O為氧分子結(jié)構(gòu)；ntot為晶胞中原子總數(shù)；nM和nO分別為晶胞中M、O原子的個數(shù)。

根據(jù)式(1)，通過第一性原理計算得到晶胞的總能以及M、O元素在其平衡態(tài)晶體結(jié)構(gòu)中的平均單個原子能量，結(jié)果列于表2。由表2可知，m- ZrO2的形成焓最低，SnO的形成焓最高，同種金屬不同類型氧化物的形成焓差值較小，均小于0.21 eV/atom。不同金屬氧化物的形成焓差值較大，Zr氧化物的形成焓比Nb氧化物的約低1.38 eV/atom，而Nb氧化物的形成焓比Sn氧化物的約低1.1 eV/atom。根據(jù)形成焓的計算結(jié)果可大致判斷出：Zr的氧化驅(qū)動力最大，Nb的次之，Sn最難氧化，其中Nb的氧化慢于Zr，與試驗結(jié)果相吻合[18]。在NIST- JANAF數(shù)據(jù)中NbO2和Nb2O5的形成焓分別為-2.73和-2.79 eV/atom，比計算值約大0.5 eV/atom，表明DFT對NbO2和Nb2O5中Nb—O化學(xué)鍵有所低估，這與著名O—O化學(xué)鍵低估問題類似[19]。

表2 由第一性原理計算得到的Zr- Sn- Nb合金氧化膜中氧化物的形成焓

2.2 氧化物穩(wěn)定性與O化學(xué)勢的關(guān)系

金屬氧化物的形成焓是反映其穩(wěn)定性的關(guān)鍵物理量，與從環(huán)境中獲取O原子的難易程度有關(guān)(即O化學(xué)勢)，隨O化學(xué)勢的變化而變化。氧化物形成焓的計算公式為：

(2)

式中μO為氧的化學(xué)勢。

圖4為氧化膜中Zr、Sn和Nb氧化物的形成焓隨O化學(xué)勢的變化?？梢钥闯?，O化學(xué)勢越高，金屬越容易從環(huán)境得到氧。比較發(fā)現(xiàn),價態(tài)相同的氧化物，如m- ZrO2、t- ZrO2、t- SnO2、t- NbO2的形成焓差值不隨化學(xué)勢變化而變化，而價態(tài)不同的氧化物形成焓的差值則隨化學(xué)勢變化而變化。

圖4 Zr- Sn- Nb合金氧化膜中Zr、Sn和Nb氧化物的形成焓隨O化學(xué)勢的變化

由于O化學(xué)勢變化不會引起m- ZrO2、t- ZrO2、t- SnO2、t- NbO2形成焓的變化，因此也不會引起Zr、Nb和Sn氧化順序的變化。由此可以推測，試驗t- ZrO2相成分的改變不是O化學(xué)勢變化引起的，而是其他原因，如壓力、晶粒大小等引起的。O化學(xué)勢變化會引起NbO2和Nb2O5的相對穩(wěn)定性改變，富氧區(qū)和貧氧區(qū)對應(yīng)的O/M界面?zhèn)群蚈/W水側(cè)的產(chǎn)物也有所變化，O/M界面附近穩(wěn)定相是NbO2，而O/W附近則為Nb2O5，這與堿性水中氧化膜的試驗結(jié)果相吻合[18]，也與氧化性相對弱的去離子水中出現(xiàn)NbO2，而氧化相對較強的堿性水常出現(xiàn)Nb2O5的結(jié)論相一致[20]。進一步推測出，O/M界面附近Sn氧化物穩(wěn)定相是SnO，而O/W附近則為SnO2。因此，根據(jù)上述結(jié)果可以定性地分析氧化膜不同位置的相種類。

2.3 溫度與壓力對O化學(xué)勢和氧化物相對穩(wěn)定性的影響

通過熱力學(xué)公式計算特定溫度和壓力下的O化學(xué)勢，得到一定溫度和壓強下氧化物的相對穩(wěn)定性，計算公式為：

(3)

(4)

對于腐蝕試驗，常用的2個水化學(xué)條件為360 °C/18.6 MPa去離子水以及400 ℃/10.3 MPa過熱水蒸氣。計算得到10.3和18.6 MPa壓力，300～500 ℃(573～773 K)溫度下O化學(xué)勢分別為-7.36～-7.22和-7.38～-7.14 eV/atom。圖5為在10.3 MPa和300～500 ℃(573～773 K)條件下，O化學(xué)勢隨溫度的變化以及氧化物形成焓隨O化學(xué)勢的變化。可以看出，NbO2和Nb2O5在O化學(xué)勢為-7.01 eV/atom時的形成焓相等，在低于-7.01 eV/atom的貧氧側(cè)，NbO2更穩(wěn)定，而在高于-7.01 eV/atom的富氧側(cè)，Nb2O5更穩(wěn)定，這與O越多，Nb2O5越容易生成是一致的。O化學(xué)勢為-7.01 eV/atom時，對應(yīng)10.3 MPa、500 ℃以上的化學(xué)勢。因此，從熱力學(xué)角度分析，在10.3 MPa、300～500 ℃(673～773 K)，Nb的氧化物穩(wěn)定相是NbO2，這與部分試驗得到的去離子水中Nb氧化物為NbO2一致[20]，但與另一部分試驗得到的Nb氧化物為Nb2O5不一致[18]，這可能是試驗條件復(fù)雜，部分條件沒有在計算中體現(xiàn)所致。SnO和SnO2在O化學(xué)勢為-7.12 eV/atom時的形成焓相等，在低于-7.12 eV/atom的貧氧側(cè)，SnO更穩(wěn)定，在高于-7.12 eV/atom的富氧側(cè)，SnO2更穩(wěn)定。O化學(xué)勢為-7.12 eV/atom時，對應(yīng)10.3 MPa、500 ℃以上的化學(xué)勢。因此，在10.3 MPa、300～500 ℃(573～773 K)，Sn的氧化物穩(wěn)定相是SnO，這能解釋Takeda等[7]的試驗結(jié)果；在18.6 MPa、300～500 ℃(573～773 K)，Nb和Sn氧化物的穩(wěn)定相仍分別是NbO2和SnO。

圖5 10.3 MPa壓力下O化學(xué)勢隨溫度的變化以及氧化物形成焓隨O化學(xué)勢的變化

圖6為在327 ℃(600 K)和10～20 MPa條件下，O化學(xué)勢隨壓力的變化。由圖6可知，600 K下，10～20 MPa壓力對應(yīng)的O化學(xué)勢為-7.262～-7.244 eV/atom。O化學(xué)勢隨壓力的增加而升高，壓力越大，金屬越容易從系統(tǒng)中獲得O原子,因此壓力越大，鋯合金越容易被氧化。NbO2和Nb2O5在O化學(xué)勢為-7.01 eV/atom時的形成焓相等，對應(yīng)的壓力大于20 MPa；SnO和SnO2在O化學(xué)勢為-7.12 eV/atom時的形成焓相等，對應(yīng)的壓力也大于20 MPa。在600 K和10～20 MPa壓力范圍內(nèi)，SnO比SnO2穩(wěn)定，NbO2比Nb2O5穩(wěn)定。因此，在此壓力范圍內(nèi)穩(wěn)定氧化物的種類不隨壓力而改變，氧化膜中Nb和Sn的氧化物穩(wěn)定相仍是NbO2和SnO。

圖6 600 K下O化學(xué)勢隨壓力的變化

3 結(jié)論

(1)Zr- Sn- Nb合金氧化膜中，m- ZrO2的形成焓最小，最穩(wěn)定，這與試驗得到的氧化膜主要物相為m- ZrO2相一致。氧化物穩(wěn)定性由強到弱的順序為Zr氧化物、Nb氧化物和Sn氧化物，氧化膜中典型析出相β- Nb的氧化速率小于Zr基體，與試驗結(jié)果相吻合。

(2)Nb和Sn氧化物的穩(wěn)定性隨O化學(xué)勢的變化而變化，貧氧區(qū)NbO2和SnO比較穩(wěn)定，富氧區(qū)Nb2O5和SnO2比較穩(wěn)定，這與試驗中NbO2常出現(xiàn)于氧化性較弱的去離子水環(huán)境中，而Nb2O5常出現(xiàn)于氧化性較強的堿性水中的結(jié)論相一致。

(3)在10.3和18.6 MPa壓力、300～500 °C(573～773 K)溫度條件下，Sn的穩(wěn)定相均為SnO，Nb的穩(wěn)定相均為NbO2，能解釋一些試驗結(jié)果。