 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 947 | 回復(fù): 0 | |||
A寶祿計算銀蟲 (初入文壇)
|
[交流]
INCAR詳解2
|
|
>SYSTEM name of System 任務(wù)的名字*** >NWRITE verbosity write-flag (how much is written) 輸出內(nèi)容詳細(xì)程度0-3 缺省2 如果是做長時間動力學(xué)計算的話最好選0或1(首末步/每步核運(yùn)動輸出) 據(jù)說也可以結(jié)合shell的tail或grep命令手動輸出 >ISTART startjob: restart選項0-3 缺省0/1 for 無/有前次計算的WAVECAR(波函數(shù)) 1-3 只計算第一/二/三晶矢方向適于slab的計算4 全部計算尤其適于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate) Default: Medium VASP4.5+采用了優(yōu)化的accurate來替代high,所以一般不推薦使用high。不過high可以確保'絕對收斂'，作為參考值有時也是必要的。同樣受推薦的是normal，作為日常計算選項，可惜的是說明文檔提供的信息不足。受PREC影響的參數(shù)有四類：ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT 如果設(shè)置了PREC，這些參數(shù)就都不需要出現(xiàn)了當(dāng)然直接設(shè)置相應(yīng)的參數(shù)也是同樣效果的，這里不展開了，隨后詳釋 >ENCUT energy cutoff in eV : default taken from POTCAR-file important! 重要到幾乎最好不要手工去設(shè)置 除非文獻(xiàn)告訴你要用多少，或者經(jīng)過結(jié)果可靠性的驗證 當(dāng)然，為了測試一下提交的任務(wù)，也不妨先設(shè)個較小的值 附加說明： 當(dāng)且僅當(dāng)POTCAR里頭沒有設(shè)置ENCUT時(其實(shí)貌似沒有才是常態(tài))，才受PREC設(shè)置影 響從POTCAR里找出相應(yīng)的ENMAX/ENMIN值來設(shè)置。 PREC= Low Medium Accurate High ENCUT= ENMIN ENMAX ENMAX 130%ENMAX 對于多個元素的POTCAR不同的ENMAX/ENMIN，都取最大值 >NGFX, NGFY, NGFZ: FFT mesh for charges 也是兩類重要的最好不要去動的參數(shù)，PREC設(shè)置將從POTCAR中自動讀取。 PREC=High,Accurate 2倍值，用來避免wrap around errors得到精確解 PREC=Low,Medium,Normal 3/4 也已經(jīng)足夠精確到 1 meV/atom >LREAL: Default= .FALSE. 贗勢的非局域部分用到的一個積分在倒格空間或者實(shí)空間都可以求值。這個選項就 是決定是在哪個空間里求。在倒格空間里，采用平面波基組求解，在實(shí)空間里，采用積 分球求解。 缺省是.FALSE，即不在實(shí)空間求。但效率會低一些。 其他選項是 O or On,A or Auto 和.True.。 On和.TRUE.的差別在于是否使用King-Smith算法優(yōu)化，Auto則自動選擇，推薦。 >ROPT: 優(yōu)化控制每個核周圍的積分球內(nèi)的格點(diǎn)數(shù)，LREAL=Auto or On 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 the real space sphere ( ROPT=1.5) For LREAL=Auto PREC= Low accuracy 1e-2 ( ROPT=0.01) PREC= Med accuracy 2e-3 ( ROPT=0.002) PREC= High accuracy 2e-4 ( ROPT=2E-4) >NELM, NELMIN and NELMDL nr. of electronic steps Default 最大電子自洽循環(huán)次數(shù) NELM = 60 最小次數(shù) NELMIN = 2 弛豫次數(shù) NELMDL = -5 if ISTART=0, INIWAV=1, and IALGO=8 -12 if ISTART=0, INIWAV=1, and IALGO=48 0 else 如果初始的波函數(shù)采取隨機(jī)賦值，即ISTART=0, INIWAV=1，那么很可能開始的值比較離 譜，那么在第一步核運(yùn)動循環(huán)之前采用NELMDL(負(fù)值)步的非自洽(保留初始的H)步計算將減少計算所需的時間。 如果NELMDL取正值，將在每次核運(yùn)動之后附加指定次數(shù)的弛豫步，目前不知道可以干嘛 >EDIFF 電子SC循環(huán)的收斂精度 缺�。�1e-4 注意，即使EDIFF=0，NELM步也會執(zhí)行 >EDIFFG 核運(yùn)動的收斂精度 缺�。篍DIFF*10 (總能量) EDIFFG<0 則在所有的力都小于EDIFFG時停止 EDIFFG=0 則在NSW步后停止 此參數(shù)不支持MD，僅用于Relax >NSW 指定核運(yùn)動步數(shù) 缺省: 0 NBLOCK and KBLOCK inner block; outer block Default NBLOCK = 1 KBLOCK = NSW >IBRION ionic relaxation: -1-Fixed 0-MD 1-quasi-New 2-CG 3-Damp 5-freq Default if NSW=0 or 1 IBRION=-1 else IBRION=0 這個參數(shù)是和ISIF;IALGO/ALGO一起決定怎么算的最重要的參數(shù) 1-3 是三種Relax的方法，受ISIF決定是否固定核位置、晶胞大小和形狀 0 是標(biāo)準(zhǔn)的ab-initio MD，不受ISIF影響，即不改變晶胞大小和形狀 5 大概是和0差不多吧？支持Hessian和Freq(僅Г點(diǎn))的計算以及部分固定的MD 詳細(xì)的要在示例中具體情況具體分析了。 >ISIF calculate stress and what to relax Default if IBRION=0 (MD) 0 else 2 ISIF│calculate │ calculate │relax │ change │ change │ force │stress tensor │ions │ cell shape │cell volume ──┼─────┼───────┼───┼──────┼────── 0 │ yes │ no │yes │ no │no 1 │ yes │ trace only │yes │ no │no 2 │ yes │ yes │yes │ no │no 3 │ yes │ yes │yes │ yes │yes 4 │ yes │ yes │yes │ yes │no 5 │ yes │ yes │no │ yes │no 6 │ yes │ yes │no │ yes │yes 7 │ yes │ yes │no │ no │yes Trace only means that only the total pressure IWAVPR prediction of wf.: 0-non 1-charg 2-wave 3-comb Default if IBRION=0 (MD) 2 if IBRION=1,2 (relaxation) 1 else (static calculation) 0 以上選項保存TMPCAR +10 則全部使用內(nèi)存，不保存此文件 IWAVPR determines how wave functions and/or charge density are extrapolated from one ionic configuration to the next configuration. >ISYM symmetry: 0-nonsym 1-usesym 是否使用對稱性 Default 1 SYMPREC determines precision of the positions in POSCAR file. Default 1e-5 LCORR Harris-correction to forces. Default .TRUE. >POTIM time-step for ion-motion (fs) Default IBRION=0 (MD) no default,必須指定，MD每步步長 IBRION=1,2,3 (relaxation) 0.5 最小化的'scaling constant',尤其是IBRION=1 >TEBEG, TEEND temperature during run (MD有效) Default: TEBEG = 0 TEEND = TEBEG 注意VASP的溫度定義與實(shí)際溫度有細(xì)微的差別，所以 TEBEG=T×(N-1)/N T為實(shí)際溫度，N為原子數(shù) SMASS 控制MD中的速度模擬方法 default -3 微正則系綜(總自由能不變) -2 保持初速度不變 -1 每NBLOCK步調(diào)整速度，來保證動能連續(xù) >=0 Nosé算法模擬正則系綜，(不懂-,-) NPACO and APACO NPACO : number of slots for pair correlation (PC) function. Default 256 APACO : maximum distance for the evaluation of PC function in A. Default 16 簡單說就是在不超過APACO的NPACO個距離上求成對相關(guān)函PCF RWIGS Wigner-Seitz半徑 DOS計算用 >NELECT 總電子數(shù) 如果系統(tǒng)不是電中性的就必須設(shè)置，所帶電荷作為均一的背景電子氣考慮 NUPDOWN default不考慮電子自旋態(tài)改變的可能 EMIN, EMAX energy-range for DOSCAR file >ISMEAR part. occupancies: -5 tet with Blochl -4-tet -1-fermi 0-gaus >0 MP 采用所謂部分占有波函數(shù)，用一個函數(shù)來平滑積分，尤其是對于金屬體系可減少k點(diǎn) Default ISMEAR = 1 如果在KPOINTS里使用了tetrahedra方法 推薦ISMEAR=5 SIGMA determines the width of the smearing in eV Default SIGMA = 0.2 >ALGO algorithm: Normal (Davidson) | Fast (mixed)| Very_Fast (RMM-DIIS) >IALGO algorithm: use only 48 (RMM-DIIS) or 38(Davidson) or 8(CG) Default IALGO = 38 for VASP4.5 算法是最重要的參數(shù)之一。一般VASP推薦使用的是以上三種算法，一般來說8/38是初期 比較快收斂，在接近平衡時采用48較快，在初期或MD時使用48可能會遇到不收斂的情況 。也可以使用ALGO參數(shù)來替代IALGO，設(shè)置Fast，VASP會先用38，再自動切換到48。 各種算法只要收斂，結(jié)果應(yīng)該一致。 另一個可能有用的選項是-1。不進(jìn)行實(shí)際的計算，只對重要的步驟做計算測試，并將測 試得到的各部分耗時輸出在OUTPUT里。 VOSKOWN use VWN interpolation 算法，default 0 不用，如果使用了PW91或者需要計 算磁性質(zhì)，不妨設(shè)為1 用 mixing tag & MAXMIX IMIX = type of mixing AMIX = linear mixing parameter AMIN = minimal mixing parameter BMIX = cutoff wave vector for Kerker mixing scheme AMIX_MAG = linear mixing parameter for magnetization BMIX_MAG = cutoff wave vector for Kerker mixing scheme for mag. WC = weight factor for each step in Broyden mixing scheme INIMIX = type of initial mixing in Broyden mixing scheme MIXPRE = type of preconditioning in Broyden mixing scheme MAXMIX = maximum number steps stored in Broyden mixer Default for US-PP non-magnetic IMIX = 4 AMIX = 0.8 BMIX = 1.0 |
回復(fù)此樓» 猜你喜歡
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有240人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有3人回復(fù)
-
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有5人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 A寶祿計算 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +13 | 帆船哥 2026-03-04 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 我嘗試在小木蟲上招計算機(jī)專業(yè)的調(diào)劑學(xué)生 +4 | littlehu66 2026-03-04 | 6/300 |
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩，總分302求調(diào)劑 +4 | 五維天空 2026-03-04 | 9/450 |
|
|
|
[考博] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-05 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-06 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +9 | Faune 2026-03-06 | 10/500 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +4 | chen_2024 2026-03-02 | 4/200 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計算機(jī)、人工智能等專業(yè)的同學(xué)報考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +5 | Seele_v5 2026-03-04 | 5/250 |
|
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[論文投稿]
通訊作者寫誰，問題是你意想不到的問題
 15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
|
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|