CN113223641B - 一种空位团簇跳出超势能谷的模拟方法 - Google Patents
一种空位团簇跳出超势能谷的模拟方法 Download PDFInfo
- Publication number
- CN113223641B CN113223641B CN202110467637.1A CN202110467637A CN113223641B CN 113223641 B CN113223641 B CN 113223641B CN 202110467637 A CN202110467637 A CN 202110467637A CN 113223641 B CN113223641 B CN 113223641B
- Authority
- CN
- China
- Prior art keywords
- transition
- energy
- vacancy
- seqtr
- cluster
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005381 potential energy Methods 0.000 title claims abstract description 24
- 238000004088 simulation Methods 0.000 title claims abstract description 15
- 230000009191 jumping Effects 0.000 title claims description 9
- 230000007704 transition Effects 0.000 claims abstract description 219
- 230000005283 ground state Effects 0.000 claims abstract description 42
- 230000004888 barrier function Effects 0.000 claims abstract description 29
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 9
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 241000970807 Thermoanaerobacterales Species 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000005182 potential energy surface Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 description 10
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 9
- 230000007547 defect Effects 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 229910052751 metal Inorganic materials 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- MKYBYDHXWVHEJW-UHFFFAOYSA-N N-[1-oxo-1-(2,4,6,7-tetrahydrotriazolo[4,5-c]pyridin-5-yl)propan-2-yl]-2-[[3-(trifluoromethoxy)phenyl]methylamino]pyrimidine-5-carboxamide Chemical compound O=C(C(C)NC(=O)C=1C=NC(=NC=1)NCC1=CC(=CC=C1)OC(F)(F)F)N1CC2=C(CC1)NN=N2 MKYBYDHXWVHEJW-UHFFFAOYSA-N 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005212 lattice dynamic Methods 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 239000011824 nuclear material Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 239000011343 solid material Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种空位团簇跳出超势能谷的模拟方法，该模拟方法主要步骤有获取空间团簇周围可能的跃迁路径能垒，建立速率表格；获取跃迁序列seqtr，将所述速率表格与所述跃迁序列seqtr进行比对，若存在循环跃迁则修改所述速率表格，禁止循环跃迁；在修改后的速率表格中选择跃迁事件并执行；对比两次跃迁，更新跃迁序列seqtr；迭代后获得体系跃迁路径、能量和结构变化。通过该模拟方法可有效访问基态附近的势能面特征，可使得空位团簇跳出含有多个亚稳态的超势能谷，有效促使体系朝着近邻基态跃迁。
Description
技术领域
本发明属于核材料辐照损伤模拟技术领域，具体涉及一种空位团簇跳出超势能谷的模拟方法。
背景技术
辐照损伤是由于中子、带电粒子或电磁波等和固体材料的点阵原子发生一系列碰撞，引起材料内部出现大量原子尺度的缺陷的过程，这个过程在很短时间内就会发生。这些缺陷经过长时间的迁移、聚集和复合等形成缺陷团簇、空位等，引起材料微观组织变化，使材料的宏观力学、热学等性能退化。因此，理解辐照缺陷团簇(如空位团簇)的行为对认识辐照损伤微结构演化至关重要。
由于含有空位团簇体系往往具有多个局域亚稳态，当体系从一个基态跃迁到近邻基态时需要历经一系列亚稳态。相对于基态而言，亚稳态由于具有较高的能量，导致空位团簇激发以后很容易返回基态或陷入基态附近由多个亚稳态构成的超势能谷中，从而阻碍体系进一步演化。目前通常是采用势能谷自动填充法(ABC)或格点动力学蒙特卡洛(LKMC)方法访问体系势能面，但是当体系中含有诸如一个亚稳态到另一个亚稳态速率快得多的事件时，这些方法往往因体系“自转”而失效，不能真正地使体系跳出势能谷。
发明内容
有鉴于此，本发明有必要提供一种空位团簇跳出超势能谷的模拟方法，该模拟方法通过记录跃迁过程中的跃迁序列和禁止循环跃迁，可有效访问基态附近的势能面特征，可使得空位团簇跳出含有多个亚稳态的超势能谷，有效促使体系朝着近邻基态跃迁。
为了实现上述目的，本发明采用以下技术方案：
本发明提供了一种空位团簇跳出超势能谷的模拟方法，包括以下步骤：
S1、获取空间团簇周围可能的跃迁路径能垒，建立速率表格；
S2、获取跃迁序列seqtr，所述跃迁序列seqtr中记录有体系内连续跃迁的初末态空位编号，将所述速率表格与所述跃迁序列seqtr进行比对，若存在循环跃迁则修改所述速率表格，禁止循环跃迁；
S3、在修改后的速率表格中选择跃迁事件并执行，其中，每个跃迁事件被选中的几率与其速率成正相关；
S4、对比两次跃迁，更新跃迁序列seqtr；
S5、判断①是否跃迁至另一基态；②质心是否发生运动；
若①②均为是，则进入步骤S6，否则返回步骤S1；
S6、S1-S5迭代结束，获得体系跃迁路径、能量和结构变化。
进一步的，步骤S1中，所述可能的跃迁路径能垒的获得，具体包括如下步骤：
获取空位团簇周围可能的跃迁路径：假设空位团簇的跃迁是空位与空位周围第一近邻和第二近邻通过原子交换而进行的，对于含有n个空位的团簇Vn，其中，某个空位Vi坐标为rvi＝[rxvi,ryvi,rzvi]，其周围原子坐标为r＝[rx,ry,rz]，则满足式(1)的原子构成空位Vi的可能跃迁末态
|rvi-r|<disn+ε (1)，
其中，disn为第一近邻距离disn1或第二近邻距离disn2，ε的范围在
获取可能的跃迁能垒：根据获得的可能的跃迁路径计算可能的跃迁能垒。
优选的，所述的根据获得的可能的跃迁路径计算可能的跃迁能垒的步骤，具体为：调用NEB方法对每一条可能的跃迁路径计算其跃迁能垒。
进一步的，步骤S1中，所述的建立速率表格的步骤，具体为：假设获得的第i条跃迁路径的跃迁能垒为Esad，则根据式(2)获得速率表格中第i项速率rate
其中，T为温度，v0为固有跃迁频率，kB为玻尔兹曼常数。
进一步的，步骤S2的具体步骤为：假设速率表格中第i项跃迁的初末态编号分别为gi、gj，比对gj是否与所述跃迁序列中的某项相同，若相同，则修改速率表格，令该项速率为0。
进一步的，步骤S3中，所述的选择跃迁事件采用蒙特卡洛方法进行选择。
进一步的，步骤S4中，所述更新跃迁序列seqtr的具体方法为：设上一次跃迁为从gi到gj，本次为从gi1到gj1，若gj＝gi1，表明体系沿着路径gi-gj-gj1连续跃迁，则将gj1添加到跃迁序列seqtr中；否则跃迁序列seqtr清空，然后将gi1、gj1添加到跃迁序列seqtr中。
进一步的，步骤S5中，所述的是否跃迁至另一基态的判断方法，具体为：
假设按照步骤S3跃迁后的团簇能量为Ei，在程序运行初始团簇基态能量为Ebase，若满足式(3)，则认为已经跃迁至另一基态
|Ei-Ebase|<ε (3)，
其中，ε的范围为0.05-0.08eV。
进一步的，步骤S5中，所述的质心是否发生运动的判断方法，具体为：
假设同时若程序运行初始团簇构型坐标为rVn0，按照步骤S3跃迁后为rVn1，每个坐标方向求和并除以n，获得质心坐标，进而得到质心位移dMC，若满足式(4)，则质心发生运动
进一步的，步骤S6中，所述体系跃迁路径、能量和结构变化通过迭代中输出的相邻两次跃迁的结构、跃迁能量曲线获得。
与现有技术相比，本发明具有以下有益效果：
本发明通过在通常格点动力学蒙特卡洛的框架中，记录含空位团簇体系跃迁过程中的跃迁序列，禁止基态附近可能重新返回基态的亚稳态之间的低能垒跃迁，避免了重复“自转”的无效跃迁，可使得体系跳出含有多个亚稳态的超势能谷。同现有的模拟方法相比，该模拟方法一方面保持了格点动力学蒙特卡洛在模拟团簇长时间跃迁行为方面的优势，另一方面在LKMC框架中抑制了循环跃迁的事件，可有效促使体系朝着近邻基态跃迁。
附图说明
图1为本发明中一些较佳实施方式中空位团簇跳出超势能谷的主要步骤流程图；
图2为图1中的算法流程图；
图3为含空位团簇体系附近势能面特点以及本发明模拟方法的原理示意图；
图4为应用实施例1中含四个空位的团簇(V4)从一个基态跃迁到另一近邻基态的过程。
具体实施方式
为了便于理解本发明，下面将结合具体的实施例对本发明进行更全面的描述。但是，本发明可以以许多不同的形式来实现，并不限于本文所描述的实施方式。相反地，提供这些实施方式的目的是使对本发明的公开内容理解的更加透彻全面。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施方式的目的，不是旨在于限制本发明。
本发明第一方面提供了一种空位团簇跳出超势能谷的模拟方法，在一些实施方式中，主要步骤如图1中所示的，包括以下步骤：
S1、获取空间团簇周围可能的跃迁路径能垒，建立速率表格；
S2、获取跃迁序列seqtr，所述跃迁序列seqtr中记录有体系内连续跃迁的初末态空位编号，将所述速率表格与所述跃迁序列seqtr进行比对，若存在循环跃迁则修改所述速率表格，禁止循环跃迁；
S3、在修改后的速率表格中选择跃迁事件并执行，其中，每个跃迁事件被选中的几率与其速率成正相关；
S4、对比两次跃迁，更新跃迁序列seqtr；
S5、判断①是否跃迁至另一基态；②质心是否发生运动；
若①②均为是，则进入步骤S6，否则返回步骤S1；
S6、S1-S5迭代结束，获得体系跃迁路径、能量和结构变化。
本发明通过记录含空位团簇体系跃迁过程中的跃迁序列seqtr，避免了重复“自转”的无效跃迁，可使得体系跳出含有多个亚稳态的超势能谷，从而有效促进体系朝着近邻基态跃迁。其中，本文中的所述的“势能面”指的是含空位团簇确定的体系的势能随其核位置改变的图形，“超势能谷”指的是在体系基态附近形成的多个亚稳态。请参阅图3，图3(a)示意的是连接体系两个基态Si和Sj中间区域势能面，说明当体系从一个基态跃迁到另一个近邻基态时须经历一系列亚稳态，而图3(b)示意的是本发明基本原理，当体系经过A-B-C-D-E系列跃迁步骤后到达F，若F与C相同，则意味着体系又回过去从C处重复演化了，通过对于E-F路径禁止跃迁，从而避免了重复“自转”的无效跃迁，促进体系朝着近邻基态跃迁。
在本发明的一些实施方式中，步骤S1中，所述可能的跃迁路径能垒的获得，具体包括如下步骤：
获取空位团簇周围可能的跃迁路径：假设空位团簇的跃迁是空位与周围第一近邻和第二近邻通过原子交换而进行的，对于含有n个空位的团簇Vn，其中，某个空位Vi坐标为rvi＝[rxvi,ryvi,rzvi]，其周围原子坐标为r＝[rx,ry,rz]，则满足式(1)的原子构成空位Vi的可能跃迁末态
|rvi-r|<disn+ε (1)，
其中，disn为第一近邻距离disn1或第二近邻距离disn2，ε为一个微小量，其选择的原则是使得空位与空位周围原子距离稍微比要考虑的近邻距离大，具体范围可为在本发明的一些实施方式中，可优选取为/>此外，本发明中第一近邻距离或第二近邻距离，根据金属微观结构的不同而有所不同，举例来说，对于体心立方金属，第一近邻距离为/>第二近邻距离为disn2＝a；对于面心立方金属，/>disn2＝a，其中，a均为晶格常数。
获取可能的跃迁能垒：根据获得的可能的跃迁路径计算可能的跃迁能垒。
在获得跃迁路径后，根据跃迁路径计算跃迁能垒，该计算方法没有特别的限定，可以是本领域中常规的计算方法，在本发明的一些具体的实施方式中，所述的根据获得的可能的跃迁路径计算可能的跃迁能垒的步骤，具体为：调用弹性带算法NEB方法对每一条可能的跃迁路径计算其跃迁能垒，由于NEB方法计算跃迁能垒的方法为本领域中的常规方法，因此，这里不再具体阐述。
在本发明的一些实施方式中，步骤S1中，所述的建立速率表格的步骤，具体为：假设获得的第i条跃迁路径的跃迁能垒为Esad，则根据式(2)获得速率表格中第i项速率rate
其中，T为温度，v0为固有跃迁频率，kB为玻尔兹曼常数。具体的说，在本发明的一些实施例中，Esad单位为eV，温度T单位为K，固有跃迁频率v0＝1012/s，玻尔兹曼常数kB＝8.617×10-5eV/K，从而获得相应的速率，并根据这些速率建立速率表格。
在本发明的一些实施方式中，通过跃迁序列记录器记录跃迁序列seqtr，步骤S2的具体步骤为：假设速率表格中第i项跃迁的初末态编号分别为gi、gj，比对gj是否与所述跃迁序列seqtr中的某项相同，若相同，说明该跃迁已经循环至之前某个态了，则修改速率表格，令对应的第i项跃迁的速率为0。
在本发明的一些实施方式中，步骤S3中，选择跃迁事件采用蒙特卡洛方法进行选择，可以理解的是，本发明中，选择跃迁事件的方法还可以采用本领域中任意的常规方法，这里没有特别的限定。跃迁事件被选择后，执行跃迁事件，则团簇中空位运动至相应位置。
进一步的，步骤S4中，所述的更新跃迁序列seqtr的在一些实施方式中，其具体方法为：设上一次跃迁为从gi到gj，本次为从gi1到gj1，若gj＝gi1，表明体系沿着路径gi-gj-gj1连续跃迁，则将gj1添加到跃迁序列seqtr中；否则跃迁序列seqtr清空，然后将gi1、gj1添加到跃迁序列seqtr中。具体的说，跃迁序列seqtr中记录的是一系列跃迁初末态，因此上一次的末态就是下一次的初态，假设跃迁序列seqtr中第一次记录的是A-B，第二次记录的若是B1-C跃迁，而且B1＝B则跃迁序列seqtr中为：A-B-C，若下一次不是从B连续跃迁则A-B清空，跃迁序列seqtr记录为B1-C。
进一步的，在本发明的一些具体实施方式中，步骤S5中，所述的是否跃迁至另一基态的判断方法，具体为：
假设跃迁后的团簇能量为Ei，初始基态能量为Ebase，若满足式(3)，则认为已经跃迁至另一基态
|Ei-Ebase|<ε (3)，
其中，ε为一个微小量，其选择与计算体系能量的系统误差有关，在本发明的一些实施方式中，取为0.05-0.08eV，优选的，取为0.05eV。
在本发明的一些实施方式中，步骤S5中，所述的质心是否发生运动的判断方法，具体为：
假设同时若初始团簇构型坐标为rVn0，跃迁后为rVn1，每个坐标方向求和并除以n，获得质心坐标，进而得到质心位移dMC，若满足式(4)，则质心发生运动
进一步的，步骤S6中，所述体系跃迁路径、能量和结构变化通过迭代中输出的相邻两次跃迁的结构、跃迁能量曲线获得。
图2中示出了本发明一实施例中模拟方法的算法流程图，具体流程如下：
S101、开始，进入S102；
S102、读入晶体构型坐标，以及空位团簇初始基态坐标，进入步骤S103；
S103、计算空位团簇周围可能的跃迁路径，并根据可能的跃迁路径计算获得对应的跃迁能垒，建立速率表格，进入步骤S104；
S104、将速率表格与跃迁序列记录器中记录的跃迁序列进行比对，判断空位是否循环跃迁，若是，则进入步骤S105，若否，则进入步骤S106；
S105、已经进入循环跃迁的空位对应速率修改为0，进入步骤S106；
S106、获得更新后的速率表格，进入步骤S107；
S107、选择跃迁事件并执行跃迁事件，进入步骤S108；
S108、更新跃迁序列，进入步骤S109；
S109、判断是否找到另一质心运动的基态，若是，则转入步骤S110，若否，则转入步骤S103。
S110、迭代结束，输出体系跃迁路径、能量和结构变化。
下面结合具体的应用实施例对本发明中空位团簇跳出超势能谷的模拟方法进行更加清楚完整的说明。
应用实施例1
本应用实施例采用本发明一较佳实施方式模拟体心立方金属铁(100)表面上空位团簇从一个基态跃迁到另一个近邻基态的过程，其中，体心立方金属铁(100)的晶格常数a为图4中示出了含四个空位的团簇(V4)从一个基态跃迁到另一近邻基态的过程，其中，图4(a)为跃迁过程中体系能量变化曲线，图4(b-f)为跃迁过程中V4结构变化过程，其中黑色球表示一个空位，黑色线连接两个互为第二近邻的空位。
具体模拟过程如下：
S1、以满足公式rvi-r<disn+ε获得可能的跃迁路径，其中，disn1＝2.473，disn2＝2.855，因在铁(100)表面空位最小距离为disn2，此处disn设置为disn2；本根据跃迁路径计算对应的跃迁能垒，根据跃迁能垒以及公式/>(设温度T＝300K，即室温)建立如下所示的速率表格：
此时，跃迁速率记录器初始化为seqtr:-1-1,gi＝-1,gj＝-1。
S2、在上述速率表格中进行跃迁事件选择，选中行号为2的项，执行后，seqtr＝311295，gi＝311,gj＝295。空位跃迁后，速率表格如下：
经比对后，发现其中行号为6的事件出现了311-295-311循环跃迁，则修改速率表格，禁止该跃迁，修改后的速率表格如下所示：
修改后，上表中0行事件被选中，执行后seqtr＝327 311,gi＝327,gj＝311；
空位跃迁后，速率表格如下：
其中7行事件被禁止，禁止后8行事件被选中执行，同时检测到体系已经跃迁到另一个基态，而且团簇初始坐标为：
| 空位号 | 空位x坐标 | 空位y坐标 | 空位z坐标 |
| 0 | 9.994 | 12.85 | 7.155 |
| 1 | 9.994 | 9.994 | 7.155 |
| 2 | 12.85 | 9.994 | 7.155 |
| 3 | 12.85 | 12.85 | 7.155 |
最终坐标为：
| 空位号 | 空位x坐标 | 空位y坐标 | 空位z坐标 |
| 0 | 9.993592 | 9.993592 | 7.154898 |
| 1 | 9.993592 | 7.138280 | 7.154898 |
| 2 | 12.848904 | 7.138280 | 7.154898 |
| 3 | 12.848904 | 9.993592 | 7.154898 |
质心坐标分别为：
MC0＝[11.4220 11.4220 7.1550]
MC1＝[11.4212 8.5659 7.1549]
计算得到质心为dMC＝2.8561，满足跃迁条件，程序停止。
通过以上实施例，可以看出，通过记录含空位团簇体系跃迁过程中的跃迁序列，禁止基态附近可能重新返回基态的亚稳态之间的低能垒跃迁，避免了重复“自转”的无效跃迁，可使得体系跳出含有多个亚稳态的超势能谷，可有效促使体系朝着近邻基态跃迁。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种空位团簇跳出超势能谷的模拟方法，其特征在于，包括以下步骤：
S1、获取空间团簇周围可能的跃迁路径及其能垒，建立速率表格；
S2、获取跃迁序列seqtr，所述跃迁序列seqtr中记录有体系内连续跃迁的初末态空位编号，将所述速率表格与所述跃迁序列seqtr进行比对，若存在循环跃迁则修改所述速率表格，禁止循环跃迁；
S3、在修改后的速率表格中选择跃迁事件并执行，其中，每个跃迁事件被选中的几率与其速率成正相关；
S4、对比两次跃迁，更新跃迁序列seqtr；
S5、判断①是否跃迁至另一基态；②质心是否发生运动；
若①②均为是，则进入步骤S6，否则返回步骤S1；
S6、S1-S5迭代结束，获得体系跃迁路径、能量和结构变化；
其中，步骤S1中，所述可能的跃迁路径能垒的获得，具体包括如下步骤：
获取空位团簇周围可能的跃迁路径：假设空位团簇的跃迁是空位与周围第一近邻和第二近邻通过原子交换而进行的，对于含有n个空位的团簇Vn，其中，某个空位Vi坐标为rvi＝[rxvi,ryvi,rzvi]，其周围原子坐标为r＝[rx,ry,rz]，则满足式(1)的原子构成空位Vi的可能跃迁末态
|rvi-r|<disn+ε (1)，
其中，disn为第一近邻距离disn1或第二近邻距离disn2，ε的范围在
获取可能的跃迁能垒：根据获得的可能的跃迁路径计算可能的跃迁能垒；
步骤S1中，所述的建立速率表格的步骤，具体为：假设获得的第i条跃迁路径的跃迁能垒为Esad，则根据式(2)获得速率表格中第i项速率rate
其中，T为温度，v0为固有跃迁频率，kB为玻尔兹曼常数。
2.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，所述的根据获得的可能的跃迁路径计算可能的跃迁能垒的步骤，具体为：调用NEB方法对每一条可能的跃迁路径计算其跃迁能垒。
3.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S2的具体步骤为：假设速率表格中第i项跃迁的初末态编号分别为gi、gj，比对gj是否与所述跃迁序列中的某项相同，若相同，则修改速率表格，令该项速率为0。
4.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S3中，所述的选择跃迁事件采用蒙特卡洛方法进行选择。
5.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S4中，所述更新跃迁序列seqtr的具体方法为：设上一次跃迁为从gi到gj，本次为从gi1到gj1，若gj＝gi1，表明体系沿着路径gi-gj-gj1连续跃迁，则将gj1添加到跃迁序列seqtr中；否则将跃迁序列seqtr清空，然后将gi1、gj1添加到跃迁序列seqtr中。
6.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S5中，所述的是否跃迁至另一基态的判断方法，具体为：
假设按照步骤S3跃迁后的团簇能量为Ei，在程序运行初始团簇基态能量为Ebase，若满足式(3)，则认为已经跃迁至另一基态
|Ei-Ebase|<ε (3)，
其中，ε的范围为0.05-0.08eV。
7.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S5中，所述的质心是否发生运动的判断方法，具体为：
假设同时若程序运行初始团簇构型坐标为rVn0，按照步骤S3跃迁后为rVn1，每个坐标方向求和并除以n，获得质心坐标，进而得到质心位移dMC，若满足式(4)，则质心发生运动
8.如权利要求1所述的空位团簇跳出超势能谷的模拟方法，其特征在于，步骤S6中，所述体系跃迁路径、能量和结构变化通过迭代中输出的相邻两次跃迁的结构、跃迁能量曲线获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110467637.1A CN113223641B (zh) | 2021-04-28 | 2021-04-28 | 一种空位团簇跳出超势能谷的模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110467637.1A CN113223641B (zh) | 2021-04-28 | 2021-04-28 | 一种空位团簇跳出超势能谷的模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113223641A CN113223641A (zh) | 2021-08-06 |
| CN113223641B true CN113223641B (zh) | 2024-02-09 |
Family
ID=77089763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110467637.1A Active CN113223641B (zh) | 2021-04-28 | 2021-04-28 | 一种空位团簇跳出超势能谷的模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113223641B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115394364B (zh) * | 2022-07-12 | 2024-02-02 | 北京应用物理与计算数学研究所 | 原子动力学蒙特卡洛模拟的原子指纹计算方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108052745A (zh) * | 2017-12-15 | 2018-05-18 | 中国科学院合肥物质科学研究院 | 一种模拟辐照缺陷与晶界协同演化的方法 |
| CN111814315A (zh) * | 2020-06-17 | 2020-10-23 | 中国科学院合肥物质科学研究院 | 一种金属材料中缺陷团簇动力学性质的计算方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100968870B1 (ko) * | 2008-08-05 | 2010-07-09 | 인하대학교 산학협력단 | 양자 광 메모리 |
-
2021
- 2021-04-28 CN CN202110467637.1A patent/CN113223641B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108052745A (zh) * | 2017-12-15 | 2018-05-18 | 中国科学院合肥物质科学研究院 | 一种模拟辐照缺陷与晶界协同演化的方法 |
| CN111814315A (zh) * | 2020-06-17 | 2020-10-23 | 中国科学院合肥物质科学研究院 | 一种金属材料中缺陷团簇动力学性质的计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| MD与KMC的耦合模拟研究与实现;李建江;魏鹏;杨少峰;贺新福;胡长军;;电子科技大学学报(04);全文 * |
| 基于空位跃迁的KMC并行实现;宋梦召;冯仰德;聂宁明;王武;;数值计算与计算机应用(02);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113223641A (zh) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Assad et al. | A hybrid harmony search and simulated annealing algorithm for continuous optimization | |
| CN111814315B (zh) | 一种金属材料中缺陷团簇动力学性质的计算方法 | |
| CN113223641B (zh) | 一种空位团簇跳出超势能谷的模拟方法 | |
| Yi et al. | A modified monkey algorithm for optimal sensor placement in structural health monitoring | |
| Pepper et al. | Solving the traveling salesman problem with annealing-based heuristics: a computational study | |
| CN109271683B (zh) | 一种面向日照约束的建筑群自动排布算法 | |
| CN108090279B (zh) | 基于杂化lkmc和okmc模拟纳米结构核材料辐照损伤的方法 | |
| Cho et al. | Suppression effect on explosive percolation | |
| Dyson et al. | Macroscopic limits of individual-based models for motile cell populations with volume exclusion | |
| CN110459269A (zh) | 一种核反应堆材料辐照损伤的多尺度耦合模拟方法 | |
| Ebeida et al. | Improving spatial coverage while preserving the blue noise of point sets | |
| Sharma et al. | Ground-state properties of electron-electron biwire systems | |
| CN107798353A (zh) | 一种间歇过程监控数据处理方法 | |
| CN116736345A (zh) | 用于mhgm多路径误差建模的半天球自适应格网划分方法 | |
| Carrasco et al. | Width and extremal height distributions of fluctuating interfaces with window boundary conditions | |
| Fang et al. | Percolation effects in the Fortuin-Kasteleyn Ising model on the complete graph | |
| CN110852605B (zh) | 一种基于信息效能的产品设计决策的确定方法及系统 | |
| Nagode et al. | An alternative perspective on the mixture estimation problem | |
| CN109684760B (zh) | 基于随机搜索算法的弹性矢量波场数值模拟方法及系统 | |
| Kelling et al. | Suppressing correlations in massively parallel simulations of lattice models | |
| Luo et al. | Multiple stiffness topology optimizations of continuum structures | |
| CN114186411A (zh) | 电动车辆的充电负荷预测方法、模型训练方法及装置 | |
| CN109492800B (zh) | 一种用于自动化仓库的车辆路径优化方法 | |
| Rossi et al. | A hybrid physics/data-driven logic to detect, classify, and predict anomalies and disruptions in tokamak plasmas | |
| Dequidt et al. | Time‐critical animation of deformable solids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |