Efficient stochastic generation of special quasirandom structures.docxVIP

特殊准随机结构的高效随机模拟方法 摘要 我们提出一种新的算法来生成特殊准随机结构(SQS)，即对处在真无序状态的每个给定原子数量的超胞进行最好的周期性超胞近似。这种方法基于蒙特卡洛模拟法模拟退火循环的过程，获得目标函数，寻求与最大值完美匹配的关联函数（而不是仅仅为了缩小SQS关联函数和无序态关联函数的差距而预设一个关联函数）。该方法优化超晶胞的形状的同时优化了原子的站位，因此能保证获得详尽的空间构型并且没有预设的超胞形状的影响。这种方法已经被应用于Alloy Theoretic Automated Toolkit (ATAT) 软件的“mcsqs”代码中，是计算多组元多亚晶格系统的一般框架，在某种程度上减少了用户需要指定输入信息的数量，提高了并形化效率。 关键词：特殊准随机结构、蒙特卡洛、合金原理、第一原理计算方法、无序态、固溶体 1.简介 因为许多具有可调谐性的技术重要材料以无序合金的形式存在，对无序晶体合金准确和有效的建模在材料科学应用广泛。一个流行、简易、通俗的方法是仅仅构建巨大的超胞和随机化每个晶格点阵的站位。然而，这种方法不是最有效率的，由于同元素原子在会占据临近的晶格点阵，影响局部关联函数，在一个有限的超胞内的随机生成结构有相当大的可能性违背了“完美”的随机性。运用SQS的概念可以系统的改进这种方法。SQS代表了在真无序状态的每个给定原子数量的超胞进行最好的周期性超胞近似。SQS是最符合标准的因为在SQS中相邻点阵站位的特定集合的关联函数可以匹配相应的真正完全无序的状态的关联函数。 SQS已经非常成功地用于获得无序材料的电学和热力学性质（见[2-8]）。SQS通常通过完全枚举法，由所有可能的原子站位的超晶胞所构成，这是一个复杂的指数阶的算法。这是唯一已知的保证获得最佳的SQS的方法，但是这只适用于计算易处理的小晶胞（当前最多可以处理25个原子，即使用效率最高的结构枚举算法[9]）。然而，电子结构方法现在可以正常处理相当大的晶胞（超过50个原子），这使得完全枚举法的使用不切实际。所以，现在一般通过随机数法来替代生成SQS[5,10,11]。 我们的成果主要有四方面。一、设计了一种新的SQS生???算法，寻找与最大值完美匹配的关联函数的SQS（而不是仅仅为了缩小SQS关联函数和无序态关联函数的差距而预设一个关联函数）。我们归纳了一个通常用于随机SQS研究的被称为“rewards”的目标函数，在误差范围内完美相关。二、我们运用了一种计算多组元多亚晶格系统的一般框架的方法[12]来研究大类技术相关合金（大多数在不同的亚晶上展现出无序性）。三、这种方法同时优化超晶胞的形状和原子的站位，因此可以获得详尽的空间结构并没有预设的超胞形状的影响。最后，我们在ATAT上实施了这个算法[12-14]，在某种程度上减少了用户需要指定输入信息的数量，提高了并行化效率。 2.聚簇展开形式 让我们先简单介绍一下多组元多亚晶格系统的聚簇展开形式[12]。σ代表变量的矢量，σi指占据晶格上i点的矢量，σi=0，…，Mi=1,如果Mi代表的元素种类只能占据点阵i。之后，我们需要一种方法来量化原子占据多重点阵如i1，i2之间的相关性。 在计算中，有着特定相关性的点阵的集合被称为聚簇，用α表示。在二元系统中（每个点自由度都为一），定义哪个阵点在聚簇中至关重要（阵点在聚簇中时αi=1，不在时αi=0）。在多组元系统之中，聚簇中的每个阵点i都有一个αi值，范围从0到Mi-1，这表示了哪个组元的自由度被考虑了（αi=0时阵点不在集团之中）。关联函数与聚簇之间的关系是 …α表示所有与α聚簇镜面对称的等价聚簇的平均值，Γα(σ)是一个聚簇函数，定义为 其中γаi，Miσi满足γ0，Miσi=1，已经以下正交关系 （二元合金之中，通常选择γ0，20=1，γ0，21=1，γ1，20=-1，γ1，21=+1，这样能直接关联到著名的伊辛模型。）多组元系统在在ATAT中通常的选择方法也已经叙述过了[12]。尽管在原理上来说所有的晶体阵点都可以算出结果，但是通常选择γ0，Miσi=1，因为这样只输出在聚簇α上的阵点的结果。 3.运算法则 我们的运算法则能这样解释：ρα（σrnd）表示完全无序状态的给定组元的关联函数，ρα（σ）表示候选的SQS相σ的关联函数。Δρασ=ρασ-ρα（σrnd）。因为在无序状态下，阵点的站位都是独立的，所以ρα（σrnd）可以很容易的算出： γai，Miσiα可以通过平均组元的阵点i计算得出。 我们的目标函数（为获得最小值）为： 其中，L表示l的最大值，A表示用户自定义的聚簇，表示用户自定义的重量。 目标函数（3）中的-ωL是我们的方法中一个重要的特色。它反映了SQS的质量通常是用来衡量完全无序状态的关联函数的值，并能够匹配完全正确。特别的，当逐