ZZPWSCF计算实例分析（5）计算Gamma点声子频率

转载自：http://old.blog.edu.cn/user1/11542/archives/2006/1577431.shtml

·PWSCF计算实例分析（5）计算Gamma点声子频率

 

 

PWSCF的ph.x模块是基于linear response理论或密度泛函微扰理论来计算声子频率。计算步骤：
1、先用pw.x进行自洽计算
2、采用ph.x计算Gamma点声子频率

一、自洽计算
&control
calculation='scf',
restart_mode='from_scratch',
prefix='si'
pseudo_dir = './',
outdir='./'
/
&system
ibrav = 2, celldm(1) =10.20, nat= 2, ntyp= 1,
ecutwfc = 18.0
/
&electrons
mixing_beta = 0.7 conv_thr = 1.0d-8
/
ATOMIC_SPECIES
Si 28.086 Si.vbc.UPF
ATOMIC_POSITIONS
Si 0.00 0.00 0.00
Si 0.25 0.25 0.25
K_POINTS (automatic)
8 8 8 0 0 0

pw.x <> si.scf.out

上面是对fcc Si进行自洽的电子结构计算，输入文件与前面对Cu的计算类似。

二、用ph.x进行Gamma点声子频率的计算
输入文件

Phonon of Si at Gamma
&inputph
tr2_ph=1.0d-14,
prefix='si',
epsil=.true.,
amass(1)=28.08,
outdir='./',
fildyn='si.dynG',
/
0.0 0.0 0.0

ph.x <> si.phG.out

上面第一行是注释行，关键词由&inputph .../给出，在'/'之后给出的q点－－Gamma点的坐标。
tr2_ph用来设置声子计算时自洽收敛的标准；
prefix与前面自洽计算中的关键词prefix一致，用来标记所计算的体系，注意它的值要与自洽计算中的一致；
epsil在q=0（也就是Gamma点时），如果体系半导体，而且epsil设置为.true.，则表示计算半导体的宏观介电常数。如果是金属性的体系，或者非Gamma的计算，则不能设置为.true.；
amass(1)用来设置体系中第一类原子的原子量；
outdir用来设置计算输出的目录；
fildyn用来设置指定要将动力学矩阵元输出到什么文件中。