中国有色金属学报 2003,(01),106-110 DOI:10.19476/j.ysxb.1004.0609.2003.01.019
合金化对非晶形成能力影响的电子理论计算
张国英 胡壮麒 张海峰
沈阳师范大学物理系,中国科学院金属研究所材料科学国家实验室,中国科学院金属研究所材料科学国家实验室 沈阳110031中国科学院金属研究所材料科学国家实验室,沈阳110016 ,沈阳110016 ,沈阳110016
摘 要:
根据分子动力学理论建立了非晶Ni64Pd16Al2 0 的结构模型 ,利用Recursion方法计算了Ni64X16Al2 0 系非晶及相应晶态的电子结构。结果表明 ,Ni64X16Al2 0 合金的非晶形成能力与过渡元素局域电子态密度的双峰位置、电子得失形成的离子键及非晶与晶态合金的总结构能差有关。局域电子态密度的结果表明 ,Al与过渡元素间存在电子耦合—共价相互作用 ,合金元素使Ni64X16Al2 合金非晶形成能力从强到弱的顺序是 :Ir,Pt,Rh ,Pd ,Au ,Ag ,Cu ;Ni64X16Al2 0 合金中离子键的存在使其非晶形成能力增强。非晶与晶态合金的总结构能差越小的过渡元素使其形成非晶的驱动力越强 ,对非晶形成越有利 ,由此得出的规律与上述共价键的作用规律一致。
关键词:
非晶结构模型 ;Recursion方法 ;电子理论 ;非晶形成能力 ;
中图分类号: TG139.8
作者简介: 张国英(1965),女,副教授,博士;
收稿日期: 2002-04-22
基金: 国家重点基础研究发展规划资助项目 (G2 0 0 0 0 672 0 1); 国家高技术研究发展规划资助项目 (2 0 0 1AA3 3 10 10 );
Electronic theory calculations of alloying effect to amorphous alloy forming ability
Abstract:
The atom structure model of Ni 64 P d 16 Al 20 amorphous alloy was set up by using molecular dynami cs theory . The electronic structure calculations of crystal and amorp hous alloy Ni 64 X 16 Al 20 by Recursion meth od reveal that the amorphous alloy forming ability is related to the location of the double-peak (or tri-peak) of local density of states (LDOS) in low energy region, electrovalent bond formed by gaining or losing electrons and the total structure energy difference of amorphous alloy Ni 64 X 16 Al 20 to crystal alloy. The common conclusion drawn from LDOS and the t otal structure energy difference of amorphous alloy Ni 64 X 16 Al 20 to crystal alloy is: when transition element Ir, Pt, Rh, Pd, Au, Ag, Cu is added into NiAl alloy, the amorphous alloy forming ability we akens in the same order. The electrovalent bond existing in Ni 64 X 16 A l 20 alloy enhances the amorphous alloy forming ability too.
Keyword:
amorphous structure model; Recursion method; electronic th eory; amorphous alloy forming ability;
Received: 2002-04-22
非晶合金具有很多奇异的性能
[1 ]
, 比如, 某些非晶合金具有很高的强度
[2 ]
, 某些合金有很好的磁性
[3 ]
等, 这些都使得非晶材料的研究成为材料科学中一个非常活跃的领域。
自从发现非晶以来, 为使其优良性能得到更好的应用, 人们一直在尝试制备大块状的非晶材料。 通过实验, 人们发现不同金属或合金形成非晶的能力相差甚远, 合适的合金化能提高非晶的形成能力。 对于大块非晶, 人们发现除了要施加很大的冷却速度、 很大的过冷度外, 在合金成分选择上也要遵循以下原则:一是要多组元, 二是原子尺寸比要大, 三是要有大的负混合焓。 比如, Nishiyama
[4 ]
研制成功的直径达72 mm 的Pd40 Cu30 Ni10 P20 非晶合金,其临界冷却速率为0.1 K/s就够。Inoue
[5 ]
, Peker
[6 ]
, 张庆生
[7 ]
等人还研制成功了大块Zr基非晶。
随着对非晶材料(尤其是大块非晶)研究的深入, 急需从理论上研究多组元非晶材料的形成能力, 以便对实验进行指导。 电子理论在材料设计中已得到了较广泛的应用
[8 ,9 ,10 ]
。 非晶材料就其结构而言, 是固化了的液体, 长程无序, 短程有序。 由于非晶不具有周期性, 目前的理论研究还不能从第一原理出发, 研究非晶的形成过程。 Recursion
[11 ]
方法是基于紧束缚近似下的电子理论方法, 它在推导过程中不需要材料具有周期性和对称性, 因此对于处理非周期性材料(如非晶、 准晶、 晶界
[9 ]
等)具有优越性。 本文作者将利用Recursion方法, 从电子理论方面研究多组元大块非晶的形成能力。
1 理论方法及非晶模型的构造
Recursion方法的基本思想是根据紧束缚近似建立哈密顿矩阵, 然后把系统的哈密顿做一次幺正变换, 变换后的哈密顿矩阵是三对角化的。 由此定义的实空间局部格林函数为
G (E )=(E -H )-1 (1)
格点态密度可由初态格林函数的虚部求出, 即
ρ 0 = ? 1 π I m ? u 0 | 1 E ? H | u 0 ? ? ? ? ( 2 )
ρ
0
=
-
1
π
Ι
m
?
u
0
|
1
E
-
Η
|
u
0
?
?
?
?
(
2
)
对于任意一个格点都可以得到相应的(2)式, 同时总态密度为各格点的态密度之和, 而局部态密度即是合金元素处的格点态密度。
由上面的推导可知, Recursion方法在计算过程中不要求系统具有周期性。 在计算过程中, 哈密顿对角矩阵元-原子轨道能取自Fischer用Hartree-Fock近似计算的结果。 我们之所以这样选是因为Hartree-Fock近似对过渡性元素给出了较好的结果, 哈密顿非对角矩阵元-原子间的跃迁积分取为Slater-Koster积分, 普适参数取自固态表
[12 ]
, 链长取为30。 Ni, Pd, Cu, Ag, Rh, Ir, Pt的价电子组态取为4s2 3d8 , 5s2 4d8 , 4s1 3d10 , 5s1 4d10 , 5s2 4d7 , 6s2 5d7 , 6s2 5d8 , 6s1 5d10 , Al的价电子组态取为3s2 3p1 。
在Recursion方法和紧束缚框架下, 体系的结构能为
U s = ∑ l U l ? ? ? ( 7 )
U
s
=
∑
l
U
l
?
?
?
(
7
)
U l = ∑ α ∫ ? ∞ E f E n α l ( E ) d E ? ? ? ( 8 )
U
l
=
∑
α
∫
-
∞
E
f
E
n
α
l
(
E
)
d
E
?
?
?
(
8
)
根据不同体系总结构能的大小, 可以讨论体系形成非晶的难易。
Recursion方法由于在构造哈密顿矩阵时不需体系具有周期性和对称性, 原则上可以选择很多原子构成原子团, 原子种类也可以选为多种。 分子动力学在模拟快速凝固形成非晶的过程方面有非常独特的效果。 本文用分子动力学方法构造非晶的原子团, 以Ni64 Pd16 Al20 的成分为参考, 由于Ni、 Pd、 Al都为面心立方结构, 原子间相互作用势选为Sutton-Chen模型
[13 ]
。 体系包含500个原子, 周期性边界条件选为立方盒子。 初始温度取为1 800 K, 冷却速度为2×1013 K/s, 每隔100 K, 进行4 000步驰豫, 记取一次结构。 模拟结果发现, 非晶的转变温度约为670 K, 取300 K时非晶结构为电子理论计算的原子团(图1所示为Z 的绝对值小于0.325 nm的所有原子在Z 方向的投影)。 以Cu, Ag, Rh, Au, Pt, Ir代替Pd计算非晶的电子结构, 计算总结构能时取中心的25个原子, 同时计算
图1 非晶态Ni64Pd16Al20合金的原子结构模型
Fig. 1 Model of atom structure of Ni64 Pd16 Al20 amorphous alloy
相应的相同体积晶态的总结构能以作比较。
2 计算结果与讨论
2.1 Ni64X16Al20非晶与晶态的总结构能差
在材料物理学中, “能量”作为热力学和动力学以及稳定和临界现象的基本参量和判据, 可以反应微观体系的相互作用, 是一个“好”的特征量
[14 ]
。 Uhlmann
[1 ]
根据固体结晶过程中晶核形成、 晶体生长以及相变动力学理论得出: 晶核形成的热力学势垒(结晶势垒)愈大, 则从液态到非晶态的驱动力越大, 即形成非晶固体就愈容易。 晶核形成的热力学势垒与单位体积液体转变成晶体后自由能下降的数值成反比。 由此可知, 单位体积液体转变成晶体后自由能下降的数值越小, 形成晶体就越难, 形成非晶的驱动力就越大, 就越容易形成非晶。 体系的自由能 F =U +TS 。 熵S 与体系的混乱程度有关, 就前面利用分子动力学建立的非晶体系而言, 用不同的过渡元素代替Pd并不增加混乱度, 因此系统的熵不变, 那么自由能的变化主要是由内能决定的, 即体系的内能(晶态相对于非晶态)降低得越少, 越容易形成非晶。 我们利用Recursion方法计算了不同过渡金属代替Pd后, NiXAl系非晶相对于晶态的能量差, 结果列于表1。
从表1可以看出, 过渡元素使得NiXAl系非晶相对于晶态的能量差按Ir, Pt, Rh, Pd, Au, Ag, Cu的次序依次升高, 按Uhlmann的理论可知其非晶形成的驱动力依次降低, 非晶形成能力依次下降。
2.2Ni64X16Al20非晶与晶态合金中过渡元素与Al间的键性
前面指出了结晶势垒与能量间的关系, 另外结晶势垒与物质的结构化学特性, 如键性、 键强、 分子大小及排列也有关
[2 ]
,无方向性的纯粹离子键和金属键物质容易形成晶体, 共价键的键长及键角不易改变, 原子不容易扩散, 有阻碍结晶的作用。 当物质中存在离子-共价混合键时, 最容易形成非晶
表1 Ni64X16Al20系非晶相对于晶态的能量差(eV)
Table 1 Energy difference of Ni64 X16 Al20 amorphous alloy compared with crystal alloy (eV)
X
Ir
Pt
Rh
Pd
Au
Ag
Cu
E (Amorphous)
-3 846.51
-3 950.60
-3 884.06
-3 995.718
-4 055.62
-4 105.84
-4 142.58
E (Crystal)
-5 082.36
-5 217.97
-5 131.98
-5 279.76
-5 358.44
-5 425.20
-5 474.71
ΔE
1 235.85
1 267.37
1 247.92
1 284.04
1 302.79
1 319.36
1 332.14
体。 另外, 共价键的键强越强, 形成非晶的能力越强。
图2、 图3给出了Ni64 X16 Al20 非晶与晶态合金中, 过渡元素及 Al的局域态密度(LDOS)。 Al在高能区有一明显的小峰, 主要是其s轨道的贡献; 较低能区为p轨道贡献。 过渡元素在高能区也有一小峰, 为s电子的贡献, 低能区的双峰或三峰结构为d电子的贡献。 过渡元素态密度低能区的双峰结构(非晶中, 部分过渡元素变为三峰结构)按Cu, Ag, Au, Pd, Rh, Pt, Ir的顺序依次向高能区转移, 峰高降低, 能区展宽, 与Al的p电子在低能区发生耦合, 形成了共价键, 且后面的元素与Al形成的共价键较强; 高能区Al的s电子与过渡元素的s电子也发生耦合, 形成了共价键。 这种共价相互作用在液态中也存在, 这样, 液体在冷却过程中,
图2、 图3中下边的元素扩散相对于上边的元素要困难, 因此当它们加入到NiAl合金中时, 下边的元素使合金容易形成非晶, 即过渡元素对NiAl系合金非晶形成能力影响从大到小的顺序是: Ir, Pt, Rh, Pd, Au, Ag, Cu。 另外, 我们还根据费米能级计算了各元素的轨道电子数, 发现Al与过渡元素间存在着电子得失(见表2), 说明Al与过渡元素之间除了有共价相互作用外, 还存在离子键。 这就使Ni64 X16 Al20 合金的非晶形成能力更强。
2.3 与以往实验及理论结果的比较
人们通过实验和结构化学理论得出原子尺寸差及元素间电负性差对非晶固体的形成有很大影响,
图2 非晶Ni64X16Al20合金中合金元素的 局域态密度
Fig.2 Local density of state of alloy elements in amorphous alloy Ni64 X16 Al20
图3 Ni64X16Al20晶态合金中合金元素的 局域态密度
Fig.3 Local density of state of alloy elements in crystal alloy Ni64 X16 Al20
表2 Ni64X16Al20合金中Al与过渡元素间的得失电子数
Table 2 Gained or lost electron numbers of Al or transition elements in Ni64 X16 Al20 alloy
X
Δn x
Δn Al
Δn Ni
Ir
-0.187
1.207
-0.410
Pt
0.134
1.165
-0.428
Rh
-0.548
1.356
-0.356
Pd
-0.083
1.298
-0.386
Au
0.749
1.003
-0.486
Ag
0.613
1.070
-0.471
Cu
0.482
1.114
-0.460
原子尺寸差越大, 电负性差越大, 非晶的形成能力越强
[1 ]
。 就Ni64 X16 Al20 系合金而言, 除Cu以外, Au, Ag与Al的原子半径差较小(见表3), 当它们加入到NiAl合金中时, 其非晶形成能力较小; 相反, Ir, Rh, Pd, Pt与Al的原子半径差较大, 当它们加入到NiAl合金中时, 非晶形成能力强。 有关的电负性理论是, 两元素间的电负性差越大, 则相互作用越强。 在同一主族中, 从上到下, 电负性减小; 同一周期中, 从左至右, 电负性增大。 那么, Ir, Pt, Rh与Al的电负性差较大, Cu, Ag与Al的电负性差最小, 可见, 由此仍然可以得出, Ir, Pt, Rh能增大NiAl合金非晶的形成能力, Cu, Ag, Au对非晶的形成能力增加不大。 显然, 这与我们用电子理论计算得出的结论是一致的, 由此可
表3 Al及不同过渡元素的原子半径[15] (nm)
Table 3 Atom radii of Al and transition elements (nm)
Ir
Pt
Rh
Pd
0.135 7
0.138 7
0.134 5
0.137 5
Au
Ag
Cu
Al
0.144 2
0.144 4
0.127 8
0.143 1
见, 我们的计算是合理的。
3 结论
根据分子动力学理论建立了非晶Ni64 Pd16 Al20 的结构模型, 用Recursion方法计算了Ni64 Pd16 Al20 系非晶及相应晶态的电子结构。 结果表明: Ni64 Pd16 Al20 合金的非晶形成能力可以由局域电子态密度(LDOS)、 Al与过渡元素间的电子转移和非晶与晶态合金的总结构能差来研究。 LDOS低能区双(或三)峰结构向高能区偏移多的过渡元素与Al形成的共价键较强, 离子键与共价键共存使Ni64 X16 Al20 合金非晶形成能力增强。 Ni64 Pd16 Al20 非晶与相应晶态的总结构能差越小的过渡元素使其形成非晶的驱动力越强, 对非晶形成越有利。 过渡元素加入到NiAl合金中, 非晶形成能力由强到弱的顺序为:Ir, Pt, Rh, Pd, Au, Ag, Cu。
参考文献
[1] 郭贻诚,王震西非晶态物理学[M]北京:科学出版社,1984.1719. GUODai cheng,WANGZhen xi.AmorphousPhysics[M ].Beijing:SciencePress,1984.1719.
[2] InoueA ,ZhangT .ImpactfractureenergyofbulkamorphousZr55Al10Cu30Ni5alloy[J].MaterTransJIM ,1996,37:17261402.
[3] GrahamCD ,EgamiT .Magneticpropertiesofamor phousalloys[J].JMagMagMater,1980(1518):13251330.
[4] NishiyamaN ,InoueA .FluxtreatedPdCuNiPamor phousalloyhavinglowcriticalcoolingrate[J].MaterTransJIM ,1997,38:464472.
[5] InoueA ,ShihataT ,ZhangT .Effectofadditionalele mentonglasstransitionbehaviorandglassformationten dencyofZrAlCuNialloys[J].MaterTransJIM ,1995,36:14201426.
[6] PekerA ,JohnsonWL .Ahighlyprocessablemetallicglass:Zr41.2Ti13.8Cu12.5Ni10.0Be22.5[J].ApplPhyLett,1993,63(17):23422347.
[7] 张庆生,吴世丁,张海峰,等淬态晶相对块状ZrAlNiCu非晶合金疲劳断裂的影响[J].金属学报,2002,38(3):265268. ZHANGQing sheng,WUShi ding,ZHANGHai feng,etal.EffectsofquenchedcrystallitesoncyclicfatiguefractureofZrAlNiCubulkamorphousalloy[J].ActaMetalSin,2002,38(3):265268.
[8] ZHANGGuo ying,LIUGui li,ZENGMei guang.Re cursivedeterminationoftheembrittlingeffectofalloyingelementontheimpurityinducedlowanglegrainboundaryembrittlement[J].JournalofIronandSteelResearch,2001,8(2):5155.
[9] 张国英,刘贵立,曾梅光,等钢中小角晶界区的电子结构及掺杂效应[J]物理学报,2000,49(7):13441347 ZHANGGuo ying,LIUGui li,ZENGMei guang,etal.Energiesofimpuritiesanddopingeffectat5.3°lowanglegrainboundaryinsteel[J].ActaPhysSin,2000,49(7):13441347.
[10] 韦文楼,郭 进,邓 文,等.LaNi5基合金电极性能与其电子结构的相关性[J]中国有色金属学报,2002,12(3):501504. WEIWen luo,GUOJin,DENGWen,etal.Correla tionbetweenelectronicstructureofLaNi5basealloysandhydrogenabsorptionproperties[J].TheChineseJournalofNonferrousMetals,2000,12(3):501504.
[11] HaydockR .SolidStatePhysics,Vol.35[M].NewYork:AcademicPress,1980.216.
[12] HarrisonWA .ElectronicStructureandthePropertiesofSolids[M].SanFrancisco:Freeman,1980.551.
[13] SuttonAP ,ChenJ.Long rangefinnis sinclairpoten tials[J].PhilosophicalMagazineLetter,1990,61(3):139146.
[14] 肖慎修,王崇愚,陈天朗密度泛函理论的离散变分方法在化学和材料物理学中的应用[M ]北京:科学出版社,1998.92. XIAOShen xiu,WANGChong yu,CHENTian lang.TheApplicationofTheDiscreteVariationalMethodinTheDensityFunctionalTheorytoChemistryandMate rialsPhysics[M].Beijing:SciencePress,1998.92.
[15] 魏全金材料电子显微分析[M ]北京:冶金工业出版社,1990.186187. WEIQuan jin.TheElectronMiroAnalysisinMaterials[M ].Beijing:MetallurgicalIndustryPress,1990.186187.