TiN多型體高壓相變的第一性原理計算-物理學報

1、物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑員園圓栽蚤暈多型體高壓相變的第一性原理計算鄢顧搖雄搖高尚鵬覮淵復旦大學材料科學系袁上海 搖圓園園源猿猿冤淵圓園員年園愿月 猿園日收到 曰圓園員年園怨月 猿日收到修改稿冤搖搖基于密度泛函理論框架下的贗勢平面波方法袁計算了 月員淵氯化鈉結構冤堯月圓淵氯化銫結構冤堯月猿淵閃鋅礦結構 冤堯月噪淵六方氮化硼結構冤堯月澡淵碳化鎢結構 冤和 月愿員淵砷化鎳結構 冤遠種 栽蚤多暈型體的晶體結構堯體積彈性模量和相對穩(wěn)定性援計算指出 袁不存在 月源淵纖鋅礦 冤結構的栽蚤暈通援過不同外壓下的晶格弛豫得到每種結構的焓袁發(fā)現(xiàn)外壓低于猿

2、源緣鄖孕葬時 栽蚤最暈穩(wěn)定的結構是 月員結構袁在 猿源緣鄖孕葬以上 月圓結構的 栽蚤最暈穩(wěn)定 援栽蚤從暈 月員結構到 月圓結構的相變伴隨著體積和彈性模量的變化 援關鍵詞 院氮化鈦 袁贗勢袁高壓相變 袁密度泛函理論孕粵悅院雜苑員員援緣醞援遭袁遠源苑援園運援原袁遠圓圓援園鑿援藻遠袁員緣援園粵援澡纖鋅礦結構轉變的壓強為原圓員園郾鄖孕葬即袁不能通過加正的外壓發(fā)生 咱圓員暫員郾引言援栽蚤的暈高壓相變需要進一步的系統(tǒng)研究 袁來澄清在高壓下是否存在比暈葬悅結造構更栽蚤是暈一種新型多功能材料袁具有高強度 堯高硬穩(wěn)定的相 袁以及各亞穩(wěn)相在不同壓強下的相對穩(wěn)定度堯耐高溫 堯耐酸堿侵蝕 堯耐磨損等一系列優(yōu)點袁在關系

3、援機械工業(yè)中被用做材料表面的耐磨堯耐腐蝕 堯耐火本文采用贗勢平面波方法袁研究多種 栽蚤異暈構材料涂層和切削工具 堯力學部件上的硬質防磨損覆體的結構和性能以及在不同壓強下的相對穩(wěn)定性援層咱員要源暫主要研究的多型體包括月員淵暈葬悅結造構 冤堯月圓淵悅澤悅造援此外 栽蚤薄暈膜還在微電子業(yè)被用來做擴散結構冤堯月猿淵閃鋅礦結構冤堯月勢壘袁在汽車和玻璃制造業(yè)充當反射材料袁在珠寶袁澡鄄業(yè)用作表面的黃金色涂 咱猿袁源暫噪淵六方氮化硼結構月暈冤月堯澡淵宰悅結構冤以及 月愿員淵暈蚤粵結澤構 冤援為了表述援栽蚤的暈熱力學穩(wěn)定相為氯化納結構 袁由于 栽蚤被暈廣泛應用做耐磨材料和切簡潔袁下文中很多地方采用雜賊則怎噪賊怎

4、符則號遭藻淵則即蚤糟澡賊削工具 袁瞬時可能承受很大的應力袁所以對 栽蚤晶暈體月員袁月圓等冤來指代各種結構 援的各種可能存在的結構在高壓環(huán)境下的相對穩(wěn)定性和性能的研究具有實際的意義援基于密度泛函理論的第一性原理計算已廣泛用于材料高壓相變和高壓下的材料性能研究咱緣要員苑暫援在計算中可以很方便地設定不同的加壓條件袁如實驗室中金剛石壓腔可以實現(xiàn)的從醞孕葬到 員園圓鄖孕葬范圍的各種壓強 袁或者大星體內部以及極端化學條件下的 栽孕葬淵員園咱緣暫猿鄖孕葬量冤級的極高壓強援栽蚤的暈結構 堯電子結構 堯力學及表面等性質的密度泛函理論計算圓郾計算方法搖搖計算采用密度泛函理論框架下的贗勢平面波方法袁交換關聯(lián)能采用孕

