CN109388892A - 一种位错分叉原子结构的建模方法 - Google Patents

Abstract

本发明公开了一种位错分叉原子结构的建模方法。该方法的主要内容包括，在给定包含晶体模型原子结构信息的文件的前提下，利用C/C++语言提取该文件中的晶体模型原子结构信息，根据拟构建位错分叉原子结构的要求自动计算出包含该要求的位错分叉原子结构的晶体模型的原子坐标，然后按分子动力学软件能识别的文件格式输出数据到文件。本发明可方便快捷地在晶体内部任意位置直接构建任意位向、任意组态的一个或多个位错分叉原子结构，为分子动力学及其它计算机仿真技术对位错反应行为的精准研究创造了有利条件。

Description

一种位错分叉原子结构的建模方法

技术领域

本发明涉及分子动力学仿真技术领域，尤其涉及一种位错分叉原子结构的建模方法。

背景技术

晶体宏观上的塑性变形实质是位错在外力作用下运动的结果。晶体的生长、形变强化、滞弹性、断裂、相变、晶体的电磁性能、晶体的光学性质、超导性以及其他许多物理、化学性质都与位错有重要的关联。因此，位错的研究无论对于科学研究还是实际应用都有着重要意义。实验上研究位错的方法有浸蚀法、缀饰法、透射电子显微分析法、X射线衍射分析法、场离子显微分析法等。这些实验技术被广泛地应用于分析研究位错的密度、分布和组态以及它们的运动和交互作用等。但是在原子尺度上的研究（比如位错芯的研究），分子动力学模拟具有重要的地位。各种位错原子结构的直接构建有利于分子动力学对位错行为更为精准的研究。本发明公开了一种位错分叉原子结构的建模方法，解决了分子动力学及其它计算机仿真研究中位错分叉原子结构直接建模的问题，本发明公开的方法可方便快捷地在晶体模型内部任何位置直接构建任意方位、任意组态的位错分叉原子结构，而且可在一个晶体模型内创建多个不同位向和组态的位错分叉原子结构。

发明内容

本发明所要解决的技术问题，就是提供一种方便快捷地构建位错分叉原子结构的方法，本方法在给定包含晶体模型的原子结构信息的文件的前提下，使用编程语言提取晶体模型的原子结构信息，根据拟构建位错分叉原子结构的要求自动计算出包含该要求的位错分叉原子结构的晶体模型的原子坐标，然后按分子动力学软件能识别的文件格式输出文件。

本发明采用的技术解决方案如下：

设拟构建的位错分叉原子结构的分叉点在P(x_p y_p z_p)处，三条位错均为刃型位错，Burgers矢量分别为[u₁ v₁ w₁]a、[u₂ v₂ w₂]a和[u₃ v₃ w₃]a，三条位错的公共滑移面为(h kl)，a为晶格常数，三条位错的Burgers矢量的关系为[u₁ v₁ w₁]a+[u₂ v₂ w₂]a=[u₃ v₃ w₃]a。

步骤一：准备包含晶体模型原子结构信息的文件。

步骤二：使用编程语言提取上述文件中的原子结构信息，将坐标系原点移动到拟构建的位错的分叉点P(x_p y_p z_p)处。旋转坐标系，使x轴在第一条位错的滑移方向[u₁ v₁w₁]上，y轴垂直于位错滑移面(h k l)，然后计算晶体模型内所有原子在新坐标系中的坐标值。

步骤三：设位错周围晶格扭曲明显的区域在x和y方向的范围为2b×2c。为构建位错分叉原子结构，晶体模型内的原子要进行相应的位移，本方法根据位错周围原子分布的特点，提出两次位移法。设在第一次位移中x方向的位移为q₁，y方向和z方向不发生位移，计算公式如下：

d₁=a(u₁ ²+v₁ ²+w₁ ²)^1/2；d₂=a(u₂ ²+v₂ ²+w₂ ²)^1/2；

当z≤0时，x’=x，

当z>0时，x’=x-z((u₁ ²+v₁ ²+w₁ ²)(u₃ ²+v₃ ²+w₃ ²)-(u₁u₃+v₁v₃+w₁w₃)²)^1/2/(u₁u₃+v₁v₃+w₁w₃)，

当-b≤x’≤b且0≤y≤c时，q₁=-(1-y/c)d₁x’/(4b)，

