最新VASP的輸入輸出文件

1、vaspvasp的輸入輸出文件的輸入輸出文件侯侯 柱柱 鋒鋒復(fù)旦大學(xué)物理系復(fù)旦大學(xué)物理系2005, 12, 3北京，宏劍hou zhufeng2內(nèi)容(1)輸入文件輸入文件 incar :計(jì)算控制參數(shù)文件計(jì)算控制參數(shù)文件 poscar :描述體系結(jié)構(gòu)的文件描述體系結(jié)構(gòu)的文件 kpoints :k點(diǎn)取樣設(shè)置文件點(diǎn)取樣設(shè)置文件 potcar :贗勢(shì)文件贗勢(shì)文件輸出文件輸出文件 outcar :最主要的輸出文件最主要的輸出文件 doscar : 電子態(tài)密度文件電子態(tài)密度文件 eigenval :本征值文件本征值文件 oszicar :每次迭代或原子遲豫每次迭代或原子遲豫(或或md)的信息的信息 chg

2、和和chgcar :電荷密度文件電荷密度文件 wavecar :波函數(shù)文件波函數(shù)文件hou zhufeng3內(nèi)容(2)輸出文件輸出文件(續(xù)續(xù)) contcar :原子遲豫或原子遲豫或md后的體系結(jié)構(gòu)文件后的體系結(jié)構(gòu)文件 ibzkpt :布里淵區(qū)中的布里淵區(qū)中的k點(diǎn)點(diǎn) pcdat :對(duì)關(guān)聯(lián)函數(shù)對(duì)關(guān)聯(lián)函數(shù) xdatcar :在在md時(shí)時(shí), 原子位置變化的跟蹤文件原子位置變化的跟蹤文件 procar和和proout : 波函數(shù)投影或分解的文件波函數(shù)投影或分解的文件 locpot :總的局域勢(shì)總的局域勢(shì) elfcar :電子局域函數(shù)電子局域函數(shù)hou zhufeng4輸入文件輸入文件hou zhufe

3、ng5incar(1) 設(shè)置計(jì)算的精度 設(shè)置是做什么計(jì)算 設(shè)置交換關(guān)聯(lián)函數(shù) 設(shè)置優(yōu)化的算法和收斂標(biāo)準(zhǔn) 設(shè)置md的步長(zhǎng)、溫度、時(shí)間 設(shè)置每個(gè)波函數(shù)或電子的占有數(shù)(也就是 smearing方法及相關(guān)的參數(shù))等等 自由格式 共100多個(gè)關(guān)鍵詞，一般都有合適的默認(rèn)值。通常需設(shè)置的有10個(gè)左右。即使incar文件無(wú)任何內(nèi)容也可以進(jìn)行計(jì)算。hou zhufeng6incar(2)：例子靜態(tài)計(jì)算靜態(tài)計(jì)算(計(jì)算總能和自洽的電荷密度計(jì)算總能和自洽的電荷密度)懶：刪掉計(jì)算目錄下的懶：刪掉計(jì)算目錄下的wavecar，incar中不作任何設(shè)置中不作任何設(shè)置推薦：自己手動(dòng)設(shè)置推薦：自己手動(dòng)設(shè)置encut, title,

4、 istart, icharg, ismear, ediff, prec等等system = siliconencut = 350istart = 0icharg = 2ismear = -5ediff = 1e-5prec = accurate設(shè)置標(biāo)題，以說(shuō)明所計(jì)算的體系設(shè)置標(biāo)題，以說(shuō)明所計(jì)算的體系設(shè)置平面波切斷動(dòng)能設(shè)置平面波切斷動(dòng)能(不采用默認(rèn)值不采用默認(rèn)值)說(shuō)明這次計(jì)算是一次全新的計(jì)算說(shuō)明這次計(jì)算是一次全新的計(jì)算按體系中的原子構(gòu)造初始的原子密度按體系中的原子構(gòu)造初始的原子密度采用四面體方法采用四面體方法電子迭代的收斂標(biāo)準(zhǔn)是電子迭代的收斂標(biāo)準(zhǔn)是1e-5精度為精度為accurate可以在一行

