[发明专利]一种蜷线位错原子结构的建模方法

说明书

本发明公开了一种蜷线位错原子结构的建模方法。该方法的主要内容包括，在给定包含晶体模型原子结构信息的文件的前提下，根据拟创建蜷线位错原子结构的Burgers矢量、位错线位置、滑移面和位错线形态的要求，使用编程语言提取文件中晶体模型的原子结构信息，自动计算出包含符合要求的蜷线位错原子结构的晶体模型的原子坐标，然后按计算机仿真技术能识别的文件格式输出数据到文件。本发明可方便快捷地在晶体内部指定位置直接创建指定方位、组态和形态的蜷线位错原子结构，为计算机仿真技术对蜷线位错的形态及行为的精准研究创造了有利条件。

技术领域

本发明涉及计算机仿真技术领域，尤其涉及一种蜷线位错原子结构的建模方法。

背景技术

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

发明内容

本发明所要解决的技术问题，就是提供一种方便快捷地创建蜷线位错原子结构的方法，本方法在给定包含晶体模型的原子结构信息的文件的前提下，根据拟创建蜷线位错原子结构的Burgers矢量、位错线位置、滑移面和位错线形态的要求，使用编程语言提取文件中晶体模型的原子结构信息，自动计算出包含符合要求的蜷线位错原子结构的晶体模型的原子坐标，然后按计算机仿真技术能识别的文件格式输出数据到文件。

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

设拟创建的蜷线位错原子结构的旋转轴线经过点P(x_p y_p z_p)，位错的Burgers矢量为[uvw]a，a为晶格常数，滑移面为(hkl)，螺距为h，阿基米德螺线系数为m。

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

步骤二：使用编程语言提取上述文件中的原子结构信息，平移坐标系使位错线经过点P，旋转坐标系使x轴的正向与[uvw]方向一致，y轴的正向与[hkl]方向一致，然后计算晶体模型内所有原子在新坐标系中的坐标值。

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

设：d=a(u²+v²+w²)^1/2，r=2πmz/h，x’=x-r cos(2πz/h)，y’=y-r sin(2πz/h)，

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

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

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

当y’c时，q=0，

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

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