加载中…PWSCF计算声子谱步骤
(2017-11-15 22:31:47)pwscf计算声子谱非常简单，第一步利用pw.x进行一步优化，第二步利用ph.x计算倒空间的动力学矩阵（下面的例子是设置的q网格，也可以计算单个点的声子谱），第三步利用q2r.x把倒空间的动力学矩阵通过傅利叶转化成实空间的，第四步利用matyyn.x输出声子模式与声子频率。最后生成的*.freq.gp文件直接导入origin，就是声子谱。pwscf的输出声子谱图的路径与castep设置是一样的，所以直接写入castep中的高对称点的字母即可。
1.# self-consistent
calculation
cat > alas.scf.in << EOF
ATOMIC_SPECIES
ATOMIC_POSITIONS (alat)
K_POINTS
EOF
$ECHO " running the scf
calculation...\c"
$PW_COMMAND < alas.scf.in > alas.scf.out
check_failure $?
$ECHO " done"
2.# phonon calculation on a
(444) uniform grid of q-points
cat > alas.ph.in << EOF
phonons of AlAs
EOF
$ECHO " running the phonon calculation
...\c"
$PH_COMMAND < alas.ph.in > alas.ph.out
check_failure $?
$ECHO " done"
3.cat > q2r.in
<<EOF
EOF
$ECHO " transforming C(q) =>
C(R)...\c"
$Q2R_COMMAND < q2r.in > q2r.out
check_failure $?
$ECHO " done"
4.cat > matdyn.in
<<EOF
EOF
$ECHO " recalculating omega(q) from
C(R)...\c"
$MATDYN_COMMAND < matdyn.in > matdyn.out
check_failure $?
$ECHO " done"
喜欢
0
赠金笔