ABINIT实战手册.doc

(),,,,,有什么结论。第一章ABINIT软件的运行./user1/11542/archives/2005/,需建立输入文件*.in,在abinit网站下载Ba、Ti、O的赝势文件,之后执行abinis程序,依次输入:!Abinit使用篇简介abinis的输入文件分为三类:定义了所有输入文件名的文件(,这个文件就是用来告诉abinit哪些文件是abinit读入参数的文件名,以及输出主要结果的文件的名称),定义了计算的控制参数的文件(比如取名为INP),赝势的文件。下面举例如何做好输入文件计算金刚石结构的Si的状态方程,由此也得到Si的晶格常数和体弹性模量。大致的步骤是在Si的晶格常数的实验值附近取11个数据点,也就是说取11个晶格常数或体积,然后计算在这些晶格常数下的总能。在计算得到总能后,采用状态方程拟合得到状态方程、平衡态时的体积(或晶格常数)和体弹性模量。本例子中采用的是LDA-HGH赝势。赝势的文件名为:。(紫色标示):#####BeginINP                   #设置关键词的文件名为INPOUT                 #主要的输出文件为OUT,该文件将被写入计算最重要的结果sii                          #赝势的文件名######ENDINP文件的内容为:#CrystallinecubicSi#ndtset 11                         #说明下面将有11组数据acell:3*          #晶格常数a=b=c,....+3*                  ##Groundstatecalculation kptopt  1                       #在k点网格取样时根据对称性来取样,并由下面的                                      #ngkpt和kptrlatt,或者nshiftk和shiftk来确定k点的数目 iscf  5                          #采用CG方法对能量进行优化,用在基态计算中。########################################################################Definitionoftheunitcell rprim              #下面三行定义了原胞的基矢,结构                     #Definitionoftheatomtypesntypat1               #定义原胞中原子的类别的数目,本例子中只有1类原子znucl14              #定义原胞中原子的核电荷数#Definitionoftheatomsnatom2             #定义原胞中原子的总个数,本例子中有2个原子typat 2*1           #定义每类原子的个数,本例子中第一类原子有2个xred                   #下面定义了原胞中原子的坐标                #Givesthenumberofband,explicitely(donottakethedefault)nband 16       #定义了要计算的能带的数目,最好按这样来设置:                        #nband>原胞中总的价电子数目/2+10#Exchange-correlationfunctionalixc1              #定义交换关联函数,本例子中,采用的是TeterPade参数化的LDA形式#Definitionoftheplanewavebasissetecut      #定义了平面波的切断动能     #Definitionofthek-pointgridngkpt888           #下面定义了k点网格取样的大小nshiftk4