張鳳春， 唐黎明， 李元喆, 吳沙麗， 李春福， 文 平， 羅 強(qiáng)

(1.成都工業(yè)學(xué)院 材料與環(huán)境工程學(xué)院， 成都 611730; 2.西南石油大學(xué) 材料科學(xué)與工程學(xué)院， 成都 610500; 3.西南石油大學(xué) 理學(xué)院， 成都 610500)

1 引 言

硫化氫酸性腐蝕一直是石油天然氣工業(yè)中的重點(diǎn)難點(diǎn)問(wèn)題，這主要是由于其往往能導(dǎo)致硫化物應(yīng)力腐蝕開(kāi)裂和氫致開(kāi)裂.其中，應(yīng)力腐蝕開(kāi)裂的危害性和破壞性極大，常常在無(wú)任何征兆的情況發(fā)生，嚴(yán)重時(shí)甚至能造成災(zāi)難性事故[1].人們圍繞其開(kāi)展了大量研究工作[2-8]，鐵素體奧氏體雙相不銹鋼的H2S應(yīng)力腐蝕開(kāi)裂研究表明，鐵素體相的應(yīng)力腐蝕開(kāi)裂敏感性較高，滲入的氫主要在鐵素體相中擴(kuò)散，而在奧氏體相中則主要存在于氫陷阱中，這說(shuō)明鐵素體相是雙相不銹鋼中的“弱相”[4,5].因此，開(kāi)展H在α-Fe中的擴(kuò)散行為研究是很有必要的.而眾所周知，合金元素與合金設(shè)計(jì)息息相關(guān)，合金元素的添加將直接影響鋼材的力學(xué)性能和化學(xué)性能，在雙相不銹鋼中，合金元素的含量更是高達(dá)30%以上，其對(duì)鋼材性能的影響更是不言而喻.因此，本文采用基于密度泛函理論的第一性原理方法，研究對(duì)比了H在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)中的擴(kuò)散行為，以期從電子結(jié)構(gòu)層次分析合金元素對(duì)α-Fe氫擴(kuò)散行為的影響作用，從而進(jìn)一步探討合金元素對(duì)H2S應(yīng)力腐蝕行為的影響，這對(duì)于完善硫化氫應(yīng)力腐蝕機(jī)制和抗H2S應(yīng)力腐蝕新材料的設(shè)計(jì)具有重要的科研價(jià)值.

2 計(jì)算模型與方法

2.1 計(jì)算模型

采用材料模擬軟件Materials Studio中的Materials Visualizer模塊首先建立了α-Fe晶體的晶胞結(jié)構(gòu).α-Fe晶胞屬于四方晶系IM-3M空間群，體心立方結(jié)構(gòu)，點(diǎn)陣參數(shù)a=b=c= 2.8664 ?，α=β=γ=90°[11]，如圖1(a)所示.并利用超晶胞方法，建立了p(1×1×2)的α-Fe晶體超晶胞模型.已有研究表明[12]，H原子在α-Fe中優(yōu)先占據(jù)四面體間隙位，因此建立了H原子位于四面體間隙位的α-Fe晶胞模型，如圖1(b)所示.而Cr、Mo、Ni、Cu、Mn等合金元素在α-Fe中的最優(yōu)占位研究明，其分別優(yōu)先占據(jù)頂角、體心、體心、頂角、體心、體心位置[13]，進(jìn)而建立了H原子位于四面體間隙位的α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞模型圖，如圖2(a)、2(b)、2(c)、2(d)、2(e)所示.

圖1 α-Fe晶胞模型圖

圖2 H原子位于α-Fe-Me-H(Me = Cr、Mo、Ni、Cu、Mn)四面體間隙的晶胞模型圖

2.2 計(jì)算方法

本文所有研究工作均采用材料模擬軟件Materials Studio中的CASTEP[14]模塊完成，探討了合金原子α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞中Me與H原子的相互作用，建立了擴(kuò)散系數(shù)與溫度間的函數(shù)關(guān)系式，并對(duì)比分析了合金元素Cr、Mo、Ni、Cu、Mn對(duì)α-Fe晶體中氫擴(kuò)散行為的影響作用.計(jì)算中電子波函數(shù)采用平面波基組展開(kāi)，交換關(guān)聯(lián)能函數(shù)采用廣義梯度近似(GGA)[15,16]下的RPBEP泛函形式，平面波截?cái)嗄苓x取為450 eV，α-Fe晶體結(jié)構(gòu)優(yōu)化的k點(diǎn)取樣采用Monkhorst-Pack[17]方法，選取為9×9×5.

3 計(jì)算結(jié)果與討論