当x’>b且0≤y≤c时，q₁=-(1-y/c)d₁/4，

当x’<-b且0≤y≤c时，q₁=(1-y/c)d₁/4，

当y>c时，q₁=0，

当-b≤x’≤b且-c≤y<0时，q₁=d₁/2+(1+y/c)d₁x’/(4b)，

当x’>b且-c≤y<0时，q₁=d₁/2+(1+y/c)d₁/4，

当x’<-b且-c≤y<0时，q₁=d₁/2-(1+y/c)d₁/4，

当y<-c时，q₁=d₁/2。

步骤四：根据上面计算得到的每个原子的第一次位移值q₁，计算晶体模型内所有原子第一次位移后的坐标值。

步骤五：旋转坐标系，使x轴在第二条位错的滑移方向[u₂ v₂ w₂]上，y轴垂直于位错滑移面(h k l)，然后计算晶体模型内所有原子在新坐标系中的坐标值。

步骤六：进行第二次位移计算，设在第二次位移中x方向的位移为q₂，y方向和z方向不发生位移，计算公式如下：

当z≤0时，x’=x，

当z>0时，x’=x+z((u₂ ²+v₂ ²+w₂ ²)(u₃ ²+v₃ ²+w₃ ²)-(u₂u₃+v₂v₃+w₂w₃)²)^1/2/(u₂u₃+v₂v₃+w₂w₃)，

当-b≤x’≤b且0≤y≤c时，q₂=-(1-y/c)d₂x’/(4b)，

当x’>b且0≤y≤c时，q₂=-(1-y/c)d₂/4，

当x’<-b且0≤y≤c时，q₂=(1-y/c)d₂/4，

当y>c时，q₂=0，

当-b≤x’≤b且-c≤y<0时，q₂=d₂/2+(1+y/c)d₂x’/(4b)，

当x’>b且-c≤y<0时，q₂=d₂/2+(1+y/c)d₂/4，

当x’<-b且-c≤y<0时，q₂=d₂/2-(1+y/c)d₂/4，

当y<-c时，q₂=d₂/2。

步骤七：根据上面计算得到的每个原子的第二次位移值q₂，计算晶体模型内所有原子第二次位移后的坐标值，由此在指定的位置构建出相应位向的位错分叉原子结构。

步骤八：按照步骤二和步骤五的逆向移动坐标系，使坐标系恢复成原来的位向。

步骤九：按分子动力学软件能识别的格式输出数据到文件。

以上内容即为本发明公开的位错分叉原子结构的建模方法的主要内容。

本发明所公开的位错分叉原子结构的建模方法，可方便快捷地在晶体内部任何位置直接构建任意方位、任意组态的位错分叉原子结构，而且可在一个晶体内直接创建多个不同位向和组态的位错分叉原子结构，有利于分子动力学及其它计算机仿真技术对位错反应行为更为精准的研究。

附图说明

图1是本发明实施实例中创建的不含位错分叉结构的B2型NiAl金属间化合物40×40×40超胞的原子结构图。

图2是本发明实施实例中超胞内原子进行两次位移后，Ovito软件显示的位错分叉原子结构图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细说明，但本发明的实施方式不限于此，在不脱离本发明上述思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的范围内。

实施例：

本实例公开了一种位错分叉原子结构的建模方法。本实例在一个B2型NiAl金属间化合物40×40×40超胞内构建一个分叉点在超胞中心P(59.64, 59.64, 59.64)处、三条刃型位错分别为[001]、[100]和[101]、三条位错的公共滑移面为(010)的位错分叉原子结构。

步骤一：采用Materials Studio创建40×40×40超胞，如图1所示，然后以car的格式输出数据文件。

步骤二：提取上述文件中的原子结构信息，将坐标原点移到超胞的中心点P(59.64, 59.64, 59.64)（原坐标系的原点在超胞的一角，x轴在晶向[100]，y轴在晶向[010]，z轴在晶向[001]），旋转坐标系，使x轴沿着晶向[001]，y轴沿着晶向[010]，计算晶体模型内所有原子在新坐标系中的坐标值。

步骤三：设定位错分叉原子结构周围晶格扭曲明显的区域在x和y方向的范围为8a×8a，a为晶格常数，计算晶体超胞内所有原子在x方向的第一次位移q1，y方向和z方向不发生位移，主要程序代码如下：

