Gupta团簇的2次融化过程及其与结构变化的关系.pdfVIP

( ) 第 30 卷 第 1 期 河北师范大学学报 自然科学版 Vol. 30 No. 1 2006 年 　1 月 Journal of Hebei Normal University (Natural Science Edition) Jan. 2006 Gupta 团簇的2 次融化过程及其与结构变化的关系 任 世 伟 (河北科技大学 理学院 ,河北 石家庄　050018) 摘 要 :利用微正则经典分子动力学方法 ,对具有 Gupta 势的 13 粒子团簇的融化过程进行了研究. 发现 在 Gupta 势的势参量 p 和 q 取值合适的情况下 ,13 粒子团簇的融化过程具有 2 步融化现象. 进一步利用原子 σ( ) 等价系数 i t 的研究还发现 ,团簇的第 1 步融化仅涉及表面粒子的异构化过程 ,而第 2 步融化中的异构化过 程则涉及团簇中所有粒子. 关键词 : 团簇; 分子动力学 ; 相对均方根键长涨落 ( ) 中图分类号:O 561. 1 　　　文献标识码 :A 　　　文章编号 2006 团簇作为一种新的物质形态 ,呈现出许多新的、奇异的物理和化学特性 ,而团簇热力学特性的变化 则是其中最基本变化之一. 对这些现象的研究可以极大地推动和促进许多相关学科 ,尤其是纳米科学与 技术方面的基础及应用研究工作. 研究表明[1～5 ] ,团簇的融化过程是非常复杂的. 人们使用各种手段和 方法对其进行了广泛的研究 ,其中包括采用各种不同作用势模型的分子动力学以及蒙特卡罗方法. 在对 13 粒子团簇的研究中 , 人们发现某些团簇具有 2 次融化过程[5～7 ] . 但截止到目前 , 人们仅仅在具有 LennardJones 等少数作用势的团簇中发现了此种现象 ,对 13 粒子团簇表现出的 2 步融化过程的认识 还很不深入. 尤其对于多体类型作用势的原子团簇 ,人们目前对其多体部分对团簇热力学过程的影响尚 缺乏较为系统而深入的研究. 本文中 ,笔者使用了微正则经典分子动力学方法 ,对具有 Gupta 势的 13 粒子等团簇的融化过程进 行了详细的研究. 通过对 Gupta 势中表征对势和多体势势范围参量的调节 ,发现在 p 和 q 取值合适的情 况下 ,13 粒子等团簇的融化过程也具有 2 步融化现象. 进一步的研究还发现 , 团簇发生这种现象是和其 结构变化紧密相联的. 1 　原子间作用势与计算过程 Gupta 势[1 ]形式为 n n n Un 1/ 2 ( ) ( ) ( ) 　　　　V = A exp [ - p r / r - 1 ] - exp [ - 2 q r / r - 1 ] . 1 2 ∑ ∑ ij 0 n ∑ ij 0 n j = 1 i ( ≠j