VASP高級教程[4]快速計算並繪製能帶圖

2021-01-18 計算凝聚態物理

本教程需要四個文件，

total-bnd.sh 

kp.py

k-points

sband.py （含自旋極化的）或者ns.py(不含自旋極化的)

 

1靜態計算後，將能帶計算的輸入文件做好（記得把CHGCAR複製到能帶計算的目錄），準備計算能帶。

2運行命令產生KPOINTS:

python2.7  kp.py >KPOINTS ，這裡根據自己需要在kp.py中修改高對稱點之間的K點個數(kp.py中「mesh=50」 這一行裡的數字就是高對沉澱之間的k點個數)，在k-popints 中修改高對稱點。


這時候會產生新的KPOINTS文件

   

這裡截圖沒有全顯示出來。

3 提交能帶計算

4 計算完成後運行total.sh提取數據。其中total.sh內容如下，如果是沒有自旋極化或者SOC的計算中，需要去掉下面標出來的一行。

5 運行total.sh:

sh  total.sh 

產生如下文件：

將這三個文件放到python運行文件目錄，打開spyder或者在terminal運行python命令


一般安裝了anaconda的都會有matplotlib庫，這裡使用matplotlib畫圖：

在anaconda中打開spyder（當然熟悉的讀者可以直接打開spyder）：



相關參數如下



運行sband.py, 這樣就完成了


在文件中找到圖片可以放到文章中使用，可以發現，能帶路徑之間的比例，小編也計算過了。

對於非自旋極化的能帶，運行ns.py就可以了

所有腳本小編已經分享在下面的百度雲盤中了，請自取

連結: https://pan.baidu.com/s/1igmN2l78ng-j3QxGInRCdw 提取碼: 9c8h 