3.1 合金原子與H原子的相互作用

為了從微觀電子層次探討合金元素與H原子的相互作用，本文分別得到了各合金體系的分波態(tài)密度圖(PDOS，Partial Density of States)，如圖3-圖7所示.從圖3(a)、4(a)、5(a)、6(a)、7(a)中可以看出，在-11.17 eV ～ -7.18 eV區(qū)域，F(xiàn)e的3d和4s軌道均與H的1s軌道之間形成明顯的共振峰，這說(shuō)明在α-Fe-Me合金體系中，H、Fe原子間存在一定的相互作用.由圖3(b)可知，在-10.68 eV～ -7.31 eV區(qū)域，Cr的3d軌道和4s軌道分別與H的1s軌道形成較為明顯的共振峰，說(shuō)明Cr、H原子間存一定的鍵合作用.由圖4(b)可知，在-11.14 eV ～ -7.78 eV區(qū)域，Mo的4d、5s軌道分別與H的1s軌道形成明顯的共振峰，說(shuō)明Mo原子的4d、5s軌道與H原子的1s軌道間存在雜化共軛作用.由圖5(b)中可知，在-10.88 eV ～ -7.45 eV區(qū)域，Ni的3d、4s軌道分別與H的1s軌道形成共振峰，這表明Ni與H原子間也存在一定鍵合作用. 由圖6(b)中可知，在-11.04 eV～-7.92 eV區(qū)域，Cu的3d、4s軌道分別與H的1s軌道形成共振峰，說(shuō)明Cu、H原子間存相互作用.由圖7(b)中可知，在-10.69 eV～-7.25 eV區(qū)域，Mn的3d、4s軌道分別與H的1s軌道形成共振峰，說(shuō)明Mn、H原子間存在一定的相互作用.

圖3 α-Fe-Cr-H晶胞體系的分波態(tài)密度圖

圖4 α-Fe-Mo-H晶胞體系的分波態(tài)密度圖

圖5 α-Fe-Ni-H晶胞體系的分波態(tài)密度圖

圖6 α-Fe-Cu-H晶胞體系的分波態(tài)密度圖

圖7 α-Fe-Mn-H晶胞體系的分波態(tài)密度圖

3.2 合金元素對(duì)氫擴(kuò)散行為的影響

3.2.1擴(kuò)散過(guò)渡態(tài)

為研究合金元素對(duì)氫擴(kuò)散行為的影響作用，本文采用“Complete LST/QST”方法分別搜索了H原子在α-Fe和α-Fe-Me(Me=Cr，Mo，Ni，Cu，Mn)兩個(gè)最近鄰T間隙的擴(kuò)散過(guò)渡態(tài)，所得擴(kuò)散勢(shì)壘曲線如圖8所示，曲線中最高點(diǎn)與最低點(diǎn)的能量差值即擴(kuò)散所需克服的最高能壘，稱為擴(kuò)散勢(shì)壘，而該能量所對(duì)應(yīng)的晶胞結(jié)構(gòu)即過(guò)渡態(tài).

過(guò)渡態(tài)搜索結(jié)果表明，H原子在α-Fe和α-Fe-Me(Me=Cr，Mo，Ni，Cu，Mn)兩個(gè)最近鄰T間隙的擴(kuò)散過(guò)渡態(tài)均為H原子處于八面體間隙位置.由圖8可知，H原子在α-Fe、α-Fe-Cr、α-Fe-Mo、α-Fe-Ni、α-Fe-Cu、α-Fe-Mn晶胞中的擴(kuò)散勢(shì)壘分別為0.2007 eV、0.0529 eV、0.5327 eV、0.0776 eV、0.1066 eV、0.0536 eV，其中，α-Fe-Cr晶胞最低，α-Fe-Mo晶胞最大.同時(shí)，對(duì)比α-Fe和α-Fe-Me(Me=Cr，Mo，Ni，Cu，Mn)晶胞的擴(kuò)散勢(shì)壘，不難發(fā)現(xiàn)，Cr、Ni、Cu、Mn元素使得H原子的擴(kuò)散勢(shì)壘降低，氫擴(kuò)散過(guò)程將更為容易.而Mo元素的固溶則使得H原子的擴(kuò)散勢(shì)壘增大，氫擴(kuò)散過(guò)程的發(fā)生將變得更為困難，這說(shuō)明Mo元素在α-Fe中的固溶對(duì)于H在α-Fe中的擴(kuò)散具有一定的阻礙作用，這也將對(duì)鐵素體組織抗H2S腐蝕產(chǎn)生一定的積極作用.

圖8 H原子在各晶胞內(nèi)的擴(kuò)散勢(shì)壘曲線

