【MS論文精讀】DFT+AIMD計算界面能和GSFE曲線,研究界面熱力學和動力學穩定性 2023年11月3日 上午9:27 ? 計算 ? 閱讀 260 研究背景 界面穩定性對用于超高溫應用的雙相Ir/Ir3X材料至關重要。盡管如此,人們對這一性質的認識和理解比較匱乏。近日,西安交通大學李燁飛、王怡然等人通過DFT和AIMD方法計算了界面能和廣義層錯能(GSFE)曲線,以研究六個Ir/Ir3X界面的熱力學和動力學穩定性。 計算方法和模型結構 作者將表面能和界面能分別在0K下用第一性原理和2000K下用AIMD方法計算,考慮到計算成本,作者僅在0K時計算GSFE。在Perdew–Burke–Ernzerhof(PBE)廣義梯度近似(GGA)的約束下,作者使用投影增強波(PAW)方法和利用CASTEP進行總能量計算,而離子和電子之間的相互作用用超軟贗勢來描述。 作者使用截止能量為600 eV的電子波函數的平面波基組,以及利用Monkhorst–Pack方案對第一布里淵區進行數值積分,相應的網格為0.03??1。同時,自洽計算的標準被設置為1meV/atom。為了獲得有限溫度下的表面能和界面能,作者在2000K下對體相、表面和界面結構進行了AIMD模擬,并選擇了400 eV的平面波截止能量,以及使晶格以1fs的時間步長弛豫2ps,使其在2000K下的NPT系綜下達到平衡。 圖1. 模型結構 作者首先對大塊Ir和六種Ir3X合金進行了幾何優化,然后將平衡態結構擴展到超胞。作者選擇低指數緊密堆積表面(111)來研究表面和界面性質,因為它是FCC(Ir)和L12(Ir3X)材料中的優先解理面。該模型在X軸和Y軸方向上采用周期性邊界條件,并且每層中具有八個原子,其中對于Ir表面為八個Ir原子,而對于Ir3X表面為六個Ir原子和2X原子。作者在Z軸方向上插入15?的真空層,以避免兩個自由面之間的相互作用。 此外,作者通過測試確定了平板模型的層數,以確保計算的準確性和效率。對于界面,作者在FCCIr(111)/L12-Ir3X(111)界面中考慮了FCC中空配位界面結構,并且將堆疊序列設置為ABCABCA(Ir側)|BCABCAB(Ir3X側)。 作者將晶格常數設為兩相的平均值,以避免由兩相之間不同的晶格常數引起的應變。圖1(a)和(b)是L12-Ir3X和FCC Ir結構,而圖1c是用于計算的界面模型,其中虛線表示界面。 結果與討論 圖2. 表面能 如圖2所示,Ir3Zr(111)和Ir3Hf(111)的表面能低于其他Ir3X表面的表面能,表明它們在0K下是最穩定的表面,其次是Ir3Ti(111)表面。此外,Ir3Ta(11 1)在0K和2000K下都具有最高的表面能,這意味著它是最活躍的界面結構。 總體而言,所有Ir3X(111)表面的表面能在溫度升高到2000 K后都會降低,這是因為振動的原子越多,在高溫下的內聚力結合行為就越弱。 圖3. 界面能 從圖3中可以看出,在0 K和2000 K時,Ir/Ir3Ti、Ir/Ir3HF和Ir/Ir3Zr界面比Ir/Ir3V、Ir/Ir3Nb和Ir/Ir3Ta界面更穩定,表明添加IVB合金元素(Ti、HF和Zr)會導致Ir中的界面結構比三種VB元素(V、Nb和Ta)中的界面結構更穩定。 圖4. 差分電荷密度 圖5. 比例關系 為了說明含有IVB和VB元素界面之間的差異,圖5中Ti、Hf和Zr元素的界面用藍色標記,而V、Nb和Ta元素的界面則用紅色標記。如圖5(a)所示,更活躍的Ir3V(111)、Ir3Nb(111)和Ir3Ta(111)表面并沒有形成更穩定的Ir/Ir3X界面。而圖5(b)中的Ir/Ir3Ti(Ir/Ir3Hf)界面表現出比Ir/Ir3V(Ir/Ir3Ta)界面更小的晶格失配,并表現出更好的界面穩定性。 這表明,由于相間邊界處的晶格畸變導致了界面的不穩定性。然而,并非所有含IVB元素的界面都具有較小的晶格失配,其中Ir/Ir3Zr和Ir/Ir3Nb就是例外。此外,作者還研究了界面電子相互作用因子,該因子表示為當兩個表面形成界面時發生的電荷再分配。六個Ir/Ir3X界面的電荷密度差(CDD)如圖4所示,粉紅色區域表示電荷損失并對應于PTCDD曲線上的正值,而藍色區域表示電荷損耗并對應于曲線上的負值。電荷再分配的增加伴隨著曲線中絕對值的增加。 圖4中的CDD和PTCDD曲線表明電荷再分配集中在兩個原子層內的界面處。通過比較界面中CDD的等效表面,可以發現含有Ti、V或Zr元素的界面比含有Nb、Hf或Ta元素的界面表現出更多的電荷再分配。 六種Ir/Ir3X結構的界面能作為電荷再分配的函數如圖5所示,兩者之間存在明顯的負相關性。在Ir和Ir3Ta之間的電荷再分配最小的情況下,形成了最不穩定的界面結構。而對于相同周期的元素,與含有V/Nb/Ta元素的界面相比,含有Ti/Hf/Zr元素的界面明顯具有更多的電荷再分配,并具有更好的界面穩定性。總體而言,較小的晶格畸變和更多的電荷再分配可以形成更穩定的界面結構,如Ir–Ir3(Ti/Hf/Zr)。 圖6. 層錯能量曲線 在圖6中,對于每個完美的晶體,當兩個界面原子平面彼此滑過時,滑移平面的單位面積能量增加,然后發現峰值位于μ/b≈0.6處,該峰值對應于不穩定的堆垛層錯能(γus)。隨后,能量降低到與本征堆垛層錯能量(γis)相對應的最小值,并且結構達到相對穩定的構型。如圖6所示。 對于最有可能的(111)滑移系統,不同界面系統的γus值的降序是Ir/Ir3Ti>Ir/Ir>Ir/Ir3Hf>Ir/Ir3V>Ir/If3Zr>Ir/Ir3Ta>Ir/Il3Nb,及大多數第二相的引入降低了銥基質的GSFE。盡管較低的GSFE值降低了界面的動力學穩定性,但同時也降低了Ir材料中由交叉滑移引起的脆性斷裂的概率。 當比較包含相同周期元素的界面時,可以發現界面GSFE值如下:Ir/Ir3Ti>Ir/Ir3V,Ir/Ir3Zr>Ir/Ir3Nb,以及Ir/Ir3Hf>Ir/Il3Ta。這表明VB元素(V、Nb和Ta)的加入提高了Ir基材料的塑性,而IVB元素(Ti、Zr和Hf)的加入改善了界面穩定性。 結論與展望 研究發現,Ir/Ir3Ti在0K和2000K時具有最低的界面能,以及最高的不穩定層錯能,表明它是最穩定的界面。此外,含有IVB元素(Ti、Hf和Zr)的界面比含有VB元素(V、Nb和Ta)的界面更穩定。 通過分析表明,Ir/Ir3X的界面能差異有助于界面晶格失配和界面處的電荷再分配,該工作為新型高溫雙相銥基合金的設計提供了理論指導。 文獻信息 Bai Xue et.al A comparative study on the stability of six Ir/Ir3X (X = Ti, V, Zr, Nb, Hf, Ta) interfaces by first-principle and AIMD calculations Applied Surface Science 2023 https://doi.org/10.1016/j.apsusc.2023.157502 點擊閱讀原文,報名計算培訓! 原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2023/11/03/68fd32e0aa/ 贊 (0) 0 生成海報 相關推薦 【DFT+實驗】浙理工紀律律/王晟、紹興文理劉準:RuX?(X = S/Se/Te)中硫族配體效應及析氫構效關系 2023年11月8日 【DFT+實驗】Chem. Mater. | 宇宙的浪漫:金基雙鈣鈦礦衍生材料 2023年11月20日 中科大江海龍/周蒙,最新Angew! 2024年3月3日 云南大學何天威副教授/鹽城師范學院孔佑超博士Small研究論文:單原子種類和配位微環境共同調控實現穩定高效多功能催化劑 2023年11月16日 【MS論文解讀】JMCA:硼烯/MoS2異質結在不同堆積條件下析氫反應的催化性能 2023年12月19日 發現38萬種新材料、17天自主合成41種新化合物,DeepMind一日兩篇論文登上Nature 2024年2月8日