ABINIT实战手册

《ABINIT实战手册》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、ABINIT-5.8.4软件的Windows版本实战守则（byH.J.Zhao）前言第一性原理计算发表文章的要点1.你算的是什么，2.你为什么要算这个东西，3.你用了怎样的算法，4.你算出来了什么，5.你怎样解释你算出来的东西，有什么结论。第一章ABINIT软件的运行http://old.blog.edu.cn/user1/11542/archives/2005/363031.shtml以BaTiO3为例，需建立输入文件*.in，在abinit网站下载Ba、Ti、O的赝势文件，之后执行abinis程序，依次输入：1.

2、bto.in2.bto.out3.btoi4.btoo5.bto.temp6.把三个赝势文件拖入程序执行框中运行即可！Abinit使用篇简介abinis的输入文件分为三类：定义了所有输入文件名的文件（比如in.files，这个文件就是用来告诉abinit哪些文件是abinit读入参数的文件名，以及输出主要结果的文件的名称），定义了计算的控制参数的文件（比如取名为INP)，赝势的文件。下面举例如何做好输入文件计算金刚石结构的Si的状态方程，由此也得到Si的晶格常数和体弹性模量。大致的步骤是在Si的晶格常数的实验值附近取

3、11个数据点，也就是说取11个晶格常数或体积，然后计算在这些晶格常数下的总能。在计算得到总能后，采用状态方程拟合得到状态方程、平衡态时的体积（或晶格常数）和体弹性模量。本例子中采用的是LDA-HGH赝势。赝势的文件名为：14si.4.hgh。in.files的内容为（紫色标示）：#####BeginINP                   #设置关键词的文件名为INPOUT                 #主要的输出文件为OUT，该文件将被写入计算最重要的结果sii                    sisi

4、14si.4.hgh      #赝势的文件名######ENDINP文件的内容为：#CrystallinecubicSi#ndtset 11                         #说明下面将有11组数据acell:3*9.8112          #晶格常数a=b=c，将从9.8112....a.u.开始增加acell+3*0.09                  #晶格常数将以0.09a.u.的间隔进行增加#Groundstatecalculation kptopt  1             

5、          #在k点网格取样时根据对称性来取样，并由下面的                                      #ngkpt和kptrlatt,或者nshiftk和shiftk来确定k点的数目 iscf  5                          #采用CG方法对能量进行优化，用在基态计算中。########################################################################Definitionoftheunitce

6、ll rprim       0.0 0.5 0.5     #下面三行定义了原胞的基矢，本例子中Si是fcc结构      0.5 0.0 0.5         0.5 0.5 0.0  #Definitionoftheatomtypesntypat1               #定义原胞中原子的类别的数目，本例子中只有1类原子znucl14              #定义原胞中原子的核电荷数#Definitionoftheatomsnatom2             #定义原胞中原子的总个数，本例子中有2

7、个原子typat 2*1           #定义每类原子的个数，本例子中第一类原子有2个xred                   #下面定义了原胞中原子的坐标  0.0 0.0 0.0        0.250.250.25    #Givesthenumberofband,explicitely(donottakethedefault)nband 16       #定义了要计算的能带的数目，最好按这样来设置:                        #nband>原胞中总的价电子数目/2+10#Ex

8、change-correlationfunctionalixc1              #定义交换关联函数，本例子中，采用的是TeterPade参数化的LDA形式#Definitionoftheplanewavebasissetecut   20.0   #定义了平面波的切断动能     #Definitionofthek-pointgri