大块金属玻璃 (BMG) 由于其特殊的结构, 具有机械性能好、抗腐蚀性能高、软磁性能优异和储氢能力高等优点, 不仅有理论研究意义, 而且有工程应用价值, 因而备受关注, 成为材料领域的研究热点之一。A.Inoue提出了BMG合金配制的三大原则:多组元、原子半径差比率大于12%以及组元间混合热负值要大[1]。多组元能够提高大块金属玻璃的玻璃形成能力, 这已被实验验证并被人们用来作为BMG合金成分配制的主要途径。近来大连理工大学利用原子半径和电子浓度的观点制备出了一系列Zr-Al-Ni-Cu大块非晶[2], 但是, 迄今还没有人研究原子半径差和化学键参数与BMG的玻璃形成能力 (GFA) 之间的具体关系。

笔者用三种方法计算了Pd基、Zr基、Mg基以及La基BMG合金原子半径差比率Δd、原子电负性差比率Δe以及化学键离子率v等三个物理参数, 以大块金属玻璃形成的临界冷却速度Rc作为合金GFA的评价参量 (见表1) [3,4,5], 通过Origin 软件处理这三个物理参量与Rc的关系, 发现四种BMG的这三个物理参量与Rc的关系不是线性关系, 而是抛物线关系。

《1 BMG合金物理参数的选取与计算》

1 BMG合金物理参数的选取与计算

BMG作为热力学上亚稳定的重要工程材料, 其形成是一个动力学问题。任何动力学过程的进行总需要克服一定的能量, 这个能量通常称为势垒或激活能。若此势垒高, 则质点活化到足以克服此势垒的概率就小, 所需要的时间就长;反之, 概率就大。结晶势垒是描述物质由非晶态转变成晶态所需要克服的能量。结晶势垒大小与物质的结构化学特性, 如键性、键强、原子大小以及原子排列几何方式等有关。化学物质的结构与性能主要由其组成的原子性质和化学键的性质所决定, 化学键表示原子间的作用力。Pauling 原子电负性是由热化学数据得出的, 因此, 原子之间的这种电负性差能反映化学键的键强和混合热的大小, 并且电负性与键性之间有定量关系, 笔者采用Hooydank关系式来计算v以表示化学键的键性[6], Δd和化学键键性的搭配能够反映原子排列的几何方式, 进而可以判定分子的空间结构。Rc是反映合金非晶形成能力

表1 BMG物理参数和临界冷却速度 (Rc) Table 1 Physical parameters and critical cooling velocity (Rc) of BMG

《表1》

