一种Pd-Ir双金属纳米团簇结构分步优化方法技术

本发明专利技术涉及一种Pd

【技术实现步骤摘要】
一种Pd
‑
Ir双金属纳米团簇结构分步优化方法


[0001]本专利技术属于纳米
，涉及一种Pd
‑
Ir双金属纳米团簇结构分步优化方法。

技术介绍

[0002]金属纳米团簇的几何结构在很大程度上决定了它的光、电、磁性、催化等物理和化学性质；因此，优化与预测金属纳米团簇的稳定结构是研究其诸多特性的基础，也是金属纳米团簇研究领域中的关键环节。目前，由于实验方法精确表征纳米团簇结构的困难,理论计算就成为研究团簇结构的主要手段。理论上，主要从两个角度来描述纳米团簇的结构。一种是电子结构的角度，例如基于量子力学的计算方法(DFT计算，从头算等)都是从电子结构的角度来描述纳米团簇结构。传统的理论计算方法如第一性原理计算，因计算过程复杂、计算量非常大，只能研究尺寸较小和某些具有特定对称性的纳米团簇，具有明显的局限性。另一种是经验势函数的角度，如通过对势(如L
‑
J势能)和多体势(如EAM势)等来描述原子的位置信息从而获得团簇的结构。因此,发展高性能的计算方法对纳米团簇的结构研究具有非常重要的理论意义。

技术实现思路

