24小时热门版块排行榜    

查看: 2144  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 化学工程321分求调剂 +11 大米饭! 2026-03-15 14/700 2026-03-17 17:11 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
请填处理意见