First-principles calculation and experimental study on the influence mechanism of diffusion activation energy of Cu atoms in current-assisted sintering
-
摘要:
从第一性原理计算与电流辅助烧结两方面出发,研究了外加电场对晶体Cu扩散激活能的影响规律。结果表明,外加电场和电流使Cu空位产生的难度降低,但是原子迁移能几乎不变,导致扩散激活能在达到电场强度(电流密度)阈值(2 V·Å‒1(307.1 A·cm‒2))后略有下降,超过阈值后会剧烈下降,最终在电场强度(电流密度)达到5 V·Å‒1(708.5 A·cm‒2)后,由于空位形成能逐渐下降到0,扩散激活能下降至临界值,扩散激活能临界值相比于阈值对应的扩散激活能降低了约60.2%。研究结果揭示扩散激活能在电场或电流作用下呈现出明显规律性的下降趋势,实验结果与第一性原理计算模拟结果呈现较好的正相关性。
Abstract:The effect of the applied electric field on the diffusion activation energy of Cu crystal was studied by the first-principles calculation and current assisted sintering. The results show that the applied electric field and current reduce the difficulty of Cu vacancy generation, but the atomic migration energy is almost unchanged, resulting in a slight decrease in diffusion activation energy after reaching the electric field intensity (current density) threshold (2 V·Å‒1 (307.1 A·cm‒2)), and a sharp decrease after exceeding the threshold; finally, after the electric field intensity (current density) reaching 5 V·Å‒1 (708.5 A·cm‒2), the vacancy formation energy gradually decreases to 0, the diffusion activation energy drops to the critical value, and the critical value of diffusion activation energy decreases by about 60.2% compared with that of the threshold value. The diffusion activation energy shows a regular decline trend under the action of electric field or current, and the experimental results show a good positive correlation with the first-principle simulation results.
-
随着现代电子科技的不断发展,电子产品已经趋向于小型化、便携化、多功能化方向发展,封装材料对于微电子器件不仅起着连接、支撑和散热的作用,还需要有高热导率,低介电损耗等性质,因此对电子封装材料的电气性能、散热性能以及稳定性能提出了更高的要求[1–3]。传统的陶瓷、金属等单一材料已经不能满足当下电子器件对于封装材料的综合要求,复合材料作为金属基封装材料的首选逐渐成为研究的热点,Mo–Cu复合材料以其优异的导热导电性能、可调的热膨胀性能和耐高温性能被应用于电子封装材料、热沉材料、电触头、电极材料和电子束靶等领域,在国防科技、航空航天领域有广泛的应用。Mo–Cu复合材料是由Mo和Cu两相组成,Mo和Cu两相不仅物理性能和晶体结构差距较大,而且Mo和Cu两相互不相溶,两者的熔点相差较大,无法使用传统的熔炼方法制备Mo–Cu复合材料。随着Mo–Cu复合材料的使用越来越广泛,对其性能的要求也逐渐提高,国内外学者对于Mo–Cu复合材料的制备和性能优化给予了持续的关注[4,5]。
目前Mo–Cu复合材料制备方法有溶渗法[6]、液相烧结法、低温活化烧结[7,8]、热压烧结法[9]、放电等离子烧结法[10]、微波烧结法[11]和激光快速成形法等[12,13]。液相烧结法作为粉末冶金常用成形方法对设备要求低,但复合材料的相对密度仍需提高。热压烧结法的加热和加压同时进行,容易获得接近理论密度的样品,但对模具和设备的要求较高,只能生产较为简单的产品,而且加工成本高、生产效率低、耗费大[14]。微波烧结法虽然烧结速度快,但依然存在产品均匀性不容易控制、加热设备复杂、微波透过材料不能烧结、出现热斑现象[15]等问题。相较于上述制备方法,熔渗法制备Mo–Cu复合材料不仅对设备的要求较低、工序简单、适合批量生产,而且制备的复合材料相对密度较高、性能好、导热和导电性能也较为理想。因此目前商用的Mo–Cu复合材料往往还是基于熔渗法制备。传统熔渗法采用纯Mo粉压制坯块,坯块骨架强度低,需要通过预烧结才能进行熔渗,预烧结会使纯Mo骨架孔隙率降低影响熔渗效率,而且Mo与Cu润湿性较差,Cu液不能充分渗入纯Mo坯块内,难以制备出高铜含量、高相对密度且综合性能优良的Mo–Cu复合材料。王新刚等[16]尝试了在W–Cu复合材料熔渗时添加诱导Cu,发现虽然由于W和Cu的密度相差过大导致在混粉过程中Cu粉分布不均匀,但基于诱导法制备的复合材料其电导率和相对密度都要高于预烧结骨架法获得的复合材料,原因是在压坯过程中添加诱导Cu,一方面增加压坯的密度和强度,减小孔隙率,从而提高熔渗时的毛细作用力,另一方面也改善了Cu液的润湿性,因此使用诱导熔渗法Cu液可以充分的填充孔隙,提高复合材料的相对密度,在工业生产中更易进行,成品率更高。韩蕊蕊等[17]进一步研究了诱导Cu粉末特性对Mo–Cu复合材料熔渗质量的影响,发现诱导Cu能够使熔渗更充分,提升复合材料成分的均匀性。但很少有文献报道诱导Cu含量对复合材料微观组织和性能的影响。
本文进一步优化制备Mo–Cu复合材料的熔渗工艺,通过诱导熔渗法改善传统熔渗技术方法,研究诱导Cu含量对于Mo–Cu复合材料微观组织及热性能的影响规律,探寻低成本制备综合性能更佳Mo–Cu复合材料的新方法,为助推Mo–Cu复合材料在高功率、高密度电子封装领域的应用提供思路和技术借鉴。
1. 实验材料及方法
1.1 实验材料
实验选用金堆城钼业股份有限公司提供的纯Mo粉,纯度为≥99.95%,平均粒径6.5 μm;北京有色金属研究院提供的电解Cu粉,纯度≥99.99%,平均粒径8.0 μm;两种粉末的微观形貌如图1所示。实验所用Cu块由洛铜集团提供,纯度为99.9%,Cu块尺寸为ϕ50 mm×1 mm。
1.2 压坯和熔渗过程
将诱导Cu粉和Mo粉,分别按照Mo–0%Cu(质量分数)、Mo–10%Cu、Mo–20%Cu、Mo–30%Cu配制成混合粉末放进钢制球磨罐中,加入大小不同的钢球(球料比1:4),使用QM-3SP04型行星球磨机混合球磨5 h,将混合后的粉末在烘箱中烘干,最后将烘好的粉末倒入钢制模具中,采用YLJ-303立式液压机进行压坯成形;将Cu块放置在坯块上表面,在真空热压炉中进行熔渗处理,熔渗温度
1350 ℃,升温速率3 ℃·min−1,保温时间4 h,氩气气氛(0.06 MPa),制成Mo–Cu复合材料。1.3 表征及性能检测
测试不同诱导Cu含量试样的相对密度、电导率、热导率以及热膨胀系数。采用S4800型场发射扫描电子显微镜(field emission scanning electron microscope,FESEM)及其配备的牛津X250型X射线能谱仪(energy disperse spectroscope,EDS)对Mo–Cu复合材料试样显微组织、断口形貌和元素分布进行表征。采用D8 ADVANCE型X射线衍射分析仪(X-ray diffraction,XRD)测定复合材料的物相组成(Cu靶,扫描速度10°/min、扫描角度为20°~80°,步长为0.02°)。
选用标准硫代硫酸钠溶液(Na2S2O3),摩尔浓度为0.04 mol·L−1,通过化学分析法测试试样的最终Cu含量[18]。使用阿基米德排水法测试试样的相对密度。采用D60K型数字金属涡流电导率测量仪测量电导率,本文中所使用的电导率表示形式为相对电导率(% IACS),是实际电导率与国际标准软铜电导率的比值。采用LFA-457型激光导热系数测量仪对熔渗试样进行测量,测试前进行喷碳处理,根据标准参照样测得试样的热扩散系数,根据公式λ=αρc计算导热系数,其中λ为热导率,单位W∙m−1∙K−1;α为热扩散系数;ρ为试样密度,单位为g∙cm−3;c为比热容,单位J∙g−1∙K−1。采用DIL402SE型热膨胀仪测试Mo–Cu复合材料熔渗试样的热膨胀曲线,升温速率5 ℃∙min−1,温度测量范围为5~500 ℃。
2. 结果与讨论
2.1 诱导铜含量对于微观组织的影响
2.1.1 微观组织分析
图2为不同质量分数诱导Cu在
1350 ℃熔渗后所得Mo–Cu复合材料的微观组织形貌,其中较光滑的相为Mo,相对粗糙相为Cu。由图2(a)可见,当不添加诱导Cu时,组织内部存在大量的细小孔隙,Mo颗粒烧结合并成形成互连网状形貌,Cu相最终的含量较少。这是由于Mo骨架无诱导Cu,Cu液润湿性差,难以充分渗入骨架中的孔隙;另一方面Mo–Mo颗粒合并引起骨架收缩,熔渗通道缩小甚至关闭,最终形成较细但是数量较多的孔隙。当诱导Cu质量分数为10%时,见图2(b),组织内部孔隙数量减少,复合材料的组织更加均匀,Cu相的相对量显著提升,但是由于诱导Cu的量较低,坯体中仍存在Mo–Mo颗粒直接接触并在熔渗过程中合并联结,导致部分区域仍存在熔渗不足的现象。当诱导Cu质量分数为20%时,见图2(c),组织内部孔隙数量进一步减小,复合材料中的Mo和Cu两相分布更加均匀。一方面得益于诱导Cu含量的增加有效减少了初始加热阶段Mo–Mo的直接接触,防止Mo骨架收缩形成的局部封闭气孔;另一方面,诱导Cu通过高温阶段自身熔化改变骨架内部孔隙尺寸以及表面状态,使熔渗阶段进入的Cu液更易渗入Mo骨架。因此Mo和Cu两相分布更加均匀,Mo或Cu偏聚现象和孔隙都显著减小,复合材料中Cu的最终相对量显著增加。但当诱导Cu质量分数上升到30%时,如图2(d),坯体中Cu容易偏聚,并导致熔渗后大量的“Cu池”出现,反而降低了对Mo颗粒的协调作用。图4为图2(c)中P点的元素点分析和X射线衍射图谱。由图4可知,复合材料中两相分别由Mo和Cu元素组成,对应的X射线衍射图谱也说明复合材料中组成相为Cu和Mo,熔渗过程中没有引入其他杂质元素。从图3中可以看出Mo和Cu两相分布均匀,Mo相呈现网络化结构,这也和图2(c)的观察结果相对应。
图5为诱导Cu含量对于Mo–Cu复合材料熔渗后最终Cu含量的影响规律和相对密度。复合材料中最终Cu含量随初始诱导Cu含量的增加呈现先升高后降低的趋势,当初始诱导Cu质量分数达到20%时,试样最终Cu质量分数达到最大值40.46%。最终Cu含量也是影响Mo–Cu复合材料热性能的一个重要因素[19]。不同质量分数诱导Cu所制备试样的相对密度变化曲线如图5所示,可以看出随着诱导Cu含量的增加Mo–Cu复合材料的相对密度先增加后减小,这个趋势也与微观形貌观察结果相对应(见图2)。当诱导Cu质量分数为20%时,试样相对密度达到最大,为98.1%。由此可见,合适的诱导Cu的存在,不但降低了Mo–Mo直接接触引起的孔隙和偏聚倾向,而且在熔渗过程中有效调节孔隙尺寸以增进Cu液渗入,可以显著提升复合材料的相对密度。当诱导Cu质量分数为30%时,相对密度降低,可归结为局部两相分离,导致Mo–Mo区域出现闭孔,Cu液无法渗入的结果。
图6为不同诱导Cu含量试样断口组织形貌图。当诱导Cu质量分数为0%时,明显看出Mo颗粒的堆积,并且存在大量的孔隙,见图6(a),断裂方式主要是Mo颗粒之间的脆性解理断裂。主要原因是由于诱导Cu质量分数为0%时,坯块内没有Cu相,并且Mo和Cu相两相的润湿性较差,Cu液难以渗入坯块的孔隙当中,导致Mo–Mo颗粒逐渐烧结合并,孔隙未能被充分填充。由图6(b)可以看出,当诱导Cu质量分数增加至10%时,Cu相呈现部分区域网状结构,韧窝断裂面积增大,Mo–Cu复合材料韧性增强,但Mo相内部仍存在大量的孔隙,可能是由于诱导Cu含量不足,局部出现Mo–Mo颗粒直接接触形成的骨架结构,浸渗时Cu液无法完全填充孔隙所致。当诱导Cu质量分数为20%时,Cu相空间网状结构更加均匀、连续,见图6(c),断口中两相均匀性明显得到改善,Mo相内部大量孔隙消失。当诱导Cu质量分数为30%时,见图6(d),断口上可以观察到Mo和Cu的偏聚现象,孔隙数量再次增多,但孔隙的尺寸明显减小。
2.2 诱导Cu含量对Mo–Cu合金性能影响
2.2.1 电导率和热导率
图7为不同Cu含量试样电导率和热导率测试结果。图7可知,当最终Cu质量分数为40.46%时,试样电导率达到最大值52.69% IACS,热导率达到最大值203.94 W·m−1·K−1。当最终Cu质量分数为33.70%时,电导率和热导率均最小,因为此时Cu含量最小,Cu含量是影响Mo–Cu复合材料电导率和热导率的最主要原因之一;另一方面当不添加诱导Cu时,复合材料组织内部存在大量的细小孔隙,且组织中的Cu相此时也没有形成连续的网状结构,存在大量的Mo颗粒聚集区,这些因素都可能导致复合材料的电导率、热导率相对不高。当坯块最终Cu质量分数达到40.46%时,此时电导率和热导率达到最大值。一方面此时Cu含量最高,电子的扩散阻力随Cu含量的增加而降低;另一方面因为当诱导Cu质量分数为20%时,此时复合材料内部大量孔隙消失,Mo和Cu两相组织均匀性更高,相对密度高[20],尤其是Cu网络结构更加均匀、连续,有利于热和电的传导,因此电导率和热导率提升。在最终Cu质量分数相近的情况下(≤40%),诱导熔渗法制备的Mo–Cu复合材料的电导率和热导率要高于缺碳预还原+氢气深脱氧方法制备的Mo–Cu复合材料(电导率49.48% IACS,热导率187.2 W·m−1·K−1)[21]。
2.2.2 热膨胀系数
图8为不同最终Cu质量分数试样的热膨胀系数变化曲线,从图中可以看出随着测试温度的升高,Mo–Cu复合材料的热膨胀系数均逐渐升高。这是因为随着温度升高,原子的热振动加剧,导致热膨胀系数增加。理论上Mo–Cu复合材料理想的组织结果应该是Mo颗粒形成连续骨架,Cu相在其间隙之中呈均匀、连续的网络分布。连续的Mo骨架是Mo–Cu复合材料保持低热膨胀系数的重要因素。但需要指出的是,尽管最终Cu质量分数为40.46%时,Mo–Cu复合材料的相对密度最高,孔隙数量明显降低,两相分布更加均匀,Mo相对Cu相的约束效应更充分,但是由于此时样品最终Cu质量分数最高,Cu的热膨胀系数远大于Mo的热膨胀系数,因此相比之下Cu质量分数为40.46%时的样品在测试温度范围内的热膨胀系数最大,室温热膨胀系数为9.4×10−6 K−1。
3. 结 论
(1)采用质量分数分别为0%、10%、20%和30%诱导Cu制备出了不同Cu含量的Mo–Cu复合材料;
(2)随着诱导Cu质量分数的增加,复合材料的Mo、Cu两相分布更加均匀,孔隙率大大降低,诱导Cu质量分数为20%时,复合材料的组织相对较优;
(3)在本研究的范围内,随着最终Cu质量分数的增加,Mo–Cu复合材料的电导率、热导率均增大,在最终Cu质量分数为40.46%时Mo–Cu复合材料电导率达到最大52.69% IACS,热导率最大为203.94 W·m−1·K−1;但是,此时复合材料的热膨胀系数也最大;通过调控诱导Cu质量分数,辅以适宜的熔渗工艺,有望制备出综合性能更佳的Mo–Cu复合材料。
-
图 15 不同内径石墨模具烧结的
$ {\text{ln}}\left[ {\left( {{\text{d}}\dot \rho /{\text{d}}T} \right) \cdot T} \right] $ 与1/T关系曲线:(a)8 mm;(b)10 mm;(c)25 mm;(d)30 mmFigure 15. Relationship curves of
$ {\text{ln}}\left[ {\left( {{\text{d}}\dot \rho /{\text{d}}T} \right) \cdot T} \right] $ and 1/T for the graphite dies with different inner diameters: (a) 8 mm; (b) 10 mm; (c) 25 mm; (d) 30 mm表 1 Cu的空位形成能
Table 1 Vacancy formation energy of Cu
eV EVa EPer EVF ‒397.78 ‒402.83 1.32 表 2 Cu扩散激活能
Table 2 Diffusion activation energy of Cu
eV 激活能 EVF ESAM EDA 计算值 1.32 0.79 2.11 实验值[23] 1.22±0.02 0.71 2.07 1.29±0.02 0.90 表 3 无电场Cu界面扩散系统各层的空位形成能
Table 3 Vacancy formation energy for each layer of the interfacial diffusion system without electric field
eV Cu界面扩散层 EVa EPer EVF 第1层 ‒195.02 ‒199.32 0.74 第2层 ‒194.61 1.15 第3层 ‒194.44 1.32 第4层 ‒194.44 1.32 -
[1] Choi J, Sung H M, Roh K B, et al. Fabrication of sintered tungsten by spark plasma sintering and investigation of thermal stability. Int J Refract Met Hard Mater, 2017, 69: 164 DOI: 10.1016/j.ijrmhm.2017.08.013
[2] Terentyev D, M Vilémová, Yin C, et al. Assessment of mechanical properties of SPS-produced tungsten including effect of neutron irradiation. Int J Refract Met Hard Mater, 2020, 89: 105207 DOI: 10.1016/j.ijrmhm.2020.105207
[3] Deng S H, Li J N, Li R D, et al. The effect of particle size on the densification kinetics of tungsten powder during spark plasma sintering. Int J Refract Met Hard Mater, 2020, 93: 105358 DOI: 10.1016/j.ijrmhm.2020.105358
[4] Li J, Ding S Y, Zhu J S, et al. Current-density dependence of activation energy in high-Tc superconductor HgBa2Ca2Cu3O8+ x. J Appl Phys, 1995, 77(12): 6398 DOI: 10.1063/1.359113
[5] Garay J E, Anselmi-Tamburini U, Munir Z A. Enhanced growth of intermetallic phases in the Ni–Ti system by current effects. Acta Mater, 2003, 51(15): 4487 DOI: 10.1016/S1359-6454(03)00284-2
[6] Kondo T, Kuramoto T, Kodera Y, et al. Enhanced growth of Mo2C formed in MoC diffusion couple by pulsed DC current. J Jpn Soc Powder Powder Metall, 2008, 55(9): 643 DOI: 10.2497/jjspm.55.643
[7] Deng S H, Li R D, Yuan T C, et al. Direct current-enhanced densification kinetics during spark plasma sintering of tungsten powder. Scr Mater, 2018, 143: 25 DOI: 10.1016/j.scriptamat.2017.09.009
[8] 娄鹤子, 王海滨, 刘雪梅, 等. WC-Co硬质合金摩擦磨损行为的分子动力学模拟. 粉末冶金技术, 2022, 40(5): 8 Lou H Z, Wang H B, Liu X M, et al. Molecular dynamics simulation on friction and wear behavior of WC–Co cemented carbides. Powder Metall Technol, 2022, 40(5): 8
[9] Mukherjee S, Banwait A, Grixti S, et al. Adsorption and diffusion of lithium and sodium on defective rhenium disulfide: A first principles study. ACS Appl Mater Interfaces, 2018, 10(6): 5373 DOI: 10.1021/acsami.7b13604
[10] Yuan Z H, Kong X G, Ma S G, et al. Adsorption and diffusion mechanism of hydrogen atom on the Li2O (111) and (110) surfaces from first principles calculations. J Nucl Mater, 2018, 513: 232
[11] Lu H J, Wu H, Zou N, et al. First-principles investigation on diffusion mechanism of alloying elements in dilute Zr alloys. Acta Mater, 2018, 154: 161 DOI: 10.1016/j.actamat.2018.05.015
[12] Lu C, Yang J, Zhao Y, et al. Influence of applied electric field on atom diffusion behavior and mechanism for W/NiFe interface in diffusion bonding of steel/NiFe interlayer/W by spark plasma sintering. Appl Surf Sci, 2020, 541: 148516
[13] Zhang L Y, Peng C T, Shi J, et al. Surface alloying of chromium/tungsten/stannum on pure nickel and theoretical analysis of strengthening mechanism. Appl Surf Sci, 2020, 532: 147477 DOI: 10.1016/j.apsusc.2020.147477
[14] Kohn W, Sham L J. Self-consistent equations including exchange and correlation effects. Phys Rev, 1965, 140(4A): A1133 DOI: 10.1103/PhysRev.140.A1133
[15] Peng H W, Perdew J P. Rehabilitation of the Perdew-Burke-Ernzerhof generalized gradient approximation for layered materials. Phys Rev B Condens Matter, 2017, 95(8): 081105 DOI: 10.1103/PhysRevB.95.081105
[16] You J S, Fang S A, Xu S Y, et al. Berry curvature dipole current in transition metal dichalcogenides family. Phys Rev, 2018, 98(12): 121109 DOI: 10.1103/PhysRevB.98.121109
[17] Hwang J, Noh S H, Han B. Design of active bifunctional electrocatalysts using single atom doped transition metal dichalcogenides. Appl Surf Sci, 2019, 471: 545 DOI: 10.1016/j.apsusc.2018.11.147
[18] Zhang X, Grabowski B, Körmann F, et al. Accurate electronic free energies of the 3 d, 4 d, and 5 d transition metals at high temperatures. Phys Rev B Condens Matter, 2017, 95(16): 165126 DOI: 10.1103/PhysRevB.95.165126
[19] Kittel C. Introduction to Solid State Physics. 8th Ed. Weenheim: Wiley, 2004
[20] wang Y, Curtarolo S, Jiang C, et al. Ab initio lattice stability in comparison with CALPHAD lattice stability. Calphad, 2004, 28(1): 79 DOI: 10.1016/j.calphad.2004.05.002
[21] Naghavi S S, Hegde V I, Wolverton C. Diffusion coefficients of transition metals in fcc cobalt. Acta Mater, 2017, 132: 467 DOI: 10.1016/j.actamat.2017.04.060
[22] Esfandiarpour A, Nasrabadi M N. Vacancy formation energy in CuNiCo equimolar alloy and CuNiCoFe high entropy alloy: ab initio based study. Calphad, 2019, 66(1): 101634
[23] 王顺花. Cu中空位的计算机模拟. 兰州铁道学院学报, 1997, 16(4): 43 Wang S H. Computer simulation of vacancies in Cu. J Lanzhou Railway Inst, 1997, 16(4): 43
[24] Makov G, Shah R, Payne M C. Periodic boundary conditions in ab initio calculations. II. Brillouin-zone sampling for aperiodic systems. Phys Rev B Condens Matter, 1996, 53(23): 15513
[25] Yu X X, Thompson G B, Weinberger C R. Influence of carbon vacancy formation on the elastic constants and hardening mechanisms in transition metal carbides. J Eur Ceram Soc, 2015, 35(1): 95 DOI: 10.1016/j.jeurceramsoc.2014.08.021
[26] Wang J, Raj R. Activation energy for the sintering of two-phase alumina/zirconia ceramics. J Am Ceram Soc, 1991, 74(8): 1959 DOI: 10.1111/j.1151-2916.1991.tb07815.x
[27] 张剑平, 左红艳, 罗军明, 等. 微波场中铜粉烧结颈形成和晶粒长大动力学. 材料热处理学报, 2018, 39(6): 7 Zhang J P, Zuo H Y, Luo J M, et al. Neek formation and grain growth kinetics of copper powder during microwave sintering. Trans Mater Heat Treat, 2018, 39(6): 7
[28] Demuynck M, Erauw J P, Biest O V D, et al. Densification of alumina by SPS and HP: A comparative study. J Eur Ceram Soc, 2012, 32(9): 1957 DOI: 10.1016/j.jeurceramsoc.2011.10.031