5、設(shè)置多個(gè)關(guān)鍵詞的值，但是每個(gè)關(guān)鍵值之間用分號(hào)可以在一行設(shè)置多個(gè)關(guān)鍵詞的值，但是每個(gè)關(guān)鍵值之間用分號(hào)(;)隔開(kāi)。隔開(kāi)。如如ismear= 0; sigma= 0.2。當(dāng)想不用當(dāng)想不用incar中某個(gè)關(guān)鍵詞的值時(shí)，在該行前面加上井號(hào)中某個(gè)關(guān)鍵詞的值時(shí)，在該行前面加上井號(hào)(#)注釋掉，注釋掉，如如#ismear=0; sigma = 0.2hou zhufeng7poscar(1) 描述了體系的結(jié)構(gòu)文件: 原(或晶)胞的基矢, 原子的位置, 原子是否移動(dòng)，原子的初始速度等) 分?jǐn)?shù)坐標(biāo)和卡笛爾坐標(biāo)cubic sic3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1dir

6、ect 0.00 0.00 0.000.25 0.25 0.25例子例子:cubic sic3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1cartesian 0.00 0.00 0.000.25 0.25 0.25設(shè)置體系的名稱設(shè)置體系的名稱晶格常數(shù)或縮放系數(shù)晶格常數(shù)或縮放系數(shù)原原(或晶或晶)胞的基矢胞的基矢每類原子的個(gè)數(shù)每類原子的個(gè)數(shù)確定按何種坐標(biāo)來(lái)寫(xiě)原子位置確定按何種坐標(biāo)來(lái)寫(xiě)原子位置第一類原子的第一個(gè)坐標(biāo)第一類原子的第一個(gè)坐標(biāo)第二類原子的第一坐標(biāo)第二類原子的第一坐標(biāo)第七行以字母第七行以字母d開(kāi)頭表示下面的是分?jǐn)?shù)坐標(biāo)，如果是以開(kāi)頭表示下面的是分?jǐn)?shù)坐標(biāo)，如果

7、是以c或或k開(kāi)頭表示下面的坐開(kāi)頭表示下面的坐標(biāo)是卡笛爾坐標(biāo)。標(biāo)是卡笛爾坐標(biāo)。hou zhufeng8poscar(2) 對(duì)原子位置進(jìn)行選擇性的遲豫(selective dynamic)cubic sic3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1selective dynamicdirect 0.000 0.000 0.000 f f f0.248 0.248 0.248 t t t例子例子:表示對(duì)體系中的原子進(jìn)行限制性的遲豫表示對(duì)體系中的原子進(jìn)行限制性的遲豫或優(yōu)化，這一行以字母或優(yōu)化，這一行以字母s開(kāi)頭。開(kāi)頭。如何寫(xiě)出具有復(fù)雜結(jié)構(gòu)的晶體的如何寫(xiě)出具有復(fù)雜

8、結(jié)構(gòu)的晶體的poscar：a) 查到該晶體的晶格常數(shù)、空間群查到該晶體的晶格常數(shù)、空間群和烏科夫和烏科夫(wyckoff)坐標(biāo)坐標(biāo); b)用用material studio中中crystal builder或其他的晶體學(xué)軟或其他的晶體學(xué)軟件畫(huà)出晶體，并得到各個(gè)原子的坐標(biāo)件畫(huà)出晶體，并得到各個(gè)原子的坐標(biāo);或或c)根據(jù)空間群國(guó)際表寫(xiě)出各等價(jià)位置根據(jù)空間群國(guó)際表寫(xiě)出各等價(jià)位置; d)晶晶體數(shù)據(jù)庫(kù)體數(shù)據(jù)庫(kù)固定第一類原子的位置，固定第一類原子的位置，對(duì)第二類原子的位置進(jìn)行優(yōu)化對(duì)第二類原子的位置進(jìn)行優(yōu)化hou zhufeng9kpoints(1) 設(shè)置布里淵區(qū)k點(diǎn)取樣或k點(diǎn)的坐標(biāo),計(jì)算能帶時(shí)k點(diǎn)的高對(duì)稱線

