| 查看: 1747 | 回复: 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 / 换了几个赝势，计算结果都问题，请问：问题出现在哪啊？我都急死了 |
回复此楼» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 单层和多层石墨电子谱和声子谱研究
已经有66人回复
- 【文献学习】Li,Na掺杂ZnO空穴的杂化密度泛函计算(PRB，2009)
已经有21人回复
- ZnO的声子计算
已经有18人回复
倔强的坚果墙
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 8437.1
- 散金: 52
- 红花: 12
- 帖子: 463
- 在线: 878.6小时
- 虫号: 1251958
- 注册: 2011-04-01
- 专业: 凝聚态物性 II ：电子结构
9楼2013-11-05 16:36:05
huai
木虫 (职业作家)
- 应助: 4 (幼儿园)
- 金币: 1048.1
- 散金: 118
- 红花: 6
- 帖子: 4522
- 在线: 554小时
- 虫号: 420388
- 注册: 2007-07-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2013-11-03 11:09:43
3楼2013-11-03 11:14:13
huai
木虫 (职业作家)
- 应助: 4 (幼儿园)
- 金币: 1048.1
- 散金: 118
- 红花: 6
- 帖子: 4522
- 在线: 554小时
- 虫号: 420388
- 注册: 2007-07-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2013-11-03 16:41:39
