VASP輸入參數含義「速查表」

2021-01-14 材料基因組高通量材料計算



共有四個,分別為:POSCAR、INCAR、KPOINTS和POTCAR。其中POSCAR為結構文件,KPOINTS為K點文件,POTCAR是贗勢文件。INCAR是VASP軟體的「C位「 輸入文件,它直接決定了做什麼計算以及如何做此次計算。


是否理解INCAR中參數含義,將直接決定:

☎  我們能否正確設置參數從而得到合理的計算結果

☎  在計算出現問題時能否向正確的方向調整參數,而不是向錯誤的方向越走越遠,最終浪費了大量的計算資源和我們最寶貴的時間。


INCAR文件中的參數眾多,對於VASP軟體使用新手來說,學習和記憶的難度都比較大。本文系統整理了INCAR中的重要參數(文章參數按排序包含了電子步、離子步參數、總體精度、動能截斷、自旋極化、範德華力修正、LDA+U、電子態密度和能帶等),涵蓋:

☎  參數作用

☎  參數值及釋義

☎  參數適用性

☎  應用參數時的注意事項等


後續我們將分享更多分子模擬學習資源,您可以關注我們的公眾號,獲得最新資源。

控制如何生成初始波函數


0:隨機生成波函數

1:從WAVECAR中讀取波函數,當讀取失敗時,會隨機生成

控制如何生成初始電荷密度


0:讀取WAVECAR生成電荷密度,當讀取失敗時,通過原子電荷密度疊加生成初始電荷密度

1:從CHGCAR中的設定讀取電荷密度

2:通過原子電荷密度的疊加方法,生成初始的猜測電荷密度

11:從CHGCAR讀取電荷密度,且在自洽循環中保持不變,適用於非自洽循環計算,如電子態密度的計算

設置自洽迭代過程優化電子波函數的算法


N:DAV算法,收斂性好,但速度慢

V:RMM算法,收斂性差,但速度快(比N快2~3倍)

F:以上算法的結合,綜合表現與V類似

Tips:結構偏離穩定結構較多的,建議用N,接近穩定結構的,可以用V。

自洽迭代循環收斂標準,單位為eV


自洽迭代過程緊鄰兩次迭代的系統能量的差與此標準相比。一般設置為1E-5~1E-6,意為兩次能量差小於設定的收斂標準,則自洽迭代結束,判定體系收斂。


自洽迭代循環的最大次數,通常設置為100


 有哪些應用注意事項?

點擊下方空白區域查看答案

過渡金屬體系、帶真空層的體系,設置了自旋極化的體系,此數值需要增加到200,300等。如果仍不能收斂,不建議繼續提高此數值,需要嘗試調整其他參數達到收斂的目標。


離子步弛豫(離子位置優化)算法


0:分子動力學模擬

1:準牛頓法

2:共軛梯度法

5:振動頻率計算

6:彈性常數計算


Tips:推薦設置值為2,如果初始結構和最終穩定結構接近,可以設置為1。

離子步弛豫收斂標準(包含能量標準和力的標準)


正值為系統能量變化(單位為eV),負值為原子上殘餘力(單位為:eV/埃)。多數情況用力的收斂標準判斷離子步弛豫是否收斂


Tips:三維結構可以收斂到-0.01~-0.03 eV/埃,低維體系可以收斂到-0.03~0.05eV/埃。

最大的離子步的數目 or 分子動力學模擬步數


當IBRION=1和2時,NSW代表最大的離子步數目

當IBRION=0時,NSW代表分子動力學模擬步數

離子步弛豫控制


參數說明如下圖:

3:全弛豫,適用於絕大多數的三維結構優化

2:為固定體積的弛豫,適用於如表面結構優化,微量摻雜體系優化。

4:固定體積,但是形狀可變的弛豫

Tips:複雜的結構優化過程,可以分步優化,即初始只優化離子(ISIF=2)→優化晶胞(ISIF=6)→優化離子和晶胞(ISIF=3)。

總體計算精度控制


受PREC影響的參數有四類(ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT)

PREC=Low | Medium | High | Normal | Accurate | Single。


