24小时热门版块排行榜    

查看: 4169  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

对马拉二胡

至尊木虫 (职业作家)

[求助] 请问谁用过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 这几个输入文件?
求使用过这个软件的虫子指导指导,万分感谢!
回复此楼
相互学习,共同进步。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)

引用回帖:
6楼: Originally posted by awdawdawdd at 2015-05-18 23:44:28
你是计算VASPWIKI上的Si_GW0的能带的例子吧
第三部你先不用把NELM的值调那么高,一般NELM=4就可以了,推荐你先用NELM=1 (G0W0)来试一试,等算完了应该会有wannier90.amn, wannier90.mmn, wannier90.eig这几个文 ...

又过了一天,我的vasprun.xml 中没有出现 <i type="logical" name="LWANNIER90" > T  </i> , 但是编译是按你说的来的,真是醉了。
相互学习,共同进步。
7楼2015-05-19 19:12:21
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:47 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:12 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见