1、.vasp參數(shù)設置詳解計算材料 2010-11-30 20:11:32 閱讀197 評論0 字號：大中小 訂閱 轉自小木蟲，略有增減 軟件主要功能： 采用周期性邊界條件（或超原胞模型）處理原子、分子、團簇、納米線（或管）、薄膜、晶體、準晶和無定性材料，以及表面體系和固體l 計算材料的結構參數(shù)（鍵長、鍵角、晶格常數(shù)、原子位置等）和構型l 計算材料的狀態(tài)方程和力學性質（體彈性模量和彈性常數(shù)）l 計算材料的電子結構（能級、電荷密度分布、能帶、電子態(tài)密度和elf）l 計算材料的光學性質l 計算材料的磁學性質l 計算材料的晶格動力學性質（聲子譜等）l 表面體系的模擬（重構、表面態(tài)和stm模擬）l 從頭分

2、子動力學模擬l 計算材料的激發(fā)態(tài)（gw準粒子修正）計算主要的四個參數(shù)文件：incar ,poscar,potcar ,kpoints,下面簡要介紹，詳細權威的請參照手冊incar文件：該文件控制vasp進行何種性質的計算，并設置了計算方法中一些重要的參數(shù)，這些參數(shù)主要包括以下幾類：精品.l 對所計算的體系進行注釋：systeml 定義如何輸入或構造初始的電荷密度和波函數(shù)：istart，icharg，iniwavl 定義電子的優(yōu)化 平面波切斷動能和綴加電荷時的切斷值：encut，enaug 電子部分優(yōu)化的方法：algo，ialgo，ldiag 電荷密度混合的方法：imix，amix，amin，b

3、mix，amix_mag，bmix_mag，wc，inimix，mixpre，maxmix 自洽迭代步數(shù)和收斂標準：nelm，nelmin，nelmdl，ediffl 定義離子或原子的優(yōu)化 原子位置優(yōu)化的方法、移動的步長和步數(shù)：ibrion，nfree，potim，nsw 分子動力學相關參數(shù)：smass，tebeg，teend，pomass，nblock，kblock，pstress 離子弛豫收斂標準：ediffgl 定義態(tài)密度積分的方法和參數(shù) smearing方法和參數(shù)：ismear，sigma 計算態(tài)密度時能量范圍和點數(shù)：emin，emax，nedos 計算分波態(tài)密度的參數(shù)：rwigs，l

4、orbitl 其它 計算精度控制：prec 磁性計算：ispin，magmom，nupdown 交換關聯(lián)函數(shù)：gga，voskown 計算elf和總的局域勢：lelf，lvtot 結構優(yōu)化參數(shù)：isif 等等。 主要參數(shù)說明如下：? system：該輸入文件所要執(zhí)行的任務的名字。取值：字符串，缺省值：system精品. ? nwrite：輸出內(nèi)容詳細程度。取值：04，缺省值：2如果是做長時間動力學計算的話，最好選0或1(首末步/每步核運動輸出)，短時運算用2，選3則會在出錯的時候給出說明信息。 ? istart：決定是否讀取wavecar文件。取值：02，缺省0/1 for無/有前次計算的wa

5、vecar（波函數(shù)）0：begin from scratch，根據(jù)iniwav初始化波函數(shù)1：restart with constant energy cut-off，從wavecar讀取波函數(shù)（重定義平面波集）2：restart with constant basis set，從wavecar讀取波函數(shù)（平面波集不變） ? icharg：決定如何建立初始電荷密度。取值：02，缺省值: if istart=0 2 else 00：由初始波函數(shù)計算電荷密度1：從chgcar文件讀取電荷密度2：使用原子電荷密度的疊加+10非自洽計算 ? ispin：是否進行spin polarized calcu

