ABINIT實战手册.docVIP

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