| 查看: 3078 | 回复: 7 | ||
[求助]
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 |
回复此楼» 猜你喜欢
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
帮忙看看fileCode
已经有5人回复
-
FileCode能看出啥？
已经有13人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
magicmonk
至尊木虫 (著名写手)
- 应助: 172 (高中生)
- 金币: 11451.3
- 红花: 24
- 帖子: 1124
- 在线: 3108.8小时
- 虫号: 1191315
- 注册: 2011-01-17
- 专业: 理论和计算化学
2楼2016-04-16 13:14:24
3楼2016-04-16 17:32:27
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
4楼2016-04-16 19:18:56
5楼2016-04-18 21:14:45
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
6楼2016-04-18 21:19:00
|
下面是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 |
7楼2016-04-18 22:30:44
8楼2016-04-22 10:00:22
100