2.3 wannier90+wanniertools

参考:

https://zhuanlan.zhihu.com/p/355317202

https://mp.weixin.qq.com/s?__biz=MzI3MzIyMjI2Nw==&mid=2247488523&idx=1&sn=90a544712f166e248993e1198eb40981&chksm=eb27cb1fdc50420921bb11d3bc8f4833c14418b0d0c0d755af6687fca024fc9d15d999ae7813&scene=21#wechat_redirect

https://blog.sciencenet.cn/blog-2909108-1154273.html

https://zhuanlan.zhihu.com/p/642116307

https://zhuanlan.zhihu.com/p/25432728

一、安装编译

wannier90安装

下载: https://wannier.org/download/

2.1版本是可以与VASP结合的最新版本。但是2.1版本默认安装与VASP接口并不好,主要是借助肖承诚博士写了一个Fortran的接口(https://github.com/Chengcheng-Xiao/VASP2WAN90_v2_fix), 需要注意的是这个接口是针对VASP 5.4.4版本的。1.2版本与VASP的接口是好的,所以1.2版本默认安装就好。对比两个版本的不同,1.2版本的主要缺点就在于不能构建向上自旋和向下自旋的能带,也就是没有自旋轨道耦合作用的铁磁和反铁磁体系。 0201

1.1 wannier90 2.1安装

首先,解压安装包

tar -zxvf wannier90-2.1.0.tar.gz

其次,进入文件夹

cd wannier90-2.1.0/

然后,准备编译文件(这里老王用的是ifort,注意要检查ifort和mpiifort执行命令)

cp config/make.inc.ifort make.inc

接着,编译

make

完成后,编译库,得到libwannier.a文件

make lib

1.2 VASP编译(注意以上接口VASP2WAN90_v2_fix只针对VASP 5.4.4)

首先拷贝VASP2WAN90_v2_fix接口文件中的mlwf.patch 到VASP代码编译目录下。然后执行如下命令。

patch -p0 < mlwf.patch

接着在VASP makefile.include 文件中加入下面两行。注意路径。

1
2
3
CPP_OPTIONS+=-DVASP2WANNIER90v2

LLIBS+=/path/to/your/wannier90_distro/libwannier.a

最后编译即可

make all

wannier-tools安装

开源软件包WannierTools,是一个用于研究新型拓扑材料的软件。此代码在紧束缚模型中工作,可以由另一个软件包Wannier90生成。

它可以通过计算Wilson loop来帮助对给定材料的拓扑相进行分类,通过角分辨光电发射(ARPES)和红外光谱技术,可以得到表面态光谱扫描隧道显微镜(STM)实验。它还可以识别Weyl/Dirac点和节点线结构的位置,计算闭合动量环周围的贝里相位和部分布里渊区(BZ)的贝里曲率。
此外,WannierTools还可以计算非磁性金属和半金属的普通磁电阻利用玻尔兹曼输运理论,计算给定磁场方向和强度下的朗道能级谱;并从超胞计算中得到展开的能谱。

wanniertools官方安装教程

http://www.wanniertools.com/installation.html

安装依赖环境

1、Intel编译器即oneapi(Fortran compiler (Gfortran or ifort),MPICH version higher than 2.1.5,Lapack and Blas library, (Intel MKL recommended))

2、Arpack-ng

安装步骤

下载wanniertools和arpack安装包并分别解压

1
2
unzip wannier_tools-master.zip
tar vzxf arpack96.tar.gz

先进入arpack,修改编译文件并编译
在ARmake.inc中修改编译路径,这里以实际安装路径为准

1
home = $(HOME)/ARPACK

修改编译器,参考VASP安装时的makefile.include

1
2
FC=f77
FFLAGS= -O -cg89

改为

1
2
3
FC         = mpiifort
FCL = mpiifort -mkl=sequential
FFLAGS = -O -cg89

然后执行make编译

make lib

直至产生libarpack_SUN4.a文件,该文件会因为不同机器不同环境而不同名称,只要确定好,然后在编译wanniertools的时候指定好就行了

在wanniertools的src文件夹的makefile文件中该默认值为

1
2
3
#ARPACK LIBRARY

ARPACK=/Users/quan/work/workplace/ARPACK/libarpack MAC.a

修改为刚刚已经生成的文件路径,这里使用的是Makefile.intel-mpi,重命名为makefile

1
cp Makefile.intel-mpi Makefile

然后编译即可

1
make

最后生成wt.x文件即为安装成功

1
cp -f wt.x../bin

二、计算过程(以石墨烯为例)

1、结构优化加自洽计算

POSCAR为

{.line-numbers}
1
2
3
4
5
6
7
8
9
10
Surf-mp-48                             
1.00000000000000
1.2344083195740001 -2.1380573715510001 0.0000000000000000
1.2344083195740001 2.1380573715510001 0.0000000000000000
0.0000000000000000 0.0000000000000000 12.0000000000000000
C
2
Direct
0.0000000000000000 0.0000000000000000 0.5000000000000000
0.3333329999917112 0.6666670000099160 0.5000000000000000

INCAR为

{.line-numbers}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
SYSTEM=Graphene
Global Parameters
ISTART = 0 (Read existing wavefunction; if there)
ISPIN = 1 (Non-Spin polarised DFT)
ICHARG = 2 (Non-self-consistent: GGA/LDA band structures)
LREAL = .FALSE. (Projection operators: automatic)
ENCUT = 500 (Cut-off energy for plane wave basis set, in eV)
PREC = A (Precision level)
ADDGRID= .TRUE. (Increase grid; helps GGA convergence)

Electronic Relaxation
ALGO = N
ISMEAR = 0 (Gaussian smearing; metals:1)
SIGMA = 0.05 (Smearing value in eV; metals:0.2)
NELM = 90 (Max electronic SCF steps)
NELMIN = 6 (Min electronic SCF steps)
EDIFF = 1E-08 (SCF energy convergence; in eV)

Ionic Relaxation
NSW = 0 (Max ionic steps)
IBRION = -1 (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
ISIF = 3 (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
EDIFFG = -2E-02 (Ionic convergence; eV/AA)

加SOC在自洽计算加上
LSORBIT=T

2、非自洽计算投影能带

在自洽计算基础上,cp -r scf band
修改INCAR

1
2
3
ISTART =  1         
ICHARG = 11
LORBIT = 10/11

KPOINTS 改为高对称点路径。

画一下能带结构,然后再做一下Fatband分析,了解费米面附近能带是有哪些成分构成的。以C的pz轨道。

3、非自洽vasp产生.win文件

在自洽计算基础上 cp -r scf wan1
修改INCAR

1
2
3
4
5
ISTART =  1         
ICHARG = 11
NBANDS = 40
LWANNIER90 = .TRUE.
NUM_WANN = 2

注意,INCAR中使用NPAR参数会出错

KPOINTS、POSCAR、POTCAR不改。

准备好wannier90.win文件,同时在INCAR里头添加一行LWANNIER90=.TRUE.和一行ICHARG=11,然后再运行一遍VASP。 顺利的话,会产生 wannier90.amn, wannier90.mmn, wannier90.eig, wannier90.chk 等4个文件。wannier90.win里头的内容也有可能会被改变。包括结构信息,k点信息。有了这三个文件就可以进一步调节参数,构造wannier函数了。请注意wannier90.win里头的num_bands要和OUTCAR里头的NBANDS一致。.mmn积分矩阵和.amn布洛赫态在轨道上的尝试投影和UNK0000x.1存储实空间单胞的布洛赫态相关信息,这些文件是运行链接wannier90的DFT计算输出的结果,而gaas.win是计算参数设置文件

这一步的wannier90.win文件只需要加入几个关键参数,其他会自动生成。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
# 解纠缠窗口
dis_win_max =
dis_win_min=
# 冻结窗口
dis_froz_max =
dis_forz_min = #(default = dis_win_min)

num_bands = #(default = nbnd - exclude_bands)
num_wann = #(num_wann <= num_bands, 等于投影总轨道数,此处为32)
exclude_bands =

begin projections ## block
...
end projections

mp_grid = 8 8 8 ## 与自洽的KPOINTS一致

4、wannier90.x计算

运行wannier90主程序。mpirun -np 40 wannier90.x wannier90 &

注意在第3步打开hr_plot = T,关闭write_hr = true;在第4步打开write_hr = true,关闭hr_plot = T

这里后面一个winner90代表wannier90.win的文件.win前面的部分。如果你的win文件是name.win运行的地方就可以是

wannier90.x name &

运行过程中当num_bands > num_wann 有个迭代过程,控制迭代步数可以通过

1
2
dis_num_iter  =  200   # 迭代1,默认200步
num_iter = 200 # 迭代2,默认200步

最后会生成wannier90.chk文件,restart计算会读取这个文件。

加入能带部分的话第一次运行就会出来,wannier90_band.dat文件,这个可以用来和DFT结果对比。
构造好Wannier函数以后,仔细检查一下wannier90.wout,
搜索”Final State”,在这里可以看到每条Wannier轨道的展宽,由此可以检查你构造的Wannier函数是否足够局域。如果有某个Wannier函数的展宽
比晶格常数大许多,就表明你的Wannier函数构造不是很理想,需要不断调节投影轨道,解纠缠窗口以及Frozen窗口等。从这些轨道的中心和展宽也能
看出轨道之间的简并性如何,由此也可以判断所构造的Wannier函数对称性如何,这个也会影响后续的分析,特别是表面态的计算。查看wannier90.wout末尾的Final State的最后一列小于3,其结果是可以的。

5、

紧接着,你就可以开始书写WannierTools的主输入文件了,
目前叫input.dat,以后会改成http://wt.in (WannierTools v2.2之后)。 这里不具体介绍每一个参数,在说明文档中有详细的介绍,中文版的说明文档正在翻译中。
input.dat文件有模板,使用的时候直接拷贝然后加以修改。很多内容都可以直接从wannier90.win和wannier90.wout中拷
贝。 Bi2Se3的主输入文件input.dat可以在附件中下载得到。我们可以通过控制输入文件中的一些参数来实现我们所需要的功能。

常见报错

  1. 首先跑接口会生成3个文件.amm .mmm .eig通过这三个文件构造wannier函数。也可以提前把.win文件准备好放进去。以上三个文件生成后运行wannier90.x会成.chk文件

  2. Restar=T是需要读取.chk文件的

Error: restart requested but wannier90.chkfile not found

原因:加了restart = plot 这个命令,但却没有wannier90.chk这个文件。Restart是需要读取wannier90.chk的。

解决办法:注释掉restart = plot,然后重新运行wannier90.x wannier90 &得到wannier90.chk.

  1. 问题:Wannier90 is running in LIBRARY MODE

Setting up k-point neighbours…

Ignoring in input file

解决办法:尽量串行,不要并行。请再次运行 $postw90 wannier90的为输出的数据

  1. Spinors与文件INCAR中的ISPIN的设置有关,同时存在或者同时去掉,否则影响.amn文件的形成。(自旋轨道耦合)

  2. 谨慎用use_bloch_phases=T一般是F

  3. 问题:param_read: mismatch in wannier90.eig

解决办法:每一步计算的核数要相等

  1. 问题:dis_windows: More states in the forze window than target WFs

解决办法:把dis_windows改大点

  1. Fewer projections defined than the number of wannier functions requested.

投影的轨道太少

解决办法:增加投影的轨道

  1. Kmesh_get_bvector: Not enongh bevectors found

解决办法:添加kemsh_tol

  1. Wanted band:1 found band:17

Wanted kpoint:2 found kpoint:1

解决办法:常见原因是num_bands设置错误,返回band计算文件夹查看NBANDS

  1. Dis_windows: Energy window contains fewer states than number sf target WFs

解决办法:把dis_windows改大点

  1. Dis_windows: more states in the frozen window than WFs

解决办法:调节frozen window, 投影太多或者这个窗口太大涉及到别的轨道

  1. Wannier 运行有两种模式一是library mode,一种是单机模式stand-alone mode

解决办法:library mode是需要进行串行。

  1. param_get_projection: too many projections defined

解决办法:begin projections random

原因:num_wann 与设置投影数不匹配。

解决办法:不加soc:num_wann 等于每种原子投影轨道数之和

soc:num_wann 等于两倍的每种原子投影轨道数之和

  1. Gnuplot的报错信息”arc_bulk.gnu”, line 26: warning: No unsable data in this plot to auto-scale axis range

  2. param_get_projection: Problem reading m state into string

解决办法:已解决Fe: s;px;py;pz;dz2;dxz;dyz;dx2-y2;dxy

不是Fe:s;p;d

  1. 能带拟合的不好可能是因为,投影的不对,还有核数和能带数的关系

解决办法:太马虎了,这样的问题肯定是K点写错了

能带布里渊区的宽度不一样有可能是wannier里规定的K点数并没有完全计算。

  1. num_iter = 2000 #wannier center 迭代次数

num_print_cycles = 40 #wannier center 迭代每40步输出一次

  1. dis_num_iter = 1200 #解纠缠的迭代步数#

  2. num_iter最小化过程中的总迭代次数。如果您希望生成投影wf而不是最大本地化的wf,则设置num_iter=0(请参阅教程中的示例8)?默认值为100

  3. dis_mix_ratio = 1.d0 # 0.5<dis_mix_ratio<1, 默认值0.5

#在解纠缠过程中,要利用混合参数进行收敛,值为0.5是“安全”的选择。使用1.0(即不混合)通常会得到更快的收敛,但可能会导致 最小化,在某些情况下是不稳定的。

dis_conv_tol = 1.0E-10 # Delta < 1.000E-10 over 3 iterations

  1. dis_win_max = 15

dis_froz_max = 14

num_wan = 12

exclude_bands : 1-11,26 #在使用这个参数时,num_bands = 总能带数减去排除在外的能带数

num_bands = 14

#如果使用能量窗口规定想要拟合的能带区间,需要同时使用这个几个参数

  1. wannier_plot_supercell =2 #在一个对应于“超胞”的网格上生成WFs

  2. wannier_plot_list =2 #列举要画出WFs的序号,与wannier center 的序号是一致的

  3. 运行接口的情况有三种,一是选择投影能带数发生了改变,二是,需要写入wannier90.win的输入文件;三是需要在INCAR添加关键参数,输出wannier.x需要的文件。(仅仅是调节能量窗口是不需要重新运行接口的)

  4. 无没有projection是肯定没有amn的

wannier90.win参数介绍

系统参数

1
2
3
4
5
6
7
8
num_wann: I, Number of WF 
num_bands: I, Number of bands passed to the code
unit_cell_cart: P, Unit cell vectors in Cartesian coordinates,不用设默认即可
atoms_cart*: P, Positions of atoms in Cartesian coordinates,会自动补全
atoms_frac*: R, Positions of atoms in fractional coordinates with respect to the lattice vectors会自动补全
mp_grid: I,Dimensions of the Monkhorst-Pack grid of k-points与自洽KPOINTS一致
spinors: L, WF are spinors考虑SOC或者不考虑,默认是关闭不考虑

任务控制

1
2
3
4
5
6
exclude_bands: I, List of bands to exclude from the calculation不考虑那几条能带
restart: S, Restart from checkpoint file
iprint: I, Output verbosity level
spin: S, Which spin channel to read
translate_home_cell: L, To translate final Wannier centres to home unit cell when writing xyz file
write_xyz: L, To write atomic positions and final centres in xyz file format

画图参数

1
2
3
4
5
6
7
8
9
10
wannier_plot: L, Plot the WF
bands_plot: L, Plot interpolated band structure
kpoint_path: P, K-point path for the interpolated band structure
bands_num_points: I, Number of points along the first section of the k-point path
bands_plot_format: S, File format in which to plot the interpolated
bands bands_plot_project: I, WF to project the band structure onto
bands_plot_mode: S, Slater-Koster type interpolation or Hamiltonian cut-off
fermi_surface_plot: L, Plot the Fermi surface
write_hr: L, Write the Hamiltonian in the WF basis
dist_cutoff: P, Cut-off for the distance between WF

解纠缠参数

1
2
3
4
5
6
dis_win_min: P, Bottom of the outer energy window 
dis_win_max: P, Top of the outer energy window
dis_froz_min: P, Bottom of the inner (frozen) energy window
dis_froz_max: P, Top of the inner (frozen) energy window
dis_num_iter: I, Number of iterations for the minimisation of ΩI

wannier化参数

1
2
3
num_iter: I, Number of iterations for the minimisation of Ω
num_print_cycles: I, Control frequency of printing
guiding_centres: L, Use guiding centres

wt.in参数介绍

可参考:https://mp.weixin.qq.com/s/HUDobUio4BdQfVLImZfYmg

先给一个模版:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
&TB_FILE
Hrfile = 'wannier90_hr.dat' !default
Package = 'VASP' !default
/


&CONTROL
AHC_calc = T
ANE_calc = F
BulkBand_calc = T
BulkFS_calc = F
BulkGap_cube_calc = F
BulkGap_plane_calc = F
SlabBand_calc = F
WireBand_calc = F
SlabSS_calc = F
SlabArc_calc = F
SlabQPI_calc = F
Z2_3D_calc = F
Chern_3D_calc = F
SlabSpintexture_calc = F
Wanniercenter_calc = T
BerryPhase_calc = T
BerryCurvature_calc = F
EffectiveMass_calc = F
/

&SYSTEM
#NSLAB = 10 ! for thin film system default=10
#NSLAB1= 4 ! nanowire system
#NSLAB2= 4 ! nanowire system
NumOccupied = 18 ! NumOccupied=num_wan
SOC = 1 ! soc on
E_FERMI = 4.4195 ! e-fermi
#Bx= 0, By= 0, Bz= 0
#surf_onsite= 0.0
/

&PARAMETERS
Eta_Arc = 0.0001 ! infinite small value, like brodening
#E_arc = -0.12 ! energy for calculate Fermi Arc
OmegaNum = 1001 ! omega number step=1meV
OmegaMin = -1 ! energy interval
OmegaMax = 1 ! energy interval
Nk1 = 151 ! number k points odd number would be better
Nk2 = 151 ! number k points odd number would be better
Nk3 = 1 ! number k points odd number would be better test
#NP = 1 ! number of principle layers default=2
#Gap_threshold = 1.0 ! threshold for GapCube output
TMin = 100     ! energy interval 注意 Tmin ≠0
TMax = 300    ! energy interval
NumT = 100  ! omega number  (300-100)/100 从100开始间隔100K输出一次
/

LATTICE
Angstrom
-2.069 -3.583614 0.000000 ! crystal lattice information
2.069 -3.583614 0.000000
0.000 2.389075 9.546667

ATOM_POSITIONS
5 ! number of atoms for projectors
Direct ! Direct or Cartisen coordinate
Bi 0.3990 0.3990 0.6970
Bi 0.6010 0.6010 0.3030
Se 0 0 0.5
Se 0.2060 0.2060 0.1180
Se 0.7940 0.7940 0.8820

PROJECTORS
3 3 3 3 3 ! number of projectors
Bi px py pz ! projectors
Bi px py pz
Se px py pz
Se px py pz
Se px py pz

SURFACE ! from version v2.2.6 on, you can specify a surface with only two lattice vectors
1 0 0
0 1 0

KPATH_BULK ! BulkBand_calc=T
4 ! number of k line only for bulk band
G 0.00000 0.00000 0.0000 Z 0.00000 0.00000 0.5000
Z 0.00000 0.00000 0.5000 F 0.50000 0.50000 0.0000
F 0.50000 0.50000 0.0000 G 0.00000 0.00000 0.0000
G 0.00000 0.00000 0.0000 L 0.50000 0.00000 0.0000

KPATH_SLAB !SlabBand_calc=T or SlabSS_calc=T
2 ! numker of k line
K 0.33 0.67 G 0.0 0.0
G 0.0 0.0 M 0.5 0.5

KPLANE_SLAB !SlabArc_calc=T, SlabSpintexture_calc= T
-0.1 -0.1 ! Original point for 2D k plane
0.2 0.0 ! The first vector to define 2D k plane
0.0 0.2 ! The second vector to define 2D k plane for arc plots

KPLANE_BULK !BulkGap_plane_calc=T, BerryCurvature_calc=T, wanniercenter_calc=T
0.00 0.00 0.50 ! Original point for 3D k plane k3=0.5, bar{a}, along k1
1.00 0.00 0.00 ! The first vector to define 3d k space plane
0.00 0.50 0.00 ! The second vector to define 3d k space plane

!> The following 5 matrices are for backup using, will not affect the main input for WannierTools
0.00 0.00 0.00 ! Original point for 3D k plane k3=0.0, bar{a}, along k1
1.00 0.00 0.00 ! The first vector to define 3d k space plane
0.00 0.50 0.00 ! The second vector to define 3d k space plane

0.00 0.50 0.00 ! Original point for 3D k plane k2=0.5, bar{c}, along ka
0.00 0.00 1.00 ! The first vector to define 3d k space plane
0.50 0.00 0.00 ! The second vector to define 3d k space plane

0.00 0.00 0.00 ! Original point for 3D k plane k2=0, bar{c}, along ka
0.00 0.00 1.00 ! The first vector to define 3d k space plane
0.50 0.00 0.00 ! The second vector to define 3d k space plane

0.50 0.00 0.00 ! Original point for 3D k plane k1=0.5, bar{c}, along kb
0.00 0.00 1.00 ! The first vector to define 3d k space plane
0.00 0.50 0.00 ! The second vector to define 3d k space plane

0.00 0.00 0.00 ! Original point for 3D k plane k1=0, bar{c}, along kb
0.00 0.00 1.00 ! The first vector to define 3d k space plane
0.00 0.50 0.00 ! The second vector to define 3d k space plane



KCUBE_BULK !BulkGap_cube_calc=T
-0.50 -0.50 -0.50 ! Original point for 3D k plane
1.00 0.00 0.00 ! The first vector to define 3d k space plane
0.00 1.00 0.00 ! The second vector to define 3d k space plane
0.00 0.00 1.00 ! The third vector to define 3d k cube


EFFECTIVE_MASS ! optional
2 ! The i'th band to be calculated
0.01 ! k step in unit of (1/Angstrom)
0.0 0.0 0.0 ! k point where the effective mass calculated.


WANNIER_CENTRES ! copy from wannier90.wout
Cartesian
-0.000040 -1.194745 6.638646
0.000038 -1.196699 6.640059
-0.000032 -1.192363 6.640243
-0.000086 -3.583414 2.908040
0.000047 -3.581457 2.906587
-0.000033 -3.585864 2.906443
-0.000001 1.194527 4.773338
0.000003 1.194538 4.773336
-0.000037 1.194536 4.773327
0.000006 -1.194384 1.130261
-0.000018 -1.216986 1.140267
0.000007 -1.172216 1.140684
0.000011 -3.583770 8.416406
-0.000002 -3.561169 8.406398
-0.000007 -3.605960 8.405979
0.000086 -1.194737 6.638626
-0.000047 -1.196693 6.640080
0.000033 -1.192286 6.640223
0.000040 -3.583406 2.908021
-0.000038 -3.581452 2.906608
0.000032 -3.585788 2.906424
0.000001 1.194548 4.773330
-0.000003 1.194537 4.773332
0.000037 1.194539 4.773340
-0.000011 -1.194381 1.130260
0.000002 -1.216981 1.140268
0.000007 -1.172191 1.140687
-0.000006 -3.583766 8.416405
0.000018 -3.561165 8.406400
-0.000007 -3.605935 8.405982

Tips:

  • 2D材料就按照bulk去算
  • 需要测试Nk1 、NK2、NK3对结果的影响
  • 必须设置的模块:&TB_FILE、&CONTROL、&SYSTEM、&PARAMETERS、LATTICE、ATOM_POSITIONS、PROJECTORS、WANNIER_CENTRES ;选择设置SURFACE、KPATH_BULK、 KPATH_SLAB、 KPLANE_SLAB、KPLANE_BULK 、KCUBE_BULK、EFFECTIVE_MASS 等。