序号合金名称Δd1Δd2Δd3Δe1Δe2Δe3v1v2v3Rc /K·s-1
1Mg75Ni15Nd100.864 330.231 190.573 810.732 670.298 610.723 460.581 40.174 630.580 56 46.1
2Mg65Ni20Nd150.809 60.238 390.527 120.619 050.343 750.622 860.454 530.207 550.463 5530.0
3Mg70Ni15Nd150.865 40.237 260.634 010.731 340.343 140.623 680.581 40.207 10.475 51178.2
4Mg65Cu25Y100.738 130.214 40.649 520.472 220.368 420.661 520.309 090.20.524 3350.0
5Mg80Ni10Nd100.914 50.230 510.691 110.831 780.339 510.695 60.7120.204 460.573 331 212.4
6Zr57Ti5Al10Cu20Ni80.687 690.196 920.619 030.487 780.240 780.610 950.322 560.136 870.478 4310.0
7Zr65Al7.5Cu17.5Ni100.699 20.206 990.611 980.498 780.243 230.634 750.332 250.138 450.489 531.5
8Zr41.2Ti13.8Cu12.5Ni10Be22.50.380 390.242 690.478 310.628 220.220 270.398 280.457 370.123 770.266 861.4
9Zr66Al8Ni260.725 410.218 490.745 680.551 720.216 220.702 410.380 950.121 210.563 1866.6
10Zr66Al8Cu7Ni190.723 640.213 440.657 490.545 020.227 770.619 010.374 590.128 520.491 6822.7
11Zr66Al8Cu12Ni140.735 120.209 830.551 580.559 960.235 810.573 90.388 850.133 660.443 489.8
12Zr66Al9Cu16Ni90.722 370.205 360.588 50.540 230.242 490.589 880.370 080.137 970.455 364.1
13Zr38.5Ti16.5Ni9.75Cu15.25Be200.373 070.233 760.444 360.571 790.221 560.357 630.400 350.124 580.233 451.4
14Zr39.88Ti15.12Ni9.58Cu13.77Be21.250.382 570.245 360.465 660.606 770.221 240.388 580.435 520.124 380.256 311.4
15Zr42.63Ti12.37Ni10Cu11.25Be23.750.384 060.247 180.493 130.654 190.219 690.412 610.486 090.123 40.279 295.0
16Zr44Ti11Ni10Cu10Be250.387 650.251 570.502 490.679 930.218 920.436 360.515 070.122 910.310 5512.5
17Zr45.38Ti9.62Ni10Cu8.75Be26.250.391 240.255 950.523 410.704 910.217 970.476 90.544 290.122 320.347 1117.5
18Zr46.25Ti8.25Ni10Cu7.5Be27.50.394 790.260 30.550 650.725 320.216 80.525 850.570 740.121 580.390 928.0
19La55Al25Ni200.820 340.281 350.571 060.632 650.319 440.275 040.462 690.190 080.169 6667.5
20La55Al25Ni10Cu100.818 180.272 730.620 890.622 450.337 840.472 430.451 850.203 250.333 9922.5
21La55Al25Cu200.816 030.264 10.563 530.612 240.355 260.255 070.441 180.2160.156 5872.3
22La66Al14Cu200.822 350.289 380.612 120.594 020.384 210.544 90.422 490.237 780.384 0937.5
23La55Al25Ni5Cu10Co50.818 180.272 730.672 530.622 450.337 840.583 460.451 850.203 250.496 0318.8
24La55Al25Ni15Cu50.819 260.277 040.724 880.627 550.328 770.575 280.457 250.196 720.430 6834.5
25La55Al25Ni5Cu150.817 10.268 420.728 450.617 350.346 670.574 090.446 490.209 680.495 835.9
26Pd77.5Cu6Si16.50.830 440.141 930.782 670.758 940.169 70.791 70.611 530.092 720.673 68100.0
27Pd77Cu6Si170.824 250.141 910.784 170.752 070.169 960.793 170.602 650.092 870.692 07125.0
28Pd40Ni40P200.791 980.167 940.417 930.446 150.169 230.373 740.287 130.092 440.2390.167
29Pd40Cu30Ni10P200.793 760.175 020.534 460.423 080.134 620.531 540.268 290.072 170.383 26 0.1
30Pd79.5Cu4Si16.50.830 710.143 290.836 840.786 740.172 950.843 610.648 440.094 660.739 63500.0

的最现实的参数, Rc越小, 表明合金非晶形成能力越强。笔者基于加权平均原理 (权重为组元的原子百分率) , 用三种不同方法分别计算了Δd, Δev:①把合金体系的组元看作两大类, 即原子半径大的和原子半径小的, 电负性大的和电负性小的;②在第一种方法的基础上消除权重的影响;③把合金中的组元两两配对进行计算, 然后平均。三种方法处理Δd的计算公式分别为

Δd1=i=1n1(xidi)-j=1n2(xjdj)i=1n1(xidi)(1)Δd2=[i=1n1(xidi)/i=1n1xi]-[j=1n2(xjdj)/j=1n2xj]i=1n1(xidi)/i=1n1xi(2)Δd3=i,j=1,ijn[xjdi-xjdj)/(xidi)]Cn2(3)

式中, Δd1, Δd2, Δd3的下标表示采用的方法, xi, din1分别为合金中大原子的原子百分含量、半径和原子个数;xj, djn2分别为合金中小原子的原子百分含量、半径和原子个数;n为合金中的组元数;Cn2为排列组合数;n=n1+n2

Δe的三种计算方法与Δd计算公式的形式完全相同。v的三种处理方法分别为:

v1=i=1n1(xiei)-j=1n2(xjej)i=1n1(xiei)+j=1n2(xjej)(4)

v2=[i=1n1(xiei)/i=1n2xi]-[j=1n2(xjej)/j=1n2xj][i=1n1(xiei)/i=1n1xi]+[j=1n2(xjej)/j=1n2xj](5)

v3=i,j=1,ijn[(xiei-xjej)/(xiei+xjej)]Cn2(6)

