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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问谁用过wannier90？急需帮助。
51/1返回列表
查看: 4599  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

对马拉二胡

至尊木虫 (职业作家)

[求助] 请问谁用过wannier90？急需帮助。 已有2人参与

想做结构的费米面，查了下wannier90这个软件可以实现。安装了wannier90 2.0版本
按照官方教程计算Si的能带。
第一步先优化，静态计算。 第二步 计算DFT virtual orbital 。 第三步结合wannier90计算

前两步都没问题，但是到了第三步vasp一直在运行（24h）PS：我用的是40核 96G的服务器。以下已经输出的文件
CHG       PCDAT        W0004.tmp  W0014.tmp  W0024.tmp  W0034.tmp  W0044.tmp  W0054.tmp  W0064.tmp
CHGCAR    POSCAR       W0005.tmp  W0015.tmp  W0025.tmp  W0035.tmp  W0045.tmp  W0055.tmp  wannier90.win
CONTCAR   POTCAR       W0006.tmp  W0016.tmp  W0026.tmp  W0036.tmp  W0046.tmp  W0056.tmp  WAVECAR
DOSCAR    PROCAR       W0007.tmp  W0017.tmp  W0027.tmp  W0037.tmp  W0047.tmp  W0057.tmp  WAVECAR.chi
EIGENVAL  quence       W0008.tmp  W0018.tmp  W0028.tmp  W0038.tmp  W0048.tmp  W0058.tmp  WAVEDER
INCAR     script.out   W0009.tmp  W0019.tmp  W0029.tmp  W0039.tmp  W0049.tmp  W0059.tmp  WFULL0001.tmp
INCAR-1   vasprun.xml  W0010.tmp  W0020.tmp  W0030.tmp  W0040.tmp  W0050.tmp  W0060.tmp  XDATCAR
KPOINTS   W0001.tmp    W0011.tmp  W0021.tmp  W0031.tmp  W0041.tmp  W0051.tmp  W0061.tmp
OSZICAR   W0002.tmp    W0012.tmp  W0022.tmp  W0032.tmp  W0042.tmp  W0052.tmp  W0062.tmp
OUTCAR    W0003.tmp    W0013.tmp  W0023.tmp  W0033.tmp  W0043.tmp  W0053.tmp  W0063.tmp

按道理应该很快就可以算完，请问这是出了什么问题呢？ 另外我还想问下，加入LWANNIER90 = .TRUE. 这个参数后，是不是会产生wannier90.amn, wannier90.mmn, wannier90.eig 这几个输入文件？
求使用过这个软件的虫子指导指导，万分感谢！
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 wannier 90 vasp+wannier90

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
相互学习，共同进步。
1楼 2015-05-16 21:25:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)
引用回帖:
4楼: Originally posted by awdawdawdd at 2015-05-17 19:03:58
没有啥文件吧。。。假定你wannier90的source code文件夹放在和VASP的source code的文件夹的同一目录的地方（也就是在VASP的文件夹用 cd  ../wannier90-1.2 可以进到wannier90的文件夹里面），然后假定你已经变异了 ...

谢谢，我按你写的重新编译了一下。但还是出现同样的问题
第一步
INCAR：
SYSTEM = Si
ENCUT = 450
ISTART =0
ICHARG=2
ISMEAR= 0
SIGMA=0.05
EDIFF=1E-6
EDIFFG=-0.01
NSW= 0
ISIF= 0
IBRION= -1
NELM=100
PREC = Accurate
PSTRESS=0
LREAL = .FALSE.
KPOINTS：
Automatic generation
0
Monhkorst-Pack
4 4 4
0.0 0.0 0.0
POSCAR：
system Si
5.430
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
2
cart
0.00 0.00 0.00
0.25 0.25 0.25
第二步
INCAR：
SYSTEM = Si
ENCUT = 450
ISTART =0
ICHARG=2
ISMEAR= 0
SIGMA=0.05
EDIFF=1E-6
EDIFFG=-0.01
NSW= 0
ISIF= 0
IBRION= -1
NELM=100
PREC = Accurate
PSTRESS=0
LREAL = .FALSE.
ALGO = Exact
NBANDS  = 64
LOPTICS=.TRUE.
NEDOS=2000
KPOINTS， POSCAR 文件与第一步相同
第三步
INCAR：
SYSTEM = Si
ENCUT = 450
ISTART =0
ICHARG=2
ISMEAR= 0
SIGMA=0.05
EDIFF=1E-6
EDIFFG=-0.01
NSW= 0
ISIF= 0
IBRION= -1
NELM=100
PREC = Accurate
PSTRESS=0
LREAL = .FALSE.
ALGO = GW0
LSPECTRAL = .TRUE.
NOMEGA = 50
NBANDS=64
LRPA = .FALSE.
LWANNIER90=.TRUE.

