lammps实例一硅的晶格常数和体弹

北京中科医院是真是假 http://m.39.net/baidianfeng/a_5154118.html

硅的晶格常数和体弹模量的计算

一、平衡晶格常数和内聚能自然条件下硅为金刚石结构(dc)。计算模拟时，我们可以假定它为各种结构，fcc,bcc,sc,dc.可以预测，模拟的dc结构的硅的体系能量最低，也即最稳定。下面我们将运用LAMMPS来对硅的各种结构进行模拟。定义晶格能量为Φ，数密度为ρ：

其中Epot为势能，N为体系总原子数，V为体系的体积。选取Stillinger-Weber(SW)，以下面命令执行lammps运算：

lmp_serialin.Silicon–logdc.log

其中，lmp_serial为lammps命令；””符号为读取符；in.Silicon为输入文件，里面包含运算所需要的各种数据和命令；-log指定输出文件的名称。

可以看到屏幕上显示出lammps运行的信息。这个计算量很小，所以很快就结束。接下来以如下命令来查看计算得到的数据：

grep^

dc.log

grep是linux中一个很重要的命令，用来搜索文本，读取匹配的行并打印出来。这里是搜索dc.log文件，将

开头的行打印出来。如下：

晶格参数为5.埃,数密度为0.，每个原子的能量为-4.eV.

下面具体来看刚才给的输入文件，in.Silicon.

dc.log文件中有原子总数的信息，

每个金刚石晶胞中有8个原子，8x3^3，所以是个原子。如下给出各种结构下的体系的原子数：

下图是计算模拟得出的各种结构下的数密度与每个原子能量的关系图。横坐标为数密度，以金刚石为例，ρ=8/5.^3=0.，也即我们直接通过grep命令得到的第二项值；纵坐标为每个原子的能量，为第三项值。

金刚石之外，还需计算其他结构。只需对in.Silicon做稍微改动：

首先，将in.Silicon复制成in.fcc:

cpin.Siliconin.fcc

然后编辑in.fcc

geditin.fcc

改动如下几项：

然后如下命令执行：

相应的，如下命令查看log文件中的数据:

以同样方法编辑in.bcc,in.sc，计算不同晶格参数时的体系能量值，并绘制下图：

下图更为细致地画出金刚石结构中，不同晶格参数所对应的内聚能。内聚能(cohesiveenergyEcoh)的定义是，最小的晶格能。由图可以得到，平衡晶格常数为a0=5.?，内聚能为Ecoh=?4.eV.

二、体弹模量

我们同时可以从晶格能曲线在最低处得到体弹模量的信息。体弹模量定义为：

V和P分别为晶胞的体积和体系的压强。我们已经得到了内聚能与晶格参数的函数关系，对于立方晶胞而言，

其中a0为平衡晶格常数，M为体积为V=a3的晶胞中的原子数目。

从多项式拟合，可以得到

此例中为3.eV/?2(a0=5.?).由上面公式，计算得到硅的体弹模量为B=.GPa.文献中的实验数据为99GPa.

可视化

每一次lammps运行后，会生成一个dump.atom文件。可以通过如下命令转换为Atomeye可读取的.cfg文件：

这时当前目录中就会生成.cfg之类的.cfg文件，然后通过Atomeye查看

Tab键可以改变观察方位

上下左右箭头键可以转动原子

PgUp和PgDn改变原子大小

Alt+1和Alt+2以及Alt+3改变原子颜色

先按一次9这个键，转动时就是每次以90度转动

文章来源于网络，如有侵权，请联系删除。

预览时标签不可点收录于话题#个上一篇下一篇

转载请注明：http://www.abachildren.com/sstx/151.html