9、 設(shè)置方式:手動(dòng)輸入所有的k點(diǎn)，按monkhorst-pack方法自動(dòng)產(chǎn)生 格式：普通，四面體，m-p網(wǎng)格，line模式 卡笛爾坐標(biāo)系和倒格子坐標(biāo)系hou zhufeng10kpoints(2) 按普通格式手動(dòng)輸入各個(gè)按普通格式手動(dòng)輸入各個(gè)k點(diǎn)點(diǎn)(推薦用在計(jì)算能帶時(shí)推薦用在計(jì)算能帶時(shí))k-points entered manually5reciprocal 0.500 0.000 0.500 1.00 0.475 0.000 0.475 1.00 0.450 0.000 0.450 1.00 0.425 0.000 0.425 1.00 0.400 0.000 0.400 1.00標(biāo)題或注釋行

10、，無(wú)特別意義標(biāo)題或注釋行，無(wú)特別意義k點(diǎn)的數(shù)目點(diǎn)的數(shù)目以字母以字母r開(kāi)頭表示開(kāi)頭表示k點(diǎn)是按倒格子坐標(biāo)系點(diǎn)是按倒格子坐標(biāo)系前三個(gè)數(shù)是前三個(gè)數(shù)是k點(diǎn)的坐標(biāo)，最后一個(gè)數(shù)是相應(yīng)點(diǎn)的坐標(biāo)，最后一個(gè)數(shù)是相應(yīng)k點(diǎn)的權(quán)重點(diǎn)的權(quán)重(下面共下面共5個(gè)個(gè)k點(diǎn)點(diǎn))如果是以卡笛爾坐標(biāo)系來(lái)寫(xiě)k點(diǎn)坐標(biāo)，則第三行以字母c開(kāi)頭。 按四面體格式手動(dòng)輸入各個(gè)按四面體格式手動(dòng)輸入各個(gè)k點(diǎn)點(diǎn)(復(fù)雜，不推薦用復(fù)雜，不推薦用)example file 4 cartesian 0.0 0.0 0.0 1. 0.0 0.0 0.5 1. 0.0 0.5 0.5 2. 0.5 0.5 0.5 4. tetrahedra 1 0.1833333

11、33333333 6 1 2 3 4 標(biāo)題或注釋行，無(wú)特別意義k點(diǎn)的數(shù)目以字母c開(kāi)頭表示k點(diǎn)是按卡笛爾坐標(biāo)系前三個(gè)數(shù)是k點(diǎn)的坐標(biāo)，最后一個(gè)數(shù)是相應(yīng)k點(diǎn)的坐標(biāo)(下面共4個(gè)k點(diǎn))當(dāng)smearing方面采用的是四面體方法時(shí)，以字母t開(kāi)頭表示下面的是描述各個(gè)四面體的信息四面體的個(gè)數(shù)，每個(gè)四面體的體積權(quán)重每個(gè)四面體頂角的權(quán)重，每個(gè)四面體頂角k點(diǎn)坐標(biāo)的序號(hào)當(dāng)當(dāng)smearing時(shí)不采用四面體方法，則時(shí)不采用四面體方法，則tetrahedra后面的內(nèi)容不起作用后面的內(nèi)容不起作用hou zhufeng11kpoints(2) 按按line模式手動(dòng)輸入各個(gè)高對(duì)稱模式手動(dòng)輸入各個(gè)高對(duì)稱k點(diǎn)點(diǎn)(用在計(jì)算能帶中，用在計(jì)

12、算能帶中，4.6以上版本支持以上版本支持)k along high symmetry lines 10 line-mode rec 0.00 0.00 0.00 ! gamma 0.50 0.50 0.00 ! x0.50 0.50 0.00 ! x 0.50 0.75 0.25 ! w0.50 0.75 0.25 ! w 0.00 0.00 0.00 ! gamma 標(biāo)題或注釋行，無(wú)特別意義標(biāo)題或注釋行，無(wú)特別意義每對(duì)高對(duì)稱點(diǎn)之間產(chǎn)生每對(duì)高對(duì)稱點(diǎn)之間產(chǎn)生10個(gè)個(gè)k點(diǎn)點(diǎn)以字母以字母l開(kāi)頭表示按開(kāi)頭表示按line模式產(chǎn)生模式產(chǎn)生k點(diǎn)點(diǎn)以字母以字母r開(kāi)頭表示開(kāi)頭表示k點(diǎn)按倒格子坐標(biāo)系點(diǎn)按倒格子坐標(biāo)

