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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 求助ZnO单层声子计算
91/1返回列表
查看: 1518  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +12 wzzz王 2026-04-09 12/600 2026-04-11 22:36 by dongdian1
[考研] 26自然地理学303分求调剂 +6 一战成硕啊啊啊� 2026-04-06 11/550 2026-04-11 21:27 by 裴宏伟
[考研] 266求调剂，一志愿哈工程电子信息，本科获多项国奖和省奖 +8 lumine1 2026-04-06 8/400 2026-04-11 18:35 by 逆水乘风
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 一志愿京区985，085401，与本科专业一致，电子信息工程， +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 环境专硕调剂 +16 会说话的肘子 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[考研] 301求调剂 +6 静静想想 2026-04-05 6/300 2026-04-10 09:15 by Delta2012
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 312求调剂 +18 gtw1 2026-04-06 20/1000 2026-04-07 18:16 by 蓝云思雨
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
信息提示
关闭
请填处理意见
关闭