3.2.2擴(kuò)散系數(shù)

為了進(jìn)一步理解合金元素Cr、Mo、Ni、Cu、Mn對(duì)α-Fe晶胞中H擴(kuò)散行為的影響作用，本文分別計(jì)算了H在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞中的擴(kuò)散系數(shù).已有研究表明，固體中雜質(zhì)原子的擴(kuò)散系數(shù)可通過(guò)式(1)進(jìn)行計(jì)算[18-24]：

(1)

式中，D為擴(kuò)散系數(shù)；L為躍遷距離；v0為雜質(zhì)原子的有效躍遷頻率，由于固體中雜質(zhì)原子的有效躍遷頻率一般為1013Hz數(shù)量級(jí)，因此本文將H原子的有效躍遷頻率值v0選為1013Hz.；NA為阿伏伽德羅常數(shù)，約6.0221 × 1023mol-1；Ea為擴(kuò)散激活能；R為理想氣體常數(shù)，約8.3141 J·mol-1·K-1；T為絕對(duì)溫度.

將各晶胞體系的v0、Ea、L值分別代入式(1)中，可得到H原子在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞中擴(kuò)散系數(shù)與溫度的關(guān)系表達(dá)式：

(2)

(3)

(4)

(5)

(6)

(7)

從上述表達(dá)式中可以看出，H原子在α-Fe晶胞中的擴(kuò)散激活能為19.4 kJ/mol，而在合金元素Cr、Mo、Ni、Cu、Mn固溶后，其擴(kuò)散激活能分別為5.1 kJ/mol、51.4 kJ/mol、7.5 kJ/mol、10.3 kJ/mol、5.2 kJ/mol，其直觀對(duì)比圖如圖9所示，由圖可知，合金元素的固溶對(duì)氫擴(kuò)散激活能的影響較大，Cr、Ni、Cu、Mn等元素使得氫擴(kuò)散激活能降低，而Mo元素卻大大提高了氫擴(kuò)散激活能，這說(shuō)明Mo元素的添加能在一定程度上阻礙氫擴(kuò)散行為，這對(duì)于抗H2S應(yīng)力腐蝕新材料的設(shè)計(jì)來(lái)說(shuō)是較為有益的.

圖9 各晶胞的氫擴(kuò)散激活能

為了更為直觀的對(duì)比H原子在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞中的擴(kuò)散行為，本文分別計(jì)算得到了在273.15 K ～ 1273.15 K溫度范圍內(nèi)，各晶胞H擴(kuò)散系數(shù)與溫度倒數(shù)關(guān)系的阿倫尼烏斯圖，如圖10所示.由圖可知，隨溫度升高，各晶胞的氫擴(kuò)散系數(shù)均呈增大趨勢(shì).同時(shí)，與α-Fe晶胞中的H擴(kuò)散系數(shù)相比，Cr、Ni、Cu、Mn元素使得H擴(kuò)散系數(shù)增大，而Mo元素的固溶使得H擴(kuò)散系數(shù)顯著降低，這說(shuō)明Mo元素的添加將使得氫擴(kuò)散激活能增大，阻礙H原子在晶格間隙位間的擴(kuò)散作用，從而防止合金鋼的氫致開(kāi)裂.因此，在鐵素體中適當(dāng)增加Mo元素的含量將有利于防止其發(fā)生開(kāi)裂行為.

4 結(jié) 論

本文采用基于密度泛函理論的第一性原理方法，研究了H原子在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)中的擴(kuò)散行為，分析了合金原子Me與H原子的相互作用，計(jì)算了H原子在α-Fe和α-Fe-Me晶胞T位間擴(kuò)散的過(guò)渡態(tài)和擴(kuò)散系數(shù)，討論了合金原子Me對(duì)H原子在α-Fe晶胞中的氫擴(kuò)散行為的影響作用，得到的主要結(jié)論如下：

(1)在α-Fe-Me晶胞中，F(xiàn)e的3d、4s軌道分別與H的1s軌道存在相互作用，而合金原子Me與H原子也存在一定的相互作用，并且這種作用主要由Me價(jià)電子的d、s軌道與H的1s軌道雜化共軛所貢獻(xiàn)；

(2)合金原子Mo的固溶增大了α-Fe晶胞的氫擴(kuò)散激活能，進(jìn)而顯著降低了氫擴(kuò)散系數(shù)，因此在鋼材中適當(dāng)添加Mo元素將有利于防止鐵素體中氫致開(kāi)裂的發(fā)生，這對(duì)于抗H2S應(yīng)力腐蝕的新材料設(shè)計(jì)具有一定的現(xiàn)實(shí)意義.