6、lation。取值：1，2（1-no，2-yes），缺省值：2 ? magmom：在icharg=2或在chgcar中未包含磁化密度（icharg=1）時，指定每個原子的初始磁化時刻。取值：實數(shù)數(shù)組，缺省值:對ispin=2 nions*1.0，對非共線型磁化體系3*nions*1.0 精品.? iniwav如何設置初始波函數(shù)，只在istart=0時使用。取值：0，1（0-最低動能的平面波，1-隨機數(shù)），缺省值：1。 ? idipol控制計算單極、偶極和四極修正。取值：14。13只計算第一/二/三晶矢方向，適于厚板（slab）的計算4 所有方向都計算，適于計算孤立分子 ? prec進動（pre

7、cession）。取值：low/medium/high/normal/ accurate/single)，缺省值: normal（vasp.4.x）；medium（vasp.5.x）vasp4.5+采用了優(yōu)化的accurate來替代high，所以一般不推薦使用high。不過high可以確?！敖^對收斂”，作為參考值有時也是必要的。同樣受推薦的是normal，適于作為日常計算選項。受prec影響的參數(shù)有四類：encut；ngx，ngy，ngz；ngxf，ngyf，ngzf；ropt。如果設置了prec，這些參數(shù)就都不需要出現(xiàn)了，當然直接設置相應的參數(shù)也有同樣效果。具體影響效果見p5354。 ? e

8、ncut平面波基組的截斷能量（ev）。取值：實數(shù)，缺省值：受prec設置影響，從potcar文件中找出相應的enmax/enmin值來設置。prec = low medium accurate normalencut = enmin enmax enmax enmax single highenmax enmax*1.3對于多個元素不同的enmax/enmin，都取最大值。精品.該參數(shù)非常重要，最好不要手工去設置，除非文獻告訴你要用多少，或者經(jīng)過結果可靠性的驗證。當然，為了測試一下提交的任務，也不妨先設個較小的值。 ? ngx，ngy，ngz：控制fft網(wǎng)格在三個晶矢方向上的格點數(shù)量。? ng

9、fx，ngfy，ngfz：控制第二次更精確的fft網(wǎng)格的格點數(shù)量。也是兩類重要的最好不要去動的參數(shù)。在未指定的情況下將根據(jù)prec的設置從potcar中自動讀取。prec=high/accurate，基組中向量的2倍值，用來避免wrap around errors，得到精確解。prec=low/medium/normal，基組中向量的3/4倍值（已足夠精確到 1 mev/atom）。 ? lreal:決定投射是在實空間還是倒易空間進行。取值：.true.（實空間）/.false.（倒易空間），缺省值：.false.用于求解贗勢的非局域部分用到的一個積分，在倒格空間里采用平面波基組求解，在實空間

10、里則采用積分球求解。其他還有兩個選項：o or on，a or auto。on和.true.的差別在于是否使用king-smith算法優(yōu)化，設為auto則進行自動選擇，推薦使用。 ? ropt:在lreal=auto or on時，優(yōu)化控制每個核周圍的積分球內(nèi)的格點數(shù)。取值：實數(shù)數(shù)組for lreal=on精品. prec=low, 700 points in the real space sphere (ropt=0.67) prec=med, 1000 points in the real space sphere (ropt=1.0) prec=high, 1500 points in

11、the real space sphere (ropt=1.5)for lreal=auto prec=low, accuracy 10-2 (ropt=0.01) prec=med, accuracy 2*10-3 (ropt=0.002) prec=high accuracy 2*10-4 (ropt=2e-4) ? nelm, nelmin and nelmdl：控制電子自洽循環(huán)步數(shù)。取值：整數(shù)nelm：電子自洽循環(huán)最大次數(shù)。缺省值：60nelmin：電子自洽循環(huán)最小次數(shù)。缺省值：2nelmdl：弛豫次數(shù)。缺省值：if istart=0, iniwav=1, and ialgo=8，-5