13、系每對(duì)高對(duì)稱點(diǎn)的坐標(biāo)。（共三對(duì)，表明沿著每對(duì)高對(duì)稱點(diǎn)的坐標(biāo)。（共三對(duì)，表明沿著gamma x w gamma計(jì)算能帶）計(jì)算能帶）如果是以卡笛爾坐標(biāo)系來(lái)寫(xiě)k點(diǎn)坐標(biāo)，則第四行以字母c開(kāi)頭。如何寫(xiě)出各種布拉菲格子對(duì)應(yīng)的倒格子中的高對(duì)稱點(diǎn)標(biāo)記和坐標(biāo)：如何寫(xiě)出各種布拉菲格子對(duì)應(yīng)的倒格子中的高對(duì)稱點(diǎn)標(biāo)記和坐標(biāo)：a)查空間群查空間群國(guó)際表國(guó)際表; b)用用xcrysden軟件畫(huà)出布里淵區(qū)，并選取高對(duì)稱點(diǎn)。軟件畫(huà)出布里淵區(qū)，并選取高對(duì)稱點(diǎn)。hou zhufeng12kpoints(3) 按按monkhorst-pack網(wǎng)格方法自動(dòng)產(chǎn)生k點(diǎn)automatic mesh0auto10 標(biāo)題或注釋行，無(wú)特別意義標(biāo)題

14、或注釋行，無(wú)特別意義0表示自動(dòng)產(chǎn)生表示自動(dòng)產(chǎn)生k點(diǎn)點(diǎn)以字母以字母a開(kāi)頭表示完全自動(dòng)以開(kāi)頭表示完全自動(dòng)以點(diǎn)為中心產(chǎn)生點(diǎn)為中心產(chǎn)生k點(diǎn)點(diǎn)長(zhǎng)度，間接確定了沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目長(zhǎng)度，間接確定了沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目 完全自動(dòng)的方式 (1)automatic mesh0gamma1010 10 0.0 0.0 0.0 完全自動(dòng)的方式 (2)標(biāo)題或注釋行，無(wú)特別意義標(biāo)題或注釋行，無(wú)特別意義0表示自動(dòng)產(chǎn)生表示自動(dòng)產(chǎn)生k點(diǎn)點(diǎn)以字母以字母g開(kāi)頭表示按開(kāi)頭表示按m-p方法以方法以點(diǎn)為中心自動(dòng)產(chǎn)生點(diǎn)為中心自動(dòng)產(chǎn)生k點(diǎn)點(diǎn)確定沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目確定沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目對(duì)所按網(wǎng)格分割產(chǎn)生的