5、月耘形式的 鄖鄖粵近似 咱圓圓暫援布里淵區(qū)采用醞燥灶噪澡燥則澤撒賊點鄄咱孕圓猿葬暫糟噪袁最大近鄰運點間距 為 園郾園源猿魡原員援選 擇 月云鄖淵雜月則燥贈鑿云藻造灶藻賊糟澡藻則咱圓源進暫行幾何優(yōu)化來確定零壓鄖燥造鑿棗雜葬澡則葬遭灶方灶法燥冤和不同外壓條件下的平衡結構和性質援本工作中所已有較多報道 咱員愿要圓園暫加外壓方式均為流體靜壓強援在給定外壓下袁各相袁但是高壓下栽蚤結暈構相變的第一性原理計算研究還很少袁我們只檢索到最近宰葬灶早的相對穩(wěn)定性可由焓勻越耘垣孕災來確定 袁焓最小的等人根據總能 鄄體積曲線推斷出栽蚤由暈 暈葬悅結造構到結構最穩(wěn)定 援所有計算工作均采用悅粵雜栽耘程孕序 咱圓緣暫鄢國家

6、自然科學基金淵批準號 院員園愿園源和園教員育愿部冤留學回國人員科研啟動基金資助的課題援覮通訊聯(lián)系人援耘鄄皂葬早蚤葬造燥院澤棗責怎岳鑿葬藻灶鑿援怎糟援灶訫圓園員中員國物理學會悅澡蚤灶孕藻澡澤贈藻澤蚤雜糟燥葬糟造蚤藻賊贈搖澡賊院賊轅轅責憎怎造蚤蚤曾責遭澡援葬贈糟援糟援灶園緣苑員搖員園圓鄄物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑員園圓完成援的總能差值和焓差值 袁兩組贗勢給出的結果差別在在贗勢平面波方法用于高壓研究時袁贗勢構造園郾園員源藻災以內 袁圖 員淵遭給冤出分別采用兩組贗勢計算過程中所采用的芯區(qū)截斷半徑則的值是一個需要的晶格常數 袁差別在 園郾園圓魡

7、以內 袁計算結果基本一糟特別注意的參量援由于所采用的壓強可能很高袁作致和預期相符 援在本論文以下計算中袁我們將統(tǒng)一為理論預測研究也不必受限于目前技術手段可能采取芯區(qū)截斷半徑較小的一組超軟贗勢援達到的壓強范圍袁在高壓下晶體的晶格常數變小袁近鄰原子距離變小袁而贗勢方法的基本概念則要求最相鄰原子如栽蚤和 暈的 則之和要等于或小于原子糟間距袁也就是近鄰原子的芯區(qū)不能重疊援程序贗勢庫所提供的贗勢的則取值一般是保證適用于常規(guī)糟壓強下的 袁所以對高壓下的適用性要進行測試援在所計算的各種結構的栽蚤中暈袁月猿栽蚤在暈同等壓強下近鄰 栽蚤原子和暈原子的間距最小援在 圓園園鄖孕葬時計算得到的近鄰的栽蚤原子與 暈原子

8、的距離則為 員郾苑圓苑魡袁小于 悅粵雜栽耘程孕序贗勢庫里的栽蚤原園的芯區(qū)半徑之和 袁則子與 暈原子超軟贗勢咱圓遠暫糟淵栽蚤垣冤則糟淵暈冤越員郾苑怨怨魡援這就意味著在此壓強條件下采用此芯區(qū)半徑的贗勢不再適用 袁必須修改 則糟淵栽蚤與冤則糟淵暈冤使袁研究壓強下的不同結構栽蚤的暈最近鄰原子距離 則都大于或等于則園淵栽蚤垣冤則淵暈冤援將 則淵栽蚤修冤糟糟糟改為 園郾怨園園魡袁則糟淵暈冤修改為園郾源圓猿魡袁即 則糟淵栽蚤垣冤則糟淵暈冤越員郾猿圓猿魡袁采用新的芯區(qū)半徑產生一組超軟贗勢援表 員的測試結果表明直到 員怨園鄖園孕葬六袁種結構的 栽蚤暈均滿足 則園躍則糟淵栽蚤垣冤則糟淵暈冤取援更小的芯區(qū)截斷半徑會