式中v1, v2, v3的下标表示采用的方法, xi, ein1分别为合金中大电负性原子的原子百分含量、电负性和原子个数;xj, ejn2分别为合金中小电负性原子的原子百分含量、电负性和个数;n为合金中的组元数;Cn2为排列组合数;n=n1+n2

《2 结果与分析》

2 结果与分析

Mg基、Zr基、La基以及Pd基等四种BMG合金的Δd, Δev与BMG合金玻璃形成Rc的关系及其拟合曲线见图1至图3, 拟合曲线方程和相关系数 (R) 见表2。

从图1至图3和表2可知, 所有拟合曲线为二

《图1》

图1 Mg基、Zr基、La基、Pd基四种BMG
的Δd与玻璃形成Rc的关系
Fig.1 The relation between Rc forming
 glass and Δd of Mg-based, Zr-based, 
 La-based and Pd-based BMG

图1 Mg基、Zr基、La基、Pd基四种BMG 的Δd与玻璃形成Rc的关系 Fig.1 The relation between Rc forming glass and Δd of Mg-based, Zr-based, La-based and Pd-based BMG  

次或三次型曲线, 不管用什么处理方法, Δe , Δdv等三个参量与BMG合金的玻璃形成临界冷却速度 (Rc) 都是在三个参量达到某个数值

《图2》

图2 Mg基、Zr基、La基、Pd基四种BMG
的Δe与玻璃形成Rc的关系
Fig.2 The relation between Rc forming
 glass and Δe of Mg-based, Zr-based, 
 La-based and Pd-based BMG

图2 Mg基、Zr基、La基、Pd基四种BMG 的Δe与玻璃形成Rc的关系 Fig.2 The relation between Rc forming glass and Δe of Mg-based, Zr-based, La-based and Pd-based BMG  

Rc有极小值, 只是对于不同的合金体系, Rc出现极小值所对应的上述三个物理参量的数值不同而已, 表明合金玻璃形成能力与这三个参数之间的关系并不是越大越好, 或越小越好。根据Pauling电负性的定义可知, 电负性差值越大, 原子混合热也是越大的, 所以, 认为原子半径差和混合热越大, 合金非晶形成能力也越大的观点值得进一步研究, 并且曲线拟合相关系数大部分在0.8以上, 可见三个参数与Rc间的极小值关系是可信的。

《图3》

图3 Mg基、Zr基、La基、Pd基四种BMG
的v与玻璃形成Rc的关系
Fig.3 The relation between Rc forming
 glass and v of Mg-based, Zr-based, 
 La-based and Pd-based BMG

图3 Mg基、Zr基、La基、Pd基四种BMG 的v与玻璃形成Rc的关系 Fig.3 The relation between Rc forming glass and v of Mg-based, Zr-based, La-based and Pd-based BMG  

从图1可知, 在五种BMG合金体系中, Rc与Δd间的关系都表现出了极小值, 这不是偶然的。根据近代物理的研究成果[7], 非晶态是具有长程无序而短程有序结构, 这种有序结构由Bernal空洞或Voronoi多面体组成, 如果原子半径差大到一定程度, 半径较小的原子将填入空洞或多面体间隙中, 而原子之间同时具有吸引力和排斥力, 空洞或多面体中填入半径小的原子后, 小原子与顶点原子之间的吸引力大于排斥力, 使得空洞或多面体变瘪, 从而使空洞或多面体之间的通道变宽, 给其他原子提供了较大的原子扩散通道, 降低了原子扩散势垒, 不利于非晶的形成;如果原子半径差小到一定程度, 原子的填充作用变得很弱, 此时虽然有很小一部分原子填入空洞或多面体中并使空洞或多面体扩大, 引起其邻近多面体或空洞之间的原子扩散通道变窄, 导致原子的扩散变得困难, 但是此时起决定作用的是原子之间所形成的空洞或多面体数量大大增加, 原子的扩散通道增加, 降低了原子的扩散势垒。由于后者起决定作用, 总的效果是原子扩散增加, 当然也不利于非晶的形成;只有当原子半径差恰当时, 才能克服上面两种情况的缺点, 有利于非晶的形成。例如Pd-Cu-Si和Mg-Ni-Nd两大体系的原子半径差在五大体系中是最大的, 并且其玻璃形成能力最差。如果把所有的物质统一来看, Pd-Cu-Si体系和Mg-Ni-Nd体系的原子半径差都处在最大位置, 因而其玻璃形成能力最差;Pd-Ni-P和Pd-Ni-Cu-P的原子半径差处于中间位置, 其玻璃形成能力最强;Zr-Al-TM (TM=Ni, Cu, 下同) 和La基两大体系的原子半径差比Pd-Ni-P和Pd-Ni-Cu-P小, 其玻璃形成能力比前两者弱。但是Zr-Ti-Be-TM体系的原子半径差最小, 其非晶形成能力比Zr-Al-TM和La基强, 这可能说明多项式关系是可信的, 因此可以预计原子半径差小于且靠近Zr-Ti-TM的合金体系, 其玻璃形成能力较强。从各自曲线的走势来看, Zr-Ti-Be-TM基、Pd基两大系统似乎已经达到了两大系统的极限, 而Zr-Al-TM基、Mg基和La基两大体系还有进一步开发的可能, 特别是La基体系表现出了三次方曲线关系, 可能存在冷却速度更低的BMG。

