〈100〉間隙型位錯(cuò)環(huán)在純鎢及含氦雜質(zhì)鎢(010)表面下運(yùn)動(dòng)行為的分子動(dòng)力學(xué)模擬*

秦夢(mèng)飛 王英敏 張紅玉 孫繼忠?

1) (大連理工大學(xué)物理學(xué)院,大連 116023)

2) (大連理工大學(xué)材料科學(xué)與工程學(xué)院,大連 116023)

1 引言

鎢具有極高的熔點(diǎn)、高熱導(dǎo)率、良好的高溫強(qiáng)度、低濺射率和氘/氚滯留等優(yōu)異性能,是未來(lái)聚變裝置中首選的面向等離子體壁材料[1-3].在聚變環(huán)境下,高通量的中子和氦、氫粒子輻照在鎢中形成大量空位和間隙原子,二者聚集形成的位錯(cuò)環(huán)是鎢壁材料中一種主要位錯(cuò)形式[4].輻照會(huì)產(chǎn)生氦泡,位錯(cuò)環(huán)的形核和穩(wěn)定與氦泡密切相關(guān)[5,6].因此,位錯(cuò)環(huán)的存在同時(shí)影響材料的力學(xué)性能[7]和氫/氦滯留[8,9].研究位錯(cuò)環(huán)的動(dòng)力學(xué)行為是理解和評(píng)價(jià)鎢壁材料服役性能的重要基礎(chǔ).

基于位錯(cuò)晶體幾何學(xué)特征,純鎢的間隙型位錯(cuò)環(huán)主要有〈100〉和1/2〈111〉兩種類(lèi)型[10-12].低能輻照實(shí)驗(yàn)表明[13],在T=773 K鎢的淺表層內(nèi)〈100〉位錯(cuò)環(huán)大約占總位錯(cuò)體積分?jǐn)?shù)的50%,在T=1073 K仍然高達(dá) 35%—50%.這表明聚變正常工作溫度范圍內(nèi)〈100〉位錯(cuò)環(huán)是一類(lèi)主要的點(diǎn)陣缺陷.同時(shí),由于氦離子主要沉積在鎢的淺表層 ≤1 μm深度范圍[14],〈100〉型位錯(cuò)環(huán)和氦離子相互作用不可避免.Yi等[13,15]對(duì)比研究了純鎢在150 keV和2 MeV鎢離子輻照后的微結(jié)構(gòu)后發(fā)現(xiàn),兩種樣品近表層中位錯(cuò)環(huán)的數(shù)目密度不同,后者幾乎是前者的兩倍.實(shí)驗(yàn)上,由于透射電子顯微觀察的薄膜樣品易于造成位錯(cuò)弛豫和組態(tài)重構(gòu),難以原位觀察近表層位錯(cuò)環(huán)在映像力作用下的運(yùn)動(dòng).前期研究主要集中于數(shù)值模擬: Groves和Bacon[16]基于連續(xù)介質(zhì)的彈性理論研究了位錯(cuò)環(huán)的旋轉(zhuǎn),計(jì)算了位錯(cuò)環(huán)的能量和線(xiàn)張力,獲得了位錯(cuò)環(huán)的穩(wěn)定晶體學(xué)取向;Fikar等[17]研究了鎢自由表面處〈100〉和1/2〈111〉型棱柱位錯(cuò)環(huán)的深度與映像力關(guān)系,發(fā)現(xiàn)尺寸較大的位錯(cuò)環(huán)較易逸出表面.這些研究主要針對(duì)近表面處位錯(cuò)環(huán)靜態(tài)的性質(zhì),例如位錯(cuò)彈性應(yīng)力場(chǎng)、映像力和能量,沒(méi)有涉及間隙型位錯(cuò)環(huán)在映像力作用下的運(yùn)動(dòng)行為.

本文針對(duì)鎢近表面〈100〉型間隙型位錯(cuò)環(huán),運(yùn)用分子動(dòng)力學(xué)方法從原子層次上研究在 (010) 表面映像力作用下,位錯(cuò)環(huán)遷移過(guò)程中取向和應(yīng)力場(chǎng)隨深度、溫度的變化,并研究了氦間隙原子對(duì)與位錯(cuò)環(huán)的相互作用.

2 計(jì)算細(xì)節(jié)

