QE画声子谱

1 q2r.x

建立q2r.in

&INPUT
fildyn = 'matdyn'
 lshift_q = .true.
   zasr = 'crystal'
  flfrc = 'ifc.dat'
/

q2r.x -in q2r.in > q2r.out

2.计算高对称点

2.1方法1：MS

1.改元胞

2.计算高对称点

2.2方法2：网站AFLOW

http://www.aflowlib.org/aflow-online/

2.2.1.输入POSCAR

2.2.2.勾选生成K点路径

2.2.3.生成

2.2.4.整理结果

reciprocal
   0.0000   0.0000   0.0000   ! \Gamma
   0.5000   0.0000   0.5000   ! X
 
   0.5000   0.0000   0.5000   ! X
   0.5000   0.2500   0.7500   ! W
 
   0.5000   0.2500   0.7500   ! W
   0.3750   0.3750   0.7500   ! K
 
   0.3750   0.3750   0.7500   ! K
   0.0000   0.0000   0.0000   ! \Gamma
 
   0.0000   0.0000   0.0000   ! \Gamma
   0.5000   0.5000   0.5000   ! L
 
   0.5000   0.5000   0.5000   ! L
   0.6250   0.2500   0.6250   ! U
 
   0.6250   0.2500   0.6250   ! U
   0.5000   0.2500   0.7500   ! W
 
   0.5000   0.2500   0.7500   ! W
   0.5000   0.5000   0.5000   ! L
 
   0.5000   0.5000   0.5000   ! L
   0.3750   0.3750   0.7500   ! K
 
   0.6250   0.2500   0.6250   ! U
   0.5000   0.0000   0.5000   ! X

2.2.4.1 把不连着的路径删掉

reciprocal
0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.0000 0.5000 ! X

0.5000 0.0000 0.5000 ! X
0.5000 0.2500 0.7500 ! W

0.5000 0.2500 0.7500 ! W
0.3750 0.3750 0.7500 ! K

0.3750 0.3750 0.7500 ! K
0.0000 0.0000 0.0000 ! \Gamma

0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.5000 0.5000 ! L

0.5000 0.5000 0.5000 ! L
0.6250 0.2500 0.6250 ! U

0.6250 0.2500 0.6250 ! U
0.5000 0.2500 0.7500 ! W

0.5000 0.2500 0.7500 ! W
0.5000 0.5000 0.5000 ! L

0.5000 0.5000 0.5000 ! L
0.3750 0.3750 0.7500 ! K

0.6250 0.2500 0.6250 ! U
0.5000 0.0000 0.5000 ! X

2.2.4.2把重复路径删掉

reciprocal
0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.0000 0.5000 ! X

0.5000 0.0000 0.5000 ! X
0.5000 0.2500 0.7500 ! W

0.5000 0.2500 0.7500 ! W
0.3750 0.3750 0.7500 ! K

0.3750 0.3750 0.7500 ! K
0.0000 0.0000 0.0000 ! \Gamma

0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.5000 0.5000 ! L

0.5000 0.5000 0.5000 ! L
0.6250 0.2500 0.6250 ! U

0.6250 0.2500 0.6250 ! U
0.5000 0.2500 0.7500 ! W

0.5000 0.2500 0.7500 ! W
0.5000 0.5000 0.5000 ! L

0.5000 0.5000 0.5000 ! L
0.3750 0.3750 0.7500 ! K

2.2.4.3 最后结果

0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.0000 0.5000 ! X
0.5000 0.2500 0.7500 ! W
0.3750 0.3750 0.7500 ! K
0.0000 0.0000 0.0000 ! \Gamma
0.5000 0.5000 0.5000 ! L
0.6250 0.2500 0.6250 ! U
0.5000 0.2500 0.7500 ! W
0.5000 0.5000 0.5000 ! L
0.3750 0.3750 0.7500 ! K

2.3方法3：网站

网址Bilbao Crystal

https://www.cryst.ehu.es/

点击Space-group symmetry
点击输入空间群数字

3.matdyn.x

matdyn.x -in matdyn.in > matdyn.out

建立matdyn.in

&input
asr = 'crystal',
flfrc = 'ifc.dat'
flfrq = 'freq'
q_in_band_form = .true.
q_in_cryst_coord = .true.
/
#上面生成的高对称点
6 #一共有几个
-0.5000    0.5000    0.5000  50 ! Z
（坐标 + 在两个高对称点之间插入几个点一般选30或50，越大越密）
 0.0000    0.0000    0.0000  50 ! G
 0.0000    0.0000    0.5000  50 ! X
 0.2500    0.2500    0.2500  50 ! P
 0.0000    0.5000    0.0000  50 ! N
 0.0000    0.0000    0.0000  50 ! G

plot “freq.gp.dat” u 1:2 w l title “matdyn.x”,
“” u 1:3 w l notitle,
“” u 1:4 w l notitle,
“freq” u 1:($2*8.06554) pt 7 title “EPW”