abinit实战手册.doc

ABINIT- 软件的 Windows 版本实战守则( by ) 前言第一性原理计算发表文章的要点 1. 你算的是什么, 2. 你为什么要算这个东西, 3. 你用了怎样的算法, 4. 你算出来了什么, 5. 你怎样解释你算出来的东西,有什么结论。第一章 ABINIT 软件的运行 ./user1/11542/archives/2005/ 以 BaTiO 3 为例, 需建立输入文件*.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 # 下面三行定义了原胞的基矢,本例子中 Si是 fcc 结构 #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# 定义了要计算的能带的数目,最好按这样来设置: # nband > 原胞中总的价电子数目/2+ 10 #Exchange-correlation functional ixc 1# 定义交换关联函数,本例子中,采用的是 Teter Pade 参数化的 LDA 形式#Definition of the planewave basis set ecut # 定义了平面波的切断动能#Definition of the k-point grid ngkpt 888# 下面定义了 k 点网格取样的大小 nshiftk 4 shiftk #Definition of the SCF procedure nstep 60# 电子自洽迭代的最大步数 diema