24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

[求助] QE计算MgB2的声子谱，声子的软化现象已有1人参与

最近在用QE计算MgB2的声子谱，但是计算后始终是这样，
想和大家请教一下这个是什么原因造成的，有可能是K点取得太小了造成的吗？目前我自己觉得可能是这样。
然后贴上我的自洽的参数，然后我的结构相关的参数我觉得应该没有错吧，难道是结构我写的不对吗。
大家帮忙看看，非常感谢了，谢谢大家！
&control
calculation='scf'
restart_mode='from_scratch',
!pseudo_dir='directory where pseudopotentials are stored/',
!outdir='directory where large files are written/'
pseudo_dir='./'
outdir='./tmp'
prefix='PHMgB2'
tstress=.true.
tprnfor=.true.
/
&system
ibrav=4,celldm(1)=5.8013,celldm(2)=5.8013,celldm(3)=1.14488,
nat=3,
ntyp=2,
ecutwfc=60,
ecutrho=480,
occupations='smearing', smearing='mp', degauss=0.01
/
&electrons
diagonalization='david'
conv_thr=1.0d-8
mixing_beta=0.7
/
ATOMIC_SPECIES
Mg  24.3050  Mg.pz-n-vbc.UPF
B 10.811  B.pz-bhs.UPF
ATOMIC_POSITIONS {crystal}
  Mg  -0.0000000050616800  0.0000000000000088  0.0000000000000476
  B 0.3333333201942207  0.6666666660819658  0.5000000000000511
  B 0.6666666548674568  0.3333333339180254  0.4999999999999015
K_POINTS {automatic}
10 10 10 0 0 0

文献的.png

我做的.png

回复此楼

» 猜你喜欢

求个博导看看已经有17人回复
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

1楼 2016-04-16 12:54:29

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

菲顶顶

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 32.6
散金: 45
红花: 10
帖子: 631
在线: 250.2小时
虫号: 1462090
注册: 2011-10-26
专业: 凝聚态物性 II ：电子结构

引用回帖:

6楼: Originally posted by y1ding at 2016-04-18 21:19:00
你把ph的输入文件也贴出来
pw的输入没啥大问题...

下面是ph文件内容：
   1 Phonons for MgB2
   2  &inputph
   3 tr2_ph=1.0d-14,
   4 amass(1)=24.3050, amass(2)=10.811,
   5 prefix='PHMgB2',
   6 ldisp=.true.,
   7 trans=.true.,
   8 outdir='./tmp',
   9 fildyn='MgB2.dyn',
   10 nq1=5,
   11 nq2=5,
   12 nq3=5,
   13  /
下面是q2r文件内容：
   1 Phonons for MgB2
   2  &input
   3 fildyn='MgB2.dyn', flfrc='PHMgB2.fc', zasr='simple'
   4  /
下面是matdyn文件内容：
   1  &input
   2          asr='simple', amass(1)=24.3050, amass(2)=10.811,
   3                flfrc='PHMgB2.fc', flfrq='PHMgB2.freq', dos=.false.
   4                      /
   5  101
   6 0.000000  0.000000  0.000000  1.00
   7 0.000000  0.025000  0.000000  1.00
   8 0.000000  0.050000  0.000000  1.00
   9 0.000000  0.075000  0.000000  1.00
   10 0.000000  0.100000  0.000000  1.00