24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1284  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

[求助] 高手指点下,声子谱

在计算ZnS B4结构的声子谱时候,发现自己计算的不对,请高手指点!

计算出来的结果,在红线的位置不分裂

文献结果

phonon.in文件如下phonons of ZnSB4
&inputph
  tr2_ph=1.0d-14,
  prefix='ZnSB4',
  epsil=.true.,
  amass(1)=32.065,
  amass(2)=65.409,
  outdir='./tmp'
  fildvscf='ZnSB4.dyn',
  ldisp=.true.
  nq1= 6,  nq2= 6,  nq3= 4,     
/
scf文件如下&control
    calculation='scf'
    restart_mode='from_scratch',
    prefix='ZnSB4',
    pseudo_dir = '~/pwscf/pseudo',
    outdir='./tmp'
    tstress = .true. ,
    tprnfor = .true. ,
    etot_conv_thr = 1.0E-6  ,
    forc_conv_thr = 1.0D-5 ,
/
&system   
    ibrav=  4,
    celldm(1) =6.023827,
    celldm(3) =1.60727173,

    nat=  4,
    ntyp= 2,
    ecutwfc = 40,
    ecutrho = 400.0 ,
/
/
&ELECTRONS
   conv_thr  =  1.0d-8  ,
   mixing_beta = 0.7 ,
/
&IONS
/
&CELL
/

ATOMIC_SPECIES
       S     32.065    S.pz-van_ak.UPF
      Zn     65.409   Zn.pz-van_ak.UPF
ATOMIC_POSITIONS (crystal)
S        0.333333333   0.666666667   0.376638686
S       -0.333333333  -0.666666667   0.876638686
Zn       0.333333333   0.666666667  -0.001638686
Zn      -0.333333333  -0.666666667   0.498361314
K_POINTS AUTOMATIC
16 16 8 0 0 0

[ Last edited by w76990480 on 2011-11-7 at 08:21 ]
回复此楼
专业文献求助
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

引用回帖:
3楼: Originally posted by gemucai at 2011-11-07 10:42:00:
gamma点的振动频率本身就不太准确,还好结构是对的,对称性没被破坏。

感觉你的结构在C方向上分得有点开,导致了在这个方向上的相互作用比较小,才会出现类似局域的即平行的声子谱。个人意见。

谢谢,有没有什么解决的方法呢?
专业文献求助
4楼2011-11-07 10:46:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

hanyanli0475

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★
w76990480(金币+5): 谢谢,一会贴出来帮我看下 2011-11-07 08:16:25
franch(金币+3): 谢谢回帖交流 2011-11-07 15:26:14
引用回帖:
1楼: Originally posted by w76990480 at 2011-11-06 20:54:35:
在计算ZnS B4结构的声子谱时候,发现自己计算的不对,请高手指点!

计算出来的结果,在红线的位置不分裂
[eimg]03/79/942733_1320584101_ ...

声子计算,那么trans也得设置为true啊,还有,自洽文件也很重要,你也应该贴出来,我自己估计,应该把k-mesh取大点,试试。
什么时候能算完!!!!
2楼2011-11-07 08:00:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

★ ★
franch(金币+2): 谢谢回帖交流 2011-11-07 15:26:31
gamma点的振动频率本身就不太准确,还好结构是对的,对称性没被破坏。

感觉你的结构在C方向上分得有点开,导致了在这个方向上的相互作用比较小,才会出现类似局域的即平行的声子谱。个人意见。
3楼2011-11-07 10:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

★ ★
贺仪(金币+2): 多谢建议 2011-11-09 08:39:13
引用回帖:
4楼: Originally posted by w76990480 at 2011-11-07 10:46:01:
谢谢,有没有什么解决的方法呢?

试着改一改晶格结构,将c方向的长度改小一点。

还有你用的是什么赝势,你算的结果总体偏小,有可能跟赝势的选择有关。
5楼2011-11-08 19:54:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +22 ORCHID1 2026-04-10 25/1250 2026-04-13 23:11 by pies112
[考研] 08工学 309分求调剂 +7 Yin DY 2026-04-08 7/350 2026-04-13 22:50 by pies112
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[材料工程] 材料调剂推荐 +9 蛋糕x2 2026-04-07 9/450 2026-04-13 09:07 by lhj2009
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 电子信息279求调剂,有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 复试调剂,一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 软件工程求调剂22软工296分求调剂,接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
请填处理意见