查看: 5973  |  回复: 16

skydark

木虫 (正式写手)

[求助] wannier90得不到hr.dat檔 已有1人参与

各位高手們:

在下希望用VASP接wannier90計算出系統的onsite energy hopping energy 等參數
據說這些參數在hr.dat檔案裏頭

我看手冊說  只要在.win檔裡頭 加上hr_plot = .true. 就可以得到hr.dat檔

但我算完之後  ?K沒有出現hr.dat檔   試了好幾天 依然不知原因何在
只好在此向各位高手請教

以下是我的操作:

我用VASP算完自洽那一步後
把CHGCAR WAVEDER等檔案copy到新資料夾
然後接著跑一次wannier90

以下是我的INCAR 與 .win檔

INCAR:

ALGO = GWO
LSPECTRAL = .TRUE.
NOMEGA = 50
NBANDS = 240
LWANNIER90 = .TRUE.


WTe2.win檔 :

num_bands         =   240
num_wann          =   24
exclude_bands     =   1-60

##dis_win_max       =   12.5
##dis_win_min       =   -7.5
##dis_froz_max      =   10
##dis_froz_min      =   0.0
##dis_num_iter      =   2000
##dis_mix_ratio     =   1.d0

num_iter          = 0
##num_print_cycles  = 1

guiding_centres=true

begin projections
W  : s ; p ; d
Te : s ; p ; d
end projections

spinors = .true.

hr_plot = .true.

!restart = wannierise
!restart = plot

#bands_plot         =   true
#bands_plot_format = xmgrace
#bands_plot_project = 1-20 # W-d

#begin kpoint_path
#X 0.5   0.0   0.0  G  0.0    0.0    0.0
#G 0.0   0.0   0.0  Y  0.0    0.5    0.0
#end kpoint_path

begin unit_cell_cart
     3.4769999999999999    0.0000000000000000    0.0000000000000000
     0.0000000000000000    6.2489999999999997    0.0000000000000000
     0.0000000000000000    0.0000000000000000   14.0180000000000007
end unit_cell_cart


begin atoms_frac
W   0.0000000000000000  0.6006199999999993  0.5000000000000000
W   0.5000000000000000  0.3993800000000007  0.0000000000000000
W   0.0000000000000000  0.0397999999999996  0.0152199999999993
W   0.5000000000000000  0.9602000000000004  0.5152199999999993
Te  0.0000000000000000  0.8576100000000011  0.6552500000000023
Te  0.5000000000000000  0.1423899999999989  0.1552500000000023
Te  0.0000000000000000  0.6463099999999997  0.1111199999999997
Te  0.5000000000000000  0.3536900000000003  0.6111199999999997
Te  0.0000000000000000  0.2984500000000025  0.8598300000000023
Te  0.5000000000000000  0.7015499999999975  0.3598300000000023
Te  0.0000000000000000  0.2072199999999995  0.4038699999999977
Te  0.5000000000000000  0.7927800000000005  0.9038699999999977
end atoms_frac


mp_grid = 8 6 4
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

golddoushi

木虫 (正式写手)

★ ★ ★ ★ ★
skydark(蛋蛋小童鞋代发): 金币+5, 谢谢交流,小木虫鼓励积极回帖 2016-09-29 21:50:00
projection 的问题是你自己定义的轨道数目和万尼尔函数的个数不一致
hr 参数则需要在wannier90.win 里面指定才会生效,vasp只读取这个,然后挂wannier作为lib,和qe并不一样

考虑soc,则是简单的万尼尔函数数目乘以二
如果是共线的自旋计算,则up,dn,分别计算就可以
10楼2016-09-28 20:31:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

漫天飘雪

版主 (知名作家)

认真做事,踏实做人

好像低版本的没有,你试试高版本的。
http://blog.sciencenet.cn/u/pfliu89
2楼2016-09-25 21:08:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

???????:
2?: Originally posted by ?????? at 2016-09-25 21:08:09
?????汾????У????????汾???

????汾???
?????l?F?????????? wannier90.x ????
???^??F??????????????}

???? ???? wannier90.x ??
???F param_get_projections: Fewer projections defined than the number of Wannier functions requested

?????????N?
3楼2016-09-26 07:53:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NickyHCQ

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
skydark: 金币+5, ★有帮助 2016-10-02 17:57:52
skydark: 金币+5, ★有帮助 2016-10-02 17:58:55
???????:
3?: Originally posted by skydark at 2016-09-26 07:53:22
????汾???
?????l?F?????????? wannier90.x ????
???^??F??????????????}

???? ???? wannier90.x ??
???F param_get_projections: Fewer projections defined than the number of Wannier functions  ...

????????????vasp.wannier   ??????????У??????????о???ò??????qe?????????????????????????????vaspWannier???м?????????????amn.mmn.wout.win.eig??

????С???Android?????
4楼2016-09-28 14:31:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NickyHCQ

木虫 (正式写手)

【答案】应助回帖

???????:
3?: Originally posted by skydark at 2016-09-26 07:53:22
????汾???
?????l?F?????????? wannier90.x ????
???^??F??????????????}

???? ???? wannier90.x ??
???F param_get_projections: Fewer projections defined than the number of Wannier functions  ...

?????????????????projection?????漰?????

????С???Android?????
5楼2016-09-28 14:33:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

引用回帖:
4楼: Originally posted by NickyHCQ at 2016-09-28 14:31:37
我试了一个月的vasp.wannier   这个对接还不行,,,好多功能感觉是得不到的,qe就轻松得到。。。。。楼主,你跑完vaspWannier是有几个输出文件??amn.mmn.wout.win.eig?
...

沒amn  其他都有
6楼2016-09-28 15:47:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NickyHCQ

木虫 (正式写手)

引用回帖:
6楼: Originally posted by skydark at 2016-09-28 15:47:43
沒amn  其他都有...

对的,我暂时认为不会出这个amn,,不管怎么弄,而且spn也不会出来,,软件有待改善

发自小木虫Android客户端
7楼2016-09-28 16:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

引用回帖:
7楼: Originally posted by NickyHCQ at 2016-09-28 16:52:50
对的,我暂时认为不会出这个amn,,不管怎么弄,而且spn也不会出来,,软件有待改善
...

projection的問題解決 現在遇到跟你一樣的問題
Error: Problem opening input file WTe2.amn
Error: examine the output/error file for details
8楼2016-09-28 16:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NickyHCQ

木虫 (正式写手)

引用回帖:
8楼: Originally posted by skydark at 2016-09-28 16:59:53
projection的問題解決 現在遇到跟你一樣的問題
Error: Problem opening input file WTe2.amn
Error: examine the output/error file for details...

这个问题我已经好久了,还有.spn缺少~~这个vaspWannier到处是问题,,

发自小木虫Android客户端
9楼2016-09-28 18:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skydark 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见