12、，if istart=0, iniwav=1, and ialgo=48，-12，else 0nelmdl可以取負值。如果初始波函數(shù)采用隨機賦值，即istart=0，iniwav=1，那么很可能開始的值比較離譜，那么在第一步核運動循環(huán)之前采用nelmdl（負值）步的非自洽（保留初始的h）步計算將減少計算所需的時間。 ? ediff：指定電子自洽循環(huán)的全局中斷條件，用于控制收斂精度。取值：實數(shù)，缺省值：10-4注意，即使ediff=0，nelm步電子自洽循環(huán)也會執(zhí)行。 ? ediffg：指定離子弛豫循環(huán)的中斷條件，用于控制核運動的收斂精度。取值，實數(shù)，缺省值：10*ediff精品.ediffg0

13、在兩個離子步的總自由能之差小于ediffg時停止ediffg=0 nos算法模擬正則系綜 ? npaco：成對相關函數(shù)的槽數(shù)。取值：整數(shù)，缺省值：256? apaco：成對相關函數(shù)求值中的最大距離（?）。取值：整數(shù)，缺省值：16簡單說就是在不超過apaco的npaco個距離上求成對相關函數(shù)pcf。 ? rwigs：給出wigner-seitz半徑，dos計算用。取值：實數(shù)數(shù)組，缺省值：從potcar文件中讀取 ? nelect：總電子數(shù)，如果系統(tǒng)不是電中性的就必須設置，所帶電荷作為均一的背景電子氣考慮。取值：實數(shù)，缺省值：-（價電子數(shù)），由poscar和potcar文件自動決定（通常不必給出）

14、。 ? nupdown上下自旋成分間的電子數(shù)之差。取值：整數(shù)，缺省值：未設置（此時將進行完全弛豫）精品. ? emin, emax：dos求值的最小/最大能量。取值：實數(shù)，缺省值：emin = - (lowest ks-eigenvalue -)，emax = - (highest ks-eigenvalue -) ? ismear：決定每個波函數(shù)的部分占位fnk如何設置。取值：-5 | -4 | -3 | -2 | 0 | n（-5-帶有blochl修正的四面體方法, -4-不帶blochl修正的四面體方法, -3-根據(jù)incar文件中提供的smearing參數(shù)執(zhí)行循環(huán)，-2-從waveca

15、r文件中讀取, -1-fermi-smearing, 0-gaussian smearing, 0-method of methfessel-paxton ordern，缺省值：1采用部分占位波函數(shù)，用一個函數(shù)來平滑積分，尤其是對于金屬體系可減少k點。 ? sigma：決定smearing的寬度（ev）。取值：實數(shù)，缺省值：0.2 ? algo：指定電子最小化算法。取值：normal(blocked davidson block iteration scheme)/veryfast(rmm-diis)/fast (前兩個算法的混合)/all(波函數(shù)的所有帶同時更新)/ damped(dampe

16、d velocity friction algorithm)，缺省值: normal ? ialgo：指定主算法（整數(shù)選擇算法）。取值：8(共軛梯度算法)/38(davidson block iteration scheme)/48(rmm-diis)，缺省值：38算法是最重要的參數(shù)之一。一般vasp推薦使用的是以上三種算法，一般來說8/38是初期比較快收斂，在接近平衡時采用48較快，在初期或md時使用48可能會遇到不收斂的情況。也可以使用algo參數(shù)來替代ialgo，設置fast，vasp會先用38，再自動切換到48。各種算法只要收斂，結果應該一致。另一個可能有用的選項是-1。不進行實際的計

17、算，只對重要的步驟做計算測試，并將測試得到的各部分耗時輸出到output中。 ? voskown：決定是否使用vwn插值算法。取值：0(不使用)/1(使用)，缺省值：0精品.如果使用了pw91泛函或需要計算磁性質時可以設為1使用。 ? mixing-tags：w imix：混合的類型，取值：整數(shù)，缺省值：4w amix：線性混合參數(shù)。取值：實數(shù)，缺省值：0.8（us-pp），0.4（paw）w amin：最小混合參數(shù)。取值：實數(shù)，缺省值：0.1w bmix：kerker混合方案的截斷波向量。取值：實數(shù)，缺省值：1.0w amix_mag：磁化過程的線性混合參數(shù)。取值：實數(shù)，缺省值：1.6w b

