第26卷第9期强激光与粒子束V o l.26,N o.9 2014年9月H I G H P OW E R L A S E R A N D P A R T I C L E B E AM S S e p.,2014飞秒激光烧蚀镍钛形状记忆合金的蚀除机理*唐一波1,陈冰2,陈志勇1,朱卫华1,李月华1,王新林1,2(1.南华大学电气工程学院,湖南衡阳421001;2.南华大学机械工程学院,湖南衡阳421001)摘要:结合双温模型的分子动力学模拟方法,研究了飞秒激光脉冲辐照B2结构镍钛合金时烧蚀阈值附近的靶材蚀除机制,数值模拟了中心波长为800n m,脉宽为100f s,能量密度为25~50m J/c m2的激光与90n m厚B2结构镍钛合金薄膜相互作用过程㊂确定了脉宽为100f s的脉冲激光与镍钛形状记忆合金相互作用的烧蚀阈值,发现烧蚀阈值条件下,靶材的蚀除机制是单纯基于应力作用的机械破碎;烧蚀阈值附近,未蚀除靶材受热影响发生无序化相变的区域较小,且随激光能量密度的降低而减小㊂提高激光功率密度,烧蚀同时呈现热机械蚀除和机械破碎机制㊂关键词:飞秒激光烧蚀; B2结构镍钛合金;双温模型;分子动力学模型中图分类号:0437; T N249文献标志码: A d o i:10.11884/H P L P B201426.091025由功能材料制备的微纳器件及结构表现出来的优越性能及其诱人的应用前景,使得功能材料微纳器件制备成为近年来超快激光微细加工领域的研究热点[1-4]㊂镍钛合金因在不同温度下可实现不同微观相结构间的可逆性转换而具备形状记忆功能,被称为形状记忆功能材料,在国防㊁军事㊁航天航空㊁生物医学以及工业等领域有着广泛的应用,镍钛形状记忆合金制备的微纳器件在医学应用等领域更是起着无可替代的作用[5-6]㊂而镍钛形状记忆合金的加工具有其特殊性,传统的加工方法产生的热效应或应力易引起加工区域产生相变,从而会改变镍钛合金的特性,影响加工零件的性能㊂加工过程带来的性能改变在微细加工中的影响尤为重要㊂因此探索新的更实用的镍钛合金微细加工方法具有重要意义㊂飞秒激光加工的超高精度㊁超小热影响区域及加工材料范围广等特性,使得采用飞秒激光进行无相变㊁微纳加工具有不可替代的优势[7-8]㊂近年来,激光技术不断的取得进展,使得飞秒激光加工有望成为镍钛形状记忆合金微纳加工的有效且先进的手段[9-12]㊂开展飞秒激光烧蚀镍钛形状记忆合金的数值模拟工作,探索不同激光参数条件下靶材的相变行为以及蚀除机制,进而为飞秒激光微纳㊁无相变加工提供理论基础和参数依据,对拓展镍钛合金的应用具有重大意义㊂L e o n i dV.Z h i g i l e i 等开展的分子动力学模拟工作总结得出,超快激光烧蚀金属时,热效应及烧蚀压力波共同作用导致靶材发生了蚀除,并提出机械破碎㊁裂散㊁液相爆炸等一系列靶材蚀除机制[13-14]㊂然而,更为详细的靶材蚀除机理以及微观蚀除现象还需进一步探讨㊂本文数值模拟并分析了飞秒激光与B2结构镍钛形状记忆合金靶材相互作用时,烧蚀阈值附近的靶材蚀除机理,给出了飞秒激光微纳㊁无相变加工镍钛形状记忆合金的参数区间㊂1数值模拟方法超快激光与金属靶材相互作用时,激光能量首先沉积到靶材电子,靶材电子被迅速加热至极高的温度,而晶格却仍处于相对 冷 的状态,导致了靶材电子与晶格之间的非热平衡㊂双温模型将靶材分为电子和晶格两个体系,分别计算电子㊁晶格体系的热传导,突破了传统热传导方程的局限,较为准确地描述了超快激光辐照下,金属靶材的非热平衡能量弛豫过程㊂而分子动力学通过求解靶材原子体系的牛顿力学方程组,追踪了靶材每一个原子的运动,详细地描述了靶材去除㊁相变等微观过程㊂结合双温模型的分子动力学方法,具备了双温模型和分子动力学的所有优点,是研究超快激光与金属材料相互作用最为常用的方法㊂描述双温模型的双温方程为C e∂T e∂t=∂∂z(k e∂T e∂z)-g(T e-T l)+S(z,t)(1)C l∂T l∂t=-∂∂z(k l∂T l∂z)+g(T e-T l)S(z,t)=I(t)(1-R)αe x p(-αz)*收稿日期:2013-12-10;修订日期:2014-04-17基金项目:国家自然科学基金项目(11174119);南华大学重点学科建设资助项目(N H X K04)作者简介:唐一波(1988 ),男,硕士研究生,主要从事超快激光与金属相互作用方面的研究;t y b19880810@s i n a.c o m.c n㊂通信作者:王新林(1970 ),男,博士㊁教授,主要从事激光与光电子技术及应用方面的研究;w x l_l y000@y a h o o.c o m.c n㊂091025-1091025-2 电子温度远小于费米温度时,C e (T e )=C e 0T e 为电子热容,其中C e 0是电子热容常量;C l 是晶格热容;k e ,k 1分别为电子㊁晶体热导率;T e ,T l 分别为电子㊁晶格温度;g 为电子晶格耦合系数;S (z ,t )表示激光光源项;R 为材料对激光的表面反射率;α为材料的吸收系数;z 为激光传输方向;I (t)表示激光光强随时间的分布㊂将飞秒脉冲激光沉积到靶材的能量以速度均衡力的形式耦合到靶材晶格体系,从而实现双温模型与分子动力学的结合㊂结合双温模型的分子动力学模拟中靶材原子体系运动方程组为m i ∂2r i ∂t 2=F i +u m i v i (2)F i =-Ñr i E t o t a l (3)式中:F i 为第i 个原子无外加载情况下靶材场受力;m i ,r i ,v i 分别为第i 个原子的质量㊁位置和速度;E t o t a l 为势能;u 为速度均衡因子,其大小主要与飞秒激光作用下金属内部电子晶格温度之差有关,均衡因子为u =ðnk =1g V N (T e -T l [])/n ˑði m i (v i )[]2(4)式中:V N 为模拟区域第N 层的体积㊂采用第二动量紧束缚势描述[15-16],势函数截断半径r c =0.42n m ㊂E t o t a l =ði ðj ʂi A αβe x p -p αβr i j d αβ-æèçöø÷éëêêùûúú1-ðj ʂi ξ2αβe x p -2q αβr i j d αβ-æèçöø÷éëêêùûúú{}1(5)式中:α和β表示原子种类;r i j 表示原子i 和j 之间的距离;参量A αβ,ξαβ,p αβ,q αβ的大小由原子间结合能㊁弹性常量㊁空穴形成能和晶格常量决定;d αβ表示原子之间的最邻近距离㊂计算B 2结构镍钛形状记忆合金热物性参数[17],双温方程采用对三角矩阵追赶法进行差分化求解,并根据冯㊃诺依曼稳定性判据原理:t F D <0.5(Δz F D )2C e (T e )/k e (T e ,T l ),取时间差分步长t F D 为0.1f s ,空间差分步Δz F D 为1n m ㊂靶材原子体系运动方程组采用V e r l e t 积分算法进行计算㊂模拟对象为体心立方结构的B 2镍钛形状记忆合金材料,靶材体系由5ˑ5ˑ90个晶胞组成,横向尺寸1.5n mˑ1.5n m ,纵向深度90n m ,总原子数为15000㊂取x O y 坐标面为横向界面,z 为纵向深度方向,激光脉冲沿z 轴入射㊂x ,y 方向施加周期性边界条件,底部施加压力传输边界条件[18]㊂采用热力学统计方法进行温度计算,基于维里理论进行压力统计[19]㊂T l =ði m i (v i )2/(3k B N )(6)p =ρk B T l +16V <ðNi =1ðj ʂi F i j ㊃r i j >(7)式中:ρ为原子密度;k B 为波耳兹曼常数;V 为体积;F i j 为第i 个原子与第j 个原子之间的相互作用力;模拟在300K 的室温下进行㊂2 结果及讨论数值模拟时,选用典型的商用化钛宝石飞秒激光器参数,脉宽为100f s ,中心波长为800n m ㊂为探寻飞秒激光与B 2镍钛合金相互作用的烧蚀阈值,在参考飞秒激光与单质金属相互作用的烧蚀机理后,选取能量密度区间为25~50m J /c m 2,以5m J /c m 2为间隔,并在阈值附近密集采样㊂分别计算了脉宽为100f s,各能量密度的激光与90n m 厚B 2结构镍钛合金薄膜相互作用的温度场㊁压力场时空分布以及原子位型㊂激光在1p s 时刻加载,整个模拟时间持续30p s ㊂2.1 烧蚀阈值烧蚀阈值是靶材开始发生蚀除的临界值,是飞秒激光加工中参数选择与工艺优化的重要判据㊂分子动力学模拟时,追踪每一个原子的运动,靶材的原子位型图包含了每一个原子的位置信息,有效地呈现了飞秒激光与B 2结构镍钛合金靶材相互作用时靶材 相变-断裂-去除 的过程㊂数值模拟时,只有60n m 以上区域发生了相变,为了清晰呈现靶材的相变过程,原子位型图只摘取60n m 以上的区域原子位置信息㊂脉宽为100f s ,能量密度分别为40,40.5和41m J /c m 2时,激光脉冲与B 2结构镍钛合金靶材相互作用时靶材辐照中心处的原子位型(蓝色点和红色点分别代表N i 原子和T i 原子)如图1(a ),(b ),(c )所示㊂如图1(a )所示,脉宽为100f s ,能量密度为40m J /c m 2的激光烧蚀时,靶材表层出现无序化,但直到30p s也未出现明显的原子分层与晶格断裂,表明无烧蚀现象发生㊂图1(b ),(c )先后在19p s 左右出现原子分层,在30p s 时分层明显,发生断裂㊂表明在该能量密度下发生烧蚀现象,但此条件下,靶材底层原子排列依然比较规强激光与粒子束091025-3F i g .1 S n a p s h o t s f r o mt h em o l e c u l a r d y n a m i c (M D )s i m u l a t i o no f l a s e r i r r a d i a t i o no naB 2N i T i a l l o y t a r ge t w i t h t h e p u l s e d u r a t i o nof 100f s a n dd i f f e r e n t f l u e n c e图1 脉宽为100f s ,不同能量密度激光脉冲辐照靶材分子动力学模拟原子位型图则㊂由图1(a ),(b ),(c )可得,脉宽为100f s 的激光与B 2结构镍钛合金靶材相互作用,能量密度约为40.5m J/c m 2时,开始发生烧蚀㊂即脉宽为100f s 的激光与B 2镍钛合金靶材相互作用的烧蚀阈值约为40.5m J /c m 2㊂2.2 烧蚀阈值附近靶材蚀除机理烧蚀阈值附近的靶材蚀除机理,可以揭示靶材开始发生蚀除的物理本质,为提高飞秒激光加工的质量及优化加工参数提供了理论指导㊂脉宽为100f s 时,飞秒激光与B 2结构镍钛合金靶材相互作用的烧蚀阈值为40.5m J /c m 2,其靶材内部温度场㊁压力场时空分布分别如图2(a ),(b )所示㊂F i g .2 T h e t e m p e r a t u r e a n d p r e s s u r e f i e l dd i s t r i b u t i o n i n t i m e a n d s pa c e f r o mt h eM Ds i m u l a t i o no f l a s e r i r r a d i a t i o no n aB 2N i T i a l l o y fi l m w i t h p u l s e d u r a t i o no f 100f s a n d f l u e n c e o f 40.5m J /c m 2图2 脉宽为100f s ,能量密度为40.5m J /c m 2的激光辐照B 2结构镍钛合金靶材分子动力学模拟温度场㊁压力场时空分布图2(a )中,红色表示熔点以上温度(1583K 以上),黄色表示热影响温度(1200~1583K ),由图2(a )可知,能量密度为40.5m J /c m 2的激光作用后,25n m 以上靶材温度达到熔点,发生熔化㊂25~50n m 的靶材温度达到热影响临界温度,晶格发生无序化相变㊂图2(b )中,黄色代表数值为正的压力,蓝色代表数值为负的拉应力,红色代表发生压力集中,深蓝色代表发生拉应力集中㊂由图2(b )可知,烧蚀时,在靶材的次表层形成压力集中,在25n m 以上相变较大的区域形成数值为负的应力场,在25n m 以下相变较小的区域形成拉应力集中㊂结合图2(a ),(b )与图1(b )可知,靶材在拉应力波叠加作用下,在温度熔点以下的35n m 处发生了蚀除㊂图2(a ),(b )中的靶材蚀除过程在图1(b )中得到了很好的呈现㊂图1(b )所示,19p s 时刻,35n m 处无序化的晶格在拉应力叠加作用下发生拉裂,产生裂隙㊂接下来裂隙持续生长,并于26p s 时刻,将靶材拉断,被拉断的靶材以约500m /s 的速度向外喷射,最终发生蚀除㊂由此可得,烧蚀阈值时,激光与B 2结构镍钛合金靶材相互作用产生的拉应力波叠加作用于受热影响发生无序化相变的晶格,导致靶材断裂,发生蚀除,蚀除机制是单纯基于应力作用的机械破碎[13]㊂烧蚀阈值附近,脉宽为100f s ,能量密度41m J /c m 2的激光与B 2结构镍钛合金靶材相互作用时,靶材内部温度场㊁压力场时空分布分别如图3(a ),(b )所示㊂由图3(a ),(b )可知,脉宽为100f s ,能量密度为41m J /c m 2的激光与B 2结构镍钛合金靶材相互作用,靶材内部温度场㊁压力场时空分布以及在热传导作用下发生熔化和无序化相变的区域与烧蚀阈值时基本相同㊂然而,41m J /c m 2的激光烧蚀时,靶材在拉应力的作用下,在温度达到熔点的24n m 处发生了蚀除㊂图3(a ),(b )中的热熔化与拉应力共同作用导致靶材发生蚀除的过程在图1(c )中得到了很好的呈现㊂图1(c )所示,18p s 唐一波等:飞秒激光烧蚀镍钛形状记忆合金的蚀除机理091025-4F i g .3 T h e t e m p e r a t u r e a n d p r e s s u r e f i e l dd i s t r i b u t i o n i n t i m e a n d s pa c e f r o mt h eM Ds i m u l a t i o no f l a s e r i r r a d i a t i o no n aB 2N i T i a l l o y fi l m w i t h t h e p u l s e d u r a t i o no f 100f s a n d f l u e n c e o f 41m J /c m 2图3 脉宽为100f s ,能量密度为41m J /c m 2的激光辐照B 2结构镍钛合金靶材分子动力学模拟温度场㊁压力场时空分布时刻,24n m 处熔融的靶材,在拉应力的作用下产生空穴㊂随即空穴持续生长,并于30p s 时刻,将靶材完全蚀断,发生蚀除㊂由此可得,脉宽为100f s ,能量密度为41m J /c m 2的激光与B 2结构镍钛合金相互作用时,热熔化与拉应力共同作用导致靶材发生蚀除,蚀除机制是基于熔化与应力作用结合的热机械蚀除[14]㊂结合烧蚀阈值附近的烧蚀现象,可以得出,脉宽为100f s ,能量密度位于烧蚀阈值附近的激光与B 2结构镍钛合金相互作用时,交替出现了机械破碎和热机械蚀除的现象㊂2.3 较高能量密度的烧蚀现象在烧蚀阈值附近,能量密度相近的激光烧蚀,靶材内部温度场和压力场时空分布基本相同的情况下,交替出现了机械破碎和热机械蚀除的现象㊂极有可能在烧蚀过程中,两种蚀除机制是同时作用的,但是在某种参数条件下,一种蚀除机制可能呈现出较强的效应㊂为了验证这一假设,我们观察了较高能量密度下的烧蚀现象㊂脉宽为100f s ,激光能量密度为50m J /c m 2的激光脉冲与B 2结构镍钛合金靶材相互作用时,靶材内部温度场㊁压力场时空分布以及原子位型演化图分别如图4(a ),(b ),(c)所示㊂F i g .4 T h e t e m p e r a t u r e a n d p r e s s u r e f i e l dd i s t r i b u t i o n i n t i m e a n d s p a c e ,a sw e l l a s t h e s n a ps h o t s o f a t o m i c p o s i t i o n s f r o mt h e M Ds i m u l a t i o no f l a s e r i r r a d i a t i o no naB 2N i T i a l l o y fi l m w i t h t h e p u l s e d u r a t i o no f 100f s a n d f l u e n c e o f 50m J /c m 2图4 脉宽为100f s ,能量密度为50m J /c m 2的激光辐照B 2结构镍钛合金靶材分子动力学模拟温度场㊁压力场时空分布以及瞬时原子位型图 图4(a )中,靶材30n m 以上区域发生了热熔化,30~55n m 区域的晶格受热影响发生了无序化相变,对比图4(a )和图2(a )可知,图4(a )中发生热熔化和无序化相变的区域比图2(a)中大㊂由此可得,相同脉宽的激光烧蚀时,能量密度越高,靶材发生热熔化和无序化相变的区域越大㊂结合4(a ),(b ),(c )可得,脉宽为100f s ,能量密度为50m J /c m 2的激光与B 2结构镍钛合金靶材相互作用时,热机械蚀除和机械破碎同时作用,分别于20n m 处和35n m 处诱使靶材发生蚀除,假设得到了验证㊂由图1(b ),(c )和图4(c )可得,脉宽为100f s 的激光与B 2结构镍钛合金靶材相互作用,能量密度为40.5m J /c m 2时,未蚀除靶材发生无序化相变区域约为14n m ;能量密度为41m J /c m 2时,未蚀除靶材发生无序化相变区域约为26n m ;能量密度为50m J /c m 2时,未蚀除靶材发生无序化相变区域约为26n m ㊂烧蚀阈值附近时,机械破碎相对于热机械蚀除而言祛除的靶材更多,未蚀除靶材受热影响发生无序化相变的区域更小㊂随着激光能量密度的增加,未蚀除靶材受热影响发生无序化相变的区域增大㊂强激光与粒子束唐一波等:飞秒激光烧蚀镍钛形状记忆合金的蚀除机理3结论采用结合双温模型的分子动力学的方法,数值模拟了脉宽为100f s,能量密度在25~50m J/c m2区间的飞秒激光与B2结构镍钛合金相互作用㊂结果发现,脉宽为100f s的激光与镍钛合金靶材相互作用时,烧蚀阈值为40.5m J/c m2;烧蚀阈值条件下,靶材蚀除机制是单纯基于应力作用的机械破碎㊂在烧蚀阈值附近,靶材蚀除的机制是热机械蚀除和机械破碎同时作用,但在特定参数条件下其中的一种蚀除机制效应更明显,从而呈现出随能量密度变化交替出现热机械蚀除和机械破碎现象;较高能量密度条件下,同时呈现热机械蚀除和机械破碎机制㊂烧蚀阈值附近,未蚀除靶材发生无序化相变的区域较小,且随激光能量密度的下降而减小,由此可见,烧蚀阈值附近是镍钛形状记忆合金微纳㊁无相变加工的理想参数区间㊂参考文献:[1] A l v a r e z-P u e b l aR,L i z-M a r z췍nL M,G a r cía d eA b a j oF J.L i g h t c o n c e n t r a t i o n a t t h e n a n o m e t e r s c a l e[J].T h e J o u r n a l o f P h y s i c a l C h e m i s-t r y L e t t e r s,2010,1(16):2428-2434.[2] L i uC h a n g,L i F e n g,M aL a i p e n g,e t a l.A d v a n c e dm a t e r i a l s f o r e n e r g y s t o r a g e[J].A d v a n c e dM a t e r i a l s,2010,22(8):E28-E62.[3] W a n g X u d o n g.P i e z o e l e c t r i c n a n o g e n e r a t o r s h a r v e s t i n g a m b i e n tm e c h a n i c a l e n e r g y a t t h e n a n o m e t e r s c a l e[J].N a n o E n e r g y,2012,1(1):13-24.[4] H uR e n z o n g,Z h uM i n,W a n g H u i,e t a l.S n b u f f e r e d b y s h a p em e m o r y e f f e c t o fN i T i a l l o y s a s h i g h-p e r f o r m a n c e a n o d e s f o r l i t h i u mi o n b a t-t e r i e s[J].A c t aM a t e r i a l i a,2012,60(12):4695-4703.[5] F uY o n g q i n g,D uH e j u n,H u a n g W e i m i n,e t a l.T i N i-b a s e d t h i n f i l m s i nM E M S a p p l i c a t i o n s:a r e v i e w[J].S e n s o r s a n dA c t u a t o r sA:P h y s-i c a l,2004,112(2):395-408.[6] B a n s i d d h iA,S a r g e a n tTD,S t u p p S I,e t a l.P o r o u sN i T i f o r b o n e i m p l a n t s:a r e v i e w[J].A c t aB i o m a t e r i a l i a,2008,4(4):773-782.[7]杨建军.飞秒激光超精细 冷 加工技术及其应用(I)[J].激光与光电子学进展,2004,41(3):42-52.(Y a n g J i a n j u n.F e m t o s e c o n d l a s e rC o l d m i c r o-m a c h i n i n g a n d i t s a d v a n c e d a p p l i c a t i o n s.L a s e r&O p t o e l e c t r o n i c sP r o g r e s s,2004,41(3):42-52)[8] G o l d m a n JR,P r y b y l a JA.U l t r a f a s t d y n a m i c s o f l a s e r-e x c i t e d e l e c t r o nd i s t r i b u t i o n s i n s i l i c o n[J].P h y s i c a l R e v i e wL e t t e r s,1994,72(9):1364-1368.[9]江伟华.高重复频率脉冲功率技术及其应用:(1)概述[J].强激光与粒子束,2012,24(1):10-15.(J i a n g W e i h u a.R e p e t i t i o nr a t e p u l s e dp o w e r t e c h n o l o g y a n d i t s a p p l i c a t i o n s:(i)I n t r o d u c t i o n.H i g hP o w e rL a s e r a n dP a r t i c l eB e a m s,2012,24(1):10-15)[10]曾曙光,张彬,孙年春,等.啁啾匹配O P C P A方案用于消除超短激光脉冲的预脉冲[J].强激光与粒子束,2012,24(2):394-398.(Z e n gS h u g u a n g,Z h a n g B i n,S u n N i a n c h u n,e ta l.C h i r p-m a t c h i n g O P C P As c h e m e f o re l i m i n a t i n gp r e-p u l s eo fu l t r a s h o r t l a s e r p u l s e.H i g h P o w e rL a s e r a n dP a r t i c l eB e a m s,2012,24(2):394-398)[11]丛培天,邱爱慈.快脉冲直线变压器气体开关技术[J].强激光与粒子束,2012,24(6):1263-1268.(C o n g P e i t i a n,Q i uA i c i.R e v i e wo n g a ss w i t c h e s d e v e l o p e d f o r f a s t l i n e a r t r a n s f o r m e r d r i v e r.H i g hP o w e rL a s e r a n dP a r t i c l eB e a m s,2012,24(6):1263-1268)[12]郭建增,刘铁根,池伟,等.大口径精密光学调整架的优化设计[J].强激光与粒子束,2013,25(2):301-304.(G u oJ i a n z e n g,L i uT i e g e n,C h iW e i,e t a l.O p t i m i z a t i o n a n d t e s t o f o p t i c a l p r e c i s i o nm o u n tw i t h l a r g e a p e r t u r e.H i g hP o w e rL a s e r a n dP a r t i c l eB e a m s,2013,25(2):301-304)[13] D e k e l E,E l i e z e r S.M e a s u r e m e n t so f l a s e r d r i v e ns p a l l a t i o n i n t i na n dz i n cu s i n g a no p t i c a l r e c o r d i n g v e l o c i t y i n t e r f e r o m e t e r s y s t e m[J].J o u r n a l o f A p p l i e dP h y s i c s,1999,86(8):4242-4248.[14] Z h i g i l e i L V,G a r r i s o nBJ.M i c r o s c o p i cm e c h a n i s m s o f l a s e r a b l a t i o no f o r g a n i c s o l i d s i n t h e t h e r m a l a n d s t r e s s c o n f i n e m e n t i r r a d i a t i o n r e-g i m e s[J].J o u r n a l o f A p p l i e dP h y s i c s,2000,88(3):1281-1298.[15] L a iW S,L i uBX.L a t t i c e s t a b i l i t y o f s o m eN i-T i a l l o yp h a s e s v e r s u s t h e i r c h e m i c a l c o m p o s i t i o na n dd i s o r d e r i n g[J]J o u r n a l o f P h y s i c s:C o n d e n s e dM a t t e r,2000,12(5):L53-L60.[16] L i JH,D a i XD,L i a n g SH,e t a l.I n t e r a t o m i c p o t e n t i a l s o f t h e b i n a r y t r a n s i t i o nm e t a l s y s t e m s a n d s o m e a p p l i c a t i o n s i nm a t e r i a l s p h y s i c s[J].P h y s i c s R e p o r t s,2008,455(1/3):1-134.[17]王新林,吴鹤.热物性参量对飞秒激光烧蚀金属影响的分子动力学模拟[J].光子学报,2009,38(12):3052-3056.(W a n g X i n l i n,W uH e.I n-f l u e n c e o f t h e r m o p h y s i c a l p a r a m e t e r sb y f e m t o s e c o n dl a s e ra b l a t i o no fm e t a l s:m o l e c u l a rd y n a m i c ss i m u l a t i o n.A c t aP h o t o n i c aS i n c a,2009,38(12):3052-3056)[18]S c häf e rC,U r b a s s e kH M,Z h i g i l e i LV,e t a l.P r e s s u r e-t r a n s m i t t i n g b o u n d a r y c o n d i t i o n s f o rm o l e c u l a r-d y n a m i c s s i m u l a t i o n s[J].C o m p u-t a t i o n a lM a t e r i a l sS c i e n c e,2002,24(4):421-429.[19] C h e n g C h a n g r u i,X uX i a n f a n.M e c h a n i s m s o f d e c o m p o s i t i o n o fm e t a l d u r i n g f e m t o s e c o n d l a s e r a b l a t i o n[J].P h y s i c a l R e v i e wB,2005,72:165415.091025-5强激光与粒子束M e c h a n i s mo f f e m t o s e c o n d l a s e r a b l a t i n g N i T i s h a p em e m o r y a l l o yT a n g Y i b o1, C h e nB i n g2, C h e nZ h i y o n g1, Z h u W e i h u a1, L iY u e h u a1, W a n g X i n l i n1,2(1.D e p a r t m e n t o f E l e c t r i c a lE n g i n e e r i n g,U n i v e r s i t y o f S o u t hC h i n a,H e n g y a n g421001,C h i n a;2.D e p a r t m e n t o f M e c h a n i c a lE n g i n e e r i n g,U n i v e r s i t y o f S o u t hC h i n a,H e n g y a n g421001,C h i n a)A b s t r a c t: T h em e c h a n i s mo f f e m t o s e c o n d(f s)l a s e r a b l a t i n g B2N i T i a l l o y h a v e b e e n i n v e s t i g a t e db y t h em o l e c u l a r d y n a m i c s i m u l a t i o n s c o m b i n e dw i t h t h e t w o-t e m p e r a t u r em o d e l.As e r i e s o f s i m u l a t i o n sh a v eb e e nc o n d u c t e d,i nw h i c h,t h e l a s e ru s e d i s 800n mi n c e n t r a lw a v e l e n g t h,100f s i n p u l s ed u r a t i o n,a n d f r o m25m J/c m2t o50m J/c m2a t f l u e n c e.T h e t a r g e t i s90n mi n d e p t h.T h e f l u e n c e t h r e s h o l d f o r100f s p u l s e d l a s e r a b l a t i n g N i T i s h a p em e m o r y a l l o y i s d e t e r m i n e d,a n d i t i s f o u n d t h a t t h e t a r-g e t i r r a d i a t e db y t h r e s h o l df l u e n c e i sa b l a t e d w h o l l y d u e t ot h e t e n s i l ea n dt h eh e a t i n d u c e d p h a s ec h a n g ea r e a i st h es m a l l e s t. W h e n t h e f l u e n c e i s i m p r o v e d,t h e p h a s e c h a n g e a r e a e x p a n d e d.K e y w o r d s:f e m t o s e c o n d l a s e r a b l a t i o n; B2N i-T i a l l o y;t w o-t e m p e r a t u r em o d e l; m o l e c u l a r d y n a m i c s s i m u l a t i o nP A C S:42.50.W k;47.11.M n091025-6。