9、使贗勢的轉移性更好 袁但是也會使贗勢更野硬冶袁也就是平面波截斷能更高 袁經收斂條件測試后袁我們取平面波截斷能為 緣緣園藻災理援論上講 袁取兩種不同的芯區(qū)半徑所構造的贗勢在其適用范圍內所得到的結果應一致袁圖 員淵葬給冤出分別采用兩組贗勢 月員栽蚤和暈 月猿栽蚤在暈 員園園鄖孕葬范圍內的各壓強下圖 員搖月員淵氯化納結構 冤和 月猿淵閃鋅礦結構 冤栽蚤在暈員園園鄖孕葬外壓范圍內的 總能差及焓差淵葬冤和晶格常數淵遭冤援叉號代表采用悅粵雜栽耘程孕序超軟贗勢庫里的贗勢淵栽蚤和 暈的芯區(qū)截斷半徑分別為 園郾怨緣圓魡和 園郾愿源苑魡冤的計算結果袁三角符號代表采用新的較小的芯區(qū)截斷半徑淵栽和蚤 暈分別為 園郾怨

10、園園魡和 園郾源圓猿魡冤產生的贗勢的計算結果搖搖搖表 員搖遠種結構的 栽蚤在暈 園鄖孕葬和 員怨園鄖園孕葬下近鄰原子的距離 淵單位 魡冤月員月圓月猿月月月愿噪澡員園鄖孕葬圓郾員猿圓圓郾圓怨猿員郾怨怨愿圓員郾員怨圓郾員遠園圓郾員緣怨員怨園鄖園孕葬員郾緣猿緣員郾遠猿遠員郾猿猿怨員緣郾圓愿員郾緣遠員員郾遠員猿猿郾計算結果和討論猿郾員郾栽蚤暈異構體的結構和性能我們采用 苑種常見的晶體結構 袁即 月員淵暈葬悅結造構冤堯月圓淵悅澤結悅構造冤堯月猿淵閃鋅礦結構 冤堯月源淵纖鋅礦結構冤堯月噪淵澡鄄月結暈構冤堯月澡淵宰悅結 構 冤以 及月愿員淵暈蚤粵結澤構冤袁作為初始模型來進行晶格弛豫計算得到不同結構栽蚤的暈平

11、衡結構參數和體積彈性模量援計算中我們發(fā)現(xiàn)以纖鋅礦結構作為起始結構進行幾何優(yōu)化后 袁內部參數 怎淵對于理想纖鋅礦結構為猿轅愿冤變?yōu)?員轅圓援也就是說 栽蚤原子和 暈原子將在同一平面上袁每一層的 栽蚤袁原暈子相間排列成六角環(huán)狀網絡袁相鄰的原子層有沿糟軸方向的栽蚤要暈鍵袁這是類似于六方氮化硼的結構淵即 月噪冤援其他遠種結構的園緣苑員搖圓園圓鄄物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑員園圓栽蚤的暈結構參數和彈性模量在表圓中給出 援算工作參考 援月圓栽蚤暈的密度最大 袁而 月猿栽蚤的暈密度最小援具有高硬度是 栽蚤獲暈得廣泛應用的一個重要實驗報 道 的 月員栽

12、蚤暈晶 格常 數為源郾圓猿愿魡咱圓苑暫袁原因袁體積彈性模量是與硬度相關的一個重要參雜賊葬皂等責棗使造用 云孕蘊粵孕宰鄄鄖計鄖算粵 的 數 值 為量袁零壓下 月員栽蚤的暈體積彈性模量是最大的 袁而 月猿咱員愿暫孕月耘近源郾圓遠魡 袁醞葬則等造采燥用贗勢平面波方法和似計算的數值為 源郾圓猿愿魡咱員怨暫袁我們的理論計算得出栽蚤的暈體積彈性模量是最小的援由圖 圓給出 月員和晶格常數為 源郾圓源愿魡袁與實驗和其他理論計算結果的月圓栽蚤的暈體積彈性模量隨壓強的變化關系可以看差別符合一般認可的 鄖鄖粵近似計算精度 援其他結構出袁月圓栽蚤暈的體積彈性模量在低壓下略小于月員的 栽蚤作暈為可能存在的亞穩(wěn)相還沒有相

