版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

【奖励】本帖被评价53次，作者gavinliu7390增加金币 42.4 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gavinliu7390

木虫 (著名写手)

[资源] Band structure of Si within GW and VASP

在 http://www.wannier.org/download.html 下载
wannier90-1.2
修改make.sys
如果装了ifort编译器
LIBS = -lmkl_intel_lp64  -lmkl_sequential -lmkl_core
然后 make all

之后修改vasp

LIB    = -L../vasp.5.lib -ldmy  \
   ../vasp.5.lib/linpack_double.o \
   ../wannier90-1.2/libwannier.a $(SCA) $(LAPACK) $(BLAS)

CPP = $(CPP_) -DMPI  -DHOST=\"LinuxIFC\" -DIFC \
   -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc  \
   -DMPI_BLOCK=8000 -Duse_collective -DscaLAPACK -DNGZhalf -DVASP2WANNIER90

之后重新make

分五步计算
1.
INCAR
ISMEAR =  0
SIGMA  =  0.05
2.
INCAR
ISMEAR =  0
SIGMA  =  0.05
ALGO = Exact
NBANDS  = 64
LOPTICS=.TRUE.
NEDOS=2000
3.
INCAR
ALGO = GW0 ; LSPECTRAL = .TRUE. ; NOMEGA = 50
NBANDS=64
LRPA = .FALSE.
LWANNIER90=.TRUE.
前三步运行vasp
4.
wannier90.win
num_wann=8
num_bands=8

# for GW uncomment
exclude_bands 9-64

Begin Projections
Si:sp3
End Projections

dis_froz_max=9
dis_num_iter=1000

guiding_centres=true

begin unit_cell_cart
   2.7150000    2.7150000    0.0000000
   0.0000000    2.7150000    2.7150000
   2.7150000    0.0000000    2.7150000
end unit_cell_cart

begin atoms_cart
Si    0.0000000    0.0000000    0.0000000
Si    1.3575000    1.3575000    1.3575000
end atoms_cart

mp_grid =    4    4    4

begin kpoints
   0.0000000    0.0000000    0.0000000
   0.2500000    0.0000000    0.0000000
   0.5000000    0.0000000    0.0000000
   0.2500000    0.2500000    0.0000000
   0.5000000    0.2500000    0.0000000
-0.2500000    0.2500000    0.0000000
   0.5000000    0.5000000    0.0000000
-0.2500000    0.5000000    0.2500000
   0.0000000    0.2500000    0.0000000
   0.0000000    0.0000000    0.2500000
-0.2500000 -0.2500000 -0.2500000
-0.2500000    0.0000000    0.0000000
   0.0000000 -0.2500000    0.0000000
   0.0000000    0.0000000 -0.2500000
   0.2500000    0.2500000    0.2500000
   0.0000000    0.5000000    0.0000000
   0.0000000    0.0000000    0.5000000
-0.5000000 -0.5000000 -0.5000000
   0.0000000    0.2500000    0.2500000
   0.2500000    0.0000000    0.2500000
-0.2500000 -0.2500000    0.0000000
-0.2500000    0.0000000 -0.2500000
   0.0000000 -0.2500000 -0.2500000
   0.0000000    0.5000000    0.2500000
   0.2500000    0.0000000    0.5000000
-0.2500000 -0.2500000    0.2500000
-0.5000000 -0.2500000 -0.5000000
   0.2500000    0.5000000    0.0000000
   0.2500000 -0.2500000 -0.2500000
-0.5000000 -0.5000000 -0.2500000
   0.0000000    0.2500000    0.5000000
-0.2500000    0.2500000 -0.2500000
-0.2500000 -0.5000000 -0.5000000
   0.5000000    0.0000000    0.2500000
-0.5000000 -0.2500000    0.0000000
   0.0000000 -0.5000000 -0.2500000
-0.2500000    0.0000000 -0.5000000
   0.2500000    0.2500000 -0.2500000
   0.5000000    0.2500000    0.5000000
-0.2500000 -0.5000000    0.0000000
-0.2500000    0.2500000    0.2500000
   0.5000000    0.5000000    0.2500000
   0.0000000 -0.2500000 -0.5000000
   0.2500000 -0.2500000    0.2500000
   0.2500000    0.5000000    0.5000000
-0.5000000    0.0000000 -0.2500000
   0.0000000 -0.2500000    0.2500000
   0.2500000    0.0000000 -0.2500000
