24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1285  |  回复: 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的回帖

gemucai

木虫 (正式写手)

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

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工329分求调剂 +18 叶zilin 2026-04-13 19/950 2026-04-13 22:29 by pies112
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 305求调剂 +6 77Qi 2026-04-07 6/300 2026-04-12 02:30 by 秋豆菜芽
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 求调剂 +6 小聂爱学习 2026-04-11 9/450 2026-04-11 21:20 by 蓝云思雨
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 复试调剂,一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
信息提示
请填处理意见