24小时热门版块排行榜    

查看: 2153  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 267一志愿南京工业大学0817化工求调剂 +8 SUICHILD 2026-03-12 8/400 2026-03-18 20:55 by winsuccess
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见