15、對(duì)所按網(wǎng)格分割產(chǎn)生的k點(diǎn)進(jìn)行平移的量（這里不平移）點(diǎn)進(jìn)行平移的量（這里不平移）第三行也可以字母m開(kāi)頭，則不以點(diǎn)為中心hou zhufeng13kpoints(4) 按按monkhorst-pack網(wǎng)格方法自動(dòng)產(chǎn)生k點(diǎn) 完全自動(dòng)的方式 (2) (續(xù))automatic mesh0monkhorst-pack1010 10 0.5 0.5 0.5 完全自動(dòng)的方式 (3)標(biāo)題或注釋行，無(wú)特別意義標(biāo)題或注釋行，無(wú)特別意義0表示自動(dòng)產(chǎn)生表示自動(dòng)產(chǎn)生k點(diǎn)點(diǎn)以字母以字母m開(kāi)頭表示按開(kāi)頭表示按m-p方法自動(dòng)產(chǎn)生方法自動(dòng)產(chǎn)生k點(diǎn)點(diǎn)確定沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目確定沿倒格子各方向上網(wǎng)格點(diǎn)的數(shù)目對(duì)所按網(wǎng)格分割產(chǎn)生

16、的對(duì)所按網(wǎng)格分割產(chǎn)生的k點(diǎn)進(jìn)行平移的量（這里不平移點(diǎn)進(jìn)行平移的量（這里不平移)如果第三行以字母g開(kāi)頭，則以點(diǎn)為中心產(chǎn)生k點(diǎn)hou zhufeng14kpoints(5) 對(duì)六角晶系的結(jié)構(gòu)，強(qiáng)烈推薦采用以對(duì)六角晶系的結(jié)構(gòu)，強(qiáng)烈推薦采用以為中心按為中心按m-p網(wǎng)格產(chǎn)生網(wǎng)格產(chǎn)生k點(diǎn)點(diǎn) k點(diǎn)網(wǎng)格大小點(diǎn)網(wǎng)格大小n1, n2, n3按倒格子基矢的長(zhǎng)度進(jìn)行合適的選取按倒格子基矢的長(zhǎng)度進(jìn)行合適的選取一般選取一般選取n1, n2和和n3為奇數(shù)，以便產(chǎn)生的為奇數(shù)，以便產(chǎn)生的k點(diǎn)包含了點(diǎn)包含了點(diǎn)點(diǎn) 設(shè)置不同的設(shè)置不同的n1,n2和和n3，看它們對(duì)總能的影響，以確定能保證一定，看它們對(duì)總能的影響，以確定能保證一定計(jì)算

17、精度的計(jì)算精度的k點(diǎn)數(shù)目點(diǎn)數(shù)目#/bin/shfor i in 3 5 7 9 11docat kpoints summarydonehou zhufeng15potcar 按按poscar文件中，每類原子的順序把各類原子相應(yīng)的贗文件中，每類原子的順序把各類原子相應(yīng)的贗勢(shì)合并到同一個(gè)文件勢(shì)合并到同一個(gè)文件(potcar)中。比如在中。比如在poscar中原中原子坐標(biāo)按子坐標(biāo)按si, o的順序?qū)懗鼋o原子的坐標(biāo)，則到贗勢(shì)庫(kù)目錄的順序?qū)懗鼋o原子的坐標(biāo)，則到贗勢(shì)庫(kù)目錄中使用下面的命令中使用下面的命令: zcat si/potcar.z potcar zcat o/potcar.z potcar 每類原

18、子的贗勢(shì)類型每類原子的贗勢(shì)類型(paw或超軟或超軟)一致一致 每類原子的贗勢(shì)類型每類原子的贗勢(shì)類型(交換關(guān)聯(lián)交換關(guān)聯(lián))要與要與incar中交換關(guān)聯(lián)設(shè)置中交換關(guān)聯(lián)設(shè)置一致一致hou zhufeng16輸出文件輸出文件hou zhufeng17outcaroutcar文件包含了vasp計(jì)算后得到的絕大部分結(jié)果，每步迭代的詳細(xì)情況。下面介紹如何從outcar取出一些有用的信息：查看所計(jì)算體系的體積，使用下面的命令grep volume outcar得到的結(jié)果如下volume/ion in a,a.u. = 32.92 222.17volume of cell : 65.84第一行給出體系的體積分別以