KPOINTS:
k-points along high symmetry lines
          64
Reciprocal
  0.00000000  0.00000000  0.00000000  1.562500e-02
  0.00000000  0.00000000  0.25000000  1.562500e-02
  0.00000000  0.00000000  0.50000000  1.562500e-02
  0.00000000  0.00000000  0.75000000  1.562500e-02
  0.00000000  0.25000000  0.00000000  1.562500e-02
  0.00000000  0.25000000  0.25000000  1.562500e-02
  0.00000000  0.25000000  0.50000000  1.562500e-02
  0.00000000  0.25000000  0.75000000  1.562500e-02
  0.00000000  0.50000000  0.00000000  1.562500e-02
  0.00000000  0.50000000  0.25000000  1.562500e-02
  0.00000000  0.50000000  0.50000000  1.562500e-02
  0.00000000  0.50000000  0.75000000  1.562500e-02
  0.00000000  0.75000000  0.00000000  1.562500e-02
  0.00000000  0.75000000  0.25000000  1.562500e-02
  0.00000000  0.75000000  0.50000000  1.562500e-02
  0.00000000  0.75000000  0.75000000  1.562500e-02
  0.25000000  0.00000000  0.00000000  1.562500e-02
  0.25000000  0.00000000  0.25000000  1.562500e-02
  0.25000000  0.00000000  0.50000000  1.562500e-02
  0.25000000  0.00000000  0.75000000  1.562500e-02
  0.25000000  0.25000000  0.00000000  1.562500e-02
  0.25000000  0.25000000  0.25000000  1.562500e-02
......
wannier.win  ：
num_bands         =   12      
num_wann          =   8

dis_win_max       = 17.0d0
dis_froz_max      =  6.4d0
dis_num_iter      =  120
dis_mix_ratio     = 1.d0

num_iter          = 50
num_print_cycles  = 10

Begin Atoms_Frac
Si  -0.25   0.75  -0.25
Si   0.00   0.00   0.00
End Atoms_Frac
   
Begin Projections     
Si : sp3
End Projections      
   
begin kpoint_path
L 0.50000  0.50000 0.5000 G 0.00000  0.00000 0.0000
G 0.00000  0.00000 0.0000 X 0.50000  0.00000 0.5000
X 0.50000 -0.50000 0.0000 K 0.37500 -0.37500 0.0000
K 0.37500 -0.37500 0.0000 G 0.00000  0.00000 0.0000
end kpoint_path


Begin Unit_Cell_Cart
-2.6988 0.0000 2.6988
0.0000 2.6988 2.6988
-2.6988 2.6988 0.0000
End Unit_Cell_Cart


mp_grid      = 4 4 4


begin kpoints
0.0000  0.0000   0.0000   
0.0000  0.2500   0.0000   
0.0000  0.5000   0.0000   
0.0000  0.7500   0.0000   
0.2500  0.0000   0.0000   
0.2500  0.2500   0.0000   
0.2500  0.5000   0.0000   
0.2500  0.7500   0.0000   
0.5000  0.0000   0.0000   
0.5000  0.2500   0.0000   
0.5000  0.5000   0.0000   
0.5000  0.7500   0.0000   
0.7500  0.0000   0.0000   
0.7500  0.2500   0.0000   
0.7500  0.5000   0.0000   
0.7500  0.7500   0.0000   
0.0000  0.0000   0.2500   
0.0000  0.2500   0.2500   
0.0000  0.5000   0.2500   
0.0000  0.7500   0.2500   
0.2500  0.0000   0.2500   
0.2500  0.2500   0.2500   
0.2500  0.5000   0.2500   
0.2500  0.7500   0.2500   
0.5000  0.0000   0.2500   
0.5000  0.2500   0.2500   
0.5000  0.5000   0.2500   
0.5000  0.7500   0.2500   
0.7500  0.0000   0.2500   
0.7500  0.2500   0.2500   
0.7500  0.5000   0.2500   
0.7500  0.7500   0.2500   
0.0000  0.0000   0.5000   
0.0000  0.2500   0.5000   
0.0000  0.5000   0.5000   
0.0000  0.7500   0.5000   
0.2500  0.0000   0.5000   
0.2500  0.2500   0.5000   
0.2500  0.5000   0.5000   
0.2500  0.7500   0.5000   
0.5000  0.0000   0.5000   
0.5000  0.2500   0.5000   
0.5000  0.5000   0.5000   
0.5000  0.7500   0.5000   
0.7500  0.0000   0.5000   
0.7500  0.2500   0.5000   
0.7500  0.5000   0.5000   
0.7500  0.7500   0.5000   
0.0000  0.0000   0.7500   
0.0000  0.2500   0.7500   
0.0000  0.5000   0.7500   
0.0000  0.7500   0.7500   
0.2500  0.0000   0.7500   
0.2500  0.2500   0.7500   
0.2500  0.5000   0.7500   
0.2500  0.7500   0.7500   
0.5000  0.0000   0.7500   
0.5000  0.2500   0.7500   
0.5000  0.5000   0.7500   
0.5000  0.7500   0.7500   
0.7500  0.0000   0.7500   
0.7500  0.2500   0.7500   
0.7500  0.5000   0.7500   
0.7500  0.7500   0.7500   
End Kpoints


