CN110120248A - 模拟纳米晶金属累积离位损伤的方法 - Google Patents
模拟纳米晶金属累积离位损伤的方法 Download PDFInfo
- Publication number
- CN110120248A CN110120248A CN201910276042.0A CN201910276042A CN110120248A CN 110120248 A CN110120248 A CN 110120248A CN 201910276042 A CN201910276042 A CN 201910276042A CN 110120248 A CN110120248 A CN 110120248A
- Authority
- CN
- China
- Prior art keywords
- defect
- rate
- event
- atom
- okmc
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Testing Or Measuring Of Semiconductors Or The Like (AREA)
Abstract
本发明公开了一种模拟纳米晶金属累积离位损伤的方法,步骤包括基本原子过程MD参数化;基本原子过程OKMC速率列表创建;OKMC调用相应速率,选择、执行事件;表征OKMC粗粒化结构;采样不同剂量下的微结构,并用MD弛豫之;将观察到的原子过程速率化;更新基本原子过程,重新运行OKMC模块。本发明优点在于通过结合原子尺度MD模拟方法和粗粒化OKMC方法,既考虑了基本原子过程,又能实现大尺度模拟,可与实验剂量率相比拟,同时考虑了基本过程累积作用下可能出现的新原子过程;同现有的模拟方法相比,避免了因MD时间尺度有限带来的辐照剂量率异常高的问题,又避免了OKMC在结构弛豫上的不足,兼顾了MD在结构弛豫和OKMC在长时间尺度模拟上的综合优势。
Description
技术领域
本发明涉及核材料辐照损伤模拟技术领域,具体涉及模拟纳米晶金属累积离位损伤的方法。
背景技术
高能粒子(如中子、离子)辐照材料时会诱导生成空位(vacancy,V)、自间隙原子(sel-interstitial atom,SIA)等辐照缺陷,引起辐照损伤。纳米结构材料通常具有较好的抗辐照损伤性能,这与材料中大量的晶界对辐照缺陷的捕获、进而促进其复合有关。目前人们研究纳米结构金属材料中辐照缺陷与界面相互作用机理时,大多关注辐照缺陷的基本原子过程,如扩散、复合过程。
而在实际的辐照环境下,纳米晶金属往往受到以一定剂量率辐照的累积辐照损伤。有部分学者采用分子动力学(molecular dynamics,MD)模拟纳米晶累积离位损伤,虽能达到几个dpa的辐照剂量,然而由于分子动力学固有的较短时间尺度限制,模拟中的剂量率往往高于实验中的剂量率数个数量级,因此,在模拟中观察的现象可能不具有物理意义。如果仅仅考虑基本的原子过程,采用粗粒化模拟技术(如实体动力学蒙特卡洛,objectkinetic Monte Carlo,OKMC)直接模拟累积离位损伤,则由于粗粒化技术中失去了原子结构信息,离位损伤累积到一定程度后可能带来新的物理过程,从而使得高剂量下的微结构演化预测失效。
发明内容
本发明的目的在于提供一种模拟纳米晶金属累积离位损伤的方法,该方法既能达到实验剂量率和辐照剂量,又能根据剂量补充、校正已有原子过程的模拟技术,以实现准确模拟纳米晶材料中累积离位损伤微结构跨时间、空间尺度的演化过程。
本发明通过以下技术方案来实现上述目的:
模拟纳米晶金属累积离位损伤的方法,步骤包括
S1:基本原子过程MD参数化
采用MD计算或者收集点缺陷及其团簇扩散、聚合、复合、溶解、发射过程的能量学、动力学和作用范围参数,并根据缺陷形成能性质确定点缺陷在晶界内部的占据位置;
S2:基本原子过程OKMC速率列表创建
根据物理图像,把MD计算的参数加以简化,计算一定温度下基本过程速率和缺陷产生速率,构建速率列表;
S3:OKMC调用相应速率,选择、执行事件
根据不同事件、过程的速率,采用轮盘赌方法选择事件,且某个事件被选中的几率与该事件速率成正比,再执行事件;
S4:表征OKMC粗粒化结构
对于OKMC模拟得到的不同辐照剂量(dose)下的微结构进行分析,得到块体区域、晶界附近和晶界内部不同类型缺陷浓度(C),建立C-dose关系;分析距晶界不同位置处缺陷类型与浓度,建立缺陷浓度与缺陷-晶界距离(d)关系:C(dose)-d;
S5:采样不同剂量下的微结构,并用MD弛豫之
将S4中得到的随辐照剂量变化的微结构离散化采样,建立合适尺寸的原子模型,并再现一定剂量下的微结构,保证粗粒化模型和原子结构中相应缺陷浓度相同、缺陷分布一致,模型建好以后,先在0K下静态弛豫模型,然后在一定温度采用MD再次弛豫模型,观察弛豫过程中体系结构的变化;
S6:将观察到的原子过程速率化
计算S5中新出现的原子过程发生的速率,确定新过程发生的临界剂量;
S7:更新基本原子过程,重新运行OKMC模块
将随辐照剂量发生的新原子过程补充进基本原子过程集合,更新事件速率列表,重新运行OKMC模块。
进一步改进在于,步骤S2中,所述简化包括将晶界附近缺陷扩散能垒设置为块体中相应值的一半;将间隙发射诱导的复合分为一定区域内的自发复合和低能垒复合两部分,自发复合区域外缺陷复合能垒正比于缺陷与复合区域的距离,收敛于块体中孤立缺陷扩散能垒;设置不动缺陷团簇尺寸上限;同时简化缺陷之间、缺陷与晶界之间作用半径。
进一步改进在于,步骤S2中,所述过程速率的计算方法为:
将一个缺陷对象视作球形,其半径r(n)设置为核心区半径(R)加上影响区壳层半径(rshell),公式如下:
式中,n为团簇Vn/SIAn中点缺陷数,a0是晶格常数;
考虑到块体中自间隙原子较大的应力场,其rshell设置为较大值;以下是块体和晶界区域缺陷半径计算公式:
一个能垒为Em的过程,其速率为:
r=v0exp[-Em/(kBT)], (6)
式中,指前因子v0取1012/s,玻尔兹曼常数kB取8.517×10-5eV/K,T为温度。
进一步改进在于,步骤S2中,所述缺陷产生速率计算方法为:
对于给定的辐照剂量率τ、晶粒尺寸L,点缺陷自间隙原子原子和空位产生速率为:
式中,计算模型体积Vmodel为L3,一个原子的体积Vatom为
进一步改进在于,步骤S3中,所述执行事件包括:
扩散过程事件:处于块体区域的缺陷对象三维随机行走,处于晶界附近的缺陷一维扩散,而处于晶界内部的缺陷在晶界面上二维随机扩散,且每次扩散,缺陷沿着一个坐标轴方向运动晶格常数一半的距离;
聚合过程事件:处于块体中的空位团簇Vn1和Vn2聚合,只考虑一个小团簇被另一个团簇捕获的情况,即限制n1<5或n2<5,自间隙团簇之间的聚合不设尺寸限制,对于弛豫晶界内部的缺陷聚合反应,只考虑点缺陷被另一个缺陷捕获的情形;聚合反应执行时,反应后的对象坐标与尺寸较大的一个一致,尺寸变为n1+n2;
复合过程事件:执行复合事件时,保留尺寸较大的一个,剩余对象尺寸为|n1-n2|;
溶解过程事件:每次只有一个V/SIA从Vn/SIAn中溶解出来,溶解后,将溶解出来的一个V/SIA放置在捕获半径外的球面上;
发射过程事件:每次只有一个V/SI从晶界发射出来,发射后的V/SIA放在晶界对V/SIA捕获区域的之外的格点处;
缺陷产生事件:在体系中随机产生一对空位和自间隙原子,产生三个[0,1]之间均分分布的随机数r1、r2、r3,则缺陷坐标为[L×r1,L×r2,L×r3]。
本发明的有益效果在于:本发明通过结合原子尺度MD模拟方法和粗粒化OKMC方法,既考虑了基本原子过程,又能实现大尺度模拟,可与实验剂量率相比拟,同时考虑了基本过程累积作用下可能出现的新原子过程;同现有的模拟方法相比,避免了因MD时间尺度有限带来的辐照剂量率异常高的问题,又避免了OKMC在结构弛豫上的不足,兼顾了MD在结构弛豫和OKMC在长时间尺度模拟上的综合优势。
附图说明
图1为本发明提出的一种模拟纳米晶金属累积离位损伤的方法流程图;
图2为本发明累积离位损伤模拟方法示意图;
图3为本发明的算法流程图;
具体实施方式
下面结合附图对本申请作进一步详细描述,有必要在此指出的是,以下具体实施方式只用于对本申请进行进一步的说明,不能理解为对本申请保护范围的限制,该领域的技术人员可以根据上述申请内容对本申请作出一些非本质的改进和调整。
结合图1和图2所示,模拟纳米晶金属累积离位损伤的方法,步骤包括
S1:基本原子过程MD参数化。本步骤中,针对辐照缺陷与晶界基本相互作用过程:V、SIA及其团簇(Vn/SIAn)的扩散、聚合、复合、溶解、发射过程,采用MD计算或者文献收集整理这些作用过程的能量学、动力学和作用范围参数,如块体区域、晶界附近和晶界内部缺陷的扩散能垒、结合能、偏聚能垒和复合能垒。同时根据缺陷形成能性质确定点缺陷在晶界内部的占据位置。
S2:基本原子过程OKMC速率列表创建。本步骤中,根据物理图像(如晶界增强的缺陷扩散、复合),把MD计算的参数加以简化。如将晶界附近缺陷扩散能垒设置为块体中相应值的一半;将间隙发射诱导的复合分为一定区域内的自发复合和低能垒复合两部分,自发复合区域外缺陷复合能垒正比于缺陷与复合区域的距离,收敛于块体中孤立缺陷扩散能垒;设置不动缺陷团簇尺寸上限。同时简化缺陷之间、缺陷与晶界之间作用半径。计算一定温度下基本过程速率和缺陷产生速率。将一个缺陷对象视作球形,其半径r(n)设置为核心区半径(R)加上影响区壳层半径(rshell),公式如下:
其中n为团簇Vn/SIAn中点缺陷数,a0是晶格常数。考虑到块体中自间隙原子较大的应力场,其rshell设置为较大值。以下是块体和晶界区域缺陷半径计算公式:
一个能垒为Em的过程,其速率为:
r=v0exp[-Em/(kBT)], (6)
其中指前因子v0通常取1012/s,玻尔兹曼常数kB取8.517×10-5eV/K,T为温度。
缺陷产生速率计算方法如下。对于给定的辐照剂量率τ、晶粒尺寸L,点缺陷自间隙原子原子和空位产生速率为:
其中计算模型体积Vmodel为L3,一个原子的体积Vatom为
S3:OKMC调用相应速率,选择、执行事件。本步骤中根据不同事件、过程的速率,采用轮盘赌方法选择事件,某个事件被选中的几率与该事件速率成正比。对于扩散事件,处于块体区域的缺陷对象三维随机行走,处于晶界附近的缺陷一维扩散,而处于晶界内部的缺陷在晶界面上二维随机扩散;每次扩散,缺陷沿着一个坐标轴方向运动晶格常数一半的距离。对于聚合事件,处于块体中的空位团簇Vn1和Vn2聚合,只考虑一个小团簇被另一个团簇捕获的情况,即限制n1<5或n2<5;考虑到自间隙团簇较高的活性,自间隙团簇之间的聚合不设尺寸限制;对于弛豫晶界内部的缺陷聚合反应,只考虑点缺陷被另一个缺陷捕获的情形;聚合反应执行时,反应后的对象坐标与尺寸较大的一个一致,尺寸变为n1+n2。对于复合事件,保留尺寸较大的一个,剩余对象尺寸为|n1-n2|。对于溶解事件:每次只有一个V/SIA从Vn/SIAn中溶解出来,溶解后,将溶解出来的一个V/SIA放置在捕获半径外的球面上。对于发射事件:每次只有一个V/SI从晶界发射出来,发射后的V/SIA放在晶界对V/SIA捕获区域的之外的格点处。
另外,对于缺陷产生事件,在体系中随机产生一对空位和自间隙原子。产生三个[0,1]之间均匀分布的随机数r1、r2、r3,则缺陷坐标为[L×r1,L×r2,L×r3]。
S4:表征OKMC粗粒化结构。此步骤中,对于OKMC模拟得到的不同辐照剂量(dose)下的微结构进行分析,得到块体区域、晶界附近和晶界内部不同类型缺陷浓度(C),建立C-dose关系。分析距晶界不同位置处缺陷类型与浓度,建立缺陷浓度与缺陷-晶界距离(d)关系:C(dose)-d。
S5:采样不同剂量下的微结构,并用MD弛豫之。本步骤中将S4中得到的随辐照剂量变化的微结构离散化采样,建立合适尺寸的原子模型,并再现一定剂量下的微结构,保证粗粒化模型和原子结构中相应缺陷浓度相同、缺陷分布一致。模型建好以后,先在0K下静态弛豫模型,然后在一定温度采用MD再次弛豫模型。观察弛豫过程中体系结构的变化。
S6:将观察到的原子过程速率化。此步骤中计算S5中新出现的原子过程发生的速率;确定新过程发生的临界剂量。
S7:更新基本原子过程,重新运行OKMC模块。本步骤中将随辐照剂量发生的新原子过程补充进基本原子过程集合,更新事件速率列表,重新运行OKMC模块。
参照图3,本发明算法流程如下:
(1)计算或收集点缺陷及其团簇扩散、偏聚、复合、聚合、溶解能垒数据,缺陷与缺陷、缺陷与晶界之间结合能数据,以及晶界作用范围数据、缺陷作用半径数据。
(2)设置模拟温度、剂量率、模型尺寸等参数;计算基本原子过程发生速率。
(3)根据缺陷所在位置,调用不同的速率,构建速率列表。
(4)采用轮盘赌的方式筛选要执行的事件。
(5)执行选择的事件。
(6)更新计时器和缺陷位置,输出剂量、缺陷对象位置、缺陷尺寸等表征量。
(7)判断辐照剂量是否超过设定值,否则转到步骤(3),是则计算不同剂量下辐照缺陷浓度及其分布。
(8)生成相应浓度的原子模型,采用有限温度MD弛豫模型。
(9)观察累积一定剂量后是否出现新过程。否则结束,是则计算新原子过程发生速率,转到步骤(3)。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。
Claims (5)
1.模拟纳米晶金属累积离位损伤的方法,其特征在于:步骤包括
S1:基本原子过程MD参数化
采用MD计算或者收集点缺陷及其团簇扩散、聚合、复合、溶解、发射过程的能量学、动力学和作用范围参数,并根据缺陷形成能性质确定点缺陷在晶界内部的占据位置;
S2:基本原子过程OKMC速率列表创建
根据物理图像,把MD计算的参数加以简化,计算一定温度下基本过程速率和缺陷产生速率,构建速率列表;
S3:OKMC调用相应速率,选择、执行事件
根据不同事件、过程的速率,采用轮盘赌方法选择事件,且某个事件被选中的几率与该事件速率成正比,再执行事件;
S4:表征OKMC粗粒化结构
对于OKMC模拟得到的不同辐照剂量(dose)下的微结构进行分析,得到块体区域、晶界附近和晶界内部不同类型缺陷浓度(C),建立C-dose关系;分析距晶界不同位置处缺陷类型与浓度,建立缺陷浓度与缺陷-晶界距离(d)关系:C(dose)-d;
S5:采样不同剂量下的微结构,并用MD弛豫之
将S4中得到的随辐照剂量变化的微结构离散化采样,建立合适尺寸的原子模型,并再现一定剂量下的微结构,保证粗粒化模型和原子结构中相应缺陷浓度相同、缺陷分布一致,模型建好以后,先在0K下静态弛豫模型,然后在一定温度采用MD再次弛豫模型,观察弛豫过程中体系结构的变化;
S6:将观察到的原子过程速率化
计算S5中新出现的原子过程发生的速率,确定新过程发生的临界剂量;
S7:更新基本原子过程,重新运行OKMC模块
将随辐照剂量发生的新原子过程补充进基本原子过程集合,更新事件速率列表,重新运行OKMC模块。
2.根据权利要求1所述的模拟纳米晶金属累积离位损伤的方法,其特征在于:步骤S2中,所述简化包括将晶界附近缺陷扩散能垒设置为块体中相应值的一半;将间隙发射诱导的复合分为一定区域内的自发复合和低能垒复合两部分,自发复合区域外缺陷复合能垒正比于缺陷与复合区域的距离,收敛于块体中孤立缺陷扩散能垒;设置不动缺陷团簇尺寸上限;同时简化缺陷之间、缺陷与晶界之间作用半径。
3.根据权利要求1所述的模拟纳米晶金属累积离位损伤的方法,其特征在于:步骤S2中,所述过程速率的计算方法为:
将一个缺陷对象视作球形,其半径r(n)设置为核心区半径(R)加上影响区壳层半径(rshell),公式如下:
式中,n为团簇Vn/SIAn中点缺陷数,a0是晶格常数;
考虑到块体中自间隙原子较大的应力场,其rshell设置为较大值;以下是块体和晶界区域缺陷半径计算公式:
一个能垒为Em的过程,其速率为:
r=v0exp[-Em/(kBT)], (6)
式中,指前因子v0取1012/s,玻尔兹曼常数kB取8.517×10-5eV/K,T为温度。
4.根据权利要求1所述的模拟纳米晶金属累积离位损伤的方法,其特征在于:步骤S2中,所述缺陷产生速率计算方法为:
对于给定的辐照剂量率τ、晶粒尺寸L,点缺陷自间隙原子原子和空位产生速率为:
式中,计算模型体积Vmodel为L3,一个原子的体积Vatom为
5.根据权利要求1所述的模拟纳米晶金属累积离位损伤的方法,其特征在于:步骤S3中,所述执行事件包括:
扩散过程事件:处于块体区域的缺陷对象三维随机行走,处于晶界附近的缺陷一维扩散,而处于晶界内部的缺陷在晶界面上二维随机扩散,且每次扩散,缺陷沿着一个坐标轴方向运动晶格常数一半的距离;
聚合过程事件:处于块体中的空位团簇Vn1和Vn2聚合,只考虑一个小团簇被另一个团簇捕获的情况,即限制n1<5或n2<5,自间隙团簇之间的聚合不设尺寸限制,对于弛豫晶界内部的缺陷聚合反应,只考虑点缺陷被另一个缺陷捕获的情形;聚合反应执行时,反应后的对象坐标与尺寸较大的一个一致,尺寸变为n1+n2;
复合过程事件:执行复合事件时,保留尺寸较大的一个,剩余对象尺寸为|n1-n2|;
溶解过程事件:每次只有一个V/SIA从Vn/SIAn中溶解出来,溶解后,将溶解出来的一个V/SIA放置在捕获半径外的球面上;
发射过程事件:每次只有一个V/SI从晶界发射出来,发射后的V/SIA放在晶界对V/SIA捕获区域的之外的格点处;
缺陷产生事件:在体系中随机产生一对空位和自间隙原子,产生三个[0,1]之间均匀分布的随机数r1、r2、r3,则缺陷坐标为[L×r1,L×r2,L×r3]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276042.0A CN110120248B (zh) | 2019-04-08 | 2019-04-08 | 模拟纳米晶金属累积离位损伤的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276042.0A CN110120248B (zh) | 2019-04-08 | 2019-04-08 | 模拟纳米晶金属累积离位损伤的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110120248A true CN110120248A (zh) | 2019-08-13 |
CN110120248B CN110120248B (zh) | 2020-12-25 |
Family
ID=67520894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910276042.0A Active CN110120248B (zh) | 2019-04-08 | 2019-04-08 | 模拟纳米晶金属累积离位损伤的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110120248B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114616362A (zh) * | 2019-10-29 | 2022-06-10 | 胜高股份有限公司 | 点缺陷模拟器、点缺陷模拟程序、点缺陷模拟方法、单晶硅的制造方法及单晶提拉装置 |
CN115394364A (zh) * | 2022-07-12 | 2022-11-25 | 北京应用物理与计算数学研究所 | 原子动力学蒙特卡洛模拟的原子指纹计算方法及装置 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853780A (zh) * | 2009-03-30 | 2010-10-06 | 索尼公司 | 离子辐射损伤预测方法和仿真器以及离子辐射设备和方法 |
CN105787227A (zh) * | 2016-05-11 | 2016-07-20 | 中国科学院近代物理研究所 | 结构材料辐照损伤的多gpu分子动力学模拟方法 |
CN107195345A (zh) * | 2017-06-09 | 2017-09-22 | 中国科学院合肥物质科学研究院 | 全能谱中子辐照损伤精确模拟系统及其算法 |
CN107530638A (zh) * | 2015-02-24 | 2018-01-02 | 渥太华大学 | 在受控击穿期间通过激光照射而在膜上定位纳米孔制造 |
CN108052745A (zh) * | 2017-12-15 | 2018-05-18 | 中国科学院合肥物质科学研究院 | 一种模拟辐照缺陷与晶界协同演化的方法 |
CN108090279A (zh) * | 2017-12-15 | 2018-05-29 | 中国科学院合肥物质科学研究院 | 基于杂化lkmc和okmc模拟纳米结构核材料辐照损伤的方法 |
CN108493108A (zh) * | 2018-05-16 | 2018-09-04 | 江苏润奥电子制造股份有限公司 | 一种高压快速软恢复二极管的制造方法 |
CN108509724A (zh) * | 2018-04-03 | 2018-09-07 | 嘉兴学院 | 一种多尺度模拟纳米颗粒多相流体特性的方法 |
CN108595842A (zh) * | 2018-04-25 | 2018-09-28 | 中国科学院合肥物质科学研究院 | 模拟粒子辐照损伤的时间步长优化方法 |
CN108866496A (zh) * | 2017-11-28 | 2018-11-23 | 中国科学院金属研究所 | 抗辐照损伤金属纳米晶/碳纳米管复合材料及其制备方法 |
CN108959709A (zh) * | 2018-06-04 | 2018-12-07 | 中国科学院合肥物质科学研究院 | 基于缺陷性质和多尺度模拟的晶界结构搜索方法 |
CN109192250A (zh) * | 2018-08-01 | 2019-01-11 | 华东理工大学 | 一种多相催化中克服表面物种快速迁移的加速模拟方法 |
-
2019
- 2019-04-08 CN CN201910276042.0A patent/CN110120248B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853780A (zh) * | 2009-03-30 | 2010-10-06 | 索尼公司 | 离子辐射损伤预测方法和仿真器以及离子辐射设备和方法 |
CN107530638A (zh) * | 2015-02-24 | 2018-01-02 | 渥太华大学 | 在受控击穿期间通过激光照射而在膜上定位纳米孔制造 |
CN105787227A (zh) * | 2016-05-11 | 2016-07-20 | 中国科学院近代物理研究所 | 结构材料辐照损伤的多gpu分子动力学模拟方法 |
CN107195345A (zh) * | 2017-06-09 | 2017-09-22 | 中国科学院合肥物质科学研究院 | 全能谱中子辐照损伤精确模拟系统及其算法 |
CN108866496A (zh) * | 2017-11-28 | 2018-11-23 | 中国科学院金属研究所 | 抗辐照损伤金属纳米晶/碳纳米管复合材料及其制备方法 |
CN108052745A (zh) * | 2017-12-15 | 2018-05-18 | 中国科学院合肥物质科学研究院 | 一种模拟辐照缺陷与晶界协同演化的方法 |
CN108090279A (zh) * | 2017-12-15 | 2018-05-29 | 中国科学院合肥物质科学研究院 | 基于杂化lkmc和okmc模拟纳米结构核材料辐照损伤的方法 |
CN108509724A (zh) * | 2018-04-03 | 2018-09-07 | 嘉兴学院 | 一种多尺度模拟纳米颗粒多相流体特性的方法 |
CN108595842A (zh) * | 2018-04-25 | 2018-09-28 | 中国科学院合肥物质科学研究院 | 模拟粒子辐照损伤的时间步长优化方法 |
CN108493108A (zh) * | 2018-05-16 | 2018-09-04 | 江苏润奥电子制造股份有限公司 | 一种高压快速软恢复二极管的制造方法 |
CN108959709A (zh) * | 2018-06-04 | 2018-12-07 | 中国科学院合肥物质科学研究院 | 基于缺陷性质和多尺度模拟的晶界结构搜索方法 |
CN109192250A (zh) * | 2018-08-01 | 2019-01-11 | 华东理工大学 | 一种多相催化中克服表面物种快速迁移的加速模拟方法 |
Non-Patent Citations (9)
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114616362A (zh) * | 2019-10-29 | 2022-06-10 | 胜高股份有限公司 | 点缺陷模拟器、点缺陷模拟程序、点缺陷模拟方法、单晶硅的制造方法及单晶提拉装置 |
CN114616362B (zh) * | 2019-10-29 | 2023-06-20 | 胜高股份有限公司 | 点缺陷模拟器、点缺陷模拟程序、点缺陷模拟方法、单晶硅的制造方法及单晶提拉装置 |
CN115394364A (zh) * | 2022-07-12 | 2022-11-25 | 北京应用物理与计算数学研究所 | 原子动力学蒙特卡洛模拟的原子指纹计算方法及装置 |
CN115394364B (zh) * | 2022-07-12 | 2024-02-02 | 北京应用物理与计算数学研究所 | 原子动力学蒙特卡洛模拟的原子指纹计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110120248B (zh) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brax et al. | Nonlinear structure formation with the environmentally dependent dilaton | |
Ichikawa | Systematic investigations of deep sub-barrier fusion reactions using an adiabatic approach | |
CN108052745B (zh) | 一种模拟辐照缺陷与晶界协同演化的方法 | |
Maxwell et al. | Atomistic simulation study of clustering and evolution of irradiation-induced defects in zirconium | |
CN110120248A (zh) | 模拟纳米晶金属累积离位损伤的方法 | |
CN108090279B (zh) | 基于杂化lkmc和okmc模拟纳米结构核材料辐照损伤的方法 | |
Dai et al. | A mechanism for basal vacancy loop formation in zirconium | |
Martin et al. | Molecular dynamics study of the bulk temperature effect on primary radiation damage in uranium dioxide | |
Mocanu et al. | Self-organized criticality in boson clouds around black holes | |
Qiu et al. | On the interactions of interstitial helium atom with helium bubble in tungsten: A molecular dynamics study | |
Warwick et al. | Dislocation density transients and saturation in irradiated zirconium | |
Van Brutzel et al. | Molecular dynamics simulations of microstructural evolution of irradiated (U, Pu) O2 studied via simulated XRD patterns | |
Sakaël et al. | Modelling of zirconium growth under irradiation and annealing conditions | |
Verma et al. | A new spatially resolved model for defects and fission gas bubbles interaction at the mesoscale | |
Kononov et al. | First-principles simulation of light-ion microscopy of graphene | |
Zheng et al. | An effective genetic algorithm to VDA with discontinuous “on-off” switches | |
Maillard et al. | Why a steady state void size distribution in irradiated UO2? A modeling approach | |
Hwa | Recognizing critical behavior amidst minijets at the large hadron collider | |
Fikar et al. | Stability of small vacancy clusters in tungsten by molecular dynamics | |
Perez-Ponce et al. | Implementing Geant4 on GPU for medical applications | |
Boev et al. | Molecular dynamic simulations of the interaction of interstitial atoms with vacancy complexes in V and V-4Ti | |
Xu et al. | Long-term behavior of vacancy defects in Pu-Ga alloy: Effects of temperature and Ga concentration | |
Torres et al. | Stochastic core–shell (SCS) approximation for accelerated atomistic modeling of irradiation-induced damage in materials | |
Viner et al. | Using a coupled dispersion model to estimate depletion of a tritium oxide plume by a forest | |
BARANI | Mechanistic modeling of fission gas behavior in conventional and advanced nuclear fuel |
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 |