ABINIT实战手册.doc

ABINIT-
(by )
前言第一性原理计算发表文章的要点
1. 你算的是什么,
2. 你为什么要算这个东西,
3. 你用了怎样的算法,
4. 你算出来了什么,
5. 你怎样解释你算出来的东西,有什么结论。
第一章 ABINIT软件的运行
./user1/11542/archives/2005/
以BaTiO3为例,需建立输入文件*.in,在abinit网站下载Ba、Ti、O的赝势文件,之后执行abinis程序,依次输入:
1.
2.
3. btoi
4. btoo
5.
6. 把三个赝势文件拖入程序执行框中
运行即可!
Abinit使用篇简介
abinis的输入文件分为三类:定义了所有输入文件名的文件(,这个文件就是用来告诉abinit哪些文件是abinit读入参数的文件名,以及输出主要结果的文件的名称),定义了计算的控制参数的文件(比如取名为INP),赝势的文件。下面举例如何做好输入文件计算金刚石结构的Si的状态方程,由此也得到Si的晶格常数和体弹性模量。
大致的步骤是在Si的晶格常数的实验值附近取11个数据点,也就是说取11个晶格常数或体积,然后计算在这些晶格常数下的总能。在计算得到总能后,采用状态方程拟合得到状态方程、平衡态时的体积(或晶格常数)和体弹性模量。
本例子中采用的是LDA-HGH赝势。赝势的文件名为:。
(紫色标示):
#####Begin
INP                    #设置关键词的文件名为INP
OUT                  #主要的输出文件为OUT,该文件将被写入计算最重要的结果
sii                    
si
si
       #赝势的文件名
######END
INP文件的内容为:
# Crystalline cubic Si
#
ndtset  11                          #说明下面将有11组数据
acell: 3*                   # 晶格常数a=b=c,....
acell+ 3*                     #
#Ground state calculation
 kptopt   1                        #在 k点网格取样时根据对称性来取样,并由下面的
                                       # ngkpt和kptrlatt, 或者nshiftk和shiftk来确定k点的数目
 iscf   5                           #采用CG方法对能量进行优化,用在基态计算中。
#######################################################################
#Definition of the unit cell 
rprim 
                 #下面三行定义了原胞的基矢,结构
              
             
#Definition of the atom types
ntypat 1                #定义原胞中原子的类别的数目,本例子中只有1类原子
znucl 14               #定义原胞中原子的核电荷数
#Definition of the atoms
natom 2              #定义原胞中原子的总个数,本例子中有2个原子
typat  2*1            #定义每类原子的个数,本例子中第一类原子有2个
xred                    #下面定义了原胞中原子的坐标
             
       
#Gives the number of band, explicitely (do not take the default)
nband  16        #定义了要计算的能带的数目,最好按这样来设置: