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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 高手指点下，声子谱
汕头大学海洋科学接受调剂
91/1返回列表
查看: 1292  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 本帖@通知

@hanyanli0475

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
专业文献求助
1楼 2011-11-06 20:54:35
已阅   回复此楼   关注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取大点，试试。
赞一下(1人) 回复此楼
什么时候能算完！！！！
2楼2011-11-07 08:00:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

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

感觉你的结构在C方向上分得有点开，导致了在这个方向上的相互作用比较小，才会出现类似局域的即平行的声子谱。个人意见。
赞一下(1人) 回复此楼
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方向的长度改小一点。

还有你用的是什么赝势，你算的结果总体偏小，有可能跟赝势的选择有关。
赞一下(1人) 回复此楼
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都可以～
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
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一下结果吧，如果好的话，别人也都可以用。如果不好的话，也不要耽误别人的时间了嘛～
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
9楼2011-11-09 14:24:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 w76990480 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010，320求调剂 +4 困困困困坤坤 2026-04-14 5/250 2026-04-14 21:29 by 哆啦A梦只是个梦
[基金申请] 山东省基金2026 +6 jerry681 2026-04-08 8/400 2026-04-14 19:28 by laobibibi
[考研] 22408 312求调剂 +15 门路摸摸 2026-04-14 16/800 2026-04-14 17:15 by zs92450
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 346分，工科0854求调剂，专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 求调剂，一志愿材料科学与工程985，365分， +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 267求调剂 +8 再忙也要吃饭啊 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