for(i=0; i<total_no_atoms; i++){

a=2.882; b=4*a; c=4*a;

u1=0; v1=0; w1=1; u3=1; v3=0; w3=1; h=0; k=1; l=0;

d1=a* sqrt(u1*u1+v1*v1+w1*w1);

x1 = atoms[i].x[0];

y1 = atoms[i].x[1];

z1 = atoms[i].x[2];

if (z1<=0) {x2=x1; }

if(z1>0){x2=x1+z1*sqrt((u1*u1+v1*v1+w1*w1)*(u3*u3+v3*v3+w3*w3))-(u1*u3+v1*v3+w1*w3)*(u1*u3+v1*v3+w1*w3);}

if ((x2>=-b) &&(x2<=b) && (y1>=0) && (y1<=c)) {q1=-d1/4*x2/b*(1-y1/c);}

if ((x2>b) && (y1<=c) && (y1>=0)) {q1=-d1/4*(1-y1/c);}

if ((x2<-b) && (y1<=c) && (y1>=0)) {q1=d1/4*(1-y1/c);}

if ((y1>c)) {q1=0;}

if ((x2>=-b) &&(x2<=b) && (y1<0) && (y1>=-c)) {q1=d1/2+d1/4*x2/b*(1+y1/c);}

if ((x2>b) && (y1>=-c) && (y1<0)) {q1=d1/2+ d1/4*(1+y1/c);}

if ((x2<-b) && (y1>=-c) && (y1<0)) {q1=d1/2-d1/4*(1+y1/c);}

if ((y1<-c)) {q1=d1/2;}}。

步骤四：根据上面计算得到的每个原子的位移值q1，计算晶体超胞内所有原子位移后的坐标值，程序代码如下：

atoms[i].x[0]+=q1。

步骤五：旋转坐标系，使x轴在位错的滑移方向[100]上，y轴垂直于位错滑移面(010)，然后计算晶体模型内所有原子在新坐标系中的坐标值。

步骤六：进行第二次位移计算，设在第二次位移中x方向的位移为q2，y方向和z方向不发生位移，主要程序代码如下：

for(i=0; i<total_no_atoms;i++){

a=2.882; b=4*a; c=4*a;

u2=1; v2=0; w2=0; u3=1; v3=0; w3=1; h=0; k=1; l=0;

d2=a* sqrt(u2*u2+v2*v2+w2*w2);

x1 = atoms[i].x[0];

y1 = atoms[i].x[1];

z1 = atoms[i].x[2];

if (z1<=0) {x2=x1; }

if(z1>0){x2=x1-z1*sqrt((u2*u2+v2*v2+w2*w2)*(u3*u3+v3*v3+w3*w3))-(u2*u3+v2*v3+w2*w3)*(u2*u3+v2*v3+w2*w3);}

if ((x2>=-b) &&(x2<=b) && (y1>=0) && (y1<=c)) {q2=-d2/4*x2/b*(1-y1/c);}

if ((x2>b) && (y1<=c) && (y1>=0)) {q2=-d2/4*(1-y1/c);}

if ((x2<-b) && (y1<=c) && (y1>=0)) {q2=d2/4*(1-y1/c);}

if ((y1>c)) {q2=0;}

if ((x2>=-b) &&(x2<=b) && (y1<0) && (y1>=-c)) {q2=d2/2+d2/4*x2/b*(1+y1/c);}

if ((x2>b) && (y1>=-c) && (y1<0)) {q2=d2/2+ d2/4*(1+y1/c);}

if ((x2<-b) && (y1>=-c) && (y1<0)) {q2=d2/2-d2/4*(1+y1/c);}

if ((y1<-c)) {q2=d2/2;}}。

步骤七：根据上面计算得到的每个原子的位移值q2，计算晶体超胞内所有原子位移后的坐标值，由此构建出相应位向的位错分叉原子结构，程序代码如下：

atoms[i].x[0]+=q2。

步骤八：按照步骤二和步骤五的逆向移动坐标系，使坐标系恢复成原来的位向。

步骤九：按分子动力学软件能识别的格式输出数据到文件。

由此，在一个B2型NiAl金属间化合物40×40×40超胞内构建一个分叉点在超胞中心、三条刃型位错分别为[001]、[100]和[101]、滑移面为(010)的位错分叉原子结构。图2为使用Ovito软件显示的位错识别图。