-0.2500000 -0.2500000 -0.5000000
   0.2500000    0.5000000    0.2500000
   0.2500000 -0.2500000    0.0000000
-0.5000000 -0.2500000 -0.2500000
   0.2500000    0.2500000    0.5000000
   0.0000000    0.2500000 -0.2500000
-0.2500000 -0.5000000 -0.2500000
   0.5000000    0.2500000    0.2500000
-0.2500000    0.0000000    0.2500000
   0.0000000    0.5000000    0.5000000
   0.5000000    0.0000000    0.5000000
   0.2500000 -0.2500000    0.5000000
   0.5000000    0.2500000 -0.2500000
-0.5000000 -0.2500000 -0.7500000
   0.2500000 -0.5000000 -0.2500000
-0.2500000    0.2500000 -0.5000000
end kpoints
运行 wannier90.x wannier90
5.
num_wann=8
num_bands=8

# for GW uncomment
exclude_bands 9-64

Begin Projections
Si:sp3
End Projections

dis_froz_max=9
dis_num_iter=1000

guiding_centres=true

# Bandstructure plot
restart       =  plot
bands_plot    =  true
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.00000 0.5000 K 0.37500 -0.37500 0.0000
K 0.37500 -0.37500 0.0000 G 0.00000  0.00000 0.0000
end kpoint_path
bands_num_points 40
bands_plot_format gnuplot xmgrace

begin unit_cell_cart
   2.7150000    2.7150000    0.0000000
   0.0000000    2.7150000    2.7150000
   2.7150000    0.0000000    2.7150000
end unit_cell_cart

begin atoms_cart
Si    0.0000000    0.0000000    0.0000000
Si    1.3575000    1.3575000    1.3575000
end atoms_cart

mp_grid =    4    4    4

begin kpoints
   0.0000000    0.0000000    0.0000000
   0.2500000    0.0000000    0.0000000
   0.5000000    0.0000000    0.0000000
   0.2500000    0.2500000    0.0000000
   0.5000000    0.2500000    0.0000000
-0.2500000    0.2500000    0.0000000
   0.5000000    0.5000000    0.0000000
-0.2500000    0.5000000    0.2500000
   0.0000000    0.2500000    0.0000000
   0.0000000    0.0000000    0.2500000
-0.2500000 -0.2500000 -0.2500000
-0.2500000    0.0000000    0.0000000
   0.0000000 -0.2500000    0.0000000
   0.0000000    0.0000000 -0.2500000
   0.2500000    0.2500000    0.2500000
   0.0000000    0.5000000    0.0000000
   0.0000000    0.0000000    0.5000000
-0.5000000 -0.5000000 -0.5000000
   0.0000000    0.2500000    0.2500000
   0.2500000    0.0000000    0.2500000
-0.2500000 -0.2500000    0.0000000
-0.2500000    0.0000000 -0.2500000
   0.0000000 -0.2500000 -0.2500000
   0.0000000    0.5000000    0.2500000
   0.2500000    0.0000000    0.5000000
-0.2500000 -0.2500000    0.2500000
-0.5000000 -0.2500000 -0.5000000
   0.2500000    0.5000000    0.0000000
   0.2500000 -0.2500000 -0.2500000
-0.5000000 -0.5000000 -0.2500000
   0.0000000    0.2500000    0.5000000
-0.2500000    0.2500000 -0.2500000
-0.2500000 -0.5000000 -0.5000000
   0.5000000    0.0000000    0.2500000
-0.5000000 -0.2500000    0.0000000
   0.0000000 -0.5000000 -0.2500000
-0.2500000    0.0000000 -0.5000000
   0.2500000    0.2500000 -0.2500000
   0.5000000    0.2500000    0.5000000
-0.2500000 -0.5000000    0.0000000
-0.2500000    0.2500000    0.2500000
   0.5000000    0.5000000    0.2500000
   0.0000000 -0.2500000 -0.5000000
   0.2500000 -0.2500000    0.2500000
   0.2500000    0.5000000    0.5000000
-0.5000000    0.0000000 -0.2500000
   0.0000000 -0.2500000    0.2500000
   0.2500000    0.0000000 -0.2500000
-0.2500000 -0.2500000 -0.5000000
   0.2500000    0.5000000    0.2500000
   0.2500000 -0.2500000    0.0000000
