qe计算电声耦合系数

1、首先进行‘scf’运算
scf.in（这一步不能用mpirun运算）
###########################################################
&CONTROL
title = 'graphene-ph' ,
calculation = 'scf' , #计算主内容
outdir = './' , #输出文件目录
pseudo_dir = '/home/txd/qe/qe-6.1/pseudo' , #赝势文件目录
prefix = 'graphene' , #输出文件抬头
tstress = .true. ,
tprnfor = .true. ,
/
&SYSTEM
ibrav = 4, #布拉维格子
celldm(1) = 4.65, #晶格矢
celldm(3) = 12.85,
nat = 2, #原子数
ntyp = 1, #原子类型
ecutwfc = 30 , #波函数的截断动能，单位Ry
ecutrho = 300 , #电荷密度和势能的截断动能，默认为4*ecutwfc
nbnd = 21, #能带数
occupations = 'smearing' , #占据态计算方法
degauss = 0.05 ,
smearing = 'marzari-vanderbilt' , #
/
&ELECTRONS
/

ATOMIC_SPECIES
C 12.00000 C.pz-kjpaw.UPF

ATOMIC_POSITIONS crystal
C 0.6670001010d0 0.3330002021d0 0.5000000000d0
C 0.3330001299d0 0.6670002598d0 0.5000000000d0

K_POINTS automatic
5 5 5 0 0 0
##################################################################################

2、用ph.x进行指定k点范围声子频率的计算
phX.in
##################################################################################
&INPUTPH
outdir = './' , #输出文件目录
prefix = 'graphene' , #抬头
fildvscf = 'graphdv' , #用于电声运算的势能变化文件
fildyn = 'graphene.dyn' #动力学矩阵文件
trans = .true.,
ldisp = .true.,
amass(1) = 12.00 , #原子质量
tr2_ph = 1e-10 , #自洽阈值
nq1=4, nq2=4, nq3=1 , #q点范围
/
（0.0 0.0 0.0） #用于gamma点计算
##################################################################################

3、用q2r.x计算力常数矩阵
q2r.in
##################################################################################
&input
zasr='simple', #如何处理‘声学支求和规则“，该规则是用在处理Born有效电荷的，要求Born有效电荷的总和是零。

#可赋的值有:
#no,表示不处理声学支求和问题
#simple, 表示通过对力常数矩阵的对角元素进行修正来考虑3支声学横模的求和处理;
#crystal,
#one-dim,
#zero-dim
fildyn='name.dyn', #读取包含了q网格点的动力学矩阵元文件
flfrc='graphene.fc', #输出力常数矩阵文件
la2F=.false. #设置是否计算出实空间中电声耦合系数
/
#################################################################################

4、用matdyn.x计算出声子色散曲线
matdyn.in
#################################################################################
&input
asr='simple',
amass(1)=12.00,
flfrc='graphene.fc', #读取力常数文件
flfrq='graphene444.freq', #输出频率本征值文件
la2F=.false.,
dos=.false. #设置计算声子态密度
（dos=.true.,fldos='phonon.dos', nk1=10, nk2=10, nk3=10, ndos=50）#计算声子态密度
/
#################################################################################

5、用plotband.x画图
plotband.in
#################################################################################
graphene444.freq #频率本征值文件
0 1700 #输出纵轴坐标
freq.plot #输出数据文件
freq.ps #输出.ps图像文件
0.0 #Fermi能
50.0 0.0 #能量度量间隔，能量平移大小


1、电声耦合常数计算参数设置 elph=.true.,5.0.2版以后改为 electron_phonon='interpolated'
2、lambda.dat解读
lambda : the electron-phonon coefficient λ
int alpha2F:α^2F(w) function
:an estimate of the critical temperature Tc????
N(Ef) :费米能级占据态？？