| 查看: 1353 | 回复: 8 | |||||||||||||||
[求助]
求助ZnO单层声子计算
|
|
刚刚接触PWSCF，用ZnO单层做尝试，计算的声子总是不对，在gamma点处的计算结果如下： omega( 1) = -4.186512 [THz] = -139.647010 [cm-1] ( 0.666757 0.000000 0.019032 0.000000 0.000000 0.000000 ) ( 0.744729 0.000000 0.021257 0.000000 0.000000 0.000000 ) omega( 2) = -4.186512 [THz] = -139.647010 [cm-1] ( 0.019032 0.000000 -0.666757 0.000000 -0.000000 0.000000 ) ( 0.021257 0.000000 -0.744729 0.000000 0.000000 0.000000 ) omega( 3) = 1.330295 [THz] = 44.373869 [cm-1] ( -0.000000 0.000000 -0.000000 0.000000 -0.815519 0.000000 ) ( 0.000000 0.000000 0.000000 0.000000 -0.578730 0.000000 ) omega( 4) = 7.850469 [THz] = 261.863452 [cm-1] ( -0.000000 0.000000 -0.000000 0.000000 -0.171074 0.000000 ) ( 0.000000 0.000000 0.000000 0.000000 0.985258 0.000000 ) omega( 5) = 17.606727 [THz] = 587.297190 [cm-1] ( 0.000819 0.000000 0.263621 0.000000 -0.000000 0.000000 ) ( -0.002996 0.000000 -0.964621 0.000000 0.000000 0.000000 ) omega( 6) = 17.606727 [THz] = 587.297190 [cm-1] ( 0.263621 0.000000 -0.000819 0.000000 -0.000000 0.000000 ) ( -0.964621 0.000000 0.002996 0.000000 0.000000 0.000000 ) 计算过程是，首先做自洽计算（原子位置已经经过严格的优化） &control calculation = 'scf' restart_mode='from_scratch', prefix='zno-001', tstress = .true. tprnfor = .true. pseudo_dir = '/home/LAB/byqu/Pseudpo-pwscf', outdir='./tmp' / &system ibrav= 4, a= 3.28 , c=15.0, nat= 2, ntyp= 2, ecutwfc = 325.0 , / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES Zn 65.39 Zn.pbe-d-hgh.UPF O 15.9994 O.pbe-hgh.UPF ATOMIC_POSITIONS {crystal} Zn 0.6666666699999979 0.3333333300000021 0.5000000000000000 O 0.3333333300000021 0.6666666699999979 0.5000000000000000 K_POINTS {automatic} 9 9 1 0 0 0 然后做声子计算 phonons for ZnO-monolayer &inputph tr2_ph = 1.0d-12 prefix = 'zno-001' outdir = './tmp' amass(1) = 65.39 amass(2) = 15.9994 fildyn = 'zno-001.dyn' trans = .true. ldisp = .true. nq1=15, nq2=15, nq3=1 / 换了几个赝势，计算结果都问题，请问：问题出现在哪啊？我都急死了 |
回复此楼» 猜你喜欢
谈谈两天一夜的“延安行”
已经有15人回复
-
EST投稿状态问题
已经有6人回复
-
职称评审没过，求安慰
已经有15人回复
-
垃圾破二本职称评审标准
已经有11人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有16人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 单层和多层石墨电子谱和声子谱研究
已经有66人回复
- 【文献学习】Li,Na掺杂ZnO空穴的杂化密度泛函计算(PRB，2009)
已经有21人回复
- ZnO的声子计算
已经有18人回复
★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:18
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:18
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-04 12:24:31
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:29
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-04 12:24:31
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:29
倔强的坚果墙
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 8437.1
- 散金: 52
- 红花: 12
- 帖子: 462
- 在线: 878小时
- 虫号: 1251958
- 注册: 2011-04-01
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-05 12:01:58
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-05 12:01:58
倔强的坚果墙
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 8437.1
- 散金: 52
- 红花: 12
- 帖子: 462
- 在线: 878小时
- 虫号: 1251958
- 注册: 2011-04-01
- 专业: 凝聚态物性 II ：电子结构