13、應的實驗結栽蚤暈隨袁著壓強的增大與月員栽蚤的暈差別變得更小 袁在構參數報道 袁我們的理論結果可供將來的實驗和計約 怨怨猿鄖孕葬以上大于月員栽蚤暈援表 圓搖遠種結構 栽蚤暈的結構參數和體彈模量結構空間群壓強 轅鄖孕葬晶格常數 轅魡密度 轅體彈模量 轅葬淵遭冤糟猿鄖孕葬早窯糟皂月員淵暈葬悅造冤云皂猿鄄皂園源援圓源愿源援圓源愿緣援猿園員圓苑員淵圓圓緣冤緣園源援園苑圓源援園苑圓遠援園愿苑源遠緣猿園園猿援遠苑怨猿援遠苑怨愿援圓緣源員圓愿源員園園園猿援圓怨苑猿援圓怨苑員員援源遠苑猿員怨愿月圓淵悅澤悅造冤孕皂猿鄄皂園圓援遠源愿圓援遠源愿緣援緣猿猿圓源苑淵圓圓員冤緣園圓援緣圓員圓援緣圓員遠援源員源源源園員園園

14、圓援源源圓圓援源源圓苑援園緣源遠員緣猿園園圓援圓遠愿圓援圓遠愿愿援愿員圓員圓源苑員園園園圓援園圓愿圓援園圓愿員圓援猿圓園猿圓園園月猿淵閃鋅礦 冤云源鄄猿皂園源援遠員猿源援遠員猿源援員愿遠員怨員淵圓員遠冤緣園源援猿源苑源援猿源苑緣援園園猿猿苑猿月孕遠猿轅皂皂糟園猿援緣猿遠源援圓園怨源援緣員園圓員猿噪淵澡鄄月暈冤淵員怨源冤緣園猿援猿猿猿源援園員苑緣援猿員苑源園員員園園猿援圓員源猿援愿怨苑緣援愿怨猿緣遠怨月孕遠鄄皂圓園圓援怨員源圓援苑園愿緣援員緣愿圓源苑澡淵宰悅冤淵員愿苑冤緣園圓援苑緣愿圓援遠園遠緣援怨愿愿源猿源員園園圓援遠緣怨圓援緣源猿遠援緣怨苑遠園圓員園園園圓援園怨愿圓援圓苑圓員員援愿遠愿圓怨源

15、園月愿孕遠猿轅皂皂糟園圓援怨愿圓緣援圓員猿緣援員員愿圓緣源員淵暈蚤粵澤冤淵員怨源冤緣園圓援愿圓苑緣援園員愿緣援怨員遠源源苑員園園圓援苑猿園源援愿怨愿遠援緣園園遠圓猿員園園園員援怨遠愿緣援園怨猿員圓援園猿猿圓遠園緣構的焓作參考值援圖 源給出了 遠種結構栽蚤的暈總能猿郾圓郾栽蚤暈各異構體外壓下的相對穩(wěn)定性和體積隨壓強的變化關系 袁從而可以分析內能和體在有外界給定壓強的條件下袁焓最低的結構最積變化分別對焓的相對大小改變的貢獻援穩(wěn)定援在壓強變化時 袁異構體的相對穩(wěn)定性可能會從圖 猿中可以看出 袁在無外壓和壓強較小的情有變化 援圖 猿給出了 遠種結構的 栽蚤在暈不同壓強下況下袁月員栽蚤焓暈最小所以是最穩(wěn)

16、定的結構袁這和 栽蚤暈的焓袁我們計算的壓強范圍從園到 員怨園鄖園孕葬為袁了最常見的相是 月員一致 援在無外壓時 袁焓值最大的是使圖線能更清楚的表達袁各壓強下的焓均取月員結月圓栽蚤暈但袁是隨著壓強的增大袁約 猿遠鄖孕葬以后 袁月圓園緣苑員搖猿園圓鄄物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑員園圓圖 圓搖月員與 月圓結構的栽蚤體暈積彈性模量隨壓強的變化圖 源搖遠種結構的栽蚤的暈總能 淵葬和冤體積 淵遭隨冤壓強的變化援為便于比較 袁圖示的能量和體積數據均對應于一個栽蚤單暈元 淵兩個原子冤圖 猿搖遠種結構的栽蚤的暈焓隨壓強的變化袁以 月員結構的焓作參考值 援

17、為便于比較袁圖示焓的數據均對應于一個栽蚤單暈元 淵兩個原子 冤栽蚤的暈焓將比 月猿結構的小 袁壓強大于 苑圓鄖孕葬時袁月圓結構的小 袁壓強大于 員遠圓鄖孕葬時袁栽蚤的暈焓將比月噪月圓栽蚤焓暈先將比 月 結構的小 袁接著在 員遠緣鄖孕葬將比澡月愿 結構的小 援壓強增加到猿源鄖緣孕葬以后 袁月圓結構的員栽蚤焓暈最小 袁所以 月圓栽蚤結暈構最穩(wěn)定 援在 猿源緣鄖孕葬左右袁栽蚤的暈月員結構和 月圓結構具有不同的體積 淵見圖猿源淵遭冤含冤有袁兩個原子的原胞體積分別為員圓園郾園魡和員員圓郾猿魡猿袁所以 栽蚤從暈 月員結構到 月圓結構的相變伴隨著體積的變化 援由于壓強孕越 原鑿耘轅鑿災為能量 鄄體積 淵耘鄄

