| 查看: 764 | 回复: 0 | ||
[求助]
求助！ QE算graphene二阶力常数， 输出文件中没有dielectric tensor 结果是F
|
|
求助！ QE算graphene二阶力常数， 输出文件中没有dielectric tensor 和Born effective charge tensor， 结果是F.求大神帮忙！！！ 以下是提交的CODE以及输出结果 （1） code: # self-consistent calculation cat > graphene.scf.in << EOF &control calculation=\'scf\' restart_mode=\'from_scratch\', tstress = .true. tprnfor = .true. prefix=\'graphene\', pseudo_dir = \'$PSEUDO_DIR/\', outdir=\'$TMP_DIR/\' / &system ibrav = 4, celldm(1)=4.61093173d0, celldm(3)=7.37704918d0, nat = 2, ntyp = 1, occupations = \'smearing\', smearing = \'mv\', degauss = 0.05d0, ecutwfc = 80.0, nbnd = 8, / &electrons conv_thr = 1.0d-8, mixing_mode = \'plain\', mixing_beta = 0.7, diagonalization = \'david\', / ATOMIC_SPECIES C 12.010700d0 C.pz-vbc.UPF ATOMIC_POSITIONS {crystal} C 0.0000000000d0 0.0000000000d0 0.5000000000d0 C 0.3333333333d0 0.6666666667d0 0.5000000000d0 K_POINTS {automatic} 28 28 1 0 0 0 EOF $ECHO " running the scf calculation...\\c" $PW_COMMAND < graphene.scf.in > graphene.scf.out check_failure $? $ECHO " done" # phonon calculation on a (14141) uniform grid of q-points cat > graphene.ph.in << EOF phonons of inas &inputph tr2_ph=1.0d-14, prefix=\'graphene\', ldisp=.true., nq1=9, nq2=9, nq3=1 amass(1)= 12.010700d0, outdir=\'$TMP_DIR/\', fildyn=\'graphene.dyn\', / EOF $ECHO " running the phonon calculation ...\\c" $PH_COMMAND < graphene.ph.in > graphene.ph.out check_failure $? $ECHO " done" cat > q2r.in <<EOF &input fildyn=\'graphene.dyn\', zasr=\'simple\', flfrc=\'graphene991.fc\' / EOF $ECHO " transforming C(q) => C(R)...\\c" $Q2R_COMMAND < q2r.in > q2r.out check_failure $? $ECHO " done" (2) output file: 1 2 4 4.6109317 0.0000000 7.3770492 0.0000000 0.0000000 0.0000000 1 \'C \' 10947.083358978085 1 1 0.0000000000 0.0000000000 3.6885245900 2 1 -0.0000000001 0.5773502692 3.6885245900 F 9 9 1 1 1 1 1 ... 输出结果是F |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : graphene991.fc
2016-11-01 07:52:59, 94.83 K
» 猜你喜欢
284求调剂
已经有10人回复
-
一志愿山东大学药学学硕求调剂
已经有4人回复
-
07化学280分求调剂
已经有4人回复
-
298-一志愿中国农业大学-求调剂
已经有12人回复
-
求材料，环境专业调剂
已经有3人回复
-
335求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
一志愿吉大化学322求调剂
已经有4人回复
-
环境学硕288求调剂
已经有8人回复
-
341求调剂(一志愿湖南大学070300)
已经有6人回复
5