qe计算电声耦合系数

1、首先进行‘scf’运算scf.in（这一步不能用mpirun运算）###########################################################&CONTROLtitle='graphene-ph',calculation='scf',#计算主内容outdir='./',#输出文件目录pseudo_dir='/home/txd/qe/qe-6.1/pseudo',#赝势文件目录prefix='graphene',#输出文件抬头tstress=.true.,tprnfor=.true.,/&SYSTEMibrav=4,#布拉维格子celldm(1)=4.65,#晶格矢celldm(3)=12.85,nat=2,#原子数ntyp=1,#原子类型ecutwfc=30,#波函数的截断动能，单位Ryecutrho=300,#电荷密度和势能的截断动能，默认为4*ecutwfcnbnd=21,#能带数occupations='smearing',#占据态计算方法degauss=0.05,smearing='marzari-vanderbilt',#/&ELECTRONS/ATOMIC_SPECIESC12.00000C.pz-kjpaw.UPFATOMIC_POSITIONScrystalC0.6670001010d00.3330002021d00.5000000000d0C0.3330001299d00.6670002598d00.5000000000d0K_POINTSautomatic555000##################################################################################2、用ph.x进行指定k点范围声子频率的计算phX.in##################################################################################&INPUTPHoutdir='./',#输出文件目录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.00.00.0）#用于gamma点计算##################################################################################3、用q2r.x计算力常数矩阵q2r.in##################################################################################&inputzasr='simple',#如何处理‘声学支求和规则“，该规则是用在处理Born有效电荷的，要求Born有效电荷的总和是零。#可赋的值有:#no,表示不处理声学支求和问题#simple,表示通过对力常数矩阵的对角元素进行修正来考虑3支声学横模的求和处理;#crystal,#one-dim,#zero-dimfildyn='name.dyn',#读取包含了q网格点的动力学矩阵元文件flfrc='graphene.fc',#输出力常数矩阵文件la2F=.false.#设置是否计算出实空间中电声耦合系数/#################################################################################4、用matdyn.x计算出声子色散曲线matdyn.in#################################################################################&inputasr='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#频率本征值文件01700#输出纵轴坐标freq.plot#输出数据文件freq.ps#输出.ps图像文件0.0#Fermi能50.00.0#能量度量间隔，能量平移大小1、电声耦合常数计算参数设置elph=.true.,5.0.2版以后改为electron_phonon='interpolated'2、lambda.dat解读lambda:theelectron-phononcoefficientλintalpha2F:α^2F(w)functionlogw:anestimateofthecriticaltemperatureTc????N(Ef):费米能级占据态？？