18、災曲線冤的負斜率 袁在某個壓強孕下兩種結構的焓相等即 孕越原淵耘圓 原耘員冤轅淵災圓 原災員冤袁所以材料在高壓下的結構相變除了通過觀察圖猿所示的焓 鄄壓強曲線外袁利用兩種結構的 耘鄄曲災線也可以判斷轉變壓強 援圖 緣是 月員栽蚤和暈 月圓栽蚤的暈 耘鄄災曲線 袁共切線斜率的負值即轉變壓強約為圓郾員愿源藻災轅魡猿袁即 約猿緣園鄖孕葬援圖 緣搖月員淵暈葬結悅構造 冤栽蚤和暈月圓淵悅澤悅結造構 冤栽蚤的暈能量 鄄體積曲線援共切線 淵虛線 冤給出轉變壓強有實驗報道在園要源緣鄖孕葬的壓強范圍內袁栽蚤暈園緣苑員搖源園圓鄄物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑

19、員園圓樣品的體積隨著壓強的增加曲線變化比較平滑袁栽蚤暈 月員躍 月躍 月愿 躍月猿躍 月 躍月圓變?yōu)楦邏合?淵員圓遠園噪員澡沒有發(fā)生相變咱圓愿圓袁怨暫園到鄖孕葬以上冤的 月圓躍月抑月愿 躍月員躍月 躍月猿援援我們的計算結果表明在澡員噪猿源緣鄖孕葬范圍內 袁月員結構都是最穩(wěn)定的袁與該實驗結果一致 援只有在 猿源鄖緣孕葬以上壓強 袁月圓栽蚤才暈是最源郾結論穩(wěn)定的 援由于從 月員結構到 月圓結構的相變可能存在的勢壘 袁在實驗上要觀察到此相變應該在猿源緣鄖孕葬利用密度泛函理論框架下的贗勢平面波方以上某個更高的壓強處援悅葬韻袁雜和則悅韻鑿韻在高壓下法袁研 究 了 栽蚤暈的 月員淵暈葬悅結造構冤堯月圓淵悅

20、澤悅結造也會發(fā)生 月員結構到 月圓結構的相變 袁實驗上觀測到構冤堯月猿淵閃鋅礦結構冤堯月淵澡原月暈結構 冤堯月淵宰悅噪的轉變壓強分別為 遠園淵依圓冤鄖孕葬悅淵葬韻冤結構 冤和 月愿澡咱猿園暫淵暈蚤粵結澤構 冤遠種異構體的平衡結構袁猿淵遠依和 怨園遠郾鄖孕葬淵悅鑿韻冤員援為了保證在極高和在不同壓強下的相對穩(wěn)定性咱猿員暫咱猿圓暫源冤鄖孕葬淵雜則韻冤袁均小于 栽蚤暈壓下近鄰原子的距離大于或等于原子贗勢芯區(qū)半的 月員到 月圓結構的轉變壓強 援量子蒙特卡羅計算預測 醞早韻從 月員結構到 月圓結構的轉變壓強要接近徑之和 袁在計算中我們使用了特別構造的具有較咱猿猿暫月員結構轉變小芯區(qū)半徑的超軟贗勢援計算結果

21、給出了六種異遠園園鄖孕葬 援云藻韻隨著壓強的增大則從構體相對穩(wěn)定性隨壓強的變化袁在高壓下 月猿結構為一種菱形結構 袁在高溫高壓條件下 淵苑園鄖孕葬以上 袁員遠園運園冤可發(fā)生 月員結構到 月愿結構的轉變 咱猿源暫和 月 結構的 栽蚤都暈是焓明顯較其他四種結構高也援噪員月猿和 月 結構的 栽蚤在暈較低壓強下比月圓栽蚤暈就是說相對較不穩(wěn)定袁而 月愿 結構和 月 結構的焓噪員澡焓更小 袁而與 月 結構以及 月愿 結構較為接近袁隨著在低壓下比月員結構的焓略大 袁隨著壓強的增大分澡員壓強的增大 袁在 猿遠鄖孕葬以后 月猿栽蚤焓暈變?yōu)樽畲?袁在別在 員園員鄖園孕葬和 員圓遠鄖園孕葬的極高壓強下焓會小苑圓鄖孕

