24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1300  |  回复: 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的回帖

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

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

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
w76990480(金币+5): 谢谢指导 2011-11-09 09:14:50
zzy870720z(金币+2): 鼓励交流 2011-11-09 12:31:48
据goldenfisher说,ak做的赝势的质量普遍不高~LZ可以试试下面这两个赝势,计算时在&system中指定input_dft为PBE或者PW91都可以~
思想重于技巧,内涵重于表象
6楼2011-11-09 09:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

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

引用回帖:
6楼: Originally posted by souledge at 2011-11-09 09:10:49:
据goldenfisher说,ak做的赝势的质量普遍不高~LZ可以试试下面这两个赝势,计算时在&system中指定input_dft为PBE或者PW91都可以~

这两个赝势我在库里面没有找到啊?还是你把名字改了?
专业文献求助
7楼2011-11-09 10:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

w76990480(金币+5): 有O的没有,在给我个!谢谢 2011-11-09 14:02:30
引用回帖:
7楼: Originally posted by w76990480 at 2011-11-09 10:03:39:
这两个赝势我在库里面没有找到啊?还是你把名字改了?

先试试结果再说~
思想重于技巧,内涵重于表象
8楼2011-11-09 13:22:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
w76990480(金币+5): 谢谢指导,我会尽快给出计算结果的,谢谢! 2011-11-09 15:06:19
贺仪(金币+2): 多谢指教! 2011-11-09 17:50:42
LZ贴一贴结果啊,事实上因为我自己的计算跟这几个元素没关系,所以我也没测试过,只是觉得理论上这几个赝势应该质量尚可。当然,跟我的计算相关的我是都测试了,效果不能说非常出众,但是感觉已经足够好了~
所以LZ就show一下结果吧,如果好的话,别人也都可以用。如果不好的话,也不要耽误别人的时间了嘛~
思想重于技巧,内涵重于表象
9楼2011-11-09 14:24:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 w76990480 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 322求调剂 +7 123安康 2026-04-12 14/700 2026-04-15 21:48 by noqvsozv
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 071000生物学调剂求助 +18 zzzzwww 2026-04-09 21/1050 2026-04-14 15:39 by zs92450
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
信息提示
请填处理意见