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.整理結果

reciprocal0.0000 0.0000 0.0000 ! \Gamma0.5000 0.0000 0.5000 ! X0.5000 0.0000 0.5000 ! X0.5000 0.2500 0.7500 ! W0.5000 0.2500 0.7500 ! W0.3750 0.3750 0.7500 ! K0.3750 0.3750 0.7500 ! K0.0000 0.0000 0.0000 ! \Gamma0.0000 0.0000 0.0000 ! \Gamma0.5000 0.5000 0.5000 ! L0.5000 0.5000 0.5000 ! L0.6250 0.2500 0.6250 ! U0.6250 0.2500 0.6250 ! U0.5000 0.2500 0.7500 ! W0.5000 0.2500 0.7500 ! W0.5000 0.5000 0.5000 ! L0.5000 0.5000 0.5000 ! L0.3750 0.3750 0.7500 ! K0.6250 0.2500 0.6250 ! U0.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 ! G0.0000 0.0000 0.5000 50 ! X0.2500 0.2500 0.2500 50 ! P0.0000 0.5000 0.0000 50 ! N0.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”

總結

以上是生活随笔為你收集整理的QE画声子谱的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。