||
优化过程省略不写了,就是选择出合适的赝势,k-点网格,Ecut等。
第一步建立*.scf.in文件, 运行pw.x <*.scf.in> *.scf.out生成后面计算所需的波函数;
此步计算后会生成一个*.scf.out文件。
第二步建立*.ph.in文件, 运行ph.x <*.ph.in> *.ph.out生成选定q网格点下的动力学矩阵,运算结束会生成 一个*.scf.out文件和 一串*.dyn.0-n文件,*.dyn.0内容是所选q网格会有多少q点,给出了具体坐标,后面的*.dyn.1-n文件分别是n个q点的具体信息文件。
第三步 建立*.q2r.in文件, 运行q2r.x <*.q2r.in> *.q2r.out 傅里叶变换至力常数矩阵,运算结束会生成*.q2r.out和*.fc两个文件。
第四步 建立*.matdyn.in文件,运行matdyn.x <*.matdyn.in > *.matdyn.out生成任意q点下的动力学矩阵,由此可以通过求解久期方程获得格波频率及色散关系。
运行结束会生成*.matdyn.out文件和*.freq文件和*.freq.pq 和*.matdyn.modes四个文件。
第五步 建立*.phdos.in文件,运行matdyn.x <*.phdos.in>*.phdos.out 计算声子态密度,计算结束后会生成*.phdos,*.phdos.out,*.matdyn.freq文件和*.matdyn.freq.pq 四个文件。
******
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 22:20
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社