18、mix_mag：磁化過程的kerker混合方案的截斷波向量。取值：實數(shù)，缺省值：1.0w wc：broyden混合方案中每步的加權因子。取值：實數(shù)，缺省值：1000.0w inimix：broyden混合方案中的初始混合類型。取值：整數(shù)，缺省值：1w mixpre：broyden混合方案中的預處理類型。取值：整數(shù)，缺省值：1w maxmix：broyden混合器中存儲的最大步數(shù)。取值：整數(shù)，缺省值：-45值得注意的是，在md或者弛豫的時候，設置maxmix（0,一般約3倍的電子sc步數(shù)）可能會大大減少核運動步數(shù)，但同時也會增加對內(nèi)存的要求。 ? lwave,lcharg：決定是否把波函數(shù)（或電

19、荷密度）寫入外部文件wavecar（或chgcar和chg）中。取值：.true./.false.，缺省值：.true. ? lvtot：決定是否把總局域勢寫入外部文件locpot中。取值：.true./.false.，缺省值：. false. ? lelf：決定是否創(chuàng)建elfcar文件。取值：.true./.false.，缺省值：. false.elfcar用于保存elf(electron localization function)。 ? lorbit：和適當?shù)膔wigs一起決定是否創(chuàng)建procar或proout文件。取值：0 | 1 | 2 | 5 | 10 | 11 | 12，缺省值：

20、0(創(chuàng)建doscar和procar文件)精品. ? npar：用于控制vasp切換到帶間并行模式。取值：整數(shù)，缺省值：總節(jié)點數(shù)。每一個節(jié)點計算一個帶，可以提高并行效率，減少通訊量，但可能會大幅增加內(nèi)存的需求。 ? nbands：給出計算中實際的總能帶數(shù)。取值：整數(shù)，缺省值：nelect/2+nions/2（非自旋），0.6*nelect + nmag（自旋）nbands的取值應使計算中包含相當數(shù)量的空帶。因為計算需要大量的空帶，至少要求1個空帶（否則vasp會給出警告）。nbands對于解決內(nèi)存需求非常重要。一般nbands在nelect/2+nions/2以上可以得到較為精確的結果，但如果內(nèi)

21、存不夠就只能減少nbands，在犧牲精度和體系大小之間平衡了。potcar文件贗勢文件，最重要的輸入文件之一?？梢岳斫鉃榉肿恿W模擬中的力場文件，但包含的信息更多。vasp將各元素優(yōu)化的incar里的參數(shù)也包含在該文件中，作為支持prec的缺省選擇。通常各元素的potcar文件已經(jīng)包含在vasp軟件包中的贗勢庫里了，用戶需要做的是確定自己具體需要哪幾種贗勢，然后按照poscar文件里的順序，將所選擇的potcar文件按順序連接起來就可以了。如以下命令：cat file1 file2 file3 potcar 贗勢庫中的贗勢文件可以分為以下幾類：l 根據(jù)產(chǎn)生方法的不同有ultra-soft贗勢（

22、uspp）和投影擴充波贗勢（paw）；l 根據(jù)交換關聯(lián)函數(shù)的不同有l(wèi)da（local density approximation）和gga（generalized gradient approximation，又可以再分為pw91和pbe）；精品.l 根據(jù)半芯態(tài)的不同元素a的贗勢文件還可以分為a，a_sv和a_pv；l 根據(jù)選取的不同截斷能量（enmax）元素a的贗勢文件還可以分為a，a_s和a_h。 vasp的贗勢文件放在目錄/vasp/potentials下，該目錄又包含五個子目錄：pot，pot_gga，potpaw，potpaw_gga，potpaw_pbe，其中每個子目錄對應一種贗勢

