Ag低指数表面非谐效应的分子动力学模拟
杨剑瑜1, 胡望宇1, 2
(1. 湖南大学 应用物理系, 长沙 410082;
2. 湖南大学 材料科学与工程学院, 长沙 410082)
摘 要: 应用分子动力学, 结合分析型嵌入原子模型研究了Ag低指数面的声子频率、 声子线宽度、 原子振动振幅和层结构因子随温度的变化关系。 结果表明, 计算值与已有的实验数据比较相符; (110)面的原子振动振幅表现出明显的各向异性; 由声子频率随温度变化关系得到的非谐常数表明(111)面的非谐效应最弱。 计算出的反映表面长程有序特征的层结构因子表明, (111)面没有表面熔化, 而(100)、 (110)面分别在低于体熔点的100K与250K时出现了表面熔化。
关键词: 表面非谐效应; 分子动力学; 嵌入原子模型 中图分类号: O736
文献标识码: A
Molecular dynamics simulation of anharmonic effects for low Miller index surfaces of Ag
YANG Jian-yu1, HU Wang-yu1, 2
(1. Department of Applied Physics, Hunan University, Changsha 410082, China;
2. College of Materials Science and Engineering, Hunan University, Changsha 410082, China)
Abstract: The anharmonic effect of temperature on phonon frequency, line width, vibrational amplitude and structure factor of atoms for the low index surface of Ag were studied by molecular dynamics simulation and the analytic embedded atom method. The results show that the calculated values agree with the experimental data. The mean square amplitude of (110) plane surface exhibits anisotropic behavior. The calculated anharmonic constants indicate that the anharmonic effect of (111) plane is the smallest. The surface melting temperature is determined by calcu-lating layer order structure factor. The (111) plane surface does not exhibit surface melting, but the (100) and (110) plane surfaces begin to melt at 100K and 250K below bulk melting temperature, respectively.
Key words: surface anharmonic effect; molecular dynamics; embeded atom method
表面非谐效应与表面扩散、 热膨胀、 相变等有重要关系, 对它的研究显得尤为重要。 晶体的各向异性, 使得表面非谐效应强弱随元素和表面不同而不同。 本文作者利用分析型嵌入原子模型(Embeded atom method, EAM)结合分子动力学计算反映Ag低指数表面非谐效应的几个物理量, 从而得到非谐效应随表面取向的关系。
1 理论模型与计算方法
从原子尺度上研究材料的结构和性能的关键是确定原子间的相互作用势函数。 嵌入原子方法在20世纪80年代初期由Daw和Baskes提出[1]。 由于其简单方便, 且物理本质较为深刻, 已成功应用于计算材料的体性质、 晶界、 原子团簇等方面[2]。 本文作者运用分析型EAM模型[3-5]研究Ag的表面非谐效应, 并运用此模型计算了Ag的表面声子谱[6]。
在某一温度下, 在等压等温(NPT)分子动力学下, 计算出Ag在某一温度下, 零压时完整晶格的晶格常数, 然后利用这个晶格常数构建表面, x, y方向运用周期性边界条件, z方向运用自由边界条件。 (111)和(100)面分别为21层各由4032和2688个原子组成, (110)面由31层2976个原子组成。 在等容等温(NVT)分子动力学下, 表面弛豫到平衡状态后, 系统在平衡状态下运行较长的时间, 统计出Ag 3个低指数面的各种非谐性质。
2 结果与讨论
2.1 表面声子态密度与线宽度
表面声子态密度可以通过计算动力学结构因子或位移自关联函数得到。 由于原子的位置难以确定, 一般情况下可以通过修正后的速度自关联函数[7]计算。 运用高斯函数光滑声子态密度, 然后由声子态密度可以计算出表面声子线宽度(表面模在半峰值时的半频宽)。
图1所示为Ag 3个低指数面在二维布里渊区边界上S1(z)模的频率、 线宽度随温度变化的关系。 利用本模型结合分子动力学(MD)与晶格动力学(LD)[6]计算出表面声子频率, 将其与实验结果(Exp)列于表1。 由表1可看出, 通过分子动力学计算出的(111)和(100)面在M-点, (110)面在Y-点的几个典型表面波频率与晶格动力学计算结果与实验结果相当吻合。 Bracco等[8]测得Ag(110)面S1(x)、 S1(z)两个表面模的线宽度都是40GHz, 计算得到两个表面模的线宽度都为110GHz, 且两者大致相符。 由图1可看出, 随着温度的增加, 频率一般逐渐减小, 而线宽度则逐渐增加, 声子频率与线宽度几乎同步变化。 (100)与(110)面的S1(z)模线宽度分别在800K与700K时突然增加, 表明在此温度下, 两个表面表现出较强的非谐效应。
通常假设表面声子频率随温度的变化在一级近似下与简谐近似下的表面声子频率成一定比例。 从以上得到的表面声子频率随温度变化关系可以计算衡量表面非谐效应强弱的一个重要参数, 即非谐常数。 它与温度无关, 能很好地描述表面非谐效应, 计算公式如下[9]:
hω(T)=hω0-Xehω0(2n0+1)(1)
式中 Xe为非谐常数; ω0为简谐近似下的表面声
图1 3个面S1(z)模的频率、 线宽度与温度的变化曲线
Fig.1 Change curves of frequency,line width and temperature of three surfaces for S1(z) mode
子频率; ω(T)是温度为T时同一个表面模的频率; n0=[exp(hω0/kT)-1]-1。 用最外层沿x、 y、 z方向极化的表面模计算了3个面的非谐常数。 (111)面沿x、 y、 z方向的非谐常数分别是0.0055, 0.0016, 0.0032; (100)面的为0.0078, 0.0078, 0.0071; (110)面的为0.0133(0.0072), 0.0064, 0.0063。 相对于(111)面, (100)和(110)面的非谐常数较大, 说明其表现出的非谐效应较强。
表1 Ag低指数面在300K时表面声子频率
Table 1 Calculated phonon frequencies (ω) of low index surface of Ag at 300K
2.2 原子振动振幅平方
原子振动振幅的平方〈u2〉可以通过统计每一层原子的位置得到[7]。 在简谐近似下, 原子振动振幅的平方与温度呈线性关系。 图2所示为本研究计算得到的Ag 3个低指数面的原子振动振幅平方与温度的关系[12]。 某些高温下原子的〈u2〉并未在图中列出, 是由于在此温度时, 表面原子扩散或蒸发, 统计出的〈u2〉已没有意义。
从图2可以看出, 在温度较低时, 〈u2〉与温度近似成线性关系; 随温度升高, 其值加速增加。 3个面在同一温度下(如600K)时, (110)面的〈u2〉最大, (100)面的次之, 而(111)面的最小。 从原子振动振幅平方的比较, (110)面表现出最强的非谐效应。 (110)面整个温度范围内最外层原子振动振幅平方都表现出各向异性, 且z方向的原子振动振幅平方小于x, y方向的原子振动振幅平方, 但第二层z方向的原子振动振幅平方大于x, y方向的原子振动振幅平方, 且大于第一层z方向的原子振动振幅平方。 这是由于(110)面开放的表面结构, 第二层原子的软通道[13]直接与真空(相当于硬壁)相联, 使得第二层原子相对第一层原子有更大的振动空间, 而平行于x方向(即[1[TX-]10]方向)的原子具有最大的振动空间, 因此x方向原子的原子振动振幅平方最大。
Busch等[14]测量了300K下Ag(110)面原子的
图2 Ag低指数面原子振动振幅平方与温度的关系
Fig.2 Relationship between atoms vibrational amplitudes square of three surfaces and temperature
原子振动振幅平方, 本研究的结果与实验结果及其他理论结果的比较列于表2。 由表2可看出, 本研究的计算值与实验值非常吻合, 实验值也表现出第一层x方向的原子振动振幅平方大于其他两个方向的原子振动振幅平方, 且第二层z方向的原子振动振幅平方与第一层相差不远, 而Rahman等[15]的计算结果与实验值相差甚远, 且未能很好地反映(110)面第一层x方向原子振动振幅平方大于其他两个方向的原子振动振幅平方。
表2 300K时(110)面原子振动振幅平方与实验值及其他理论值的比较
Table 2 Comparison of atoms mean-square vibrational amplitudes measured by Busch with calculated results and other theory results for (110) surface at 300K (10-2nm2)
2.3 层结构因子
长程有序可以通过统计结构因子(structure factor)[7]S(kl)来衡量, k是二维倒格矢, S(kl)反映了第l层沿k方向的原子有序化程度。 图3所示为3个低指数面第一层的S(kl)与温度的关系, k分别沿平行于表面的x、 y方向。 由图3可看出, S(kl)随温度的升高不断减小, 表明无序度增加。 (111)面在整个计算温度范围内都没有出现表面熔化(一般认为当结构因子接近于0时, 表面熔化)。 当接近于体熔点(本模型计算的完整晶格熔点为1150K)时, (001)面出现了表面熔化, (110)面发生表面熔化温度比熔点低250K左右[16]。 计算结果还表明, 表面结构因子小于体内结构因子, 说明随温度的升高, 晶格从外向内逐层失去有序性。
图3 Ag 3个低指数面第一层沿x、 y方向的结构因子与温度的关系
Fig.3 Relationship between structure factors along x, y direction of first layer for three low index surfaces of Ag and temperature
3 结论
1) 在整个计算温度范围内, 表面声子频率随温度升高逐渐软化, 线宽度逐渐增加。
2) 由表面声子频率计算的非谐常数结果表明, (110)和(100) 面的非谐性较大。
3) 原子振动振幅平方的计算结果说明(110)面表现出最强的非谐效应。
4) 计算出的层间结构因子说明, (111)面表面未熔化, 而(100)和(110)则在低于体熔点时出现了表面熔化。
5) (110)面表现出最强的非谐效应, 而密排的(111)面最稳定。
REFERENCES
[1]Daw M S, Baske M I. Embedded-atom method: derivation and application to impurities, surfaces, and other defects in metals[J]. Phys Rev B, 1984, 29(12): 6443-6453.
[2]Daw M S, Foiles S M, Baske M I. The embedded-atom method—a review of theory and applications[J]. Mater Sci Rep, 1993, 9(7-8): 251-310.
[3]HU Wang-yu, ZHANG Bang-wei, HUANG Bai-yun. et al. Analytic modified embedded atom potentials for HCP metals[J]. J Phys Conden Matter, 2001, 13(6): 1193-1213.
[4]Hu W Y, Fukumoto M. The application of the analy-tic embedded atom potentials to alkali metals[J]. Modelling Simul Mater Sci Eng, 2002, 10(6): 707-726.
[5]HU Wang-yu, DENG Hui-qiu, YUAN Xiao-jian, et al. Point-defect properties in HCP rare earth metals with analytic modified embedded atom potentials[J]. Eur Phys J B, 2003, 34(2): 429-449.
[6]杨剑瑜, 邓辉球, 胡望宇. Ag(110)表面声子谱的分析型EAM模型计算[J]. 物理学报, 2004, 53(6): 1946-1951.
YANG Jian-yu, DENG Hui-qiu, HU Wang-yu. Surface phonon on Ag(110): a study by means of analytic embedded-atom method[J]. Acta Phys Sin, 2004, 53(6): 1946-1951.
[7]AL-Rawi A N, Kara A, Rahman T S. Anharmonic effects on Ag(111): a molecular dynamics study[J]. Surf Sci, 2000, 446(1-2): 17-30.
[8]Bracco G, Bruschi L, Tatarek R, et al. Anomalous linewidth behaviour of the S3 surface resonance on Ag(110)[J]. Euro Phys Lett, 1996, 34(9): 687-692.
[9]AL-Rawi A N, Kara A, Rahman T S. Comparative study of anharmonicity: Ni(111), Cu(111), and Ag(111)[J]. Phy Rev B, 2002, 66(16): 165439.
[10]Ponjee M W G, Flipse C F J, Denier A W. Experimental observation of vibrational modes on Ag(111) along Γ[TX-]M[TX-] and Γ[TX-]K[TX-] [J]. Phys Rev B, 2003, 67(17): 174301.
[11]Tatarek R, Bracco G, Tommasini F, et al. Surface resonant phonons of Ag(110)[J]. Surf Sci, 1989, 211-212(2): 314-322.
[12]Botez C E, Elliott W C, Miceli P F. Thermal expansion of the Ag(111) surface measured by X-ray scattering[J]. Phys Rev B, 2001, 63(15): 113404.
[13]Marari N, Vanderbilt De.Vita A D, et al. Thermal contraction and disordering of the Al(110) surface[J]. Phys Rev Lett, 1999, 82(16): 3296-3299.
[14]Busch W B, Gustafsson T. Temperature dependence structure of clean Ag(110) studied by medium energy ion scattering[J]. Surf Sci, 1998, 407(1): 7-15.
[15]Rahman T S, Tian Z J, Black J E. Surface disordering roughening and premelting of Ag(110)[J]. Surface Sci, 1997, 374(1-3): 9-16.
[16]Molenbroek A M, Frenken J W M. Anharmonicity but absence of surface melting on Al(001)[J]. Phys Rev B, 1994, 50(15): 11132-11141.
(编辑李艳红)
基金项目: 国家自然科学基金资助项目(50371026)
收稿日期: 2005-07-15; 修订日期: 2005-08-20
作者简介: 杨剑瑜(1971-), 男, 讲师, 博士研究生
通讯作者: 胡望宇, 教授, 博士; 电话: 0731-8823971; E-mail: Wangyuhu2001cn@yahoo.com.cn