22、葬以后 月 栽蚤的暈焓僅次于月猿栽蚤而暈大于其他于 月員結構 袁但會大于高壓下的月圓結構 袁即不存在噪源種結構的 栽蚤暈從援圖 源淵遭冤可以看出 袁月猿結構和 月焓最小的壓強范圍援月員栽蚤暈直至 猿源緣鄖孕葬都具有噪結構的 栽蚤的暈體積一直是最大的兩種袁內能隨壓強最小的焓 袁也就是結構最穩(wěn)定袁在 猿源緣鄖孕葬以上袁的迅速增加是導致它們焓超過月圓栽蚤的暈主要原因 援月圓栽蚤焓暈最小 援從 月員結構到 月圓結構的相變是體由圖 猿可以明顯看到袁在較高壓強下月猿栽蚤暈和 月積不連續(xù)相變 援月圓栽蚤的暈彈性模量在低壓下略小噪于同等壓強下的月員栽蚤暈而袁在約 怨怨猿鄖孕葬以上略栽蚤的暈焓明顯高于其他四種相

23、袁因此屬于很難在高大于 月員栽蚤暈我援們的計算結果表明栽蚤在暈承受較壓下存在的相 援月 與 月愿 結構的 栽蚤暈的焓在研究壓強范圍內大應力的切削等應用中袁當等效壓強小于 猿源緣鄖孕葬澡員十分接近 袁在外壓較小時 袁略大于 月員栽蚤暈隨援著壓強的范圍內 袁應不需擔心結構相變問題援本論文對于增大袁在大約 員園員鄖園孕葬時袁月各種 栽蚤亞暈穩(wěn)相在不同壓強下的結構和性能的研澡 栽蚤的暈焓將比 月員結構究結果為將來的高壓理論和實驗工作提供有益的的小袁在大約 員圓遠鄖園孕葬時袁月愿月員結員 栽蚤的暈焓將比參考 援構的小 援遠種結構的相對穩(wěn)定性排序由無外壓時的咱員暫 雜糟澡憎葬運則員扎怨愿悅苑則蚤砸賊藻援增

24、雜援燥造雜蚤賊鑿葬醞賊葬藻賊藻雜則糟援蚤員援猿圓員員咱員園暫鄖怎葬孕灶云袁宰葬灶早悅再袁再怎栽圓園園悅愿澡蚤灶孕援澡贈澤月援員苑猿園源園咱圓暫 砸蚤糟噪藻閱則雜遭袁贈月怎則灶孕藻允賊員賊怨愿雜苑怎則悅棗燥援葬賊栽援藻糟澡灶猿燥猿員造怨援員咱員員暫在澡怎允袁再怎允載袁宰葬灶再早 允袁悅澡藻載灶砸袁允蚤灶云早匝圓園園悅愿澡蚤灶援咱猿暫 災葬造增燥災鑿員葬怨怨雜遠怎則悅棗燥援葬栽賊藻援糟澡灶愿燥園遠造員援孕澡贈澤月援員苑圓圓員遠咱源暫 再葬灶孕載袁宰怎在鄖袁載怎允宰袁在澡葬灶再早允袁蘊蚤載袁在澡葬灶宰早宰咱員圓暫云藻灶勻早允袁蘊蚤云怎醞 圓園園悅怨澡蚤灶孕援澡贈澤月援員愿員緣苑源圓園園允源援雜贈灶賊

25、悅澡則援贈澤猿賊猿怨援苑淵源蚤灶悅澡蚤灶藻咱澤閆藻鵬冤勛 堯吳志國 堯咱員猿暫栽葬灶蘊暈袁勻怎悅耘袁再怎月砸袁悅澡藻載灶砸圓園園悅苑澡蚤灶孕援澡贈澤員援遠徐建偉 堯張玉娟 堯李搖鑫堯張偉偉 圓園園人源工晶體學報 猿猿怨苑源暫猿苑苑圓咱緣暫 孕蚤糟噪悅葬允則袁鑿暈藻藻鑿砸澤允圓園員暈園葬賊醞援葬賊藻怨則遠援圓源咱員源暫悅澡藻閱灶袁悅澡藻允灶閱袁在澡葬蘊燥勻袁宰葬灶悅早蘊袁再怎月勻袁雜澡閱蚤勻咱遠暫 雜蚤造孕葬袁澤再葬賊允藻砸澤袁勻葬贈灶孕藻閱澤圓園園孕愿澡贈澤砸援藻增月援苑愿員苑源員園員圓園園悅怨澡蚤灶孕援澡贈澤月援員愿園苑猿愿咱苑暫 再怎砸袁在澡葬匝灶袁閱藻允燥灶早蘊澡悅藻圓園園粵苑灶早藻悅憎

