Mobius反演收敛性分析.docVIP

院级本科生科技创新项目 研 究 报 告 项 目 名 称 M?bius反演收敛性分析 立 项 时 间 2011年4月 计划完成时间 2012年4月 项目负责人 丁以民 项 目 成 员 丁以民、王越、李子扬、陈丹 学院与班级 数理学院物理0901 北 京 科 技 大 学 教 务 处 摘 要 我国稀土和锕系资源储量均居世界前列，稀土和锕系金属作为我国特色资源，不仅为特种功能材料所必需，也是我国今后大力发展的核材料的工业原料。多年来我国稀土行业一直因低技术附加值出口而饱受诟病，其深层次原因是我国缺乏相应的关键技术。针对我国资源优势，经济发展需求和科研需要，运用晶格反演原子间相互作用势对稀土和锕系金属间化合物各种结构、各种体系的研究将为开展相关工作奠定基础。根据实验上实际存在的不同材料，以贴近对象结构的晶格体系，利用DFT+U的方法进行系列性第一原理计算，运用新推导的不同结构反演系数计算出不同相、不同体系的Chen-M?bius晶格反演原子间相互作用势，完成针对不同体系的原子间相互作用势库。探讨各种结构性质、力学、热力学及其它功能性质。利用原子间相互作用势通过分子力学及分子动力学方法获得各种复杂材料的驰豫结构，进一步作第一原理计算，分析各种电子性质。并进一步探讨同一种材料不同温度下不同单一相之间原子间相互作用势的差别。对具有应用潜力的稀土与锕系金属间化合物以及锕系氧化物体系进行深入研究。 近年来从数论基础上发展起来的陈氏三维晶格反演进行了概述, 并介绍了其理论上的严密性以及反演结果的确定性，对它在材料的结构计算、动力学特性的研究、稀土-过渡金属金属间化合物新材料模拟设计方面的应用进行了简单的叙述。根据半群结构模型和数论反演方法推算晶体和界面系统原子间结合能进而获得复杂难处理体系原子间势是得到国际认可的独创性方法，已在诸多材料体系中获得应用。 近十年来，陈难先院士研究组先后发展了单重结合能反演、多重结合能反演和界面粘结能反演等三类原子相互作用势反演方法，用于建立稀土磁性金属间化合物等原子势库，针对国内外合作研究单位提供势库服务，在新型功能材料研究中发挥了作用，在该领域获得了国际同行的好评。 通过“M?bius反演收敛性分析”的有关工作，可以较深入地学习和掌握陈氏晶格反演方法，并有可能为确定和拓宽该方法的使用范围打下一些基础，进而为国家基础研究发展计划（973计划）项目中的课题“复杂材料相互作用势的创新与拓展”项目贡献一点力量。 目 录 1 绪论 11 1.1 项目背景条件 11 1.2 项目意义及主要工作 12 2 文献综述 13 3 基本理论 13 3.1原子间相互作用势理论 13 3.2 陈氏三维晶格反演定理 16 4．三维晶体反演系数的编程计算 18 5．Fe的体心立方(BCC)反演收敛性和原子间相互作用势的计算 19 6．Cu的面心立方(FCC)反演收敛性和原子间相互作用势的计算 22 7.结论 32 参考文献 32 附录：计算反演系数的matlab编程 33 致 谢 46 绪论 项目背景条件 随着近代科学技术的发展，及计算机软硬件的日新月异，使计算模拟在材料科学的研究中已深入到了原子尺度。从某种程度上来说，实际实验研究往往告诉我们发生了什么，而基于计算的数值实验却让我们可以分析发生这种现象的原因，并预测在其它相似的条件下会发生什么。因此，从局部的点缺陷、面缺陷到大规模的裂纹乃至宏观的疲劳断裂，原子级的计算模拟都提供了实验手段难以预测的微观物理图像。而原子间相互作用势是原子级材料计算和设计的基础。从本世纪20年代开始，人们就提出了各种各样的描述原子间相互作用的函数形式，如早期著名的Lennard-Jones势、Morse势、Buckingham势、Born-Mayer势等原子间相互作用对势形式。近廿年来，基于第一性原理计算的原子间相互作用势，由于其强大的理论背景而日益受到人们的重视。人们在不断地探索如何基于第一性原理计算，准确而又简洁地获得原子间相互作用势的方法。 1990年，陈难先先生创造性地将数论中的M?bius变换应用到凝聚态物理的研究，尤其是在三维晶格反演中的应用，给出了一条从第一性原理结合能曲线得到原子间相互作用对势的解析途径。这种自主创新的数论反演方法，在1990 年就得到Nature的评述，并被Physical Review 等重要杂志称为陈氏反演公式。十多年来已运用晶格反演方法获得的原子间相互作用势对体材料作出了大量应用，所涉及的材料有一系列金属和