下圖總結了PREC的設置值與其控制的其他參數值的關係。


上圖如何解讀?

以ENCUT為例,當INCAR文件中未設置ENCUT參數值,此時如果將PREC設置為Normal,則ENCUT將被自動設置為POTCAR(贗勢文件)中最大的ENMAX對應的數值。

平面波截斷能,控制用多少平面波展開波函數


Tips:推薦設置為POTCAR(贗勢文件)中參數ENMAX值的1.0-1.3倍,從而兼顧計算的精度與效率。

設置自洽迭代循環中新舊電荷密度的混合比例,默認值為0.4


自洽迭代不收斂該如何調整此參數?

點擊下方空白區域查看答案

自洽迭代不收斂時,可以降低此數值,即降低新的電荷密度混入的比例。


設置展寬方法,決定如何進行電子軌道分數佔據


默認值為1

值為:-5 -4 -3 -2 0 N(正值)

0:高斯方法(高斯展寬),適用於導體、半導體和絕緣體,或不清楚體系的導電性質時使用。展寬由SIGMA確定。

-5:四面體方法,適用於半導體和絕緣體


Tips:明確的半導體結構使用-5,其他絕大部分情況使用0。注意設置-5時K點網格個數要大於5。

展寬的寬度,與ISMEAR聯用,單位為eV,默認值為0.2


是否考慮電子自旋極化。默認值為1


1:非自旋極化,即每個軌道上自旋向上和自旋向下的電子數量相等。適用於非磁性體系。

2:自旋極化。適用於磁性體系,如鐵磁,反鐵磁材料。

對於磁性材料,設置初始原子磁矩


默認值是每個原子的磁矩均為1

對於含有d、f軌道的原子,需根據原子核外電子排布情況設置相應的數值

可以設置每個原子初始磁矩方向,正值為自旋向上,負值為自旋向下。

Hubbard U模型,描述電子強關聯能


參數組為:LDAU、LDAUL、LDAUU、LDAUJ

LDAU=.T.   使用LDA+U功能

LDAUL=-1/1/2/3 分別對應不加U和p、d、f軌道加U

LDAUU、LDAUJ 分別設置U和J值


Tips:LDA、GGA等交換相關泛函通常會低估帶隙,通過設置此參數可以得到更接近於真實值的帶隙計算結果。

範德華力計算方法,在DFT能量計算基礎上增加範德華力修正


10:DFT-D2方法

11:DFT-D3方法

推薦首選更新的DFT-D3方法。

 哪些體系需要使用此參數呢?

點擊下方空白區域查看答案

在計算表面吸附(物理吸附);弱相互作用佔體系能量比例較大的體系,如分子晶體,層狀結構體系時,要使用此參數。

範德華力計算方法

與IVDW的計算方法不同,此種方法將修改DFT內的能量計算。包含vdW-DF,optPBE-vdW, optB88-vdW, optB88b-vdW, vdw-DF2。


指定體系總電子數目


通常不用修改,默認值結合POSCAR和POTCAR自動計算體系價電子數。

若要實現體系帶電,則可控制此參數實現,在總電子數目基礎上減少n個電子,則體系帶n個電荷正電,反之是帶負電。


Tips:無法將電荷指定在某個特定原子上

控制是否輸出波函數和電荷密度


.T.: 寫波函數和電荷密度

.F.: 不寫波函數和電荷密度

態密度數據點個數,通常設置1000-3000


總態密度投影


10:投影到s、p、d、f軌道

11:投影到s、px、py、pz.....軌道

原子Wigner-Seitz半徑,控制分態密度強度


總能帶數

通常計算無需修改參數。

Tips:特別注意,光學性質計算設置為默認值的3~4倍



這個VASP參數分享節約了我一個月的時間快快分享給我的師妹哇咔咔,2020年6月30日,開始學習vasp我只是要發個彈幕,我是彈幕我驕傲果斷加關注,連刷有VASP學習問題組個群一起討論啊良心up主貌似這些VASP的參數我都看懂了呢期待更多資源謝謝分享

