二元Morse团簇的启发式优化-计算机技术专业论文.docxVIP

独创性声明 本人声明所呈交的学位论文是我个人在导师指导下进行的研究工作及取得的 研究成果。尽我所知，除文中已经标明引用的内容外，本论文不包含任何其他个人 或集体已经发表或撰写过的研究成果。对本文的研究做出贡献的个人和集体，均已 在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。 学位论文作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，即：学校有权 保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借 阅。本人授权华中科技大学可以将本学位论文的全部或部分内容编入有关数据库进 行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 保密□， 在 年解密后适用本授权书。 不保密□。 （请在以上方框内打“√”） 学位论文作者签名： 指导教师签名： 日期： 年 月 日 日期： 年 月 日 摘 要 原子分子团簇的优化主要是对其几何构型的优化。由于团簇会根据其原子或分 子的分布位置不同而呈现出不同的物理化学性质，并且这些性质还会随着所包含原 子、分子的数目变化而产生变化，因此，优化并求解团簇的基态稳定构型是团簇研 究的核心工作。 团簇理论研究方法很多，人们通常采用构建经验势能模型的方法来描述原子或 分子间的相互作用。构建经验势能模型是用一组含有特定参数的表达式来描述原子 间的相互作用力。我们通过求解相应的势能函数的全局最小值来求得团簇的最优构 型。然而团簇的空间构型随着团簇规模增大而呈现指数增长，已被证明是 NP-hard 问题。而对于 NP-hard 问题，以往研究表明不大可能存在多项式时间复杂度的准确 求解算法。因此，我们的研究重点就是设计出高效的启发式算法来解决团簇构型优 化问题。 Morse 团簇势能是一种简单的对势能，只有一个参数，这个参数决定了势能井 的宽度，可以用来对多种原子材料进行模型构建。它同时可以用来作为一个相当严 格的评价全局优化算法优劣性的测试系统。二元 Morse 团簇与一元 Morse 团簇不同 的地方在于，它由两种类型的原子构成，因此在构型的构成上比一元 Morse 更为多 样化，而其理论研究难度也相对高很多。对于一元 Morse 团簇，我们使用改进的动 态格子算法。该算法主要是先初始化一个团簇格局，然后用拟牛顿算法快速收敛得 到一个局部最优值，然后再使用传统的动态格子搜索算法跳出这个局部值寻找更低 的势能构型。在算法的最后一个阶段，我们对动态格子算法作了改进，增加了一个 表面操作算法，这是在动态格子算法基础上做的一个更精细化的全局扰动。实验数 据表明，表面操作算法在动态格子算法的基础上往往还能进一步降低团簇势能，得 到更优构型。经过计算，该算法求得的在材料系数为 3，原子个数为 5-80 及 147 时 的团簇势能值均能快速达到剑桥数据库上面记录的当前最低值。这证明了该算法求 解 Morse 团簇的全局最优构型是高效的并且稳定的。 二元 Morse 团簇由两类原子构成，由于同一个二元团簇存在同分异构体，所以 其结构优化问题涉及到连续优化以及离散优化。连续优化的目的是为了找到整个团 簇的最优构型，而离散优化的目的是在连续优化的基础上找到两类原子的最佳分布 位置。针对二元 Morse 团簇的特点，我们在一元 Morse 团簇优化算法的基础上增加 了一个交换算法，即在得到团簇的最优构型后，我们将两种原子中若干能量较高的 原子进行位置交换。实验数据表明，在使用了交换算法后，团簇势能往往能进一步 降低，得到更优构型。我们用改进的动态格子算法外加交换算法求得了两种原子材 料系数分别为 5 和 7，总原子个数为 38 时，团簇所对应的最优构型。而当前国内外 各学术文献对二元 Morse 团簇的最低势能均无记载。 关键词：团簇，动态格子，启发式算法，交换算法 ABSTRACT The main purpose of the optimization of atomic and molecular clusters is to optimize the geometry. The cluster will present the different physical and chemical properties according to the different locations of the atoms or molecules. And these properties also change along with the number of different atoms and molecules in the cluster. Therefore, it is of great importance to