两款费米面绘制软件ifermi和fermisurfer

pymatgen for parsing DFT calculation outputs.BoltzTrap2 for band structure interpolation.trimesh for manipulating isosurfaces.matplotlib, mayavi, and plotly for three-dimensional plotting.

pip install ifermi

sudo aptitude install libwxgtk3.0-dev

sudo aptitude install libwxgtk3.0-gtk3-dev

./configuremakesudo make install

ifermi info

Fermi Surface Summary=====================
 # surfaces: 5 Area: 32.745 ??2
Isosurfaces~~~~~~~~~~~
 Band Area [??2] Dimensionality Orientation ------ ------------ ---------------- ------------- 6 1.944 2D (0, 0, 1) 7 4.370 1D (0, 0, 1) 7 2.961 2D (0, 0, 1) 8 3.549 1D (0, 0, 1) 8 3.549 1D (0, 0, 1)

ifermi info --property velocity

Fermi Surface Summary=====================
 # surfaces: 5 Area: 32.75 ??2 Avg velocity: 9.131e+05 m/s
Isosurfaces~~~~~~~~~~~
 Band Area [??2] Velocity avg [m/s] Dimensionality Orientation ------ ------------ -------------------- ---------------- ------------- 6 1.944 7.178e+05 2D (0, 0, 1) 7 4.370 9.092e+05 quasi-2D (0, 0, 1) 7 2.961 5.880e+05 2D (0, 0, 1) 8 3.549 1.105e+06 quasi-2D (0, 0, 1) 8 3.549 1.105e+06 quasi-2D (0, 0, 1)

ifermi plot #默认绘图ifermi plot --filename my_vasprun.xml #指定文件ifermi plot -f my_vasprun.xml #指定文件ifermi plot --interpolation-factor 10 #插值 10为插值因子，默认是8ifermi plot -i 10 #插值 10为插值因子，默认是8ifermi plot --mu 1 #指定费米面，mu为化学势，默认值是0ifermi plot --property velocity #计算费米面群速度，单位m/sifermi plot --property spin #自旋轨道耦合或者非线性磁性费米面上的自旋，计算要求ISYM=-1ifermi plot --property spin --projection-axis 0 0 1 #表示自旋在001方向的投影ifermi plot --no-wigner #不显示Wigner–Seitz胞，默认显示Wigner–Seitz胞ifermi plot --type matplotlib #使用matplotlib绘制，可以用poltly或者mayavi具体设置见官方文档ifermi plot -t matplotlib #同上--type（-t）

ifermi plot --spin up #绘制上自旋 ifermi plot --spin down #绘制下自旋

ifermi plot --o fermi-surface.jpg #同上--output(-t)ifermi plot --output fermi-surface.html #输出.html格式

ifermi plot --azimuth 120 --elevation 5 #改变视角，方位角azimuth和仰角elevationifermi plot -a 120 -e 5 #同上，为简写形式

ifermi plot --property velocity --property-colormap viridis

ifermi plot --property velocity --cmin 0 --cmax 5 #调节color的最大最小值

ifermi plot --property velocity --projection-axis 0 0 1 --property-colormap RdBu

ifermi plot --property velocity --projection-axis 0 0 1 --property-colormap RdBu \ --vector-property --vector-colormap RdBu --hide-surface --hide-labels

注：--vnorm 可用于调节箭头大小 --vector-spacing 可用于调节箭头密度

ifermi plot --slice 0 0 1 0 # 表示001方向，化学势mu为0

ifermi plot --mu -0.85 -i 10 --slice 0 0 1 0 --property spin --hide-cell \ --hide-labels --projection-axis 0 1 0 --property-colormap RdBu

ifermi plot --mu -0.85 -i 10 --slice 0 0 1 0 --property spin --hide-cell \ --hide-labels --projection-axis 0 1 0 --property-colormap RdBu \ --vector-property --vector-colormap RdBu --vnorm 5 --vector-spacing 0.025

Quantum ESPRESSOSuperconducting-Toolkit (SCTK)OpenMXRESPACKVASPkitc2xWannier Berri

40 40 36 (1)0 (2)3 (3)1.0000000 0.57735026 -0.0000000 (4)0.0000000 1.1547005 0.0000000 (5)0.0000000 -0.0000000 0.87206507 (6)2.91340202E-02 (7)2.93242838E-022.98905596E-023.08193434E-02::0.143937960.128004880.0000000 (8)0.362698170.716756941.05351131.3644149::-26.409407-19.318560-10.315671

Background color #背景颜色，RGBLine width #布里渊区边界额nodal line线的宽度Line color #线的颜色Band #能带选择Brillouin zone #布里渊区形式，first BZ还是Primitive BZ，设置后需要更新Number of Brillouin zone #用于形成布里渊区阵列Color bar # 费米面的color barNodal line #输入文件中性质数据为零的线Fermi energy #移动费米面，需要更新Tetrahedron #绘图方案，需要更新

Color scale mode #需要更新

Equator # 需要更新

plot "fermi_line.dat" w l, "bz_line.dat" w l

View point