24小时热门版块排行榜    

查看: 2162  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

soiu2000

银虫 (初入文坛)

[求助] ZnO的声子计算

最近新学pwscf,想计算ZnO Gamma点声子。但是怎么算都和文献上差很远,我有的文献上看到别人可以算的很接近。不知道是什么原因。还请各位高手指点一二。
附上输入文件。

omega(  1 -  1) =        -15.8  [cm-1]   --> A_1             I+R
omega(  2 -  3) =         41.6  [cm-1]   --> E_1             I+R
omega(  4 -  5) =        103.5  [cm-1]   --> E_2             R  
omega(  6 -  6) =        234.2  [cm-1]   --> B_2               
omega(  7 -  8) =        315.2  [cm-1]   --> E_1             I+R
omega(  9 -  9) =        323.1  [cm-1]   --> A_1             I+R
omega( 10 - 11) =        369.1  [cm-1]   --> E_2             R  
omega( 12 - 12) =        493.8  [cm-1]   --> B_2               

而文献上
omega LO(A1)=574 [cm-1]
omega LO(E1)=583 [cm-1]
omega TO(A1)=380 [cm-1]
omega TO(E1)=407 [cm-1]
omega      (E2)=437 [cm-1]
omega      (E2)=101 [cm-1]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-10-31 21:27:21
首先,我不明白你的pw的in文件中ecutwfc要设置这么高,K点网格又相对比较低?ph的in文件中你只算Gamma点,为什么还要写q网格?还有,你的结果中,有一点点虚频的存在,前面的relax完全了吗?

最后,我想问你下,你用pwscf 4.3.1的并行效率这么样?我自己最近也刚装4.3.1 怎么感觉并行效率很低啊。CPU使用率都只有40%左右。你的呢?
3楼2011-10-28 18:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

ridge8888

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-10-31 21:27:34
&electrons
    conv_thr =  1.0d-10
    mixing_beta = 0.7
    electron_maxstep=200
/
ATOMIC_SPECIES
O  15.999  O.pz-van_ak.UPF
Zn 65.39  Zn.pz-van_ak.UPF
ATOMIC_POSITIONS {crystal}
Zn       0.000000000   0.000000000  -0.000969838
Zn       0.666666667   0.333333333   0.499030162
O        0.000000000   0.000000000   0.378963856
O        0.666666667   0.333333333   0.878963856



ATOMIC_SPECIES 和 ATOMIC_POSITIONS {crystal}中的原子的顺序,颠倒了。是否可以颠倒呢?
2楼2011-10-28 14:43:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-10-31 21:27:57
引用回帖:
1楼: Originally posted by soiu2000 at 2011-10-28 13:59:38:
最近新学pwscf,想计算ZnO Gamma点声子。但是怎么算都和文献上差很远,我有的文献上看到别人可以算的很接近。不知道是什么原因。还请各位高手指点一二。
附上输入文件。

omega(  1 -  1) =        -15.8  [cm ...

感觉你的 ecutwfc = 500.0,   ecutrho = 1200.0设置太高了,一般来说要做个收敛测试,获取适当的截断。就你的情况,以我的经验,你分别设置成30,350试试。
4楼2011-10-31 19:24:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-10-31 21:28:08
引用回帖:
3楼: Originally posted by fzuvivi at 2011-10-28 18:54:28:
首先,我不明白你的pw的in文件中ecutwfc要设置这么高,K点网格又相对比较低?ph的in文件中你只算Gamma点,为什么还要写q网格?还有,你的结果中,有一点点虚频的存在,前面的relax完全了吗?

最后,我想问你下 ...

显然是你安装的问题。就我的情况看,CPU利用率一般在97~99%。
5楼2011-10-31 19:25:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 286分人工智能专业请求调剂愿意跨考! +3 lemonzzn 2026-03-17 3/150 2026-03-19 23:58 by 23Postgrad
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +5 吃吃吃才有意义 2026-03-19 5/250 2026-03-19 16:18 by 30660438
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 材料专硕274一志愿陕西师范大学求调剂 +8 薛云鹏 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
信息提示
请填处理意见