从图2和图3可知, Rc与Δe的关系都表现出了极小值的关系, 这也有它的必然性。由于两原子靠近时, 一般会形成化学键, 这种化学键通常同时包含共价键和离子键成分。共价键具有方向性和饱和性[8], 作用范围小, 键长和键角不易改变, 原子不易扩散, 有阻碍结晶的作用;离子键没有方向

表2 Mg基、Zr基、La基、Pd基四种BMG的三个参数 (Δd, Δe, v) 与其临界冷却速度 (Rc) 的拟合曲线方程及其平方相关系数 (R) Table 2 The fit curve equations and squared-coefficients (R) between three related parameters (Δd, Δe, v) and the Rc of Mg-based, Zr-based, La-based and Pd-based BMG

《表2》

BMG
参量
Mg 基Zr-Al-TM 基Zr-Ti-TM 基La 基Pd 基
ΔdΔevΔdΔevΔdΔevΔdΔevΔdΔev
A59.16E376.75E230.27E249.00E315.14E389.90E224.23E26.13E23.24E21.46E9-4.01E6-1.31E64.32E57.80E33.38E3
A3.201E61.225E661.44E212.12E350.16E240.33E232.16E252.27E341.77E38.86E51.47E58.49E43.84E37.58E46.33E4
A34.26E35.066E6-1.28E6807.6119.66E211.71E2461.9036.4568.9452.35E32.14E21.33E22.07E31.26E36.78E2
B1-14.92E4-26.71E3-14.69E3-13.94E4-58.96E3-51.93E3-6.05E3-20.36E2-14.53E2-5.36E91.96E78.87E6-1.07E6-2.86E4-1.79E4
B1-2.72E7-6.89E6-65.77E3-11.92E4-41.23E3-58.27E3-27.01E3-47.10E4-66.93E4-9.58E6-1.26E6-1.20E6-4.48E4-1.01E6-1.56E6
B1-11.87E4-2.28E77.53E6-27.50E2-6.54E3-50.74E2-20.94E2-2.90E2-170.2-6.99E3-7.38E2-4.79E27.65E3-5.03E3-3.71E3
B293.74E322.62E316.92E399.19E357.31E375.78E33.78E316.93E216.39E26.54E9-3.19E7-2.00E76.65E52.44E42.07E4
B25.79E79.68E617.66E429.28E484.79E321.07E456.74E31.06E62.68E63.45E73.60E65.64E61.31E53.33E69.51E6
B210.23E43.41E7-1.47E72.36E35.46E35.53E32.37E35.24E25.62E25.24E37.30E25.47E26.71E34.67E34.42E3
B3000000000-2.66E91.73E71.51E7000
B3000000000-4.15E7-3.42E6-8.80E6000
B30-1.70E79.58E6000000000000
R0.904010.897340.935550.550480.867550.865960.991250.991680.954650.936020.960260.955810.548410.853420.92255
R0.504151.000000.58350.982970.987030.986760.987580.98190.982370.934130.944070.950810.976140.993030.9928
R0.862120.980280.891920.979080.880150.925720.966050.966060.979780.976010.875420.848110.743190.687960.70413

