版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

a_lin916

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 2
在线: 1.8小时
虫号: 3272405
注册: 2014-06-13
专业: 凝聚态物性 II ：电子结构

[求助] 求助！ 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,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有5人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有9人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有9人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有8人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有11人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有12人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有9人回复

1楼 2016-11-01 08:26:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 a_lin916 的主题更新