【正文】 原子的總個(gè)數(shù)，實(shí)際上它的默認(rèn)值是natom（即原子的總個(gè)數(shù)）。 tolwfr2 shiftk2 其他的參數(shù)是用來(lái)設(shè)置晶格結(jié)構(gòu)以及原子的種類和坐標(biāo)，平面波切斷動(dòng)能，自洽收斂參數(shù)。 下面介紹采用 ABINIT 程序電子態(tài)密度時(shí)的步驟，解釋 需要輸入的關(guān)鍵詞，以及計(jì)算結(jié)果的處理說(shuō)明。) open(11,file=output, status=39。 nstep 60 自洽迭代時(shí)的最大步數(shù) 計(jì)算得到的本征值文件 zroo_DS2_EIG 內(nèi)容為： Eigenvalues (hartree) for nkpt= 57 k points: kpt 1, nband= 30, wtk= , kpt= (reduced coord) kpt 2, nband= 30, wtk= , kpt= (reduced coord) 。 計(jì)算完后對(duì) OUT 進(jìn)行分析， grep 39。Etotal39。在計(jì)算得到總能后，采用狀態(tài)方程擬合得到狀態(tài)方程、平衡態(tài)時(shí)的體積（或晶格常數(shù)）和體彈性模量。ecut39。下面采用 ndset 這個(gè)關(guān)鍵詞，在輸入文件中輸入這兩步的控制參數(shù)，一次性計(jì)算完成得到能帶結(jié)構(gòu)。number of bands to be plotted39。 編譯后得到 ，運(yùn)行 按提示輸入相應(yīng)的值； 第一次提示要輸入的量是你計(jì)算了多少條能帶 ； 第二次提示要輸入的量是你要畫出多少條能帶。當(dāng)設(shè)置 prtdos=2時(shí)，表示采用四面體方法來(lái)計(jì)算總態(tài)密度，此時(shí)不需要設(shè)置 occopt 和 tsmear了，所采用的 k 點(diǎn)需要由參數(shù) ngkpt 或 kptrlatt 參數(shù)來(lái)確定。 第五章 ABINIT 計(jì)算 基態(tài)幾何結(jié)構(gòu) 實(shí)例： 計(jì)算 AlAs 的基態(tài)幾何結(jié)構(gòu)和性質(zhì) 計(jì)算晶體 AlAs 的總能 指定進(jìn)行基態(tài)計(jì)算 kptopt 1 程序自動(dòng)產(chǎn)生 k點(diǎn)，并考慮對(duì)稱性 tolvrs 自洽循環(huán)計(jì)算停止的標(biāo)準(zhǔn) iscf 5 表示采用算法 5進(jìn)行自洽循環(huán)計(jì)算 定義原胞 acell 3* 這等價(jià)于定義 rprim 這些同 acell一起定義了基矢 ，無(wú)默認(rèn)值 定義原胞中原子的類型 ntypat 2 表示原胞中有 2類原子 znucl 13 33 關(guān)鍵詞 znucl用來(lái)指明每類原子的原子序數(shù)，其順序必須與在上 面的 ’files’文件中列出的贗勢(shì)文件相一致，這里表明第一類原子 是 Al，第二類是 As 定義原胞中的原子 natom 2 表明原胞中有 2個(gè)原子 typat 1 2 第一個(gè)原子是屬于第一類的 (Al)，第二個(gè)原子是輸入第二類的 (As) xred 這個(gè)關(guān)鍵詞按下面的格式用來(lái)指明原子的位置，每個(gè)原子有三個(gè)數(shù) 這三個(gè)數(shù)給出了第 1個(gè)原子的分?jǐn)?shù)坐標(biāo) 這三個(gè)數(shù)給出了第 2個(gè)原子 的分?jǐn)?shù)坐標(biāo) 明確給出所要計(jì)算的能帶數(shù)（最好不采用默認(rèn)值） nband 4 對(duì)絕緣體（且要是 DFT 能正確計(jì)算出真正是絕緣體），則在響應(yīng) 函數(shù)的計(jì)算中沒(méi)有必要考慮計(jì)算導(dǎo)帶 定義交換關(guān)聯(lián)函數(shù) ixc 1 采用 Teter Pade 參數(shù)化的 LDA 定義平面波基矢 ecut 平面波切斷動(dòng)能為 Hartree 定義 k點(diǎn)網(wǎng)格的設(shè)置 kptrlatt 4 4 4 按 Cartesian 坐標(biāo)來(lái)設(shè)置 k點(diǎn)網(wǎng)格，這個(gè)網(wǎng)格是立方的而且它對(duì) 4 4 4 應(yīng)的是所謂的 8 x 8 x 8 MonkhorstPack 網(wǎng)格。 在計(jì)算出來(lái)的 _DOS 文件中給出了總態(tài)密度的值，其格式如下： ABINIT package : DOS file nsppol = 1, nkpt = 145, nband(1)= 12 Tetrahedron method For identification : eigen(1:3)= Fermi energy : The DOS (in electrons/Hartree/cell) and integrated DOS (in electrons/cell) are puted, at 38001 energies (in Hartree) covering the interval between and Hartree by steps of Hartree. index energy(Ha) DOS integrated DOS 0 1 2 .......... index給出的是能量點(diǎn)數(shù)的序號(hào)，第二列給出的 是能量的數(shù)值（單位是 Hatree），在畫圖是一般以eV，并減去前面給出的 Fermi energy。這里以采用 版本計(jì)算 MgO 電子態(tài)密度（包括總態(tài)密度和 spd 分波態(tài)密度等 )為例子。new39。 第一行告訴你輸出的本征值的單位是 Hatree，共計(jì)算了 27個(gè)特殊 k 點(diǎn)。Total energy(eV)39。 下面以采用 ABINIT 計(jì)算立方的 ZrO2 晶體為例： 采用的贗勢(shì)是 和 (它們均是 LDA的 TM 贗勢(shì) )。 第一章 ABINIT 軟件的運(yùn)行 以 BaTiO3為例，需建立輸入文件 *.in，在 abinit 網(wǎng)站下載 Ba、 Ti、 O 的贗勢(shì)文件，之后執(zhí)行 abinis 程序，依次輸入： 1. 2. 3. btoi 4. btoo 5. 6. 把三個(gè)贗勢(shì)文件 拖入程序執(zhí)行框中 運(yùn)行即可！ Abinit 使用篇 簡(jiǎn)介 abinis 的輸入文件分為三類：定義了所有輸入文件名的文件（比如 ，這個(gè)文件就是用來(lái)告訴 abinit 哪些文件是 abinit 讀入?yún)?shù)的文件名，以及輸出主要