利用大型分子動(dòng)力學(xué)模擬軟件Lammps開(kāi)展模擬工作.在鎢的BCC晶格中填充 [010] 型啞鈴間隙原子對(duì)構(gòu)建位錯(cuò)環(huán).考察了廣泛采用的兩個(gè)勢(shì)函數(shù): W-He[18]和W勢(shì)函數(shù)[19],計(jì)算得到[010]型啞鈴間隙原子對(duì)的形成能分別為10.29 eV和12.72 eV.由W-He勢(shì)函數(shù)得出的形成能介于參考值9.88 eV[20]和11.63 eV[21]之間,更適合于本文[010]間隙位錯(cuò)環(huán)的模擬計(jì)算.該勢(shì)函數(shù)為W-He體系的嵌入型相互作用勢(shì),包含Ackland的W-W相互作用勢(shì)、Beck的He-He相互作用勢(shì)以及Juslin的W-He相互作用勢(shì)三部分.

模擬單元中選取 [100],[010],[001] 分別為x,y和z方向,其中x和z方向設(shè)置為周期性邊界條件,y方向設(shè)置為固定邊界條件.模擬單元尺寸為16 nm ×18 nm ×16 nm.為構(gòu)建自由表面,將y方向上y＞16 nm的部分設(shè)置為真空,同時(shí)將y＜1 nm部分的鎢原子固定,防止粒子漂移.位錯(cuò)環(huán)的產(chǎn)生通過(guò)刪除y方向上中心處近表面的兩層原子,并將其替換為 [010] 啞鈴型間隙原子對(duì)來(lái)構(gòu)建(b⊥慣習(xí)面).采用最快速下降法進(jìn)行能量最小化,獲得力學(xué)平衡態(tài).在NVT系綜下,設(shè)定時(shí)間步長(zhǎng)為1 fs,分別在400,800 和1200 K下模擬半徑1 nm位錯(cuò)環(huán)的遷移(本模擬實(shí)驗(yàn)的位錯(cuò)環(huán)半徑都為1 nm),整個(gè)模擬時(shí)長(zhǎng)為5 ns.當(dāng)位錯(cuò)環(huán)遷移出自由表面后,則停止模擬.指定位錯(cuò)環(huán)環(huán)心代表環(huán)的位置,b⊥n代表柏氏矢量垂直于表面法線(xiàn)方向的位錯(cuò)環(huán),b//n代表柏氏矢量平行于表面法線(xiàn)方向的位錯(cuò)環(huán).

不同深度位錯(cuò)環(huán)的形成能可用分子動(dòng)力學(xué)模擬求出,通過(guò)(1)式計(jì)算映像力:

式中,σ是單位長(zhǎng)度位錯(cuò)環(huán)的映像力;b和d分別是位錯(cuò)環(huán)的柏氏矢量大小和直徑;E(h)=Ef(h)-E0,其中Ef(h) 是h深度位錯(cuò)環(huán)的形成能,E0是位錯(cuò)環(huán)位于模擬單元中心的形成能,這里Ef(h)=Eloop-Ep,Eloop是體系包含位錯(cuò)環(huán)條件下的能量,Ep是完美晶體的能量.

本文基于最簡(jiǎn)單的聚集方式——氦原子對(duì),分析氦與位錯(cuò)環(huán)的相互作用.位錯(cuò)環(huán)和氦原子對(duì)的結(jié)合能為

其中EB是位錯(cuò)環(huán)與 [112] 方向的氦間隙原子對(duì)的結(jié)合能,Etotal是體系中包含位錯(cuò)環(huán)和氦間隙氦原子對(duì)的總能量(假設(shè)位錯(cuò)環(huán)和氦原子距離足夠遠(yuǎn)),Eloop+2He是位錯(cuò)環(huán)與隙氦原子對(duì)發(fā)生相互作用時(shí)體系的總能量.

3 結(jié)果與討論

3.1 間隙位錯(cuò)環(huán)取向的影響

