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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 求助ZnO单层声子计算
91/1返回列表
查看: 1346  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