[0003]有鉴于此，本专利技术的目的在于提供一种Pd
‑
Ir双金属纳米团簇结构分步优化方法。
[0004]为达到上述目的，本专利技术提供如下技术方案：
[0005]一种Pd
‑
Ir双金属纳米团簇结构分步优化方法，该方法包括以下步骤：
[0006]S1：设置待优化的双金属团簇尺寸和构成金属纳米团簇的金属原子种类；
[0007]S2：以系统时间戳为随机数种子，在半径为R0的球类或者边长为R的正方体内，随机生成原子总数为N的双金属纳米团簇群，所产生的金属团簇种群有N个具有对称性的金属原子构成，两种金属团簇的原子个数之比由S1给出；
[0008]S3：考虑采用单金属最优结构作为初始结构的部分候选解，一部分利用S2中随机生成部分解；
[0009]S4：通过BasinHopping算法，简称BH算法，对随机生成的初始结构进行局部优化，采用Gupta势函数计算双金属能量；结合蒙特卡洛随机扰动算法，初步筛选局部最优解；
[0010]S5：通过遗传算法对局部最优结构进行有序处理，再进行选择，交叉，变异操作；
[0011]S6：计算所有个体局部最优解的势能函数值，并按能量高低进行排序，然后采用轮盘赌策略选择两个个体用于后面的交叉、变异操作，在进行后面的操作之前，需要将选中的两个个体按照目前所有种群中能量最低的个体的方向进行调整，与PSO算法中鸟群觅食的方式相同；
[0012]S7：采用球切的方式进行交叉；对选择阶段选出的两个个体P1和P2，通过固定其坐标位置，找出每个个体的质心，然后以质心为球心，以相同的半径对两个个体进行球切，将每个个体分成球内和球外两个部分；最后将个体P1的球内部分与P2的球外部分结合形成两个新个体，或者最后将个体P1的球外部分与P2的球内部分结合形成两个新个体；
[0013]S8：以较低的概率对交叉操作产生的新个体进行变异操作；具体为对个体中的所有原子以0.02的概率进行坐标位置的随机变化，原子的变化范围与其当前坐标有关；
[0014]S9：将新产生的个体与当前局部最优个体进行比较，若新个体优于当前局部最优个体，则用新个体替换当前种群中适应度最差的个体；
[0015]S10：利用BHGA算法对合金团簇的结构进行初步优化后，采用DFT计算对一系列候选解优化，获得稳定构型。
[0016]可选的，所述S2包括以下步骤：
[0017]S21：限定金属原子生成的区域；
[0018]S22：指定生成的金属团簇尺寸及两种原子的数量个数；
[0019]S23：在指定半径球内随机生成3*N个原子坐标；
[0020]S24：初始结构的限制条件如下式所示：
[0021][0022][0023]r
ci
是原子i对应的共价半径。
[0024]可选的，所述S3包括以下步骤：
[0025]S31：在算法的初始化阶段，通过采用结构随机初始化方法和引入单金属团簇最优结构作为初始解；
[0026]S32：将初始结构的个体数量30个进行划分，1～10个和10～20个分别为单金属的最优结构，20～30个则采用随机生成结构。
[0027]可选的，所述S4包括以下步骤：
[0028]S41：采用Gupta势函数计算能量，并利用最陡下降法进行优化，在多个局部最小值中用BH算法进行跳跃，实现全局的搜索；
[0029]S42：通过设置势函数梯度小于10
‑
6和迭代次数作为最陡下降法的限制条件；
[0030]S43：所采用Gupta表达式如下式所示：
[0031][0032]其中r
ij
表示原子i与原子j之间的距离，r0表示两原子间的特征长度，用块状材料的第一近邻距离表示，a0代表晶格参数；A
ij
表示原子间排斥强度系数，ξ
ij
表示跳跃积分，其取值由原子类型决定。
[0033]可选的，所述S7包括以下步骤：
[0034]S71：将父代个体和扰动个体的质心放置到坐标原点并选择一个经过原点的平面；这个平面将父代个体和扰动个体各切成两部分，若被切成的两部分恰好各拥有为L和N
‑
L个原子，则交换两个具有相同原子数的部分结构，从而拼接成新的团簇；否则，对选择的平面进行旋转，直到能够把两个父代个体切出的两部分拥有相同的原子数；
[0035]S72：假设两个父代个体P1，P2，球切交叉法通过计算各个父代个体的球心和半径，选出一个切球面，利用切球面代替切平面来实现父代个体的剪切和拼接；切球面将父代个体分成球内和球外两部分，如果两个父代个体的球内部分恰好具有相同的原子数，则交换球内部分产生新的子代个体又称为试验个体，即子代个体是由P1的球内部分和P2的球外部分组成。
[0036]可选的，所述S10包括以下步骤：
[0037]S101：在DFT计算方法中，采用的软件包是分子动力学模拟软件包VASP，是基于DFT并利用平面波赝势方法进行从头算的软件；在VASP中，对价电子与离子实之间的相互作用采用投影缀加投影波方法PAW来描述；
[0038]S102：通过将团簇放在一个足够大的超晶胞中，确保有的真空分离；所有候选构型都采用DFT方法局部优化，构型中的所有原子都做驰豫，直到原子间的相互作用力小于计算结束的条件是体系的总能量收敛到10
‑6eV，即达到电子基态的标准。
[0039]本专利技术的有益效果在于：本专利技术从智能计算的角度出发，针对物理、化学等领域亟需解决的纳米团簇结构优化问题，研究专利技术一种BHGA结合DFT的双金属纳米图簇结构优化方法，兼容准确性搞，收敛性好，计算小。
[0040]本专利技术的其他优点、目标和特征在某种程度上将在随后的说明书中进本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种Pd
‑
Ir双金属纳米团簇结构分步优化方法，其特征在于：该方法包括以下步骤：S1：设置待优化的双金属团簇尺寸和构成金属纳米团簇的金属原子种类；S2：以系统时间戳为随机数种子，在半径为R0的球类或者边长为R的正方体内，随机生成原子总数为N的双金属纳米团簇群，所产生的金属团簇种群有N个具有对称性的金属原子构成，两种金属团簇的原子个数之比由S1给出；S3：考虑采用单金属最优结构作为初始结构的部分候选解，一部分利用S2中随机生成部分解；S4：通过Basin Hopping算法，简称BH算法，对随机生成的初始结构进行局部优化，采用Gupta势函数计算双金属能量；结合蒙特卡洛随机扰动算法，初步筛选局部最优解；S5：通过遗传算法对局部最优结构进行有序处理，再进行选择，交叉，变异操作；S6：计算所有个体局部最优解的势能函数值，并按能量高低进行排序，然后采用轮盘赌策略选择两个个体用于后面的交叉、变异操作，在进行后面的操作之前，需要将选中的两个个体按照目前所有种群中能量最低的个体的方向进行调整，与粒子群优化PSO算法中鸟群觅食的方式相同；S7：采用球切的方式进行交叉；对选择阶段选出的两个个体P1和P2，通过固定其坐标位置，找出每个个体的质心，然后以质心为球心，以相同的半径对两个个体进行球切，将每个个体分成球内和球外两个部分；最后将个体P1的球内部分与P2的球外部分结合形成两个新个体，或者最后将个体P1的球外部分与P2的球内部分结合形成两个新个体；S8：以较低的概率对交叉操作产生的新个体进行变异操作；具体为对个体中的所有原子以0.02的概率进行坐标位置的随机变化，原子的变化范围与其当前坐标有关；S9：将新产生的个体与当前局部最优个体进行比较，若新个体优于当前局部最优个体，则用新个体替换当前种群中适应度最差的个体；S10：利用盆地跳变遗传算法BHGA对合金团簇的结构进行初步优化后，采用密度泛函理论DFT计算对一系列候选解优化，获得稳定构型。2.根据权利要求1所述的一种Pd
‑
Ir双金属纳米团簇结构分步优化方法，其特征在于：所述S2包括以下步骤：S21：限定金属原子生成的区域；S22：指定生成的金属团簇尺寸及两种原子的数量个数；S23：在指定半径球内随机生成3*N个原子坐标；S24：初始结构的限制条件如下式所示：S24：初始结构的限制条件如下式所示：r
ci
是原子i对应的共价半径，M表示的是结构初始化过程中初始个体的数量。3.根据权利要求2所述的一种Pd
‑
Ir双金属纳米团簇结构分步优化方法，其特征在于：所述S3包括以下步骤：S31：在BHGA的初始...

【专利技术属性】
技术研发人员：范天娥，蒋为，张艳如，易守银，杜培宁，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：