隨機(jī)追蹤一個(gè)b//n型位錯(cuò)環(huán)(T=800 K,h=5 nm)的動(dòng)態(tài)演化過(guò)程(圖1).t＜0.170 是位錯(cuò)環(huán)運(yùn)動(dòng)的第1階段,發(fā)生慣習(xí)面翻轉(zhuǎn) (t為約化時(shí)間,t=tsim/t0,tsim為模擬時(shí)間,t0=2500 ps).具體地,初始時(shí)刻(t0=0.001)慣習(xí)面為 (010) 面(圖1(a)),且b⊥慣習(xí)面,位錯(cuò)環(huán)開(kāi)始在映像力作用下移動(dòng),向表面逃逸.在t=0.018 時(shí)慣習(xí)面變?yōu)?011)面(圖1(b)),隨后位錯(cuò)環(huán)下半部分開(kāi)始加快運(yùn)動(dòng)速度,下半部分與上半部分位置逐漸持平,慣習(xí)面又變回 (010) 面(圖1(c)和圖1(d)).t＞0.170是位錯(cuò)環(huán)運(yùn)動(dòng)的第2階段,即快速遷移階段.當(dāng)t=0.239 時(shí)慣習(xí)面轉(zhuǎn)變?yōu)?(0ˉ11) 面(圖1(e)),此時(shí)位錯(cuò)環(huán)距離自由表面最近處約1 nm.t=0.245 時(shí),位錯(cuò)環(huán)的一端先與自由表面接觸,位錯(cuò)環(huán)中的〈100〉片段開(kāi)始分解為 1/2[1ˉ1ˉ1],1/2[ˉ11ˉ1],1/2[ˉ1ˉ11]和 1/2[111] 位錯(cuò),形成位錯(cuò)網(wǎng)絡(luò)(圖1(g)),最終t=0.247時(shí),位錯(cuò)環(huán)被自由表面完全吸收(圖1(h)).

圖1 b//n 位錯(cuò)環(huán)的取向變化和分解(T=800 K,初始深度 h=5 nm) (a) t=0.001;(b) t=0.018;(c) t=0.115;(d) t=0.170;(e) t=0.239;(f) t=0.245;(g) t=0.246;(h) t=0.247;圖中粉紅顏色代表 〈100〉位錯(cuò),綠顏色代表 1/2〈111〉位錯(cuò),藍(lán)色箭頭表示柏氏矢量方向Fig.1.Orientational change and decomposition of a b//n dislocation loop at T=800 K and the initial depth of h=5 nm:(a) t=0.001;(b) t=0.018;(c) t=0.115;(d) t=0.170;(e) t=0.239;(f) t=0.245;(g) t=0.246;(h) t=0.247 . The pink loop stands for the 〈100〉 type dislocation,the green one for 1/2〈111〉.The blue arrows denote the direction of Burgers vector.

位錯(cuò)環(huán)向表面遷移過(guò)程中應(yīng)力變化情況如圖2所示.在表面映像力作用下,位錯(cuò)環(huán)各部分應(yīng)力場(chǎng)不均勻,導(dǎo)致其出現(xiàn)慣習(xí)面翻轉(zhuǎn)的現(xiàn)象.在高溫(高于鎢的德拜溫度約400 K)T=800 K,激活大量空位和原子交換,造成間隙原子向低能量的(110) 密堆原子面偏聚,〈100〉位錯(cuò)環(huán)分解為1/2[ˉ11ˉ1]與 1/2[111].位錯(cuò)環(huán)完全遷移出表面后,在表面形成局部凸起結(jié)構(gòu)(圖3).位錯(cuò)環(huán)與自由表面反應(yīng)前,它以刃位錯(cuò)的形式存在(圖4(a));而在二者反應(yīng)時(shí),位錯(cuò)以刃/螺的混合位錯(cuò)形式出現(xiàn)(見(jiàn)圖4(b)).位錯(cuò)環(huán)向自由表面移動(dòng)的過(guò)程中伴隨著內(nèi)應(yīng)力釋放(圖2(a)—(f)).體系釋放內(nèi)應(yīng)力有多種方式,例如 Gao等[22]提出位錯(cuò)環(huán)分解達(dá)到釋放內(nèi)部應(yīng)力的目的.本模擬發(fā)現(xiàn)在鎢近表層處,b//n間隙位錯(cuò)環(huán)通過(guò)取向變化和向外移動(dòng)釋放應(yīng)力.

圖2 b//n 位錯(cuò)環(huán)移動(dòng)過(guò)程中的應(yīng)力變化圖(T=800 K,初始深度 h=5 nm) (a)—(e) 位錯(cuò)環(huán)向表面移動(dòng)過(guò)程中的應(yīng)力圖,其中(a) t=0.014,(b) t=0.018,(c) t=0.066,(d) t=0.115,(e) t=0.239;(f) 位錯(cuò)環(huán)完全移出表面后的應(yīng)力圖,t=0.247Fig.2.Stress variation of a b//n dislocation loop at T=800 K and initial depth of h=5 nm in motion: (a)-(e) The stress of the dislocation loop as it moves toward the surface ((a) t=0.014,(b) t=0.018,(c) t=0.066,(d) t=0.115,(e) t=0.239);(f) the stress after the dislocation loop has been completely removed from the surface,t=0.247.

