24小时热门版块排行榜    

查看: 1476  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
/
换了几个赝势,计算结果都问题,请问:问题出现在哪啊?我都急死了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 306求调剂 +4 来好运来来来 2026-03-22 4/200 2026-03-22 12:43 by 素颜倾城1988
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
请填处理意见