搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

W-In体系溶质晶界偏聚行为的第一性原理计算

王奇 唐法威 侯超 吕皓 宋晓艳

W-In体系溶质晶界偏聚行为的第一性原理计算

王奇, 唐法威, 侯超, 吕皓, 宋晓艳
PDF
HTML
导出引用
导出核心图
  • 基于第一性原理构建了钨基合金体系的溶质偏聚模型, 以W-In体系为例研究了不同浓度下溶质的晶界偏聚行为和成键特征, 从电子结构层面揭示了W-In体系的键合作用, 预测了W-In体系界面稳定性随溶质浓度的变化规律. 结合键布居、电荷密度、差分电荷密度和态密度等电子结构分析, 发现了W-In体系中溶质原子在偏聚过程中的键性转变特征, 阐明了W-In键由晶粒内部的离子键过渡为晶界区域强共价键的微观机理. 模型计算首次得到了W-In体系中溶质本征偏聚能随In浓度的非单调变化规律, 结合键合作用和能量分析揭示了溶质浓度对本征偏聚能的影响机制. 计算预测了W-In体系达到高热稳定性所需的最佳溶质浓度范围和应避开的溶质浓度范围. 本研究为具有高温稳定性的钨基合金材料的设计与制备提供了理论基础和定量化指导.
      通信作者: 宋晓艳, xysong@bjut.edu.cn
    • 基金项目: 国家重点研发计划(批准号: 2018YFB0703902, 2016YFB0700503)和国家自然科学基金(批准号: 51631002, 51425101)资助的课题.
    [1]

    Zhou X Q, Li S K, Liu J X, Wang Y C, Wang X 2010 Mater. Sci. Eng. A 527 4881

    [2]

    Scapin M 2015 Int. J. Refract. Met. Hard Mater. 50 258

    [3]

    Nguyen Manh D, Muzyk M, Kurzydlowski K J, Baluc N L, Rieth M, Dudarev S L 2011 Key Eng. Mater. 465 15

    [4]

    Tschopp M A, Murdoch H A, Kecskes L J, Darling K A 2014 JOM 66 1000

    [5]

    Posthill J B, Hogwood M C, Edmonds D V 1986 Powder Metall. 29 45

    [6]

    Gul H, Uysal M, Çetinkaya T, Guler M O, Alp A, Akbulut H 2014 Int. J. Hydrogen Energ. 39 21414

    [7]

    Millett P C, Selvam R P, Saxena A 2007 Acta Mater. 55 2329

    [8]

    Hirouchi T, Takaki T, Tomita Y 2010 Int. J. Mech. Sci. 52 309

    [9]

    Song X, Zhang J, Li L, Yang K, Liu G 2006 Acta Mater. 54 5541

    [10]

    Liu F, Kirchheim R 2004 Scr. Mater. 51 521

    [11]

    Liu F, Kirchheim R 2004 J. Cryst. Growth 264 385

    [12]

    Liu F, Yang G, Kirchheim R 2004 J. Cryst. Growth 264 392

    [13]

    Liu F, Kirchheim R 2004 Thin Solid Films 466 108

    [14]

    Darling K A, Vanleeuwen B K, Koch C C, Scattergood R O 2010 Mater. Sci. Eng. A 527 3572

    [15]

    Chookajorn T, Murdoch H A, Schuh C A 2012 Science 337 951

    [16]

    Kawazoe Y 2001 Mater. Design 22 61

    [17]

    Bond A D, Solanko K A, Jacco V D S, Neumann M A 2011 CrystEngComm 13 1768

    [18]

    Braithwaite J S, Rez P 2005 Acta Mater. 53 2715

    [19]

    Yamaguchi M, Kaburaki H, Shiga M 2004 J. Phys.:Condens. Matter 16 3933

    [20]

    Reza M, Laws K J, Nikki S, Michael F 2018 Acta Mater. 158 257

    [21]

    Wu X, You Y W, Kong X S, Chen J L, Luo G N, Lu G H, Liu C S, Wang Z 2016 Acta Mater. 120 315

    [22]

    孟凡顺, 李久会, 赵星 2014 物理学报 23 237102

    Meng F, Li J H, Zhao X 2014 Acta Phys. Sin. 23 237102

    [23]

    Tang F, Liu X, Wang H, Hou C, Lu H, Nie Z, Song X 2019 Nanoscale 11 1813

    [24]

    Scheiber D, Pippan R, Puschnig P, Ruban A, Romaner L 2016 Int. J. Refract. Met. Hard Mater. 60 75

    [25]

    Segall M D, Lindan P J D, Probert M J 2002 J. Phys.:Condens. Matter 14 2717

    [26]

    Vanderbilt D 1990 Phys. Rev. B 41 7892

    [27]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [28]

    Ceperley D M, Alder B J 1980 Phys. Rev. Lett. 45 566

    [29]

    Pfrommer B G, Cote M, Louie S G, Cohen M L 1997 J. Comput. Phys. 131 233

    [30]

    Scheiber D, Razumovskiy V I, Puschnig P, Pippan R, Romaner L 2015 Acta Mater. 88 180

    [31]

    Zdanuk E J, Krock R H 1969 US Patent 3 423 203

    [32]

    Chelikowsky J R, Cohen M L 1976 Phys. Rev. B 14 556

    [33]

    Trelewicz J R, Schuh C A 2009 Phys. Rev. B 79 094112

    [34]

    Asta M, Wolverton C, Ozoliņš V 2004 Phys. Rev. B 69 144109

  • 图 1  $ {\varSigma 3} $ {111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $晶界模型和溶质偏聚位点示意图

    Fig. 1.  Diagram of $ {\varSigma 3} ${111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $ grain boundary model and solute segregation sites.

    图 2  不同溶质浓度下晶内和晶界区域的W—In键Mulliken布居分析 (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

    Fig. 2.  Mulliken population analysis of W—In bond at grain interior and grain boundary at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122.

    图 3  不同溶质浓度下W-In界面特性的电荷密度和差分电荷密度 (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

    Fig. 3.  Charge density and charge density difference of W-In interface characteristics at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

    图 4  不同溶质浓度晶内和晶界区域W和In原子的PDOS (a) xA = 0.0488; (b) xA = 0.122

    Fig. 4.  PDOS of W and In atoms at grain interior and grain boundary at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.122

    图 5  不同位点下偏聚能随溶质浓度的变化关系

    Fig. 5.  Segregation energy as a function of solute concentration corresponding to different sites.

    图 6  溶质浓度为0.0976时溶质偏聚前后体系晶内和晶界区域W和In原子PDOS

    Fig. 6.  PDOS of W and In atoms at grain interior and grain boundary before and after segregation at solute concentration of 0.0976.

    图 7  W-In体系溶质浓度选取范围示意图

    Fig. 7.  Diagram of solute concentration selection range in W-In system.

    表 1  不同溶质浓度下的W-In体系弹性常数计算结果

    Table 1.  Calculation results of elastic constants of W-In system at different solute concentrations. GPa

    溶质浓度 C11 C12 C44 C11C12 C11 + 2C12
    0 501.5 203.4 127.2 298.1 908.3
    0.0625 488.0 201.7 140.0 246.3 2128.5
    0.125 364.5 222.3 143.9 220.6 809.1
    0.25 208.0 245.9 138.6 –37.9 699.8
    下载: 导出CSV
  • [1]

    Zhou X Q, Li S K, Liu J X, Wang Y C, Wang X 2010 Mater. Sci. Eng. A 527 4881

    [2]

    Scapin M 2015 Int. J. Refract. Met. Hard Mater. 50 258

    [3]

    Nguyen Manh D, Muzyk M, Kurzydlowski K J, Baluc N L, Rieth M, Dudarev S L 2011 Key Eng. Mater. 465 15

    [4]

    Tschopp M A, Murdoch H A, Kecskes L J, Darling K A 2014 JOM 66 1000

    [5]

    Posthill J B, Hogwood M C, Edmonds D V 1986 Powder Metall. 29 45

    [6]

    Gul H, Uysal M, Çetinkaya T, Guler M O, Alp A, Akbulut H 2014 Int. J. Hydrogen Energ. 39 21414

    [7]

    Millett P C, Selvam R P, Saxena A 2007 Acta Mater. 55 2329

    [8]

    Hirouchi T, Takaki T, Tomita Y 2010 Int. J. Mech. Sci. 52 309

    [9]

    Song X, Zhang J, Li L, Yang K, Liu G 2006 Acta Mater. 54 5541

    [10]

    Liu F, Kirchheim R 2004 Scr. Mater. 51 521

    [11]

    Liu F, Kirchheim R 2004 J. Cryst. Growth 264 385

    [12]

    Liu F, Yang G, Kirchheim R 2004 J. Cryst. Growth 264 392

    [13]

    Liu F, Kirchheim R 2004 Thin Solid Films 466 108

    [14]

    Darling K A, Vanleeuwen B K, Koch C C, Scattergood R O 2010 Mater. Sci. Eng. A 527 3572

    [15]

    Chookajorn T, Murdoch H A, Schuh C A 2012 Science 337 951

    [16]

    Kawazoe Y 2001 Mater. Design 22 61

    [17]

    Bond A D, Solanko K A, Jacco V D S, Neumann M A 2011 CrystEngComm 13 1768

    [18]

    Braithwaite J S, Rez P 2005 Acta Mater. 53 2715

    [19]

    Yamaguchi M, Kaburaki H, Shiga M 2004 J. Phys.:Condens. Matter 16 3933

    [20]

    Reza M, Laws K J, Nikki S, Michael F 2018 Acta Mater. 158 257

    [21]

    Wu X, You Y W, Kong X S, Chen J L, Luo G N, Lu G H, Liu C S, Wang Z 2016 Acta Mater. 120 315

    [22]

    孟凡顺, 李久会, 赵星 2014 物理学报 23 237102

    Meng F, Li J H, Zhao X 2014 Acta Phys. Sin. 23 237102

    [23]

    Tang F, Liu X, Wang H, Hou C, Lu H, Nie Z, Song X 2019 Nanoscale 11 1813

    [24]

    Scheiber D, Pippan R, Puschnig P, Ruban A, Romaner L 2016 Int. J. Refract. Met. Hard Mater. 60 75

    [25]

    Segall M D, Lindan P J D, Probert M J 2002 J. Phys.:Condens. Matter 14 2717

    [26]

    Vanderbilt D 1990 Phys. Rev. B 41 7892

    [27]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [28]

    Ceperley D M, Alder B J 1980 Phys. Rev. Lett. 45 566

    [29]

    Pfrommer B G, Cote M, Louie S G, Cohen M L 1997 J. Comput. Phys. 131 233

    [30]

    Scheiber D, Razumovskiy V I, Puschnig P, Pippan R, Romaner L 2015 Acta Mater. 88 180

    [31]

    Zdanuk E J, Krock R H 1969 US Patent 3 423 203

    [32]

    Chelikowsky J R, Cohen M L 1976 Phys. Rev. B 14 556

    [33]

    Trelewicz J R, Schuh C A 2009 Phys. Rev. B 79 094112

    [34]

    Asta M, Wolverton C, Ozoliņš V 2004 Phys. Rev. B 69 144109

  • [1] 张杨, 宋晓艳, 徐文武, 张哲旭. SmCo7纳米晶合金晶粒组织热稳定性的热力学分析与计算机模拟. 物理学报, 2012, 61(1): 016102. doi: 10.7498/aps.61.016102
    [2] 白静, 王晓书, 俎启睿, 赵骧, 左良. Ni-X-In(X=Mn,Fe和Co)合金的缺陷稳定性和磁性能的第一性原理研究. 物理学报, 2016, 65(9): 096103. doi: 10.7498/aps.65.096103
    [3] 张章, 熊贤仲, 乙姣姣, 李金富. Al-Ni-RE非晶合金的晶化行为和热稳定性. 物理学报, 2013, 62(13): 136401. doi: 10.7498/aps.62.136401
    [4] 沈 祥, 聂秋华, 徐铁峰, 高 媛. Er3+/Yb3+共掺碲钨酸盐玻璃的光谱性质和热稳定性的研究. 物理学报, 2005, 54(5): 2379-2384. doi: 10.7498/aps.54.2379
    [5] 王艳, 曹仟慧, 胡翠娥, 曾召益. Ce-La-Th合金高压相变的第一性原理计算. 物理学报, 2019, 68(8): 086401. doi: 10.7498/aps.68.20182128
    [6] 汪志刚, 张杨, 文玉华, 朱梓忠. ZnO原子链的结构稳定性和电子性质的第一性原理研究. 物理学报, 2010, 59(3): 2051-2056. doi: 10.7498/aps.59.2051
    [7] 陈献, 程梅娟, 吴顺情, 朱梓忠. 石墨炔衍生物结构稳定性和电子结构的第一性原理研究. 物理学报, 2017, 66(10): 107102. doi: 10.7498/aps.66.107102
    [8] 宫长伟, 王轶农, 杨大智. NiTi形状记忆合金马氏体相变的第一性原理研究. 物理学报, 2006, 55(6): 2877-2881. doi: 10.7498/aps.55.2877
    [9] 罗明海, 黎明锴, 朱家昆, 黄忠兵, 杨辉, 何云斌. CdxZn1-xO合金热力学性质的第一性原理研究. 物理学报, 2016, 65(15): 157303. doi: 10.7498/aps.65.157303
    [10] 裴海林, 康建海, 姜恩永, 宋庆功, 郭 英. 插层化合物LixTiS2中Li离子-空位二维有序结构稳定性的第一性原理研究. 物理学报, 2007, 56(8): 4817-4822. doi: 10.7498/aps.56.4817
    [11] 刘琪, 管鹏飞. La65X35(X=Ni,Al)非晶合金原子结构的第一性原理研究. 物理学报, 2018, 67(17): 178101. doi: 10.7498/aps.67.20180992
    [12] 陈家华, 刘恩克, 李勇, 祁欣, 刘国栋, 罗鸿志, 王文洪, 吴光恒. Ga2基Heusler合金Ga2XCr(X = Mn, Fe, Co, Ni, Cu)的四方畸变、电子结构、磁性及声子谱的第一性原理计算. 物理学报, 2015, 64(7): 077104. doi: 10.7498/aps.64.077104
    [13] 胡前库, 侯一鸣, 吴庆华, 秦双红, 王李波, 周爱国. 过渡金属硼碳化物TM3B3C和TM4B3C2稳定性和性能的理论计算. 物理学报, 2019, 68(9): 096201. doi: 10.7498/aps.68.20190158
    [14] 张凯旺, 孟利军, 李 俊, 刘文亮, 唐 翌, 钟建新. 碳纳米管内金纳米线的结构与热稳定性. 物理学报, 2008, 57(7): 4347-4355. doi: 10.7498/aps.57.4347
    [15] 卢顺顺, 张晋敏, 郭笑天, 高廷红, 田泽安, 何帆, 贺晓金, 吴宏仙, 谢泉. 碳纳米管包裹的硅纳米线复合结构的热稳定性研究. 物理学报, 2016, 65(11): 116501. doi: 10.7498/aps.65.116501
    [16] 刘乐, 汤建, 王琴琴, 时东霞, 张广宇. 石墨烯封装单层二硫化钼的热稳定性研究. 物理学报, 2018, 67(22): 226501. doi: 10.7498/aps.67.20181255
    [17] 杨慎东, 宁兆元, 黄峰, 程珊华, 叶超. a-C:F薄膜的热稳定性与光学带隙的关联. 物理学报, 2002, 51(6): 1321-1325. doi: 10.7498/aps.51.1321
    [18] 熊小涛, 朱逢吾, 滕蛟, 蔡建旺, 赖武彦. NiFe/FeMn双层膜交换偏置的形成及热稳定性研究. 物理学报, 2004, 53(1): 272-275. doi: 10.7498/aps.53.272
    [19] 陈美娜, 张蕾, 高慧颖, 宣言, 任俊峰, 林子敬. Sm3+,Sr2+共掺杂对CeO2基电解质性能影响的密度泛函理论+U计算. 物理学报, 2018, 67(8): 088202. doi: 10.7498/aps.67.20172748
    [20] 王啸天, 代学芳, 贾红英, 王立英, 刘然, 李勇, 刘笑闯, 张小明, 王文洪, 吴光恒, 刘国栋. Heusler型X2RuPb (X=Lu, Y)合金的反带结构和拓扑绝缘性. 物理学报, 2014, 63(2): 023101. doi: 10.7498/aps.63.023101
  • 引用本文:
    Citation:
计量
  • 文章访问数:  256
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-10
  • 修回日期:  2019-02-19
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-04-01

W-In体系溶质晶界偏聚行为的第一性原理计算

  • 北京工业大学材料科学与工程学院, 新型功能材料教育部重点实验室, 北京 100124
  • 通信作者: 宋晓艳, xysong@bjut.edu.cn
    基金项目: 国家重点研发计划(批准号: 2018YFB0703902, 2016YFB0700503)和国家自然科学基金(批准号: 51631002, 51425101)资助的课题.

摘要: 基于第一性原理构建了钨基合金体系的溶质偏聚模型, 以W-In体系为例研究了不同浓度下溶质的晶界偏聚行为和成键特征, 从电子结构层面揭示了W-In体系的键合作用, 预测了W-In体系界面稳定性随溶质浓度的变化规律. 结合键布居、电荷密度、差分电荷密度和态密度等电子结构分析, 发现了W-In体系中溶质原子在偏聚过程中的键性转变特征, 阐明了W-In键由晶粒内部的离子键过渡为晶界区域强共价键的微观机理. 模型计算首次得到了W-In体系中溶质本征偏聚能随In浓度的非单调变化规律, 结合键合作用和能量分析揭示了溶质浓度对本征偏聚能的影响机制. 计算预测了W-In体系达到高热稳定性所需的最佳溶质浓度范围和应避开的溶质浓度范围. 本研究为具有高温稳定性的钨基合金材料的设计与制备提供了理论基础和定量化指导.

English Abstract

    • 钨基合金因其优异的力学性能被广泛应用于国防军工、航空航天和机械加工等领域, 有着不可替代的地位[1-3]. 研究表明: 如果将钨基合金制备成纳米晶材料, 与同成分的粗晶材料相比, 通常具有更加优异的力学性能[4]. 例如Hogwood等[5]观测到纳米晶钨基高比重合金的最高拉伸强度可达1700 MPa, 相当于传统钨基合金的两倍. 但是, 纳米晶材料因其晶界所占比例较高, 通常稳定性较差, 易发生热失稳导致晶粒快速长大, 致使纳米晶材料的优异性能大幅度削弱甚至丧失[6]. 因此, 纳米晶钨基合金的稳定性研究具有重要的研究价值[7-9].

      围绕纳米晶合金体系热稳定性较弱的关键问题, 研究人员进行了大量的实验研究和理论计算. Liu和Kirchheim[10-13]在多类合金体系中采用基于平衡态溶质偏聚理论的热力学计算和实验, 证实部分溶质元素通过晶界偏聚作用可有效增强纳米晶合金体系的热稳定性, 抑制中高温时纳米晶粒快速长大现象. Darling等[14]在Fe-Zr体系中的热力学计算和分子动力学模拟表明溶质元素的晶界偏聚效应可有效降低晶界形成能, 从而稳定纳米晶结构. Schuh等[15]基于偏聚热力学和蒙特卡罗模拟, 获得了不同温度下钨基纳米晶体系的稳定化元素优选图. 以W-Ti体系为例进行了实验验证, 结果表明在1100℃高温下退火1周仍能保持纳米晶平均晶粒尺寸约22 nm, 显示出W-Ti纳米晶体系优异的热稳定性. 综上, 选择合适的溶质元素, 使其产生强的偏聚作用, 能有效降低晶界形成能, 从而显著提高体系的热稳定性, 稳定纳米晶组织.

      第一性原理计算作为微观尺度计算中的重要研究方法, 在合金材料的研究中得到了广泛的应用[16]. 根据第一性原理的预测结果, 优选不同合金体系中适合的掺杂元素, 可有效减少实验的盲目性, 提高研发效率, 同时揭示元素作用的微观机理, 为新材料设计提供理论支撑[17]. 进一步地, 对于纳米合金偏聚体系, 第一性原理可以从体系电子结构和能量的角度对体系的偏聚行为进行量化分析, 且所需经验参数极少, 因此非常适合用于界面特征分析和强偏聚元素筛选[18,19]. 例如, Reza等[20]在铝合金中采用第一性原理, 研究了界面的电子结构和溶质原子与界面原子的键合方式, 通过电荷密度和轨道相互作用分析揭示了不同族溶质元素在界面处成键的键强和键性的变化规律. Wu等[21]通过第一性原理对多类钨晶界类型和多种钨基合金掺杂元素体系进行了偏聚能、断裂能和强化能的计算, 研究了多种溶质元素对界面结合强度的影响机理. 孟凡顺等[22]通过第一性原理界面拉伸模拟, 分析了Zn掺杂对铜基材料特征晶界的成键能力和力学性能的影响机制. 综上, 基于第一性原理可有效探索不同溶质元素对钨基合金晶界性质的影响作用. 进一步地, 研究表明强偏聚元素通常有助于增强钨基纳米晶合金的热稳定性, 抑制高温下发生晶粒快速长大, 钨基合金体系中的强偏聚掺杂元素包括In, Ru, Sc和Th等[15] .

      目前, 第一性原理较多的集中于溶质元素的优选, 缺乏溶质浓度对合金体系中溶质原子在界面扩散、偏聚、析出等行为的影响机制研究, 致使其预测结果与实际材料体系的实验结果仍存在较大偏差. 基于以上原因, 本文以目前鲜见报道的W-In体系强偏聚体系为例, 基于第一性原理晶界偏聚模型, 引入溶质浓度效应, 系统研究了W-In体系在不同溶质浓度和界面偏聚位点下的溶质偏聚行为和键合特征, 从电子结构层面揭示了W-In体系的键合方式和界面稳定化机制, 预测高热稳定性W-In体系的最佳溶质浓度添加范围.

    • 由于具有较高的晶界占据比例和较低的晶界能量, $ {\varSigma 3}${111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $晶界模型在钨基合金的研究中具有较高的代表性[21,23]. 其中, Wu等[21]利用$ {\varSigma 3} $ {111}$ \left\langle {\; 1\; \bar1 \; 0\; } \right\rangle$等界面研究了过渡金属偏聚对钨晶界的强化作用; Scheiber等[24]$ {\varSigma 3} $ {111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle$等界面对过渡元素掺杂对晶界内聚力影响进行研究. $ {\varSigma 3} $ {111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $晶界沿公共晶面构成镜面对称的位向关系, 属于低能晶界. 本文选取$ {\varSigma 3} $ {111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $晶界模型作为研究对象, 构建了W-In体系的$ {\varSigma 3} $ {111}$ \left\langle {\; 1\; \bar1 \; 0\; } \right\rangle$晶界模型. 晶界模型层数对偏聚模型的计算精度具有一定的影响, 本文研究重点是偏聚行为随溶质浓度和偏聚位点的变化规律, 综合考虑计算精度和计算成本, 模型以晶界中心平面为界, 左右各20层原子, 共计41个原子. 图1的第一性原理晶界模型包括晶界区域和晶内区域两部分, 其中晶界区域为: 以晶界中心平面为界上下各包含7层原子. 溶质原子偏聚到晶界区域时, 考虑了3类不同的偏聚位点.

      图  1  $ {\varSigma 3} $ {111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $晶界模型和溶质偏聚位点示意图

      Figure 1.  Diagram of $ {\varSigma 3} ${111}$\left\langle {\; 1\; \bar1 \; 0\; } \right\rangle $ grain boundary model and solute segregation sites.

    • 本文所有的计算工作, 均采用Materials Studio软件中的CASTEP[25]软件包完成. 本文采用了超软赝势(ultrasoft pseudopotential, Usp)[26]进行模型计算, W和In的价电子结构分别为5s25p65d46s2和4d105s25p1, 泛函形式采取GGA-PBE (Perdew-Burke-Ernzerhof generalized gradient approximation)方法[27,28]. 对晶体结构和界面结构的优化计算采用BFGS方法[29]搜索最优结构, 平面波截断能为400 eV, k点设置为为5 × 5 × 1, 各项收敛参数为能量小于2 × 10–6 eV/atom, 原子间的相互作用力小于0.05 eV/Å, 最大应力偏差小于0.1 GPa, 原子的最大位移变化小于0.002 Å.

      对溶质晶界偏聚行为的定量描述, 可以直接计算从晶界较远处偏聚到晶界区域的能量变化, 从而获得偏聚能. W-In体系中偏聚能计算公式[30]可表示为

      $\Delta {E^{{\rm{seg}}}}\left( {{x_{\rm{A}}}} \right) = E_{{\rm{GB}}}^{{\rm{seg}}}\left( {{x_{\rm{A}}}} \right) - E_{{\rm{GB}}}^{{\rm{bulk}}}\left( {{x_{\rm{A}}}} \right), $

      式中$E_{{\rm{GB}}}^{{\rm{seg}}}\left( {{x_{\rm{A}}}} \right) $为偏聚后溶质浓度为xA时体系总能量; $E_{{\rm{GB}}}^{{\rm{bulk}}}\left( {{x_{\rm{A}}}} \right) $为偏聚前溶质浓度为xA时体系总能量. 二者差值即为溶质原子偏聚过程引起的能量变化. 偏聚能数值越小代表体系的偏聚能力越强. 一方面, 模型构建中采用不同数目的溶质原子替换溶剂原子以此引入溶质浓度效应, 并在某一浓度下构建多种溶质原子分布模式进行能量计算, 从中选取能量最低的分布模式进行该浓度点下的偏聚行为计算. 另一方面, W和In为不互溶体系[31], 关于W和In的化合物鲜见报道, 因此排除二者在偏聚过程中存在其他相或化合物生成的可能性.

      此外, 对于不同溶质浓度下W-In体系的弹性常数计算, 采用W晶胞的2 × 2 × 2超胞进行模型构建, 总计包含16个原子, 计算中最大应变振幅设置为0.003, 每步应变步数为6.

    • 布居分析可得到电子在晶体内的分布情况, 量化键合作用的类型和强度, 同时能够为体系基态性质提供更多细节信息. 通常认为键布居数越接近0, 键的离子性越强; 反之越偏离0, 共价性越强. W-In体系晶粒内部和晶界区域的部分Mulliken键布居计算结果如图2所示.

      图  2  不同溶质浓度下晶内和晶界区域的W—In键Mulliken布居分析 (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

      Figure 2.  Mulliken population analysis of W—In bond at grain interior and grain boundary at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122.

      通过图2对比晶内区域和晶界区域的W— In键布居值可以发现, 在不同溶质浓度下, In原子由晶粒内部偏聚到晶界区域的过程中, 均发生了明显的键性变化: W—In键由离子键转变为共价键. 为了进一步探究W-In体系的键性转变现象, 接下来通过电荷密度图、差分电荷密度图和态密度图对键性转变现象进行系统的分析.

      图3展示了不同溶质浓度下W-In体系的电荷密度和差分电荷密度. 如图3左图所示, 不同溶质浓度下界面处W—In键之间的电荷密度均小于晶内区域W—In键之间的电荷密度, 体现出与键布居结果相悖的电荷密度分布特征. 上述分析表明, 电荷密度没有反映出晶界W—In键的共价键特征. 其主要原因是: 在晶内区域, 原子排布有序, In原子掺杂之后, 与周围W原子的平均成键距离略短于界面区域W—In键的平均成键距离, 不同区域W—In键键长差异, 使得界面区域W—In键之间的电荷密度分布相比晶内区域W—In键的电荷密度分布可能略低, 二者对比分析导致界面W—In键的共价特征不明显. 为了排除电荷密度分析中晶内和晶界区域W—In键键长对键性分析的干扰, 进一步采用差分电荷密度对键性进行了分析, 以避免单纯通过电荷密度分析体系键性所导致的结果误判.

      图  3  不同溶质浓度下W-In界面特性的电荷密度和差分电荷密度 (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

      Figure 3.  Charge density and charge density difference of W-In interface characteristics at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.0732; (c) xA = 0.122

      不同溶质浓度下的差分电荷密度图如图3右侧图所示. 差分电荷密度相比电荷密度, 反映的不是电荷聚集的程度, 而是电荷转移情况, 因此即便键长变化导致电荷密度变化, 但共价键(电子共用)和离子键(电子转移)的本质差异仍可以体现. 对应于电荷密度, 对不同的界面模型进行了差分电荷密度进行了计算, 结果如图3右侧图所示, 其中蓝色代表缺失电子, 红色代表富集电子, 白色代表不变.

      通过观察晶内区域和晶界区域的差分电荷密度图可以发现, 在不同溶质浓度下, 晶内区域W—In键对应的W和In原子之间为白色, 表明二者之间不存在电子共用情况, 即不体现共价键的特性; 而在晶界区域, 构成W—In键的W原子和In原子之间呈现红色, 表明二者之间之间存在共用电子现象, 体现共价键特性. 通过上述分析, 采用差分电荷密度图可以明显区分晶内区域和晶界区域的W—In键的键性差异.

      为了进一步研究键性变化, 对成键原子的态密度进行分析, 结果如图4所示. 我们可以发现晶界区域W原子和In原子偏态密度(PDOS)的共振强度明显高于晶粒内部, 表明在晶界区域形成共价键. 共价键的形成使得更容易产生强的偏聚效果, 使晶界能大大降低, 提高稳定性.

      图  4  不同溶质浓度晶内和晶界区域W和In原子的PDOS (a) xA = 0.0488; (b) xA = 0.122

      Figure 4.  PDOS of W and In atoms at grain interior and grain boundary at different solute concentrations: (a) xA = 0.0488; (b) xA = 0.122

      溶质元素偏聚前后和基体的键合方式发生突变, 这在合金界面分析中较为少见, 结合W-In体系的电子结构特点, 现对体系偏聚过程中的键性转变现象作如下分析: W的价电子结构为 5s25p65d46s2, In的价电子结构为 4d105s25p1. 晶内区域W和In原子排布规则, 其中W原子提供4—6个价电子, In原子则提供3个价电子, 此时对应的W—In键中, 原子之间电荷分布的中心向W+4~+6离子实偏移, 具有离子键的成分[32]. 而在晶界区域, 由于结构较为松散, 空间较大, W周围存在大量悬链键, 未成键电子较多的聚集在W原子附近, 导致W原子的价态降低, 氧化性减弱, In原子周围的电子不再发生明显的转移, 而更多地出现在W和In之间, 导致键性从离子性转变为共价性.

    • 基于不同溶质浓度下的偏聚模型, 得到了多类界面偏聚位点下偏聚能随溶质浓度的变化关系, 如图5所示.

      图  5  不同位点下偏聚能随溶质浓度的变化关系

      Figure 5.  Segregation energy as a function of solute concentration corresponding to different sites.

      图5可以发现, 随着溶质浓度的增加, 偏聚能呈非单调变化的规律, 表明溶质偏聚能力随溶质浓度的增加先减弱再增强, 存在极值规律. 在极值点之前, 随溶质浓度的增加, 偏聚能趋于零, 表明其偏聚能力快速衰减; 但是随着溶质浓度的进一步增加, 超过极值点对应的浓度时, 体系偏聚能力显著提升. 根据经典理论, 通常认为偏聚能低于–0.5 eV的体系为强偏聚体系[33], 由图5可知, 当溶质浓度大于0.0976时, W-In体系再次成为典型的强偏聚体系.

      为了分析偏聚能随溶质浓度变化的极值规律, 进一步对偏聚前后W-In体系的电子结构进行了分析. 如图6所示, 极值点所对应的溶质浓度为0.0976时, 偏聚前后体系中W和In原子的PDOS存在类似的共振关系, 说明此时In原子在晶内和晶界处与W原子的成键差异较小, In原子并不具有较强的晶界偏聚趋势. 进一步地, PDOS同时反映了W-In体系偏聚前后体系的形成能. 较小的PDOS差异说明W-In体系在该浓度点下偏聚前后的形成能差异较小. 形成能[34]可描述为

      图  6  溶质浓度为0.0976时溶质偏聚前后体系晶内和晶界区域W和In原子PDOS

      Figure 6.  PDOS of W and In atoms at grain interior and grain boundary before and after segregation at solute concentration of 0.0976.

      ${E_{\rm{f}}} = E\left( {{{\rm{W}}_p}{\rm{I}}{{\rm{n}}_{\rm{q}}}} \right) - pE\left( {\rm{W}} \right) - qE\left( {{\rm{In}}} \right), $

      式中pq分别为W和In在超晶胞中的原子数; pE(W)和qE(In)为其单质中的单原子能量; $E\left( {{{\rm{W}}_p}{\rm{I}}{{\rm{n}}_q}} \right)$${{\rm{W}}_p}{\rm{I}}{{\rm{n}}_q}$超晶胞的总能量. 结合(1)式和(2)式可知, 偏聚能与形成能之间存在如下关联:

      $\Delta {E^{{\rm{seg}}}}\left( {{x_{\rm{A}}}} \right) = {E_{{\rm{f}}2}} - {E_{{\rm{f}}1}}, $

      其中Ef1Ef2分别为偏聚前后体系形成能. 可见偏聚前后的体系的形成能差异决定了体系偏聚能的大小. 因此, 一方面, 极值点下PDOS的较小差异导致形成能差异较小, 最终使偏聚能趋近0, 偏聚能力下降. 另一方面, 由图4可以看出, 对比极值点较远处的浓度点, 浓度0.0488和0.122的PDOS偏聚前后共振峰差异显著, 偏聚后晶界处形成的强共价键使偏聚后的W-In体系形成能明显降低, 因此, 此时体系偏聚前后形成能差异较大, 最终使体系的偏聚能较负, 更容易产生偏聚行为.

      不同溶质浓度下偏聚前后体系在晶内和晶界区域的键强和电子结构差异, 导致了不同溶质浓度下体系整体稳定性的差异, 并反映在偏聚前后的形成能差异中, 最终导致偏聚能随溶质浓度变化的非单调变化规律.

    • 为了确定W-In体系达到强偏聚所需的最佳溶质浓度范围, 进一步对其力学稳定性进行了分析. W-In属于立方晶系, 具有3个独立的弹性矩阵元(C11, C12, C44), 通过基于第一性原理的弹性常数计算, 得出结果见表1. 根据体心立方力学稳定性判据:

      溶质浓度 C11 C12 C44 C11C12 C11 + 2C12
      0 501.5 203.4 127.2 298.1 908.3
      0.0625 488.0 201.7 140.0 246.3 2128.5
      0.125 364.5 222.3 143.9 220.6 809.1
      0.25 208.0 245.9 138.6 –37.9 699.8

      表 1  不同溶质浓度下的W-In体系弹性常数计算结果

      Table 1.  Calculation results of elastic constants of W-In system at different solute concentrations. GPa

      $ {C_{11}} - {C_{12}} > {\rm{ }}{0}, {C_{11}} > {\rm{ }}{0}, {C_{44}} > {\rm{ }}{0}, {C_{11}} + {\rm{ }}{2}{C_{12}} > {\rm{ }}{0}. $

      表1可以发现, 当溶质浓度为0.25时, C11-C12为负值, 表明W-In体系出现力学失稳. 因此, 保守估计W-In体系的力学稳定浓度范围可确定在溶质浓度0.125之内. 当偏聚能低于–0.5 eV后, W-In体系成为典型的强偏聚体系. 结合体系稳定性和偏聚能两方面的分析, 最终得到了W-In体系的溶质浓度添加范围预测图, 如图7所示. 当溶质浓度小于0.0632和溶质浓度大于0.106时, W-In体系属于强偏聚体系; 当溶质浓度为0.0632—0.106时, 体系为弱偏聚体系, 甚至在溶质浓度约为0.09时, 偏聚能接近0, 此时通过偏聚效应无法起到降低晶界能的效果, 说明在溶质添加过程中应尽量避免该浓度范围. 具体来看, 当溶质浓度小于0.0632时, 体系的溶质偏聚能力虽然较强, 但同时受限于较低的溶质浓度, 对于晶界能量的降低会起到一定的作用, 但可能无法发挥到最佳水平, 因此该浓度为可选范围, 但并非最适浓度范围. 进一步地, 当溶质浓度为0.106—0.125时, W-In体系的溶质偏聚作用较强, 且具有较高的溶质浓度, 体系通过晶界偏聚作用对晶界能的降低起到显著作用, 因为该浓度范围是W-In体系最适添加浓度范围. 当溶质浓度大于0.125时, W-In体系可能出现力学失稳, 导致相变和相分解等, 因此也不利于获得稳定的纳米晶体系.

      图  7  W-In体系溶质浓度选取范围示意图

      Figure 7.  Diagram of solute concentration selection range in W-In system.

      综上, 结合W-In体系的电子结构和能量计算, 系统地揭示了W-In体系晶界偏聚过程中的溶质元素键性转变机制, 定量描述了体系的溶质偏聚行为, 预测了W-In体系热稳定性随溶质浓度的变化规律. 以上计算为钨基体系溶质元素优选提供了理论支撑, 同时为钨基纳米晶合金材料设计研发提供了新的思路.

    • 本文基于第一性原理溶质偏聚模型, 研究了W-In体系不同溶质浓度下溶质晶界偏聚行为和成键特征, 得到如下结论:

      1) 基于第一性原理, 首次建立了考虑溶质浓度效应的W-In体系晶界偏聚模型, 预测了该体系溶质原子偏聚过程中W—In键从离子性向共价性的转变规律. 结合几何结构和差分电荷密度分析, 揭示了键性变化的内在机理: 晶界结构和晶内结构的差异, 导致晶界区域内W原子价态降低, 氧化性减弱, 最终导致W—In键性转变.

      2) 通过不同溶质浓度偏聚前后体系稳定性和电子结构, 揭示了偏聚能随溶质浓度非单调变化规律的内在机理, 进一步结合W—In键的PDOS分析, 获得了W-In体系偏聚能随溶质浓度变化规律的电子结构机理: 溶质浓度较大地影响了W-In键偏聚前后的键强度, 导致溶质浓度趋于0.0976时偏聚能力明显下降, 最终形成偏聚能随溶质浓度的极值规律.

      3) 结合物相力学稳定性和溶质晶界偏聚分析, 在不考虑空位浓度的情况下, 模型预测其纳米晶体系中具有高热稳定性的最佳溶质浓度添加范围: 0.106—0.125; 同时也得到了需要规避的溶质浓度范围: 0.0632—0.106. 该模型对于W-In纳米晶合金材料的设计研发尤其是浓度调控起到了理论指导作用.

参考文献 (34)

目录

    /

    返回文章
    返回