| 查看: 838 | 回复: 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
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
20