Au团簇在金表面沉积过程的分子动力学模拟.pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号:061120 Au团簇在金表面沉积过程的分子动力学 模拟· 高虹，赵良举，曾丹答，高丽娟 重庆大学动力工程学院，重庆400044 联系电话:023E-mail:gaohong1021@ 摘要 在冷喷涂过程中，喷涂拉子被超音速气流加速到较高的速度，在低于喷涂材料熔点的温度下撞 击基体，发生剧烈的塑性变形而沉积形成涂层。但是由于高速粒子碰掩变形的瞬时特点，不能对拉子 变形沉积过程进行直接观察，本文通过对灿 团簇在彻基体上的沉积过程的分子动力学模拟，可以观 察到团簇撞击基体并在基体上沉积的过程，以及团簇和基体的形貌变化:本文还探讨了荃体的温度变 化，发现在披击过程中，基体的局部区域有熔化现象:通过计算团簇原子进入墓体表面层的数f及团 簇与基体间的最终接触面积，探讨了影响团簇沉积过程的主要因素。 关键词 冷喷涂，分子动力学棋拟，团簇，沉积过程 1.引言 冷气动力喷涂技术(CGDS:ColdGasDynamicSpray，简称冷喷涂CS:ColdSpray) 是近年来发展起来的一项新兴的材料改性技术。在冷喷涂过程中，喷涂粒子被超音速气 流加速到较高的速度 (300-1200mls)，在远低于喷涂材料熔点的温度下撞击基体，发生 剧烈的塑性变形而沉积形成涂层。与目前应用广泛的热喷涂技术相比，基体及涂层材料 没有受到高温熔化、气化和氧化，因此可以克服热喷涂工艺中因高温熔化所造成的诸多 缺陷，特别适合非耐高温工件的材料表面改性1[11M. 由于高速粒子碰撞变形的瞬时特点，不能对粒子变形沉积过程进行直接观察，所以 现有的研究冷喷涂方法多为宏观状态的数值模拟方法。李文亚等阔采用有限元数值计算 方法研究了冷喷涂过程中to粒子与Cu基体的碰撞变形行为，并探讨了粒子速度、温 度等对沉积过程的影响;Dykhuizen等4(1采用数值模拟方法在研究。粒子撞击不锈钢基 体表面时发现粒子与基体两者都产生了巨大的形变，并研究了不同撞击速度对基体表面 的影响。 随着近年来计算机技术的飞速发展，计算材料科学己成为材料学家、凝聚态物理学 家和化学家研究的热点，其中蒙特卡洛方法和分子动力学方法等已成为人们利用计算机 研究金属特性的有力工具。分子动力学 (MD)模拟是对于原子核和电子所构成的多体 系统，用计算机模拟原子核的运动过程，从而计算系统的结构和性质的方法。国内外的 学者们己将MD广泛应用于电离的原子簇束沉积((ICBD)过程[51-[41，用于了解ICED的 微观机制、所生长的薄膜质量与镀膜环境 (如轰击原子簇能量、基体温度等)的关系。 但是，ICBD镀膜所使用的粒子的能量一般为几eV到几十eV，相对于冷喷涂的入射能 量要大得多，所以IC31D过程的研究结论不能适用于冷喷涂。国内外将D 应用于冷喷 .基金资助:国家自然科学基金青年基金资助项目(编号 涂沉积过程的文献还非常少见，本文将分子动力学模拟方法应用于研究冷喷涂中沉积的 微观过程，探讨了团簇与基体之间的相互作用，以及撞击速度、团簇大小对沉积过程的 影响。 么计算模型 本文的计算模型如图1所示。 .目口.盯口囚称门. 图1计算棋型海体和团簇都已经过sops的等温弛邃) 模拟中将基体表面取为 ((301)面，坐标原点选在基体的中心，基体体积为 16aox16aoxSao(ao=0.40704=，为Au的晶格常数)，包含8192个Au原子.基体为 FCC结构，共16层，每层包含512个原子。与 (001)面平行的XY方向采用周期性边 界条件，用于模拟腹无限靶。为研究碰撞中基体的温度变化，在模拟中表面8层原子采 用NV 系综，下面8层原子采用Berendsen热浴法进行温度修正，模拟NVT系综。基 体温度为300K. 本文的模拟中有3种团簇，Au152,Au276,Au420,Au604，分别包含152,276, 420,604个Au原子。初始时刻按FCC分布。团簇的入