26、澡藻皂隕援灶賊耘援鑿蚤源賊遠咱員緣暫允蚤在勻袁在藻灶載早勻袁勻怎再允袁栽葬醞灶 匝圓園園粵愿糟賊孕葬澡贈雜澤蚤援灶緣援苑員員猿遠猿苑緣淵猿蚤悅灶澡蚤灶藻咱澤季藻正冤華 堯曾祥華 堯胡永金 堯譚明秋 圓園園物愿咱愿暫 粵造責賊藻雜噪袁閱蚤怎灶則葬灶鑿醞怎圓則園鑿園怎雜怨燥造雜蚤賊鑿葬悅賊燥藻皂皂怎員灶源援怨猿源緣理學報 緣苑猿苑緣猿暫咱怨暫 醞葬載鄖袁蘊蚤葬孕灶袁早醞蚤葬蘊燥袁月蚤雜藻宰袁在澡葬灶悅早運袁載怎蘊袁允蚤葬灶早 咱員遠暫允蚤鄖云袁在澡葬灶再早蘊袁悅怎勻蚤蘊袁蘊蚤載云袁在澡葬云燥袁醞藻灶悅早醞袁雜燥灶早允允圓園園孕怨澡贈澤雜援賊葬賊雜怎燥澤造月蚤圓鑿源蚤圓遠員猿圓在云圓園園粵怨糟賊孕葬

27、澡贈澤雜援蚤灶緣援愿源員園淵猿蚤悅灶澡蚤灶藻咱澤姬藻廣冤福 堯張艷園緣苑員搖緣園圓鄄物 理 學 報搖粵糟賊孕葬澡贈雜澤蚤援灶搖援災燥造遠援園袁暈燥緣援淵圓園員冤搖員園緣苑員園圓麗堯崔紅玲 堯李曉鳳 堯趙搖 峰堯孟川民 堯宋振飛圓園園物怨理學報緣愿源員園猿暫咱圓緣暫悅造葬雜則允噪袁雜藻早葬醞造閱造袁孕蚤糟噪悅葬允則袁鑿勻葬澤灶孕蚤允責袁孕則燥遭醞藻允則袁賊砸藻棗澤運燥袁灶孕葬贈醞灶藻悅圓園園在緣援運則蚤澤賊葬圓造圓造園緣燥遠早苑則援咱員苑暫蘊俟醞 再袁悅澡藻在灶宰袁蘊蚤蘊載袁蘊蚤砸怎孕袁宰葬灶早宰 運圓園園粵遠糟賊葬咱圓遠暫災葬灶鑿藻則閱遭員蚤怨造怨賊孕園澡贈澤砸援藻增月援源員苑愿怨圓孕澡贈澤雜

28、援蚤灶緣援緣猿緣苑淵遠蚤悅灶澡蚤灶藻澤呂藻夢冤雅咱堯陳洲文 堯李立新 堯劉咱圓苑暫雜糟澡觟藻灶暈遭員藻怨則緣早粵源糟賊悅葬澡藻皂雜援糟葬灶愿鑿圓援員猿日平 堯王文魁 圓園園物遠理學報 緣緣猿緣苑遠暫咱圓愿暫鄖怎燥在栽袁孕藻灶云早袁悅澡藻勻灶勻圓園員允園援雜燥怎賊澡哉憎灶藻蚤澤暈增賊葬援賊蚤燥灶葬造援咱員愿暫雜賊葬皂悅責袁棗醞造葬灶灶澤賊宰葬袁鑿粵賊澤葬澡砸蚤袁云則藻藻皂粵葬允灶圓園園孕員澡贈澤援淵暈葬賊雜援糟蚤耘援鑿蚤冤賊猿援遠員源緣郭咱振堂 堯彭搖放堯陳海花 圓園員西園南砸藻增月援遠猿員緣緣員園遠民族大學學報 淵自然科學版 冤猿遠員源緣暫咱員怨暫醞葬則醞造袁燥醞蚤造皂災葬圓灶園園孕園澡贈澤砸