23、形式：pot = pp, lda；pot_gga = pp, gga；potpaw = paw, lda；potpaw_gga = paw, gga, pw91；potpaw_pbe = paw, gga, pbe。每個目錄中每種元素根據(jù)截斷能量和半芯態(tài)的不同還會有多個對應的贗勢文件存在。在具體選取時可以參考各版本同目錄下的v_rhfin和psctr文件，這兩個文件說明了該版本的贗勢是如何生成的。 選好所要使用的贗勢之后，進入對應的目錄，會看到里面有四個文件：potcar.z，psctr.z，v_rhfin.z和ws_ftp.log?，F(xiàn)在需要用到的是第一個，把各元素的該文件解壓（zcat po

24、tcar.z file），然后用cp或mv命令把這些文件移到工作目錄里，再用cat命令把它們合并到potcar中，就得到了我們需要的potcar。注意要記住這里元素的排列順序，以后在poscar文件中各元素的排列就是按照這個順序。 poscar文件位置文件。描述所計算體系的晶胞參數(shù)、原子個數(shù)及晶胞中原子的位置，以及分子動力學計算時原子的初始速度。poscar文件示例：si-fcc5.430.0 0.50 0.500.50 0.00 0.50精品.0.50 0.50 0.002direct0.0 0.00 0.000.25 0.25 0.25 poscar文件的結構：第1行：任意文字注釋。第2行

25、：晶格常數(shù)（?），也是晶矢的縮放系數(shù)，后面所有長度值得自原值除以此值。a=b=c時取a即可，否則一般取三者最大值，若取負值，則為晶胞體積（?3）。第3-5行：定義晶矢。第6行：每種元素的原子個數(shù)，特別注意順序,要與下面的坐標順序以及potcar中的順序一致。第7行：可省略，無需空行。做動力學時，用于指定是否需要固定部分原子的坐標。若是，此行以s或者s作為首字母即可。第8行開始為原子的坐標位置，格式為option linecoordinate1 of element1coordinate2 of element1.coordinaten of element1option linecoordin

26、ate1 of element2coordinate2 of element2.coordinatem of element2.精品.其中，option line指定輸入坐標的格式，除了第一個以外，如果后面的輸入格式同前，則都可以無空行省略。option line可指定的輸入坐標格式有兩種：dordfor direct modecorcorkorkfor cartesian mode顧名思義，前者是定義在三個晶矢方向上的坐標：r=r1x+ r2y +r3z，r1,r2,r3為前面的晶矢，x,y,z為輸入的三個坐標，r為原子坐標位矢，是相對于晶矢給出的，而后者則是以笛卡爾坐標系來給出原子的絕對坐

27、標，實際就是將直角坐標除以前面第二行定義的晶胞常數(shù)后得到的坐標值。如果第7行設定為s(selective dynamic)，則可以用以下形式定義各坐標是否可以移動：selective dynamicscartesian0.00 0.00 0.00 t t f0.25 0.25 0.25 f f fkpoints文件設置布里淵區(qū)k點網(wǎng)格取樣大小或能帶結構計算時沿高對稱方向的k點。一般來說，k點越密越多，計算精度也就越高，但計算成本也會越高。因此，對于原子數(shù)較多的體系的計算，需要謹慎的嘗試k點數(shù)目，在避免或者預先評估wrap-around error的前提下盡量減少k點數(shù)目。 kpoints文件示

28、例：k-points along high symmetry lines11reciprocal0.00 0.00 0.00 1.000.05 0.00 0.05 1.00精品.0.05 0.00 0.50 1.00 kpoints文件結構：第1行：注釋行。第2行：k點總數(shù)，或者0（自動生成網(wǎng)格（automatic k-mesh generation）。如果是前者，給出k點總數(shù)，又分為兩種情況：m：全手動（entering all k-points explicitly），手動輸入即自定義各個k點的坐標和權重，推薦只在能帶計算時使用。此時：第3行：輸入格式標識。直角坐標（cartesian）或