以上是所有文件，实在找不出
赞一下 回复此楼
相互学习，共同进步。
5楼2015-05-18 17:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
对马拉二胡: 金币+20, ★★★很有帮助 2015-05-17 08:54:11
你把你的INCAR贴上来
另外，你重新编译了VASP没？

编译：
Mind: VASP needs to be compiled with the following additional precompiler flag:
-DVASP2WANNIER90
and the variable LIB in the makefile must contain an entry that points to libwannier.a.
For instance:
LIB     = -L../vasp.5.lib -ldmy  \
     ../vasp.5.lib/linpack_double.o $(SCA) $(LAPACK) $(BLAS)
might be changed to
LIB     = -L../vasp.5.lib -ldmy  \
     ../vasp.5.lib/linpack_double.o ../wannier90-1.2/libwannier.a $(SCA) $(LAPACK) $(BLAS)
depending on where you have installed WANNIER90, obviously.

ref：http://cms.mpi.univie.ac.at/wiki/index.php/LWANNIER90
赞一下 回复此楼
2楼2015-05-17 04:47:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)
引用回帖:
2楼: Originally posted by awdawdawdd at 2015-05-17 04:47:34
你把你的INCAR贴上来
另外，你重新编译了VASP没？

编译：
Mind: VASP needs to be compiled with the following additional precompiler flag:
-DVASP2WANNIER90
and the variable LIB in the makefile must ...

谢谢，我想应该是因为没有编译的原因。我用的是wannier90 2.0版本，是不是不能和VASP结合用啊？
另外，能看看你编译的LWANNIER90文件吗？
赞一下 回复此楼
相互学习，共同进步。
3楼2015-05-17 08:54:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
对马拉二胡: 金币+10, ★有帮助 2015-05-18 11:40:23
引用回帖:
3楼: Originally posted by 对马拉二胡 at 2015-05-17 08:54:03
谢谢，我想应该是因为没有编译的原因。我用的是wannier90 2.0版本，是不是不能和VASP结合用啊？
另外，能看看你编译的LWANNIER90文件吗？...

没有啥文件吧。。。假定你wannier90的source code文件夹放在和VASP的source code的文件夹的同一目录的地方（也就是在VASP的文件夹用 cd  ../wannier90-1.2 可以进到wannier90的文件夹里面），然后假定你已经变异了wannier90的library（得到了libwannier.a）

然后
1.在vasp的source code的Makefile里面找到类似于
LIB     = -L../vasp.5.lib -ldmy  \
     ../vasp.5.lib/linpack_double.o $(SCA) $(LAPACK) $(BLAS)

改成（或者说加上“../wannier90-1.2/libwannier.a”）使之变成类似于
LIB     = -L../vasp.5.lib -ldmy  \
     ../vasp.5.lib/linpack_double.o ../wannier90-1.2/libwannier.a $(SCA) $(LAPACK) $(BLAS)

2.在VASP的source code的Makefile里面找到你的编译的选项（Pre-compiler flags）加上-DVASP2WANNIER90，类似于这样：
CPP    = $(CPP_) -DMPI  -DHOST=\"LinuxIFC\" -DIFC \
     -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
     -DMPI_BLOCK=8000 -DRPROMU_DGEMV  -DRACCMU_DGEMV -DVASP2WANNIER90

3.然后再VASP的source code里面先make clean，然后make，等待编译完成就可以了，如果卡在了最后面，有可能是你VASP和wannier90的编译器用的不一样，一般把wannier90的编译器换成和VASP一样就可以了咯
赞一下 回复此楼
4楼2015-05-17 19:03:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 20:36 by xxniao123
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +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
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