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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 高手指点下，声子谱
91/1返回列表
查看: 1563  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
关闭
请填处理意见
关闭