Claims (1)

1.一种位错分叉原子结构的建模方法，其特征在于该方法在给定包含晶体模型的原子结构信息的文件的前提下，使用编程语言提取晶体模型的原子结构信息，根据拟构建位错分叉原子结构的要求自动计算出包含该要求的位错分叉原子结构的晶体模型的原子坐标，然后按分子动力学软件能识别的文件格式输出文件，主要步骤如下：

步骤一：准备包含晶体模型原子结构信息的文件；

步骤二：使用编程语言提取上述文件中的原子结构信息，将坐标系原点移动到拟构建的位错的分叉点P(x_p y_p z_p)处；旋转坐标系，使x轴在第一条位错的滑移方向[u₁ v₁ w₁]上，y轴垂直于位错滑移面(h k l)，然后计算晶体模型内所有原子在新坐标系中的坐标值；

步骤三：设位错周围晶格扭曲明显的区域在x和y方向的范围为2b×2c，a为晶格常数，第三条位错的滑移方向为[u₃ v₃ w₃]；本方法根据位错周围原子分布的特点，提出两次位移法；设在第一次位移中x方向的位移为q₁，y方向和z方向不发生位移，计算公式如下：

d₁=a(u₁ ²+v₁ ²+w₁ ²)^1/2；d₂=a(u₂ ²+v₂ ²+w₂ ²)^1/2，

当z≤0时，x’=x，

当z>0时，x’=x-z((u₁ ²+v₁ ²+w₁ ²)(u₃ ²+v₃ ²+w₃ ²)-(u₁u₃+v₁v₃+w₁w₃)²)^1/2/(u₁u₃+v₁v₃+w₁w₃)，

当-b≤x’≤b且0≤y≤c时，q₁=-(1-y/c)d₁x’/(4b)，

当x’>b且0≤y≤c时，q₁=-(1-y/c)d₁/4，

当x’<-b且0≤y≤c时，q₁=(1-y/c)d₁/4，

当y>c时，q₁=0，

当-b≤x’≤b且-c≤y<0时，q₁=d₁/2+(1+y/c)d₁x’/(4b)，

当x’>b且-c≤y<0时，q₁=d₁/2+(1+y/c)d₁/4，

当x’<-b且-c≤y<0时，q₁=d₁/2-(1+y/c)d₁/4，

当y<-c时，q₁=d₁/2；

步骤四：根据上面计算得到的每个原子的第一次位移值q₁，计算晶体模型内所有原子第一次位移后的坐标值；

步骤五：旋转坐标系，使x轴在第二条位错的滑移方向[u₂ v₂ w₂]上，y轴垂直于位错滑移面(h k l)，然后计算晶体模型内所有原子在新坐标系中的坐标值；

步骤六：进行第二次位移计算，设在第二次位移中x方向的位移为q₂，y方向和z方向不发生位移，计算公式如下：

当z≤0时，x’=x，

当z>0时，x’=x+z((u₂ ²+v₂ ²+w₂ ²)(u₃ ²+v₃ ²+w₃ ²)-(u₂u₃+v₂v₃+w₂w₃)²)^1/2/(u₂u₃+v₂v₃+w₂w₃)，

当-b≤x’≤b且0≤y≤c时，q₂=-(1-y/c)d₂x’/(4b)，

当x’>b且0≤y≤c时，q₂=-(1-y/c)d₂/4，

当x’<-b且0≤y≤c时，q₂=(1-y/c)d₂/4，

当y>c时，q₂=0，

当-b≤x’≤b且-c≤y<0时，q₂=d₂/2+(1+y/c)d₂x’/(4b)，

当x’>b且-c≤y<0时，q₂=d₂/2+(1+y/c)d₂/4，

当x’<-b且-c≤y<0时，q₂=d₂/2-(1+y/c)d₂/4，

当y<-c时，q₂=d₂/2；

步骤七：根据上面计算得到的每个原子的第二次位移值q₂，计算晶体模型内所有原子第二次位移后的坐标值，由此在指定的位置构建出相应位向的位错分叉原子结构；

步骤八：按照步骤二和步骤五的逆向移动坐标系，使坐标系恢复成原来的位向；

步骤九：按分子动力学软件能识别的格式输出数据到文件。