相關焦點

  • 教程| 初學者如何選擇合適的機器學習算法(附速查表)
    機器學習算法速查表機器學習算法速查表可幫助你從大量算法之中篩選出解決你的特定問題的算法,同時本文也將介紹如何使用該速查表。由於該速查表專門針對數據科學和機器學習的初學者,所以在探討這些算法之時,我們做了一些簡化的假設。本文中所有推薦的算法均來自於程序編譯反饋以及若干個數據科學家、機器學習專家和開發者的建議。
  • 前端開發相關速查表Cheatsheets整理集合
    當然寫多了以後自然而然的查 API 文檔的次數就會越來越少。但是當剛進一個領域再回來後發現自己又是什麼都不懂。聰明的前端開發工程在乎的是寫程式的思維,而不API文檔。因為記住所有API幾乎是不可能的。速查表(Cheatsheets)裡頭存放著濃縮過後的精華,當你是個有經驗的人,使用 速查表(Cheatsheets)能夠幫自己快速的複習,並且用最短的時間讓自己恢復到開發狀態,因此常常會去搜集各個前端開發相關速查表(Cheatsheets)有利於我們快速開發,節省一些不必要的時間浪費。
  • VASP主要參數詳解
    # 輸出電子局域函數LVTOT =     # 是否把總局域勢寫入LOCPOTVASP程序的每個參數都有多個選項,針對具體的計算體系以及使用情景需要給出合理的設置。比較方便的查詢方式如下:https://www.vasp.at/wiki/index.php/The_VASP_Manual。歡迎轉發!
  • vasp入門到精通[23]計算團簇時候參數怎麼設置
    使用vasp計算團簇小技巧1建立模型,如果是單元素的模型,要看多少原子,原子多了,就有很多種對稱性的可能,vasp無法做到搜索出基態
  • 史上最簡單的VASP安裝教程-非虛擬機
    /vasp.5.4.4-ok-6-201901.tar.bz2所在的目錄使用root執行:tar -xvf intel-all.tar.gz -C /opt/mkdir /opt/softwaretar -xvf vasp.5.4.4-ok-6-201901
  • VASP實用教程-從入門到進階-第一彈:VASP簡介及計算前期準備工作
    下面先簡單的介紹一下VASP軟體的基本原理以及一些重要參數的含義。VASP是一款基於第一性原理的軟體包,基本原理就是近似求解薛丁格方程,進而得到體系的電子態和能量。VASP是基於贗勢平面波基組的密度泛函理論(DensityFunctional Theory,DFT)的程序,通過VASP進行模擬計算的過程就是在DFT框架內求解Kohn-Sham方程。
  • 【2016最新】中國部分省市光伏電站最佳安裝傾角及發電量速查表
    【2016最新】中國部分省市光伏電站最佳安裝傾角及發電量速查表北極星太陽能光伏網訊:說明:(1)、速查表中發電量的計算已考慮79%的系統效率;(2)、速查表已根據當地經緯度換算出組件的最佳安裝傾角;(3)、速算表中的每瓦年發電量與電站實際裝機容量的乘積就是該電站的年發電量;1.北京市太陽能發電量統計表(晶矽組件)
  • VASP教程之Phonopy聲子譜計算(Pymatflow篇)
    使用matflow,我們可以方便地用一行命令完成整個計算任務的輸入文件的生成提交。下面我們將以NaCl固體聲子譜的計算為例,展示如何使用matflow完成此項工作。為了完成上述計算任務,只需使用Pymatflow一條命令即可:~$ mflow vasp -r 7 --cif
  • 常用螺紋尺寸參數簡明速查表
    一個常用螺紋尺寸參數簡明速查表,把外螺紋大徑和內螺紋小徑用一個圖表來表示,清楚,明了,簡單易懂,方便查詢。常見螺絲:<一>頭部剖溝,螺絲刀咀型:a. Slotted: 一字螺絲種類2b.***常用螺紋尺寸參數簡明速查表:
  • Learn VASP The Hard Way (Ex-0)
    );3.3) Atom:https://atom.io/ (這個Atom是文本編輯工具,不是原子的意思,Windows,Linux均可使用,要求會輸入文件,保存文件。4.2) p4vasp 網址:http://www.p4vasp.atLinux,Windows均可安裝,Linux下功能更加強大,Linux下安裝(Ubuntu系統): sudo apt-get install p4vasp  (一行命令搞定!!!)
  • 伺服驅動器一線維修速查手冊
    適讀人群 :電氣自動化工程技術人員,尤其是電氣設備維修保養人員等《伺服驅動器一線維修速查手冊
  • VASP5.4.1及+VTST編譯安裝
    官網:https://www.vasp.atVTST是VASP的過渡態工具,下載地址:http://theory.cm.utexas.edu/vtsttools/index.htmlVASP及VTST安裝均無需root權限,普通用戶在自己的目錄下也可以安裝使用
  • S參數的含義
    一般地,對於一個網絡有Y、Z和S參數可用來測量和分析,Y稱為導納參數,Z稱為阻抗參數,S稱為散射參數;前兩個參數主要用於集總電路,Z和Y參數對於集總參數電路分析非常有效,各參數可以很方便的測試;但是在微波系統中,由於確定非TEM波電壓、電流非常困難,而且在微波頻率測量電壓和電流也存在實際困難。
  • 常用貼片電阻阻值速查表
    可是很多人對貼片電阻的標識數據不是很了解,電阻小且不好測量,為解決部分人員對貼片電阻標識的不解,也為大家以後方便速查,本人通過各種電子書籍參考,特製作出該速查文檔用於電子愛好者速查貼片電阻阻值。希望大家能喜歡,本word文檔只在《數碼之家》論壇發表,轉載請說明出處。  下面列出了常用的5%和1%精度貼片電阻的標稱值和換算值,僅供大家使用時參考。
  • 使用vasp進行分子動力學模擬
    vasp做分子動力學的好處,由於vasp是近些年開發的比較成熟的軟體,在做電子scf速度方面有較好的優勢。缺點:可選系綜太少。
  • 分析LED顯示屏控制軟體三種主要參數
    慧聰LED屏網報導       根據生產廠家的不同,其LED顯示屏的參數項數量也會有所不同,經過分析歸類,可分為基本參數、輔助參數以及核心參數三種。  第一、基本參數  基本參數是顯示屏的基礎參數,如果設置不正確,則不能通信或不顯示以及顯示不正常。基本參數包括顯示屏寬高度、控制卡地址、波特率、IP位址、埠號、MAC地址、子網掩碼、網關、刷新頻率以及移位時鐘頻率共10項。
  • 肥料施用量計算公式及使用問題速查表!
    具體到化肥品種,實物化肥用量則要通過下式推算:計劃作物施肥量(kg)=實物化肥用量(kg)×有效成分含量(%)   其中:   1.計劃產量所需養分總量kg)=(計劃產量/100)×每形成100kg產量所需養分數量(參照不同作物形成100kg經濟產量所需養分大致數量表)。計劃產量則是當地作物3年平均產品產量再增加10%-15%。
  • VASP態密度與分子軌道
    對舊版本vasp,建議分兩步進行:1.啟用對稱的自洽計算(ISYM = 2)2.在對稱關閉的情況下重新計算部分電荷密度(ISYM = 0)在步驟2,VASP編寫的相位因子通常只能用作定性測量軌道向原子球體的投影,具體的計算方式並不十分合理。除定性計算外,不應再使用LORBIT = 12。
  • VASP6.1.1電聲耦合計算方法筆記
    https://www.vasp.at/wiki/index.php/Electron-phonon_interactions_from_Monte-Carlo_sampling#Input 方法主要分為兩種,1. 蒙特卡洛抽樣2.
  • VASP計算非線性磁矩和磁各向異性能(自旋軌道耦合)小結
    注意:VASP的輸入輸出的磁矩和類自旋量都會按照這個SAXIS方向,包括INCAR中的 MAGMOM行,OUTCAR和PROCAR.文件中的總磁矩和局域磁矩,WA.VECAR中的類自旋軌道和CHGCAR中的磁性密度。