ABINIT实战手册.doc

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.in2. bto.out3. btoi4. btoo5. bto.temp6. 把三个赝势文件拖入程序执行框中运行即可！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 calculationkptopt 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 functionalixc 1 #定义交换关联函数，本例子中，采用的是Teter Pade 参数化的LDA形式#Definition of the planewave basis setecut 20.0 #定义了平面波的切断动能 #Definition of the k-point gridngkpt 8 8 8 #下面定义了k点网格取样的大小nshiftk 4 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5#Definition of the SCF procedurenstep 60 #电子自洽迭代的最大步数diemac 12.0 #介电常数设置tolvrs 1.0d-20 #电子自洽收敛的标准- END-计算完后，得到OUT文件，用下面的命令：grep volume OUT得到如下的内容：Unit cell volume ucvol= 2.3610688E+02 bohr3Unit cell volume ucvol= 2.4266422E+02 bohr3Unit cell volume ucvol= 2.4934185E+02 bohr3Unit cell volume ucvol= 2.5614088E+02 bohr3Unit cell volume ucvol= 2.6306239E+02 bohr3Unit cell volume ucvol= 2.7010748E+02 bohr3Unit cell volume ucvol= 2.7727725E+02 bohr3Unit cell volume ucvol= 2.8457279E+02 bohr3Unit cell volume ucvol= 2.9199518E+02 bohr3Unit cell volume ucvol= 2.9954553E+02 bohr3Unit cell volume ucvol= 3.0722493E+02 bohr3然后用下面的命令：grep Etotal OUT得到如下的内容： Etotal= -7.92750029752797E+00 Etotal= -7.92997465524506E+00 Etotal= -7.93167675973445E+00 Etotal= -7.93266612552653E+00 Etotal= -7.93299797094926E+00 Etotal= -7.93272412167288E+00 Etotal= -7.93189304516315E+00 Etotal= -7.93055314505404E+00 Etotal= -7.92874706916830E+00 Etotal= -7.92651686675655E+00 Etotal= -7.92390242137627E+00因此，Volume 和Etotal对应的关系为：2.3610688E+02 -7.92750029752797E+002.4266422E+02 -7.92997465524506E+002.4934185E+02 -7.93167675973445E+002.5614088E+02 -7.93266612552653E+002.6306239E+02 -7.93299797094926E+002.7010748E+02 -7.93272412167288E+002.7727725E+02 -7.93189304516315E+002.8457279E+02 -7.93055314505404E+002.9199518E+02 -7.92874706916830E+002.9954553E+02 -7.92651686675655E+003.0722493E+02 -7.92390242137627E+00下面就可以用 Birch-Murnaghan 3阶状态方程进行 (Birch F, Phys. Rev. 71, p809 (1947)拟合得到体弹性模量和平衡状态下的体积：V0 = 263.276940709097 a.u.3B0 = 95.497(GPa)第二章 ABINIT参数设定与收敛测试/u2/85704/showart_1935274.html ABINIT计算晶体时主要参数的确定(切断动能和k点网格)。 在采用平面波赝势法进行固体的电子结构计算时，为了确保计算的精度和计算结果的可靠性，在计算晶体的物理性质之前，要进行几个重要参数的测试，以保证这些参数的选取使得计算结果有很好的收敛性，这些参数是平面波的切断动能和k点网格取样的大小。 下面以采用ABINIT计算立方的ZrO2晶体为例： 采用的赝势是40zr.psp_mod和8o.psp_mod(它们均是LDA的TM赝势)。 输入文件，in.files的内容为： INP OUT zroi zroo zro ./40zr.psp_mod ./8o.psp_mod 在测试平面波切断动能的收敛性时，我们通过设置平面波切断动能从20 Ha开始，以2 Ha递增，直到58 Ha，其他参数不变计算ZrO2的总能。在ABINIT的输入文件中很方便的通过ndtset来设置。 输入文件INP如下（紫色标示）： # Crystalline ZrO2-cubic ndtset 20 #表示有20组数据 ecut: 20.0 #平面波切断动能从20 Ha开始 ecut+ 2 #以2 Ha递增，也就是Ecut为20.0 + i*2.0, (i从1到20) #Definition of the unit cell acell 3*9.65285 #设置晶格常数a=b=c为9.65285 a.u. rprim 0.0 0.5 0.5 #同上面的acell确定了原胞的基矢 0.5 0.0 0.5 0.5 0.5 0.0 #Definition of the atom types ntypat 2 #设置原胞的原子种类数，这里有2类原子 znucl 40 8 #每类原子的核电荷数Z natom 3 #原胞中总的原子数目，这里原胞总共有3个原子 typat 1 2*2 #第一类原子有1个，第二类原子有2个 xred #下面的按分数坐标给原胞中原子的坐标位置 0.00 0.000 0.000 0.25 0.25 0.25 0.75 0.75 0.75 #Definition of the k-point grid kptopt 1 #设置生成k点的方法，这里表明有ngkpt和nshiftk来确定k点网格的大小 ngkpt 8 8 8 #设置对布里渊区进行8x8x8网格的划分 nshiftk 4 #对划分得到的k点按下面的偏移量进行平移 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 #Definition of the SCF procedure #ecut 40.0 iscf 5 #自洽迭代中的算法，这里采用cg方法优化。 toldfe 1.0d-10 #自洽迭代迭过程总能收敛的标准 diemac 3.0 #体系宏观的介电常数 nstep 60 #最大自洽迭代步数 #Definition of the outfile information prtwf 0 #不输出波函数文件。 运行下面的命令 abins OUT 对OUT采用进行分析，用下面的命令取出一些数据: grep ecut OUT summary grep Total energy(eV) OUT summary 得到summary类似如下数据: 2.00000000E+01 -2.16713704994688E+03 2.20000000E+01 -2.17244914205186E+03 2.40000000E+01 -2.17629539439455E+03 2.60000000E+01 -2.17888687489711E+03 2.80000000E+01 -2.18045824630388E+03 3.00000000E+01 -2.18139622512928E+03 3.20000000E+01 -2.18192168325676E+03 3.40000000E+01 -2.18219516090481E+03 3.60000000E+01 -2.18232706572898E+03 3.80000000E+01 -2.18237668482667E+03 4.00000000E+01 -2.18239923089573E+03 4.20000000E+01 -2.18241955706857E+03 4.40000000E+01 -2.18244775901834E+03 4.60000000E+01 -2.18248714583389E+03 4.80000000E+01 -2.18253519502611E+03 5.00000000E+01 -2.18258901301065E+03 5.20000000E+01 -2.18264245486986E+03 5.40000000E+01 -2.18269485143783E+03 5.60000000E+01 -2.18274313138312E+03 5.80000000E+01 -2.18278585826504E+03 画图可以看到当Ecut取40 Ha时，体系的总能有很好的收敛了。 下面测试对k点网格的，in.files的输入文件同上。我们这里把k点网格从6x6x6开始增加，以2x2x2递增。INP的输入文件（紫色标示）： # Crystalline ZrO2-cubic ndtset 10 ngkpt: 6 6 6 ngkpt+ 2 2 2 #Definition of the unit cell acell 3*9.65285 rprim 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 #Definition of the atom types ntypat 2 znucl 40 8 natom 3 typat 1 2*2 xred 0.00 0.000 0.000 0.25 0.25 0.25 0.75 0.75 0.75 #Definition of the k-point grid kptopt 1 #ngkpt 8 8 8 nshiftk 4 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 #Definition of the SCF procedure ecut 40.0 iscf 5 toldfe 1.0d-10 diemac 3.0 nstep 60 #Definition of the outfile information prtwf 0 运行 abinis OUT& 计算完后对OUT进行分析， grep ngkpt OUT comment grep Total energy(eV) OUT comment 取出得到的数据如下： 28.0000 -2182.416564 60.0000 -2182.416514 110.0000 -2182.416506 182.0000 -2182.416526 280.0000 -2182.416499 408.0000 -2182.416512 570.0000 -2182.416506 770.0000 -2182.416514 1012.0000 -2182.416515 1300.0000 -2182.416504 画图可以看到8x8x8的k点网格就能保证体系的总能很好的收敛了。 第三章 ABINIT计算晶体的能带结构/user1/11542/archives/2006/1273410.shtml采用第一原理的电子结构计算方法来计算晶体的能带结构一般来说，要进行两个步骤。这不论是采用VASP、PWSCF还是ABINIT这些程序。步骤为：先进行自洽的电子结构迭代得到自洽计算得到的电荷密度，然后读入这个自洽得到的电荷密度，进行非自洽的计算得到体系的本征值。 下面以立方的ZrO2晶体为例，采用ABINIT来计算。下面采用ndset这个关键词，在输入文件中输入这两步的控制参数，一次性计算完成得到能带结构。并对相关的相关的输入参数进行解释： 采用的赝势是:40zr.psp_mod和8o.psp_mod，输入文件 in.files的内容为 INP OUT zroi zroo zro ./40zr.psp_mod ./8o.psp_mod 此in.files文件的内容设置了主要输入文件的名称和赝势文件所在的目录。 主要输入文件INP的内容为（紫色标示）： # Crystalline ZrO2-cubic ndtset 2 #表示有两组控制参数：第一组控制参数用来设置自洽计算，第二组是非自洽的本征值计算 #Dataset 1 : usual self-consistent calculation kptopt1 1 #自洽计算中设置k点网格取样的方法，表明采用ngkpt和shfitk来设置k点网格。 nshiftk1 4 #使生成的k点进行平移。 shiftk1 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 ngkpt1 8 8 8 #K点网格取样，网格划分的分割数 prtden1 1 #表明输出电荷密度文件。 tolvrs1 1.0d-20 # 自洽计算收敛的标准 iscf1 5 #自洽迭代计算时，采用 CG方法来优化有效势 #Dataset 2 : #用来设置本征值计算时的参数 iscf2 -2 #表明非自洽计算 getden2 -1 #读入上一组数据进行自洽计算得到的电荷密度文件 kptopt2 -5 #负数，表示下面计算能带计算时，有5段特殊线(由6个特殊k点来确定) ndivk2 10 12 18 8 8 #每段特殊线上分几等份 kptbounds2 #特殊k点的坐标 0.5 0.25 0.75 # W point 0.5 0.0 0.0 # L point 0.0 0.0 0.0 # Gamma point 0.5 0.0 0.5 # X point 0.5 0.25 0.75 # W point 0.0 0.0 0.0 # Gamma enunit2 0 #输出本征值时，本征值的单位，这里为0表示是以Hatree为单位给出 prteig2 1 #表明输出本征值到文件中 #Definition of the unit cell acell 3*9.48196 #设置晶常数 rprim 0.0 0.5 0.5 #设置计算原胞的基矢(同上面的acell一起构成原胞的基矢) 0.5 0.0 0.5 0.5 0.5 0.0 #Definition of the atom types ntypat 2 #原胞的原子种类数目 znucl 40 8 #每类原子的核电荷数 natom 3 #原胞中总的原子数目 typat 1 2*2 #表示第一个原子是第一类的，后面两个原子是第二类的。 xred #以分数坐标给出原子的位置 0.00 0.000 0.000 0.25 0.25 0.25 0.75 0.75 0.75 #Definition of the SCF procedure nband 30 #在计算中考虑多少条能带 ecut 40.0 #平面波切断动能 diemac 3.0 #体系的宏观介电常数，给一个近似值就可以了。 nstep 60 #自洽迭代时的最大步数 计算得到的本征值文件zroo_DS2_EIG内容为： Eigenvalues (hartree) for nkpt= 57 k points: kpt# 1, nband= 30, wtk= 1.00000, kpt= 0.5000 0.2500 0.7500 (reduced coord) -1.41786 -0.62938 -0.62938 -0.61620 -0.26378 -0.26378 0.14567 0.14721 0.23910 0.23910 0.27711 0.32275 0.49670 0.50302 0.57458 0.57458 0.61760 0.78769 0.90930 0.92430 0.92430 1.05159 1.05159 1.10749 1.20089 1.28891 1.28891 1.35618 1.42279 1.46219 kpt# 2, nband= 30, wtk= 1.00000, kpt= 0.5000 0.2250 0.6750 (reduced coord) -1.41788 -0.63312 -0.62545 -0.61596 -0.26959 -0.25851 0.13111 0.16454 0.22812 0.24784 0.27846 0.32146 0.49126 0.50549 0.56954 0.58345 0.61790 0.78670 0.89599 0.91521 0.94726 1.02556 1.08306 1.08822 1.21602 1.28607 1.29698 1.37003 1.41244 1.48297 。 把它们转换为k点的绝对值和本征值后，画图如下：附1 处理abinit计算得到的本征值文件/user1/11542/archives/2006/1275248.shtml在采用abinit计算固体的能带结构时，abinit输出的本征值文件格式为： Eigenvalues (hartree) for nkpt= 27 k points: kpt# 1, nband= 20, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord) -0.37264 -0.18371 -0.15298 -0.15298 -0.11442 -0.01609 0.01563 0.02497 0.03905 0.03905 0.04309 0.05183 0.05183 0.07649 0.07991 0.07991 0.08411 0.08411 0.08768 0.11715 kpt# 2, nband= 20, wtk= 1.00000, kpt= 0.0000 0.0000 0.0192 (reduced coord) -0.37255 -0.18459 -0.15289 -0.15289 -0.11343 -0.01601 0.01572 0.02505 0.03913 0.03913 0.04317 0.04982 0.04982 0.07658 0.08000 0.08000 0.08627 0.08627 0.08776 0.11724 。 第一行告诉你输出的本征值的单位是Hatree，共计算了27个特殊k点。 下面是分别给出每个k点的的坐标，以及所计算的本征值。 但是要把这些本征值和k点的坐标对应起来画图，怎么办呢？这里提供一个我根据pwscf程序包中的一个代码改写后的小工具plotband_abinit.f90，它用来把每个k点与前一个k点之间的间隔一一加起来得到的距离和本征值对应起来。 plotband_abinit.f90的内容如下： ! for ABINIT ! program prog real, allocatable : e(:,:) real, allocatable : k(:,:) real, dimension(3) :k0,a character(len=32): input, output character(len=32): xx, yy write(6,*) number of calculated bands read(5,*) nbands write(6,*) number of bands to be plotted read(5,*) nbands2 write(6,*) number of k-points read(5,*) nk write(6,*) fermi level (eV) read(5,*) ef write(6,*) name bande.in e bande.out read(5,*) input,output allocate(e(nk,nbands) allocate(k(nk,3) open(10,file=input, status=old) open(11,file=output, status=new) read(10,*) xx, xx,xx,xx, nkpt, xx, xx do i=1,nk read(10,*)yy, ik, yy, nb, yy, wtk, yy, (k(i,j),j=1,3),yy write(6,*)(k(i,j),j=1,3) read(10,*) (e(i,n),n=1,nbands) write(13,9030) (e(i,n),n=1,nbands) 9030 format (8f9.5) enddo do j=1,nbands2 dk=0 do i=1,nk if (i.eq.1) then k0=k(i,:) endif a=k(i,:)-k0 dk=dk+sqrt(dot_product(a,a) write(11,*)dk,(e(i,j)-ef)*27.2116 k0=k(i,:) enddo write(11,*) enddo stop end program prog 采用任何一种fortran90的编译器(比如pgf90和ifc，ifort)编译后。 pgf90 plotband_abinit.f90 它的输入文件为计算完后的本征值。注意这里是针对abinit计算输出的本征值单位为Hatree。要补充的是：abinit再输出本征值文件时有个小小的bug，在以eV单位输出本征值时，输出的本征值文件中会漏掉第51个的本征值。因此建议在采用abinit计算输出本征值时以Hatree为单位。 编译后得到a.out，运行a.out按提示输入相应的值； 第一次提示要输入的量是你计算了多少条能带；第二次提示要输入的量是你要画出多少条能带。 第三次提示要输入的量是你计算了多少个k点；第四次提示要输入的是体系的费米能级。(注意不要从自洽计算得到的 OUT中取出费米能级，abinit给出的有错误。而是要从计算态密度时的得到的态密度文件中取出。) ；第五此提示要输入的是包含了计算得到的本征值的文件名，以及转换后输出的文件名。 然后就可以对转换后的文件进行画图了。附2 不同晶格的简约布里渊区形状及特殊K点坐标1. FCC晶格的简约布里渊区形状及特殊K点坐标2. BCC晶格的简约布里渊区形状及特殊K点坐标3. SC晶格的简约布里渊区形状及特殊K点坐标4. 六角晶格的简约布里渊区形状及特殊K点坐标5. 正交晶格的简约布里渊区行状及特殊k点坐标以上是采用xcrysden软件画出了各种晶格的简约布里渊区形状，并标出了几个特殊K点的位置、标记及其坐标（以倒格子基矢为坐标系）。第四章 ABINIT计算电子态密度/user1/11542/archives/2007/1664470.shtml电子态密度(Density of States, DOS)定义为单位能量范围内所允许的能级数。在电子结构计算程序中进行计算时，一般是按如下的步骤来进行的：先采用小的k点网格进行自洽计算，然后采用较密的k点网格进行非自洽计算。下面介绍采用ABINIT程序电子态密度时的步骤，解释需要输入的关键词，以及计算结果的处理说明。这里以采用ABINIT5.2.4版本计算MgO电子态密度（包括总态密度和spd分波态密度等)为例子。在例子中MgO的晶格常数(它是面心立方晶体，NaCl结构，a=7.921 a.u.)1取为实验值，平面波切断动能和k点网格大小分别取为24 Ha，10x10x10（自洽计算中）和16x16x16（态密度计算中）。1 O.L. Anderson and P. Andreatch, J. Am. Ceram. Soc. 49, 404 (1966).1. 计算总态密度 通过abinit中的ndtset设置两组数据进行计算，先是自洽计算，后是态密度计算的设置。在计算态密度时，需将prtdos设置为1、2或3。当prtdos设置1时，采用离散(smearing)的方法来确定态密度，即由本征值、k点权重(参数wtk的值)、每个本征值的占有数（由occopt设置的参数所选择的方法来确定占有数)以及离散展宽参数(tsmear参数设置的值)来计算，此时需要注意或设置wtk、occopt和tsmear参数。一般对一维体系或者在k点较少的情况采用prtdos=1的方法来计算态密度，以得到较为光滑的DOS曲线。当设置prtdos=2时，表示采用四面体方法来计算总态密度，此时不需要设置occopt和tsmear了，所采用的k点需要由参数ngkpt或kptrlatt参数来确定。当prtdos=3时，表示采用四面体方法计算分波态密度(angular-momentum projected)，同时需设置原子球半径。以计算面心立方MgO的总态密度为例子。i). in.files输入文件的内容:dos.outmgoimgoomgo12mg.pspnc8o.pspncii)主要输入文件tdos.in的内容（紫色标注）：ndtset 2kptopt 1 nshiftk1 4shiftk1 0.5 0.5 0.50.5 0.0 0.00.0 0.5 0.00.0 0.0 0.5ngkpt1 8 8 8prtden1 1 toldfe1 1.0d-6iscf1 5 #Dataset 2 : DOS calculationiscf2 -3getden2 1getwfk2 1prtdos2 2dosdeltae 0.00005 ngkpt2 16 16 16tolwfr2 1.0d-16shiftk20.0 0.0 0.00.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0#Definition of the unit cellacell 3*7.921rprim0.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0ntypat 2znucl 12 8natom 2typat 1 2 xred0.0 0.0 0.00.5 0.5 0.5nband 12ixc 1 ecut 24.0nstep 45diemac 9.0 这里利用ndtset设置两组数据，第一组是进行自洽计算以得到电荷密度，第二组是非自洽计算以得到总态密度。在第一组数据中:kptopt 1 #设置通过ngkpt由程序自动生成k点nshiftk1 4 #由于是fcc结构，为了使得的产生的k点对计算更有效，这里将所产生的k点进行移动shiftk1 #设置了移动的量0.5 0.5 0.50.5 0.0 0.00.0 0.5 0.00.0 0.0 0.5ngkpt1 8 8 8 #设置k点网格的大小。由于是在自洽计算中，此时设置的数较小。prtden1 1 #表示要将电荷密度输出到文件中，以供第二组中的计算用到。toldfe 1.0d-6 #自洽迭代计算中总能变化收敛的标准iscf1 5 #表示采用CG方法来优化总能以得到基态在第二组数据中：iscf2 -3 #表示进行非自洽计算，一般用在态密度和STM计算以及线性响应计算中对ddk微扰的计算。getden2 1 #表示从上一步中得到的电荷密度作为初始的电荷密度读入getwfk2 1 #表示读入上一步的波函数prtdos2 2 #表示要输出总态密度（或进行总态密度计算)dosdeltae 0.00005 #态密度计算时能量轴上的能量刻度ngkpt2 16 16 16 #在态密度计算时，k点网格应该设置的密一些。tolwfr2 1.0d-16shiftk20.0 0.0 0.00.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0其他的参数是用来设置晶格结构以及原子的种类和坐标，平面波切断动能，自洽收敛参数。在计算出来的_DOS文件中给出了总态密度的值，其格式如下：ABINIT package : DOS filensppol = 1, nkpt = 145, nband(1)= 12Tetrahedron methodFor identification : eigen(1:3)= -0.429 0.203 0.203Fermi energy : 0.20340927The DOS (in electrons/Hartree/cell) and integrated DOS (in electrons/cell) are computed,at 38001 energies (in Hartree) covering the intervalbetween -0.5000 and 1.4000 Hartree by steps of 0.00005 Hartree.index energy(Ha) DOS integrated DOS0 -0.50000 0.0000 0.00001 -0.49995 0.0000 0.00002 -0.49990 0.0000 0.0000.index给出的是能量点数的序号，第二列给出的是能量的数值（单位是Hatree），在画图是一般以eV，并减去前面给出的Fermi energy。第三列是总态密度值，第四列是电子数。2. 计算各个原子的分波态密度（紫色标注）ndtset 2kptopt 1nshiftk1 4shiftk10.5 0.5 0.50.5 0.0 0.00.0 0.5 0.00.0 0.0 0.5ngkpt1 8 8 8prtden1 1toldfe1 1.0d-6iscf1 5#Dataset 2 : DOS c