圖3 b//n 位錯(cuò)環(huán)與表面相互作用后的表面形貌Fig.3.Surface morphology after the interaction between b//n dislocation loop and surface has ceased to occur.

圖4 b//n 位錯(cuò)環(huán)演化過(guò)程中位錯(cuò)類(lèi)型變化圖 (a) 位錯(cuò)環(huán)未與表面發(fā)生作用時(shí)的狀態(tài);(b) 位錯(cuò)環(huán)與表面發(fā)生強(qiáng)烈作用時(shí)的狀態(tài).圖中藍(lán)色為刃位錯(cuò),紅色為螺位錯(cuò)Fig.4.Variation of the dislocation geometry of a b//n dislocation loop in motion: (a) The geometry of the dislocation loop away from the surface;(b) the geometry change of the dislocation loop while interacting with the surface.Blue represents the edge dislocation,red represents a screw dislocation.

圖5是b⊥n型位錯(cuò)環(huán)(T=800 K;h=2,3,5 nm)在近表面處的演化過(guò)程.圖5(a)和圖5(b)顯示,t=2 時(shí)位錯(cuò)環(huán)(h=3,5 nm)仍處于原來(lái)深度,并未發(fā)生位置變化.當(dāng)位錯(cuò)環(huán)距離自由表面更近時(shí),發(fā)生分解.圖5(c) 所示h=2 nm處的位錯(cuò)環(huán)在t=1.824時(shí)一部分轉(zhuǎn)變?yōu)?/2〈111〉片 段(圖5(d)),t=1.836 時(shí)位錯(cuò)環(huán)開(kāi)始被自由表面所吸收,形成復(fù)雜的位錯(cuò)網(wǎng)絡(luò)(圖5(e)),到t=2 時(shí)位錯(cuò)環(huán)完全逃出自由表面(圖5(f)).與b//n取向的情形類(lèi)似,b⊥n位錯(cuò)環(huán)移出表面的過(guò)程中也發(fā)生應(yīng)力釋放和刃螺位錯(cuò)轉(zhuǎn)變.

圖5 b⊥n型位錯(cuò)環(huán)的取向變化和分解(T=800 K) (a) t=0.018;(b) t=0.115;(c) t=0.170;(d) t=0.245;(e) t=0.246;(f) t=0.247;圖中粉紅顏色的為 〈100〉 型位錯(cuò),綠顏色的為 1/2〈111〉 型位錯(cuò),底部物理量代表位錯(cuò)環(huán)初始深度Fig.5.Orientational change and decomposition of a b⊥n dislocation loop at T=800 K: (a) t=0.018;(b) t=0.115;(c) t=0.170;(d) t=0.245;(e) t=0.246;(f) t=0.247.The pink loop stands for the 〈100〉 dislocation,the green part for 1/2〈111〉,the physical quantity at the bottom represents the initial depth of the dislocation loop.

對(duì)比5 nm深度范圍內(nèi)b//n和b⊥n型位錯(cuò)環(huán)的運(yùn)動(dòng),發(fā)現(xiàn)前者更容易移出表面,后者只有在距離表面小于 2 nm才發(fā)生逃逸.因?yàn)榍罢呦虮砻娴倪w移可以通過(guò)位錯(cuò)滑移來(lái)實(shí)現(xiàn),而后者只能通過(guò)位錯(cuò)攀移來(lái)實(shí)現(xiàn),攀移比滑移更加困難.

3.2 溫度和深度的影響