29、倒易坐標（reciprocal）。同樣的，c、c、k、kfor cartesian，其他首字母則自動切換到reciprocal。第4n行:逐個k點的坐標描述。格式為“x y z w”。其中“x y z”是該點的三個坐標，w是權重。所有k點的權重相互之間的比例對了就行，vasp會自動歸一的。注意c坐標和r坐標的定義：c：k=(2/a)(x y z)r：k=x*b1+y*b2+z*b3（b1-b3為倒易晶矢）這里“x y z”只是代表了坐標的順序，與坐標軸無關。 一般如非必要，可以先用自動模式生成k點，vasp會自動生成一個簡約化后的k點矩陣，存于ibzkpt文件，可以直接復制里面的數(shù)據(jù)到kpoi

30、nts文件中使用，這也是該輸入法的主要用途，可以減少重復自動生成格點的時間。另一個用途是為了做精確的dos（density of status）計算，由于這類計算所需的k點數(shù)極大，通過全手動盡可能的優(yōu)化k點也就必需了。 l：半手動/線形模式（line-mode）：在計算能帶結構時使用，此時需要精確地選取k點，在指定的高對稱性方向上生成指定數(shù)目的k點。此時：第2行：指定兩點間生成的k點數(shù)，不同于全自動的總k點數(shù)。第2.5行：l表示是線形模式。精品.第3行：輸入格式標識，同前。第4n行：每行描述一個點，格式為“x y z”。每兩行的點連成一條直線，在兩點間生成指定數(shù)目的k點。每隔兩行間以空行區(qū)分，

31、例如：k-points along high symmetry lines10line-modereciprocal0 0 00.5 0.5 0 0.5 0.5 00.5 0.75 0.25 最常用的方法是讓程序自動生成k點網(wǎng)格，用戶只需定義網(wǎng)格取樣大小即可。此時：第2行：0第3行：生成k點的方法?！癮”：全自動模式，可以看作以點為圓心以l為半徑做圓，當然各晶矢不同時，相應的圓就被拉成了橢圓，以確保三個倒格晶矢方向上覆蓋的k點數(shù)都為l“g”：/gamma方法“m”：monkhorst-pack方法其它首字母：高級模式，就是用c坐標或r坐標直接輸入新的晶矢，一般不使用第4行：在各個晶矢方向上分割

32、各晶矢的點數(shù)。一般為奇數(shù)，使得產(chǎn)生的k點是以gamma點為中心的。第5行：是否移動網(wǎng)格點以及移動多少，也即偏移原點的位矢，一般設成“0 0 0”，表示不移動。例如：automatic generation0精品.monhkorst-pack9 9 90.0 0.0 0.0 vasp的輸出文件主要包括：l outcar：包含運行vasp作業(yè)后得到的絕大部分計算結果以及每步迭代的詳細情況；l chg, chgcar：給出體系的電荷密度。兩個文件內(nèi)容相同，但前者的數(shù)據(jù)精度比后者略低一些；l doscar, eigenval：分別給出所計算體系的電子態(tài)密度和本征值。兩個文件中的能量值都是絕對的；l w

33、avecar：給出所計算體系的電子波函數(shù)，二進制文件，不可編輯；l oszicar：包含每次迭代或離子移動情況的簡單匯總；l contcar：給出離子進行弛豫時，每次移動后體系的晶格參數(shù)，與poscar內(nèi)容相同；l pcdat, xdatcar：給出有關分子動力學模擬中的一些結果，如配對相關函數(shù)；0人 | 分享到： 閱讀(197)|評論(0)|引用(0)|舉報精品.交流【活動or資源】vasp晶體結構優(yōu)化 youzhizhe(金幣+10): 活動獎勵。 2011-05-26 14:34:11youzhizhe:標題高亮 活動貼。 2011-05-26 14:34youzhizhe:版主推薦 2