29、援藻增月援遠圓圓愿怨怨咱圓怨暫悅澡藻勻灶勻袁孕藻灶云早袁醞葬燥勻運袁雜澡藻鄖灶再袁蘊蚤藻則皂勻葬孕灶袁灶蘊蚤在袁咱圓園暫蘊蚤蘊怎醞袁宰葬灶雜早匝袁再藻勻匝圓園園允緣援孕澡贈院澤悅援燥灶鑿藻醞灶葬澤賊援賊藻則援雜澡怎允云圓園員允園援粵責責孕造澡援贈澤員援園苑員員猿緣園猿員苑緣猿猿緣咱猿園暫允藻葬灶砸造袁燥粵扎澡則藻栽灶袁澤醞葬燥勻運袁月藻造孕造醞員怨苑雜怨糟蚤藻圓灶園糟愿遠藻圓怨咱圓員暫宰葬灶粵早允袁雜澡葬灶雜早蘊袁閱怎再袁運燥灶再早袁在澡葬灶蘊早允袁悅澡藻蘊灶袁咱猿員暫雜葬賊再燥袁允藻葬灶砸造員燥怨扎愿允員援鄖藻燥責澡砸贈藻澤澤愿援援遠月員員苑苑猿在澡葬閱燥閱袁蘊蚤在怎運圓園員悅園燥皂責怎醞賊

30、葬援賊藻雜則糟援蚤源援愿苑園緣咱猿圓暫蘊蚤勻怎袁醞葬燥勻袁醞葬鑿鑿怎醞則袁贈閱蚤灶再早袁醞藻灶再早袁勻覿怎澤藻則皂閱葬灶灶咱圓圓暫孕藻則鑿允藻孕憎袁月怎則噪運藻袁耘則灶扎藻醞則澡員燥怨棗怨孕遠澡贈澤砸援藻增蘊援藻賊苑賊苑援圓園園孕源澡贈澤砸援藻增月援苑園園怨源員員源猿愿遠緣咱猿猿暫粵造棗閱侉袁粵造棗則藻鑿醞澤袁澤月燥則灶燥鑿澡允燥袁鄖造蚤賊造造醞葬允灶袁栽燥憎造醞藻閱則袁咱圓猿暫醞燥灶噪澡燥勻則允澤袁孕賊葬糟允噪閱員怨苑孕遠澡贈澤砸援藻增月援員猿緣員愿愿暈藻藻鑿砸澤允圓園園孕緣澡贈澤砸援藻增月援苑圓園員源員員源咱圓源暫孕棗則燥皂月皂藻鄖則袁悅燥賊醞藻袁蘊燥怎蚤雜藻鄖袁悅燥澡藻醞灶蘊員怨怨允苑援

31、咱猿源暫醞怎則葬噪葬醞皂袁蚤勻蚤則燥運澤袁藻韻灶燥雜袁栽澤怎糟澡栽蚤袁贈隕葬澤澤澡醞蚤袁噪蚤悅燥皂責怎孕賊澡援贈員澤猿援員猿猿宰葬賊葬灶栽怎圓噪園蚤園孕源澡贈澤耘援葬則賊孕澡造援葬灶隕藻灶賊賊援藻員則源援圓遠苑猿粵遭蚤灶蚤糟賊葬蚤造燥糟怎造燥葬棗責賊則蚤藻燥澤灶澤怎則藻鄄責蚤澡灶葬鑿澤怎藻糟藻鑿賊則葬灶澤燥蚤棗栽賊蚤責暈燥燥灶造贈賊鄢贈責藻澤覮鄖怎載蚤燥灶鄖早葬搖燥雜澡葬灶早鄄孕藻灶早淵閱藻責葬則賊燥皂棗醞藻葬灶賊賊藻則雜蚤糟葬蚤造藻袁澤灶云糟怎藻鑿葬哉灶灶蚤增藻袁則雜澤澡蚤葬賊灶贈早搖澡圓葬園蚤園源猿悅猿澡袁蚤冤灶葬淵砸藻糟藻蚤猿增園粵藻怎鑿早怎圓澤園賊員園則曰藻增蚤皂澤葬藻灶鑿怎澤糟則則藻蚤糟責藻賊蚤猿增雜藻鑿責賊藻圓皂園遭員藻園則冤粵遭澤賊則葬糟賊月葬澤藻燥