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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问谁用过wannier90？急需帮助。
51/1返回列表
查看: 4235  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +5 @taotao 2026-03-19 5/250 2026-03-24 21:07 by greychen00
[考研] 0703化学调剂，求导师收 +7 天天好运来上岸� 2026-03-24 7/350 2026-03-24 20:26 by peike
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境） +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
信息提示
关闭
请填处理意见
关闭