圖6描述的是位錯(cuò)環(huán)形成能Ef(h) 與距離表面深度h的關(guān)系.位錯(cuò)環(huán)(b//n和b⊥n)的形成能均隨深度增大呈增大趨勢(shì),在h＞5 nm時(shí)趨于恒定值172.5 eV.在h＜5 nm的范圍內(nèi),b//n位錯(cuò)環(huán)的形成能為160.0—172.5 eV,而b⊥n位錯(cuò)環(huán)的形成能為 169.0—172.5 eV.前者形成能對(duì)深度的負(fù)梯度遠(yuǎn)大于后者,后者更加穩(wěn)定.圖7是位錯(cuò)環(huán)映像力σ隨深度h的變化情況.映像力隨深度減小而增大,在h＞5 nm時(shí),兩種位錯(cuò)環(huán)的映像力都很小,可以忽略;當(dāng)h＜5 nm時(shí),b//n位錯(cuò)環(huán)映像力的變化范圍是0—790 MPa,b⊥n位錯(cuò)環(huán)映像力的變化范圍是0—120 MPa.對(duì)比可知b//n位錯(cuò)環(huán)更易受到表面效應(yīng)的影響,而b⊥n位錯(cuò)環(huán)在材料中更穩(wěn)定,其他取向位錯(cuò)環(huán)的穩(wěn)定性介于它們之間.

圖6 位錯(cuò)環(huán)(b//n和b⊥n)形成能隨距離表面深度的變化Fig.6.Variation of the formation energies of the dislocation loops (b//n and b⊥n) with increasing depth.

圖7 位錯(cuò)環(huán)(b//n和b⊥n)映像力隨距離表面深度的變化Fig.7.Image forces of the dislocation loops (b//n and b⊥n) with increasing depth.

針對(duì)溫度因素,在T=400,800,1200 K時(shí)各生成30個(gè)含有位錯(cuò)環(huán)(b//n)的鎢樣本,模擬得到相應(yīng)的位錯(cuò)環(huán)在鎢內(nèi)部的滯留時(shí)間(圖8),t400K＞t800K＞t1200K.由T=可知,系統(tǒng)溫度越高位錯(cuò)環(huán)遷移速度越大,越容易克服其束縛勢(shì)壘遷移出表面,進(jìn)而使得位錯(cuò)環(huán)在鎢中的滯留時(shí)間變短.這里mi是原子質(zhì)量,vi是原子速度,N是原子數(shù),kB是玻耳茲曼常數(shù)

圖8 不同溫度和深度下 b//n 位錯(cuò)環(huán)的滯留時(shí)間變化,圖中編號(hào)1—10,11—20 和21—30分別是初始深度h=3,5,7nm的位錯(cuò)環(huán)Fig.8.Retention time variation of the b//n dislocation loops with temperature and depth.The numbers 1-10,11-20,and 21-30 are used to specify the dislocation loops located at the initial depth of h=3,5,7 nm,respectively.

圖8還給出不同深度(h=3 nm (1—10號(hào)),h=5 nm (11—20號(hào)),h=7 nm (21—30號(hào)))b//n位錯(cuò)環(huán)在材料中的滯留時(shí)間,滯留時(shí)間隨深度增大呈增加趨勢(shì).計(jì)算時(shí),當(dāng)模擬時(shí)間達(dá)到2即5 ns時(shí)停止模擬,這樣所有滯留時(shí)間超過(guò)5 ns的位錯(cuò)環(huán)都標(biāo)記在tsim/t0=2 處.對(duì)比圖7,在h=3 nm時(shí)位錯(cuò)環(huán)映像力大于h=5,7 nm情形.由此推知離自由表面越近,位錯(cuò)環(huán)所受映像力越大,越容易脫離鎢點(diǎn)陣的束縛,環(huán)滯留時(shí)間變短.當(dāng)h=5 nm時(shí),位錯(cuò)環(huán)受到的映像力很小,但發(fā)生漲落,使得它在初始位置附近微小移動(dòng),導(dǎo)致位錯(cuò)環(huán)以一定概率移動(dòng)到映像力較大的作用區(qū)域.此外,演化過(guò)程中環(huán)形狀會(huì)發(fā)生改變,產(chǎn)生額外的力.這些因素最終使部分位錯(cuò)環(huán)也能成功向表面處移動(dòng).對(duì)于b⊥n型位錯(cuò)環(huán)而言,它很難逃逸出表面,但是溫度升高易于發(fā)生位錯(cuò)類(lèi)型的轉(zhuǎn)變(詳見(jiàn)附錄A).

3.3 氦原子對(duì)位錯(cuò)環(huán)運(yùn)動(dòng)的影響

為分析氦與位錯(cuò)環(huán)之間的相互作用,計(jì)算了位錯(cuò)環(huán)(b//n和b⊥n)與方向Ⅰ和方向Ⅱ上氦間隙原子對(duì)的結(jié)合能(圖9).當(dāng)b//n時(shí),沿方向Ⅰ的結(jié)合能由4.5 eV逐漸減小至約0 eV(圖9(a)).可見(jiàn)距離位錯(cuò)環(huán)越近,氦間隙原子對(duì)越容易被吸收,吸收范圍約是3a＜r＜ 6a(r指氦原子距位錯(cuò)環(huán)中心的距離,a為晶格常數(shù)).方向Ⅱ上的結(jié)合能隨深度增加而減小,其中表面第1層原子下的值最大(圖9(b)),在此方向上b//n位錯(cuò)環(huán)排斥氦間隙原子對(duì).對(duì)于b⊥n型位錯(cuò)環(huán)而言,得到的結(jié)果與前面不同.沿方向Ⅰ,結(jié)合能從 -3 eV升至約0 eV(圖9(c)),在此方向表現(xiàn)為排斥氦原子,排斥范圍約為0 ＜r＜ 4a.沿方向Ⅱ,b⊥n位錯(cuò)環(huán)對(duì)氦間隙原子對(duì)呈現(xiàn)吸收作用,距離表面最近處結(jié)合能最大(圖9(d)).

圖9 位錯(cuò)環(huán)與氦間隙原子對(duì)的結(jié)合能 (a),(b) b//n 位錯(cuò)環(huán)分別在方向Ⅰ和方向Ⅱ上的結(jié)和能;(c),(d) b⊥n 位錯(cuò)環(huán)分別在方向Ⅰ和方向Ⅱ上的結(jié)和能,loop1—loop5 分別對(duì)應(yīng) h=1.5,2,2.5,3,3.5 nm的位錯(cuò)環(huán).(a),(c) 橫坐標(biāo)是氦原子對(duì)距離位錯(cuò)環(huán)的距離;(b),(d) 橫坐標(biāo)是氦原子對(duì)距表面的深度Fig.9.Binding energy of the dislocation loop with a pair of helium atoms: (a),(b) The binding energies in direction Ⅰ and Ⅱ of b//n dislocation loops;(c),(d) the binding energies for the b⊥n dislocation loop.loop1-loop5 correspond to the dislocation loops located at h=1.5,2,2.5,3,3.5 nm.The abscissa of (a) and (c) represents the distances between the pair of helium atoms and the dislocation loop,and the ordinate of (b) and (d) is for the depth.

圖10 給出了不同條件下b//n位錯(cuò)環(huán)在鎢近表層的滯留時(shí)間統(tǒng)計(jì)結(jié)果.在T=800 K,h=3 nm和T=1200 K,h=3,5 nm下分別有7/10,10/10,6/10的位錯(cuò)環(huán)移出自由表面,其他均滯留于材料內(nèi)部.與圖8相比,當(dāng)有氦間隙原子對(duì)存在時(shí),位錯(cuò)環(huán)在鎢內(nèi)的滯留時(shí)間明顯增加,同時(shí)能遷移出表面的數(shù)量也明顯變少.

圖10 不同溫度下 b//n 位錯(cuò)環(huán)的滯留時(shí)間圖(體系中含氦原子對(duì)) 圖中編號(hào)1 — 10,11—20和21—30分別是初始深度h=3,5,7 nm的位錯(cuò)環(huán)Fig.10.Variation of the retention time with temperature and depth for the b//n dislocation loops existing in the W lattice cell containing a pair of helium atoms.The numbers 1-10,11-20,and 21-30 are used to specify the dislocation loops located at the initial depth of h=3,5,7 nm,respectively.

在鎢近表層,b//n比b⊥n位錯(cuò)環(huán)更易向自由表面逃逸,且逃逸深度更大,后者傾向于釘扎在鎢材料內(nèi)部.無(wú)論哪種位錯(cuò)環(huán)都會(huì)在有限的區(qū)域內(nèi)吸附表層的氦原子,而由于b⊥n位錯(cuò)環(huán)穩(wěn)定存在于材料內(nèi)部,更易為氦原子提供吸附位點(diǎn),導(dǎo)致材料淺層區(qū)域氦濃度增大,為以后生成高密度的氦團(tuán)簇提供條件[23],氦團(tuán)簇生長(zhǎng)為氦泡,向外遷移并破裂會(huì)產(chǎn)生絨毛[24],進(jìn)一步影響偏濾器的使用壽命.

4 結(jié)論

本文運(yùn)用分子動(dòng)力學(xué)方法,從原子尺度上模擬研究鎢近表面處納米位錯(cuò)環(huán)的動(dòng)態(tài)演化.模擬結(jié)果發(fā)現(xiàn):

1)對(duì)于真空自由表面映像力作用下的b//n和b⊥n位錯(cuò)環(huán),前者受力不均,易向表面逃逸,后者傾向滯留于鎢內(nèi)部.位錯(cuò)環(huán)遷移過(guò)程中發(fā)生慣習(xí)面翻轉(zhuǎn),伴隨內(nèi)應(yīng)力釋放和刃螺位錯(cuò)轉(zhuǎn)變.位錯(cuò)環(huán)完全逸出后在鎢表面形成局部凸起結(jié)構(gòu).

2)當(dāng)溫度接近于0 K且深度h＜5 nm時(shí),b//n和b⊥n位錯(cuò)環(huán)受到的映像力分別為0—790 MPa和0—120 MPa.在有限溫度下,b//n位錯(cuò)環(huán)在T≥800 K且h≤5 nm范圍逸出表面的概率超過(guò)90%,在其他情況下滯留率約達(dá)100%.b⊥n位錯(cuò)環(huán)在T=800 K且h≤2 nm范圍也發(fā)生遷移,深度增大位錯(cuò)環(huán)不易移動(dòng),但T＞800 K時(shí)部分〈100〉環(huán)片段會(huì)分解為 1/2〈111〉位錯(cuò).

3)當(dāng)鎢中存在氦間隙原子對(duì)時(shí),在h＞1.5 nm深度的位錯(cuò)環(huán)平面內(nèi)對(duì)氦原子對(duì)的最大吸收范圍是3a＜r＜ 6a,垂直于位錯(cuò)環(huán)平面方向上的最大排斥范圍是0 ＜r＜ 4a.對(duì)比無(wú)氦原子存在的情形,位錯(cuò)環(huán)在材料內(nèi)的滯留概率增大.其中T=800,1200 K且h=5 nm時(shí),位錯(cuò)環(huán)在材料內(nèi)滯留的概率分別增加70%和40%.b⊥n位錯(cuò)環(huán)受到氦間隙原子對(duì)制約也更難遷移出表面.另外,位錯(cuò)環(huán)和氦原子之間的相互作用造成鎢材料內(nèi)氦分布不均勻.

附錄A 溫度對(duì)b⊥n位錯(cuò)環(huán)運(yùn)動(dòng)的影響

圖A1給出了深度h=5 nm條件下,溫度對(duì)b⊥n位錯(cuò)環(huán)運(yùn)動(dòng)的影響.圖 A1(a)和圖A1(b) 表明在400 K,位錯(cuò)環(huán)類(lèi)型以及位置幾乎未改變,只有形狀發(fā)生了一點(diǎn)變化.當(dāng)T=800 K,雖然位錯(cuò)環(huán)的位置不變,但是在t=0.980 時(shí)位錯(cuò)環(huán)中有一小部分發(fā)生了〈100〉向1/2〈111〉位錯(cuò)的轉(zhuǎn)變(見(jiàn)圖A1(c)和圖A1(d)).當(dāng)T=1200 K時(shí),位錯(cuò)環(huán)中有更多部分向 1/2〈111〉位錯(cuò)轉(zhuǎn)變,并且轉(zhuǎn)換頻率比800 K時(shí)更高(見(jiàn)圖A1(e)—(h)).雖然b⊥n型位錯(cuò)環(huán)很難逃逸出表面,但是溫度升高易于發(fā)生位錯(cuò)類(lèi)型的轉(zhuǎn)變.

圖A1 不同溫度 b⊥n 位錯(cuò)環(huán)的運(yùn)動(dòng)演化過(guò)程(初始深度h=5 nm) (a) t=0.980;(b) t=2.000;(c) t=0.980;(d) t=2.000;(e) t=0.540;(f) t=0.960;(g) t=1.600;(h) t=2.000;圖中粉紅顏色的為 〈100〉 型位錯(cuò),綠顏色的為 1/2〈111〉 型位錯(cuò)Fig.A1.Orientational change and decomposition of a b⊥n dislocation loop at the initial depth of h=5 nm: (a) t=0.980;(b) t=2.000;(c) t=0.980;(d) t=2.000;(e) t=0.540;(f) t=0.960;(g) t=1.600;(h) t=2.000 . The pink loop stands for the 〈100〉 type dislocation,the green one for 1/2〈111〉.