注:A, B1, B2, B3分别表示拟合多项式的零次、一次、二次以及三次方项的常系数;①、②和③分别表示数据处理方法

性和饱和性, 离子晶体倾向于紧密堆积, 离子键占优势的物质容易形成晶体, 不易形成非晶体。在化合物中, 电负性相差大的元素都以离子键结合, 电负性小的元素则以共价键结合, 居于两者之间的是以离子-共价混合键结合。由离子-共价混合键组成的物质, 既有离子晶体容易变更键角、易造成无对称变形的趋势, 又有共价键不易变更键长和键角的趋势, 最容易成为玻璃。前者造成玻璃结构的长程无序, 后者造成结构的近程有序。根据量子化学理论[9], 当电负性差大到一定程度时, 离子率大, 共价键所占的比重就少。按照原子团簇理论可知, 此时所形成的原子团簇大而松散, 虽然熔体的粘度提高了, 但是这种松散的结构使得原子扩散变得容易;根据高温熔体原子分子共存理论[10], 化合物离子性越大, 它在高温熔体中稳定存在的概率就越大, 如果冷却速度不是足够大, 就会从熔体中以晶体形式析出。综合这两点可以得出, 电负性差大到一定程度时不利于非晶的形成;电负性差小到一定程度时离子率小, 共价键所占的比重大, 所形成的原子团簇小而数量多, 虽然原子在团簇内部不易扩散, 并且熔体的粘度可能有所增加, 但是, 团簇之间通道数的增加, 又有利于原子的扩散, 这种小分子数量越多, 其在高温熔体中存在的概率就越大, 熔体在凝固过程中结晶的概率就大增。因此, 总体来说, 不利于非晶的形成。只有电负性差适中、离子率即离子键含量适当时, 所形成团簇的尺寸、松紧程度和数量适中, 熔体的粘度增加, 原子的扩散不很大, 才是有利于非晶形成的最佳状态。例如, Pd-Cu-Si和Mg-Ni-Nd这两大体系就是一个在电负性和离子率两方面走极端的很好例子。如果把所有BMG合金作为一个物质系统, Pd-Cu-Si体系和Mg-Ni-Nd体系的Δev都同时处在最大处, 其玻璃形成能力差;Zr-Ti-TM系处于中间, 其玻璃形成能力比前两者大;Pd40Ni40P20和Pd40Ni10Cu30P20处于最小的位置, 玻璃形成能力最强。这又一次表明上述的多项式关系是可信的, 从而可以预计Δev比这两个大并接近它们的合金, 其玻璃形成能力是比较好的。综合ΔevRc关系的走势可知, Pd基和Zr-Ti-Be-TM基两大系统也已经达到了极限;而Mg基、La基和Zr-Al-TM基还有希望获得新的形成能力比现有的好的BMG。

总之, 要正确地判定BMG合金的玻璃形成能力, 必须综合Δd, Δev三个参数来考察, 例如Pd-Cu-Si和Mg-Ni-Nd两大体系的三个物理参数都同时走极端, 因而其玻璃形成能力在五大体系中是最差的;而Zr-Al-TM和La基两大体系的上述三个参数都比前面两大体系小, 其玻璃形成能力比前两大体系强;Pd-Ni-Cu-P和Pd-Ni-P三个物理参量处在五大体系的中间位置, 其玻璃形成能力最强。Zr-Ti-Be-TM体系三个物理参量在五大体系中是最小的, 虽然其玻璃形成能力较强, 但比Pd-Ni-Cu-P和Pd-Ni-P差。这又一次表明, 玻璃形成能力强的合金体系, 其各个组元的原子百分含量、组元原子之间的半径差以及电负性差这三者之间的搭配是和谐的, 这样形成的团簇分子数量、分子构型以及小原子填充效应等方面的搭配是最佳, 置换元素的最佳添加量规则可能与此规律有关。

《3 结论》

3 结论

1) Δd, Δev与合金的玻璃形成能力都出现极小值现象, 对于不同的合金体系, 极小值出现的位置不同。

2) 玻璃形成能力强的合金体系, 其各个组元的原子百分含量、组元原子之间的半径差以及电负性差这三者之间的搭配是最和谐的。

3) 金属玻璃合金的最优成分范围以及各种添加元素的最佳添加量现象, 可能与Δd, Δev出现极小值现象有关。