应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 求助ZnO单层声子计算
91/1返回列表
查看: 1515  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

quby7503

银虫 (著名写手)

[求助] 求助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
/
换了几个赝势，计算结果都问题，请问：问题出现在哪啊？我都急死了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-03 10:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)
怎么不对了？没看出问题。
赞一下 回复此楼
2楼2013-11-03 11:09:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)
引用回帖:
2楼: Originally posted by huai at 2013-11-03 11:09:43
怎么不对了？没看出问题。

前三个频率对应着声学支，在gamma点时，频率应该为0，即使算上数值误差的话，也应该是个很小很小的数值才是。但我计算的频率都太大了
赞一下 回复此楼
3楼2013-11-03 11:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)
★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:18
引用回帖:
3楼: Originally posted by quby7503 at 2013-11-03 11:14:13
前三个频率对应着声学支，在gamma点时，频率应该为0，即使算上数值误差的话，也应该是个很小很小的数值才是。但我计算的频率都太大了...

我没有太好的建议，我认为这是数值误差。前三个比起后三个要小很多。我只是做过对角化动力学矩阵得到声子谱。对pwscf不了解，你问一下他们用过的人，看虚频是多大的量级。
赞一下 回复此楼
4楼2013-11-03 16:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-04 12:24:31
franch: 金币+2, 谢谢回帖交流 2013-11-04 18:19:29
计算表面声子比较复杂，你是否考虑了不同真空层厚度。
赞一下 回复此楼
未来就是现在
5楼2013-11-04 08:28:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)
引用回帖:
5楼: Originally posted by future_wl at 2013-11-04 08:28:04
计算表面声子比较复杂，你是否考虑了不同真空层厚度。

谢谢，真空层是经过严格优化的。
赞一下 回复此楼
6楼2013-11-04 12:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倔强的坚果墙

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
quby7503: 金币+3 2013-11-05 12:01:58
没有计算过ZnO，不过我建议你从以下几个方面入手：1.设置encut的值（需要测试）2.检查ecutwfc设置是否足够（印象中graphene的设置达到500Ry）3.K点是否测试达到收敛 4.结构是否充分弛豫
赞一下 回复此楼
像小强一样活着
7楼2013-11-04 21:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)
引用回帖:
7楼: Originally posted by 倔强的坚果墙 at 2013-11-04 21:05:45
没有计算过ZnO，不过我建议你从以下几个方面入手：1.设置encut的值（需要测试）2.检查ecutwfc设置是否足够（印象中graphene的设置达到500Ry）3.K点是否测试达到收敛 4.结构是否充分弛豫

谢谢，都是经过测试的
赞一下 回复此楼
8楼2013-11-05 12:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倔强的坚果墙

铁杆木虫 (正式写手)
引用回帖:
8楼: Originally posted by quby7503 at 2013-11-05 12:01:50
谢谢，都是经过测试的...

那运用声学支求和规则画声子谱看看，假如虚频不是很大，也是可以接受的
赞一下 回复此楼
像小强一样活着
9楼2013-11-05 16:36:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 quby7503 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学专硕调剂 +7 ? 一路生?花? 2026-04-10 8/400 2026-04-11 16:44 by 加一梨梨
[考研] 求调剂 +13 张番茄不炒蛋 2026-04-10 14/700 2026-04-11 16:02 by vgtyfty
[考研] 316求调剂 +4 想读研究生( ?∵ 2026-04-07 4/200 2026-04-11 11:43 by zhq0425
[考研] 296求调剂 +6 汪！？！ 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 0854调剂 +11 长弓傲 2026-04-09 12/600 2026-04-11 11:16 by zhq0425
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 材料与化工调剂 10+11 下一站上岸@ 2026-04-10 36/1800 2026-04-11 10:26 by 89436494
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 346，工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 085600材料与化工301分求调剂院校 +33 刺痛jk 2026-04-06 34/1700 2026-04-09 18:31 by hy861222
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 071000生物学，一志愿深圳大学296分，求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 工科370求调剂 +3 溏心煎鸡蛋 2026-04-05 3/150 2026-04-06 10:55 by 这是一个无聊的�
信息提示
关闭
请填处理意见
关闭