19、3/atom, a.u.3/atom為單位給出的。第二行給出體系的體積是以3/unit cell為單位給出的。查看所計(jì)算體系的總能，使用下面的命令當(dāng)ismear = -5時(shí)，free energy toten是與energy without entropy是相等，則用grep toten outcar得到結(jié)果如下free energy toten = -7.910804 ev當(dāng)ismear等于其他的值時(shí)，free energy toten是與energy without entropy是不相等，則用grep entropy= outcar得到結(jié)果如下energy without entropy

20、= -7.910804 energy(sigma-0) = -7.910804在計(jì)算體系的結(jié)合能時(shí)，體系的總能取為energy without entropy后面的值。hou zhufeng18outcar查看所計(jì)算體系的費(fèi)米能級(jí)，使用下面的命令grep fermi outcar | tail -1得到的結(jié)果為bzints: fermi energy: 6.171330; 20.000000 electrons上一行中第一個(gè)數(shù)就是體系的費(fèi)米能級(jí)，第二個(gè)數(shù)就是體系的總價(jià)電子數(shù)。注釋：對(duì)半導(dǎo)體的體系，vasp取價(jià)帶頂作為費(fèi)米能級(jí)。對(duì)呈現(xiàn)金屬性的體系，費(fèi)米能級(jí)就是該體系的真實(shí)（具有物理意義的）費(fèi)米能

21、級(jí)。查看所計(jì)算體系的倒格子基矢在采用vi對(duì)outcar編輯時(shí)，用下面的命令來(lái)查找g/reciprocal lattice vectors 或 g/recip查看所計(jì)算體系中原子的受力情況在采用vi對(duì)outcar編輯時(shí)，用下面的命令來(lái)查找g/total-force原子所受的力的單位是ev/angstrom。hou zhufeng19doscar對(duì)非自旋極化的，設(shè)置了對(duì)非自旋極化的，設(shè)置了rwigs和和lorbit=1或或2，11或或12(只對(duì)只對(duì)paw)2 2 1 0 0.1891678e+02 0.3768172e-09 0.3768172e-09 0.3768172e-09 0.500000

22、0e-15 1.00000000000000d-004 car lan 7.84401537 -16.51100092 301 1.85674053 1.00000000 -16.511 0.0000e+00 0.0000e+00 -16.430 0.0000e+00 0.0000e+00 .eminemax在在emin和和emax之間的點(diǎn)的數(shù)目之間的點(diǎn)的數(shù)目費(fèi)米能級(jí)費(fèi)米能級(jí)能量能量態(tài)密度態(tài)密度態(tài)密度積分值態(tài)密度積分值先是總的態(tài)密度先是總的態(tài)密度301行，接著是各個(gè)原子的球諧函數(shù)投影或分解的態(tài)密度行，接著是各個(gè)原子的球諧函數(shù)投影或分解的態(tài)密度(每個(gè)原子的也是每個(gè)原子的也是301行行).即即en

23、ergy s py pz px dxy dyz dz2 dxz dx2 tot對(duì)自旋極化的情況，設(shè)置了對(duì)自旋極化的情況，設(shè)置了rwigs和和lorbit，可以同類推：，可以同類推：energy s(up) s(down), 這些可以數(shù)據(jù)可以按各個(gè)原子取出，畫(huà)分波態(tài)密度，或者直接通過(guò)這些可以數(shù)據(jù)可以按各個(gè)原子取出，畫(huà)分波態(tài)密度，或者直接通過(guò)p4vasp來(lái)顯示。來(lái)顯示。hou zhufeng20doscardoscar給出了體系的電子態(tài)密度給出了體系的電子態(tài)密度(單位為單位為sates/ev.unit cell)和態(tài)密度和態(tài)密度積分的值積分的值(也就是電子數(shù)也就是電子數(shù))。在計(jì)算時(shí)，需在在計(jì)算時(shí)，

