| 查看: 795 | 回复: 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文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
2026博士还有哪些学校有名额
已经有4人回复
-
上海大学实验技术岗位非升即走
已经有5人回复
-
窗边初夏的小雨
已经有11人回复
-
西南大学考核制博士
已经有7人回复
-
裁员滚滚，退居二线
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
有带发论文的吗
已经有4人回复
-
风把牡丹吹跑了
已经有9人回复
-
青C资助名额大幅增加！
已经有15人回复