-0.5000000 -0.2500000 -0.2500000
   0.2500000    0.2500000    0.5000000
   0.0000000    0.2500000 -0.2500000
-0.2500000 -0.5000000 -0.2500000
   0.5000000    0.2500000    0.2500000
-0.2500000    0.0000000    0.2500000
   0.0000000    0.5000000    0.5000000
   0.5000000    0.0000000    0.5000000
   0.2500000 -0.2500000    0.5000000
   0.5000000    0.2500000 -0.2500000
-0.5000000 -0.2500000 -0.7500000
   0.2500000 -0.5000000 -0.2500000
-0.2500000    0.2500000 -0.5000000
end kpoints
再次运行 wannier90.x wannier90
能带即为wannier90_band.dat
如果你有gnuplot，可以直接执行wannier90_band.gnu。

GW-band-Si.png

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : BAND_STRUCTURE_GW_VASP.pdf

2013-11-05 13:15:10, 193.1 K

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

学习资料（理工类）	VASP	第一性原理相关文档	VASP
VASP软件从入门到高深	第一性原理计算经验	VASP算例	wannier 90
vasp+wannier90

» 猜你喜欢

OpenAI又攻克数学难题了，离AI做有机合成还有多远？已经有12人回复
极限拉扯，如何抉择已经有13人回复
课题组会帮你通过考核么？已经有7人回复
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍已经有6人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有9人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有8人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有10人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有9人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有7人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

VASP光学性质计算步骤已经有10人回复
vasp5.2如何计算DOS和band structure？已经有10人回复
vasp中利用GW算法算能带结构已经有9人回复
深圳超算中心vasp编译及任务提交问题已经有25人回复
ABINIT和VASP相比，计算GW，哪个更快啊？已经有13人回复
[求助成功]VASP 计算 graphene 的错误已经有15人回复
vasp5.2的编译问题非常感谢！！已经有14人回复
【求助】GW0-Band structure erro 已经有10人回复
【求助】VASP计算光学性质时的GW方法（RPA, no local field effects) 已经有5人回复
【求助】Hse06计算时KPOINTS出错已经有30人回复
【求助】vasp 编译错误已经有5人回复
【其他】关于GW（准粒子计算能带结构等）方法的工作已经有6人回复
【求助】PBE0 计算Si、Ge能带结构已经有9人回复
【其他】免费软件学习资料大征集（abinit, pwscf, lammps, CPMD, XMD 等等免费软件）已经有21人回复

1楼 2013-11-05 13:16:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

1024247900

木虫 (正式写手)

应助: 15 (小学生)
金币: 3000.5
帖子: 314
在线: 202.3小时
虫号: 2325629

★★★★★ 五星级,优秀推荐

楼主V5!

[ 发自小木虫客户端 ]

回复此楼

24楼2014-04-26 23:44:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 64 个回答

简单回复

4901051793楼

2013-11-05 17:17 回复

五星好评顶一下，感谢分享！

漫天飘雪5楼

2013-11-05 18:45 回复

五星好评顶一下，感谢分享！

tiaozaoququ7楼

2013-11-05 19:30 回复

五星好评顶一下，感谢分享！

bearpow8楼

2013-11-06 00:07 回复

五星好评顶一下，感谢分享！

乐石||9楼

2013-11-06 10:27 回复

五星好评顶一下，感谢分享！

tangsw91111楼

2013-11-06 19:11 回复

五星好评顶一下，感谢分享！

tangchao200812楼

2013-11-06 23:36 回复

五星好评顶一下，感谢分享！

Quan.13楼

2013-11-07 05:46 回复

五星好评顶一下，感谢分享！

xmc21214楼

2013-11-07 09:14 回复

五星好评顶一下，感谢分享！

查看全部 64 个回答

☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级

普通表情龙兔虎猫高级回复 (可上传附件)

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8	wangyl_123 2026-10-10	12/600	2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯，如何抉择 +9	Isak0320 2026-10-12	13/650	2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么？ +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一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5	qR2AdCWTVtkw 2026-10-09	8/400	2026-10-11 22:26 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	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文章，我:8O.5.5.1O.54,科目全,可十急 +5	qR2AdCWTVtkw 2026-10-09	8/400	2026-10-11 13:51 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 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4	qR2AdCWTVtkw 2026-10-09	4/200	2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6	zapen 2026-10-08	8/400	2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5	jing816y 2026-10-08	9/450	2026-10-10 09:05 by zapen