24、需在incar文件中設(shè)置了文件中設(shè)置了rwigs以及以及l(fā)orbit(可選的可選的)只設(shè)置了只設(shè)置了rwigs的值的值: 對(duì)非自旋極化對(duì)非自旋極化2 2 1 0 0.1891678e+02 0.3768172e-09 0.3768172e-09 0.3768172e-09 0.5000000e-15 1.00000000000000d-004 car lan 7.84401537 -16.51100092 301 1.85674053 1.00000000 -16.511 0.0000e+00 0.0000e+00 -16.430 0.0000e+00 0.0000e+00 .eminemax

25、在在emin和和emax之間的點(diǎn)的數(shù)目之間的點(diǎn)的數(shù)目費(fèi)米能級(jí)費(fèi)米能級(jí)能量能量態(tài)密度態(tài)密度態(tài)密度積分值態(tài)密度積分值先是總的態(tài)密度先是總的態(tài)密度301行，接著是各個(gè)原子的行，接著是各個(gè)原子的s, p, d投影或分解的態(tài)密度投影或分解的態(tài)密度(每個(gè)原子的也是每個(gè)原子的也是301行行).即即energy s-dos p-dos d-dos對(duì)自旋極化的情況，只設(shè)置了對(duì)自旋極化的情況，只設(shè)置了rwigs：energy dos(up) dos(dwn) integrated dos(up) integrated dos(dwn)energy s-dos(up) s-dos(down) p-dos(up) p

26、-dos(dwn) d-dos(up) d-dos(dwn)hou zhufeng21chg和chgcar電荷密度文件，兩者的格式、內(nèi)容一樣。給出了原胞基矢和原子位置，以及電荷密度文件，兩者的格式、內(nèi)容一樣。給出了原胞基矢和原子位置，以及電荷密度電荷密度(r)*v的值等的值等. 當(dāng)當(dāng)lcharg=.false.時(shí)，時(shí)，chg*文件中無(wú)數(shù)據(jù)文件中無(wú)數(shù)據(jù)lan 5.32900000000000 0.000000 0.500000 0.500000 0.500000 0.000000 0.500000 0.500000 0.500000 0.000000 1 1direct 0.000000 0.0

27、00000 0.000000 0.500000 0.500000 0.500000 56 56 56 0.48995073152e+02 0.49141986100e+02 0.49542734834e+02 0.50083226401e+02 0.50592000474e+02 0.50865422427e+02 0.50697616547e+02 0.49910660506e+02 0.48380553212e+02 0.46055231054e+02 0.42962328993e+02 0.39206114637e+02 0.34954657491e+02 0.30419549263e+

28、02 0.25831206134e+02這一部分的數(shù)據(jù)和這一部分的數(shù)據(jù)和poscar的相同的相同fft 3維網(wǎng)格的值維網(wǎng)格的值電荷密度的值，按下面的格式寫(xiě)出來(lái)的：電荷密度的值，按下面的格式寫(xiě)出來(lái)的：用用vaspview軟件可以直接顯示它軟件可以直接顯示它hou zhufeng22eigenvalue靜態(tài)或遲豫計(jì)算時(shí)，在計(jì)算最后輸出的所有k點(diǎn)的kohn-sham本征值在md計(jì)算時(shí)，輸出的是與contcar中的原子結(jié)構(gòu)一致的預(yù)測(cè)波函數(shù)。 2 2 1 1 0.1891678e+02 0.3768172e-09 0.3768172e-09 0.3768172e-09 0.5000000e-15 1.00000000000000d-004 car lan 14 91 11 0.5000000e+00 0.0000000e+00 0.5000000e+00 0.1098901e-01 1 -13.8002 2 -12.9600 3 -12.9600 4 -10.1