wannier90計算貝里曲率、反常霍爾電導

前言

假設你已知曉 如何介面VASP或QE， 使用 wannier90 計算得到幾個基本的輸入輸出檔案：

wannier90。win

wannier90。eig

wannier90。mmn

wannier90。wout

wannier90。amn

我們看一下如何在wannier90的部分計算貝里曲率以及反常霍爾電導

增加的引數

在wannier90。win中加入以下引數

berry = trueberry_task = ahcberry_kmesh = 40 40 70berry_curv_adpt_kmesh = 5 5 5fermi_energy = -4。9269bands_plot = Tberry_curv_unit = bohr2iprint = 2kpath = truekpath_task = bands+curv#kpath_bands_color = spinkpath_num_points=1000kslice = truekslice_task = curv+fermi_lineskslice_2dkmesh = 200kslice_corner = 0。0 0。0 0。0kslice_b1 = -0。6667 1。3333 0。0kslice_b2 = 1。0 0。0 0。0begin kpoint_pathG 0 0 0 M 0。5 0 0 M 0。5 0 0 K 0。3333333333333333 0。3333333333333333 0。000K 0。3333333333333333 0。3333333333333333 0。000 G 0 0 0end kpoint_path

postw90

計算完成後，我們需要使用postw90。x 進行計算

結果

生成如下相關的輸出檔案：

wannier90-kslice-curv_z+fermi_lines。pywannier90-kslice-curv_y+fermi_lines。pywannier90-kslice-curv_x+fermi_lines。pywannier90-kslice-curv。datwannier90-kslice-bands。datwannier90-kslice-coord。datwannier90-bands+curvwannier90-curv。dat

使用對應的python指令碼繪圖：

上面的案例是我使用的一個很隨意的，自己研究的例子。後面會系統跟進wannier90系列，希望讀者可以多多點贊、分享。