34、011-05-26 15:03youzhizhe:單行區(qū)域置頂 2011-05-26 15:07賀儀:關閉主題 修改中！ 2011-05-26 15:19賀儀:打開主題 編輯，美化完畢！ 2011-05-26 15:49youzhizhe:版主推薦(有效期40天) 2011-05-26 16:40youzhizhe:退還金幣50個 2011-07-18 14:56精品.結構優(yōu)勢是我們做計算的的第一步，往往也是最重要的一步，一個不合理的結夠優(yōu)化可能導致錯誤的結果。在論壇里面也有很多蟲子們問關于結構優(yōu)化的問題，貌似大家對結構優(yōu)化不是很有把握。我在這里開貼只是拋磚引玉，希望大家積極交流。下面是我的結

35、構優(yōu)化的一些經(jīng)驗，希望對大家有點用第一步:可以用vesta打開cif文件直接保存為vasp格式就可以，三樓youzhizhe另一種方法：建模一般都是ms建好之后 修改*.cell文件得到poscar具體的方法參考站內(nèi)的信息（沒找到最開始的網(wǎng)頁，誰找到了貼上來）這是*.cell文件的內(nèi)容%block lattice_cart 7.365858463230930 0.816888331419153 1.518567634398170 0.000000000000000 7.749147189276010 0.1149693472816510.000000000000000 0.0000000000

36、00000 16.835000015679299%endblock lattice_cart%block positions_frac h 0.1045034131293120 0.0506482790444403 0.1503325580937660 h-0.2037023292911190 0.2015583952382200 0.1604502537744340 h-0.4269999607791240 0.3837162987940400 0.0980194481736503 h-0.3453870956286900 0.4791801646100701-0.0460018173909

37、388 h 0.2784615598004869 0.2515795170279700-0.1346530187895070 h 0.2619963944132780-0.5890935655390530 0.3725661174512691 h 0.2595685601772050-0.5519954246490220 0.1897443712411160 h 0.0767881789054980-0.6048204220967270 0.2666625899874510 h 0.2771409424864049-0.7522385903899170 0.2411634679797081 h

38、 0.5634931211089379-0.7084903644666670 0.2828144510356441 h 0.5814851393159790-0.5287740945641601 0.3433449582781520 h 0.5793331415735910-0.5003094920852541 0.2386888180011851 h 0.4946819829083929-0.2212541120652400 0.3755202742234630 %endblock positions_frac%block kpoints_list 0.0000000000000000 0.

39、0000000000000000 0.2500000000000000 1.000000000000000 %endblock kpoints_list%block cell_constraints 1 2 3把第一行%block lattice_cart改成system=grphene（你體系的名稱 可以隨便寫）然后后面填一行1再把%endblock lattice_cart%block positions_frac改為 h c n si(原子的類型)5242 2 2（原子的數(shù)目）direct （坐標表示的方式，對于這種方法，只能用d）然后把后面的原子符號h c 等去掉h-0.0000000

40、000000000-0.0000000000000001 0.7500000000000000c 0.1666666666666668 0.0833333333333332 0.7500000000000000把%endblock positions_frac以后的全部刪掉就行了保存為poscar 就可以了此時要注意保存的格式使用unix格式，ultredit可以選擇保存格式，如果用記事本編輯的話,所有的行前不能有空格（從*.cell里面復制過來的空格，自己加的沒關系）否則提示找不到原子，論壇里面有很多這樣的問題。第二步：開始結構優(yōu)化了，我們一般使用腳本來實現(xiàn)，這里使用的是侯老師所說的方法，體積與能量的曲線然后擬合去能量最小的體積進行進一步優(yōu)化。先講能量與體積的曲線怎么得來的。incar最主要