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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 求助ZnO单层声子计算
51/1返回列表
查看: 1538  |  回复: 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的回帖

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

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的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料相关专业344求调剂双非工科学校或课题组 +23 hualkop 2026-04-12 25/1250 2026-04-16 22:12 by SUSE_CL
[考研] 化学070300 求调剂 +28 哈哈哈^_^ 2026-04-12 28/1400 2026-04-16 21:36 by 大力水手力大无�
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 271求调剂 +35 2261744733 2026-04-11 37/1850 2026-04-16 13:16 by jiahl2024
[考研] 药学专硕调剂 +9 ? 一路生?花? 2026-04-10 11/550 2026-04-16 10:49 by noqvsozv
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 人工智能320调剂08工类还有机会吗 +18 振—TZ 2026-04-10 19/950 2026-04-14 10:34 by screening
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 0859，337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 一志愿京区985，085401，与本科专业一致，电子信息工程， +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
信息提示
关闭
请填处理意见
关闭