相關焦點

  • VASP高級教程[6]計算並繪製三維能帶圖腳本分享
    VASP計算並繪製三維能帶圖腳本分享有時候大家需要就算某一能量區間的三維能帶圖，來了解某一區域的能帶色散情況，比如Dirac cone
  • VASP計算雜化能帶詳細步驟教程
    cp PRIMCELL.vasp POSCAR後，用 VASP-PBE優化結構。 KPATH.in能帶只針對於 PRIMITIVE CELL，缺少這一步，你可能得到錯誤的結果。如果有必要，比較 KPATH.in文件中的能帶路徑是否與在線能帶路徑產生工具 SeeK-Path產生的一致，包括比較 PRIMCELL.vasp和 HIGH_SYMMETRY_POINTS文件。
  • VASP實用教程-從入門到進階-第一彈:VASP簡介及計算前期準備工作
    本教程主要是為了幫助各位剛剛接觸VASP並且還沒有入門的朋友快速入門，少走彎路，筆者在學習VASP的時候比較艱難，初期沒有人帶，主要是從網上查找一些教程，但是網上的教程比較零散，而且每個人寫的教程都不相同，對於初學者來講很難分辨對錯，這樣的結果就導致嘗試了很多都沒有得到好的結果，容易產生負面情緒，喪失信心。
  • 深入分析能帶結構(十)-雜化泛函能帶
    深入分析能帶結構(二)-VASPKIT能帶圖計算深入分析能帶結構(三)-Origin畫能帶圖深入分析能帶結構(四)-高通量生成能帶圖pymatgen深入分析能帶結構(五)-魔角石墨烯的構建與其能帶結構淺析深入分析能帶結構(六)-對贗原子做投影，電子化合物(下)深入分析能帶結構(七)-VASPKIT
  • 深入分析能帶結構(三)-Origin畫能帶圖
    往期回顧： 深入分析能帶結構(一)，討論了能帶的成因，軌道作用的方式、強弱、維度對能帶結構的影響。深入分析能帶結構(二)-VASPKIT能帶圖計算，討論VASPKIT自動生成K點路徑的方式，和使用VASPKIT+VASP計算能帶結構圖的方法。
  • 深入分析能帶結構(十一)-自旋軌道耦合
    能帶圖計算深入分析能帶結構(三)-Origin畫能帶圖深入分析能帶結構(四)-高通量生成能帶圖pymatgen深入分析能帶結構(五)-魔角石墨烯的構建與其能帶結構淺析深入分析能帶結構(六)-對贗原子做投影，電子化合物(下)深入分析能帶結構(七)-VASPKIT計算有效質量深入分析能帶結構
  • VASP計算二維材料的載流子遷移率
    ，由eV換算為原子單位制下的能量CBM/27.21然後選取Γ-X方向上以Γ開始的4-8個點的數據畫能帶圖，如下:The time execution %T" >> time.infompirun -np 16 /opt/soft/strainvasp5.4.4/vasp.5.4.4/build/std/vasp > logdate "+02 Today's date is: %D. The time finish %T" >> time.infocp .
  • 史上最簡單的VASP安裝教程-非虛擬機
    本文是針對vasp初學者的安裝教程編譯器以及VASP都已編譯號直接解壓到系統中即可用，故不用配置其它的庫文件以及環境；本教程適用於任意平臺安裝centos7的伺服器以及pc機（若在其它linux發行版本中安裝請諮詢小編）。
  • VASP入門到精通:[10]VASP計算能帶(2)
    syml文件的寫入和能帶的計算這裡繼續上一節內容，就是關於syml的內容2第一行數字4，
  • VASP態密度、能帶、DFT+U實例分析教程
    今天主要記錄強關聯電子體系+U、DOS和能帶計算過程中參數的設置，相信大家在讀一些計算類文獻時會看到文章中計算部分有The GGA + U calculations were performed using the model proposed by Dudarev et al這樣一段話
  • VASP入門教程三:輸入輸出文件介紹
    VASP入門教程三：輸入輸出文件介紹1.
  • vasp入門到精通[23]計算團簇時候參數怎麼設置
    使用vasp計算團簇小技巧1建立模型，如果是單元素的模型，要看多少原子，原子多了，就有很多種對稱性的可能，vasp無法做到搜索出基態
  • 使用腳本繪製吸附勢能面
    下面的文獻裡DOI：10.1063/1.4901055有一張圖非常生動，講的分別是Li在Li(001)以及Mg在Mg(0001)表面的吸附勢能分布，我們不僅可以知道哪些位置有利於吸附，還可以根據吸附能分布，描繪出分子在表面的擴散路徑，這有助於我們使用NEB方法算擴散能壘。其實思路很簡單，就是將表面網格化，算分子在格點的吸附能，再畫出等高線圖。
  • VASP教程之Phonopy聲子譜計算(Pymatflow篇)
    計算固體聲子譜常用的方法包括有限位移、DFPT方法等，本文將為您帶來使用matflow簡化VASP結合Phonopy計算固體聲子譜的內容。
  • VASP5.4.1及+VTST編譯安裝
    VASP採用平面波贗勢（或綴加投影波）方法進行從頭計算分子動力學模擬的軟體包。VASP.5.4.1編譯安裝：1.解壓縮VASP源碼包：tar xvf vasp.5.4.1.tar.bz2，解壓縮後的目錄為vasp.5.4.1。
  • Learn VASP The Hard Way (Ex-0)
    主動去做本節的練習或者重複其他人的計算。>:4.1) Material Studio 自行下載（百度裡面一大堆）,安裝（教程也是一大堆）。4.2) p4vasp 網址：http://www.p4vasp.atLinux,Windows均可安裝,Linux下功能更加強大,Linux下安裝(Ubuntu系統): sudo apt-get install p4vasp  (一行命令搞定!!!)
  • 特此推薦【免費】VASP零基礎特訓營
    ****************【免費】VASP零基礎特訓營，科研人必備的DFT理論計算10大專題75講直播課https://mp.weixin.qq.com/s/v9aYaRA0Ja1qhfi8UGb7Zg華算科技、呂梁雲計算中心、源資科技、計算材料學、催化開天地強強聯合，推出VASP零基礎直播課程，10大專題75講直播課，帶您快速完成VASP軟體計算入門
  • VASP主要參數詳解
    LORBIT = 10   # 是否創建PROCAR或PROOUT文件ICHARG = 11  # EMIN = XX     # 計算DOS的最小能量EMAX =  XX    # 計算DOS的最大能量NEDOS = 1000   # 計算DOS的點數NBANDS = XXX   #
  • 深入分析能帶結構(十二)-躍遷偶極矩計算
    能帶圖計算深入分析能帶結構(三)-Origin畫能帶圖深入分析能帶結構(四)-高通量生成能帶圖pymatgen深入分析能帶結構(五)-魔角石墨烯的構建與其能帶結構淺析深入分析能帶結構(六)-對贗原子做投影，電子化合物(下)深入分析能帶結構(七)-VASPKIT計算有效質量深入分析能帶結構