| 查看: 752 | 回复: 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
» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复