非晶合金结构演变影响玻璃形成能力的分子动力学研究.pdf
《非晶合金结构演变影响玻璃形成能力的分子动力学研究.pdf》由会员分享,可在线阅读,更多相关《非晶合金结构演变影响玻璃形成能力的分子动力学研究.pdf(9页珍藏版)》请在咨信网上搜索。
1、第41卷第1期2024年2 月J.At.Mol.Phys.,2024,41:016007(9pp)原子与分子物理学报JOURNAL OF ATOMIC AND MOLECULAR PHYSICS非晶合金结构演变影响玻璃形成能力的分子动力学研究Vol.41 No.1Feb.2024史天姿,龙志林,张新瑞(湘潭大学土木工程学院,湘潭41110 5)摘要:通过分子动力学方法模拟了Tizs,Al2s、Ni s o Zr s o 和CusoZrso非晶合金的玻璃转变过程,得到并分析了平均原子体积和双体分布函数等结构参数,并应用Voronoi多面体指数分析法统计了玻璃转变过程中二十面体及类二十面体团的数量
2、,通过分析团簇在玻璃转变过程中种类和数量的涨落趋势,研究了非晶原子由短程序连接至中程序再至长程无序的动力学演化过程,结果表明,非晶合金的玻璃形成能力以及塑性变形能力与动力学演化过程中Voronoi团簇的种类和数量密切相关在玻璃转变过程中局部五次对称性高的团簇倾向于连接在一起形成链,从而密排整个空间,降低系统的动力学行为从而提高玻璃形成能力塑性形变倾向于发生在局部五次对称性较低的区域在玻璃转变温度附近团簇种类和数量的突变反映出非晶合金的自组织临界行为,蕴含着丰富的非线性动力学现象。关键词:非晶合金;原子团簇;玻璃形成能力;分子动力学中图分类号:0 46 9The influence of str
3、ucture evolution of metallic glass onglass forming ability by molecular dynamics study文献标识码:AD01:10.19855/j.1000-0364.2024.016007SHI Tian-Zi,LONG Zhi-Lin,ZHANG Xin-Rui(Department of Physics,Xiangtan University,Xiangtan 411105,China)Abstract:The glass transformation processes of Tizs Al2s,NisoZrso an
4、d Cuso Zrso were simulated by molecular dy-namics method.The average atomic volume and radial pair distribution function were obtained.Besides,thenumber of icosahedral and icosahedral-like clusters were calculated by Voronoi polyhedron index analysis.E-ventually,the dynamic evolution of amorphous at
5、oms from short order to medium order and then to long rangedisorder was studied by analyzing the type and number of clusters during glass transformation.The results showthat the glass forming ability is closely related to the type and number of Voronoi clusters.Besides,the clusterswith high local qu
6、intuple symmetry tend to join together in order to form chains during the glass transformationprocess,which densely fills the whole space and reduces the dynamic behavior of the system.Furthermore,theplastic deformation tends to occur at the area with low local quintuple symmetry.It is intriguing to
7、 note that theabrupt changes of the types and number of clusters near the glass transformation temperature reflect the self-or-ganized criticality behavior of metallic glass.Key words:Metallic glass;Atomic clusters;Glass forming ability;Molecular dynamics收稿日期:2 0 2 2-0 8-30基金项目:国家自然科学基金(5197 118 8)作
8、者简介:史天姿(1998 一),女,陕西省咸阳市人,硕士生,主要研究领域为非晶合金力学行为,通讯作者:龙志林.E-mail:l o n g z l x t u.e d u.c n016007-1第41卷1引 言非晶合金又称金属玻璃,由于没有晶体合金中诸如位错等结构缺陷从而表现出超高强度、韧性、弹性极限等优异的力学性能因此受到广泛关注。然而,玻璃转变的本质(包括玻璃结构和玻璃转变行为)是长期以来困扰凝聚态物理和材料科学领域的一个重大难题,非晶合金作为从原子尺度研究该类问题的极佳物理模型,肩上背负着沉重的使命2 宏观无序的非晶结构中存在有序特征,非晶合金中已被证明有短程序(SRO)和中程序(MRO
9、)且与非晶的诸多性质如玻璃形成能力(CFA)有重要联系3-8 然而,由于非晶制备的困难以及实验观测手段的限制,非晶在玻璃转变过程中的微观结构演化过程至今没有一个清晰的图像,分子动力学(MD)模拟是弥补这些不足的重要手段.采用MD方法分析液态金属和金属玻璃的微观结构和玻璃转变行为迄今已取得诸多重要进展9-15 如Park等 在非晶体系微观结构与宏观性能关系的研究中发现,原子堆积密度与非晶固体的弹性模量、强度和塑性密切相关,由Voronoi多面体定义的各种原子团簇构成了非晶的SRO,并对非晶固体的力学性能有决定性影响Cheng等10 利用CuZr体系证明了二十面体团簇是金属过冷液体非Arrheni
10、us动力学降速的微观起源,从结构动力学角度解释了非晶依赖成分的力学性质以及其强度与玻璃转变温度之间的关系:Sha等11 在模拟CuZr体系的玻璃转变过程中发现大量以Cu原子为中心的类二十面体团簇可使体系的堆积方式更为致密,从而提高体系的CFA.Li等12 在MRO尺度通过XRa y 衍射实验揭示了CuZr非晶以二十面体团簇构成弦状主干并以类液结构填充剩余空间的堆积方式并与经典MD 模拟结果一致Wei等13 在CuZrAl三元体系的MD模拟中发现Al元素的含量对体系团簇种类和数量的分布有显著影响,从而在该体系GFA中扮演关键角色.然而,上述研究主要关注的二十面体团簇随Cu含量变化较大,在其他不含
11、Cu元素的体系中几乎监测不到,表明该参数用于描述非晶结构特征并不具有广泛代表性,需要引进新的结构参数.且上述研究大都致力于多元非晶体系的模拟,但合金中含有多个组元使得对非晶GFA、玻璃转变等基础问题的研究和计算机模拟分析变得十分复原子与分子物理学报杂,因此,开发和研究二元块体非晶合金具有重要意义:作为典型的二元非晶态材料,NiZr 系可在相当宽的成份范围(2 0-7 0 at%Zr)内形成非晶态从而引起人们很大兴趣:最近的研究发现,Ni-Zr非晶合金的屈服强度随材料尺寸的降低而增加,当薄膜厚度低于36 0 nm时,拉伸过程中不出现剪切带,具有非常大的塑性,甚至能达到15%的断裂应变15,是研究
12、尺寸效应和非晶增韧的绝佳模型TiA l合金具有低密度、高比强度、抗高温、抗氧化、耐腐蚀等特殊性能,是航空航天和高温应用领域具有吸引力的材料,且该合金具有很强的结晶性,是研究非晶晶化行为的良好模型16,17 ,本文利用MD手段,模拟了 TizsAl2s、Ni s o Zr s o和CusoZrso三种典型并极具对比性的二元非晶合金模型的玻璃转变过程,并引进更具广泛性的新结构参数局部五次对称性(LFFS)研究了非晶由SRO连接至MRO再至长程无序的动力学演化过程,并由此解释了三种非晶合金宏观力学性质迥异的原子结构层面原因。2模拟计算的条件和方法利用开源的大规模原子分子并行程序LAMMPS模拟了Ti
13、zsAl2s和NisoZrso非晶合金的玻璃转变过程表1列举了模拟过程中两种非晶的不同参数设置,考虑到原子间的多体相互作用,TirsAl2s合金原子间相互作用势选用Mishin18拟合第一性原理数据得到的金属合金镶嵌原子势(EA M),Ni s o Zr s o 合金原子间相互作用势选用Sheng等19 拟合得到的EAM势,如图1所示,TirsAl2s(Ni s o Zr s o)模拟在含有40 0 0(32 0 0 0)个原子505050A(7 2 7 2 7 2 A 3)的模拟盒子中进行,其中Ti(Ni)和Al(Zr)的原子数分别为30 0 0(16 0 0 0)和10 0 0(16 0
14、0 0)个冷却过程中采用等压等温系综并利用Nose-Hoover算法控制温度和压强,边界条件采用立方周期性边界条件,时间步长设定为2 fs.首先,根据最小势能原理,系统处于稳态时的能量最小,如图2 可见,让系统在2 0 0 0 K温度下弛豫10 0 ps(2 0 p s)后体系原子处于稳定状态,即得到热力学平衡的非晶合金熔融液态结构,可以开始模拟玻璃转变随后,分别以Q1=0.51013 Ks-1,Q2=0.5 1012 Ks-(Q;=1.0 1013Ks-1,Q4=1.0101Ks)的冷却速率快速降温016007-2第1期第41卷至30 0 K,并保温2 0 ps得到系统的微观结构,每隔10
15、0 0 步记录一次系统的位置信息最后,运用Table 1 Molecular dynamics simulation parameters of two metallic glasses.研究对象EAM势Ti7s Al2sMishin YNisoZrsoSheng H W史天姿,等:非晶合金结构演变影响玻璃形成能力的分子动力学研究表1两种非晶合金的分子动力学模拟参数设置原子总数模拟盒子NA340005050503200072 72 72第1期微观结构可视化分析程序Ovito分析其微观结构特征。驰豫时间ps10020冷却速率Ks-1Q1=0.5 10 13Q2=0.5 1012Q;=1.0 10
16、13Q4=1.0 10124.05(a)4.10F-4.15(a)图1分子动力学结构模型示意图(a)T i r s A l 2 5,红色原子为Ti,蓝色原子为Al;(b)Ni s o Zr s o,紫色原子为Ni,绿色原子为ZrFig.1TThe schematics of molecular dynamics structuremodels:(a)T i r s A l 2 s,t h e r e d a t o m s a r e T i,b l u eatoms are Al;(b)NisoZrso,the purple atoms areNi,green atoms are Zr.3模
17、拟计算结果与讨论在金属液体冷却过程中,通常采用平均原子体积(体系总体积与所含原子数N的比)描述体系密度的变化情况图3显示出TirsAls和NisoZrso合金在不同冷却速率下平均原子体积随温度的变化情况由图3(a)可以发现,TirsAl2s合金在Q.冷却速率下,平均原子体积随着温度的降低缓慢且连续地减少;而在Q2冷却速率下,体系平均原子体积在10 0 0 K左右发生突变,表明此时TizsAl2s合金结构发生重大变化,即由无序的液态结构形成有序的晶态结构根据Abraham对形成玻璃系统的模拟研究表明:冷却过程中,如果平均原子体积与温度曲线的斜率发生突变,说明形成晶体结构结合图3(a)可以确定Ti
18、sAl2s合金发生玻璃转变的温度大约为10 0 0 K.同理,由图3(b)可知,NisoZrso合金在Q,和Q4两种冷却速率下均形成了非晶结构Q1与Q2、Q,与Q4均只相差一个数量级,但两种合金却表现出完全不同的玻璃转变行为,这表明NisoZrso合金的GFA胜于TizsAl2s合金至少一个量级,TirsAls合金十分容易发生晶化,目前的实验手段很难获得稳定的TirsAls非晶合金:(b)-4.204.2504.6-4.8-5.0-5.2-5.4-5.6图2平均原子势能随驰豫时间的变化图(a)T i r sAl2s;(b)NisoZrsoFig.2 Plots of average atomi
19、c potential energy with re-laxation time:(a)T i z s,A l 2 s;(b)Ni s o Zr s o-双体分布函数(RDF)g(r)描述的是在原点处存在粒子的情况下在距离处找到粒子的概率,它本质上是粒子间距离的直方图,通过粒子的数密度(即粒子总数除以模拟盒子体积)进行归一化,是MD模拟中表示微观结构的一个重要方法.和原子的偏双体分布函数定义如下NNaNego(r)4mrpN.Nei=1i=1其中p表示N个原子体系的原子数密度;N。和N016007-320Relaxation Time(ps)05Relaxation Time(ps)40106
20、01580100(b)20(1)Calculated第41卷20.0(a)19.519.018.518.0F17.517.0018.5(b)18.0F17.517.0图3平均原子体积在不同冷却速率下随温度的变化关系图(a)T i z s A l 2 s;(b)Ni s o Zr s oFig,3 Plots of average atomic volume with temperatureat different cooling rates:(a)Tizs Al2s;(b)NisoZtso分别表示类和类原子的数目;(类)和原子j(类)的距离图4为本文模拟的TizsAl2s合金在2 0 0 0
21、K下的双体分布函数与参考文献值的对比。从图中可以看出,在2 0 0 0 K下,双体分布函数的第一峰位和峰高与Pei 等1计算结果符合良好,表明本文的模拟过程具有足够精度由于非晶的短程序,即有确定的最近邻及次近邻配位层,其RDF曲线有清晰的第一峰及第二峰图5(a)为TirsAl2s合金在Q,冷却速率下2000400K温度范围每隔2 0 0 K的双体分布函数演化曲线集可以看到,随着温度的降低,g(r)曲线的第一峰(r3A 处)有右移趋势,宽度随温度降低而变窄且峰值逐渐增大,这是由于冷却过程中原子不断被冻结,临近中心原子第一峰位处的原子数逐渐增多所导致的,表明中心原子的配位数增加,短程有序性增强当温
22、度高于10 0 0 K时,g(r)曲线变化不明显,表明体系进入熔融态后内能和结构均变化不大值得注意的是,当温原子与分子物理学报度低于10 0 0 K后,第二峰位置(r5A 处)几乎没有变化,但却发生明显的劈裂现象,表明非晶结构形成,这与根据平均原子体积随温度变化情况(如图3)分析玻璃转变温度的结果一致图5(b)为NisoZrso合金在Qs冷却速率下2 0 0 0 40 0 K温度Q范围每隔2 0 0 K的双体分布函数演化曲线集:可以看出,随温度的降低g(r)曲线的第一峰同样出Q2现右移,且峰宽变窄峰高增加,当r6A后g(r)曲线的上下波动趋于平缓,类似于空间分布完全5001000Tempera
23、ture(K)QQ4400800Temperature(K)第1期1500200012001600表示原子i金样品的偏双体分布函数无论TizsAl2s(图6 a)还是NisoZrso(图6 b)的偏双体分布函数均展现出了第二峰劈裂的特性,表明MD模拟所得样品均处于非晶态TizsAl2s非晶的三个偏双体分布函数中,Ti-Al和 Ti-Ti 曲线基本重合且第一峰强度最大,表明Ti原子吸引Ti原子和Al原子的能力均较强并倾向于排列在第一近邻的位置上;铝原子之间吸引力较弱,也排在第一近邻位置上然而,NisoZrso非晶的三个偏双体分布函数的第一峰016007-4无序的理想气体的双体分布函数图,表明Ni
24、soZrso非晶在中程尺度(5Ar20A)的无序性更强.不同的是,当温度低于8 0 0 K后,NisoZrso合金g(r)曲线的第一峰开始出现明显的劈裂现象,表明NisoZrso非晶在短程尺度(r5A)的有序性更强当温度降至40 0 K时第二峰开始发生劈裂劈裂峰数量的增加以及玻璃转变温度的降低均表明NisoZrso合金的非晶形成能力胜于TizsAl2s合金,这也与平均原子体积随温度变化情况的分析结果一致3.020002.52.0F31.5601.00.5F0.00图4本文计算得到的TirsAl2s合金在2 0 0 0 K时的RDF与参考文献17 的对比曲线Fig.4 The calculate
25、d RDF of Tizs Al2s at temperature2000 K with the available data.图6 给出了30 0 K下MD模拟得到的非晶合OReferenceod8024r(A)68第41卷(a)400K60600K800K1000K1200K1400K1600K1800K2000K02(b)400K600K800K1000K1200K1400K1600K1800K2000K02图52000400K温度范围每隔2 0 0 K的双体分布函数曲线集(a)Q i 冷却速率下的TirsAl2s合金;(b)Q 冷却速率下的NisoZrso合金Fig.5The radi
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 合金 结构 演变 影响 玻璃 形成 能力 分子 动力学 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。