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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4337)
>基金申请 (773)
>休闲灌水 (686)
>招聘信息布告栏 (607)
>导师招生 (282)
>文献求助 (231)
>虫友互识 (211)
>论文投稿 (155)
>有奖起名 (79)
>考博 (74)
>硕博家园 (64)
>博后之家 (55)
>学术会议 (51)
>教师之家 (35)
>能源 (30)
>找工作 (29)

返回列表

【奖励】本帖被评价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

» 猜你喜欢

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

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的回帖

scuedu2

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 818.3
帖子: 179
在线: 280.5小时
虫号: 2382371

谢谢

发自小木虫Android客户端

回复此楼

46楼2016-01-25 13:19:23

已阅回复此楼关注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 个回答

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

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

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

[文学芳草园] 梦想 +5	myrtle 2026-08-26	7/350	2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5	ssyjh 2026-08-27	8/400	2026-09-01 12:58 by ssyjh
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13	Kittylucky 2026-08-27	14/700	2026-09-01 11:06 by feng6531
[论文投稿] 小白求助投论文要求的highlights应该如何写 5+3	l1963982152 2026-08-29	4/200	2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7	gdd2018 2026-08-28	12/600	2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7	雨打竹帘 2026-08-30	11/550	2026-08-31 23:16 by hittle2008
[基金申请] 投票: 有多少人是今天查系统知道结果的？ +17	爱看书的可乐 2026-08-26	19/950	2026-08-31 21:30 by xiangy672
[基金申请] 能否申诉？ +7	echo8914667 2026-08-30	8/400	2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7	349506619 2026-08-28	7/350	2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中，共勉 +12	eulota 2026-08-26	12/600	2026-08-31 08:42 by ZJTJZ
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6	gouxfjh 2026-08-28	11/550	2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13	wangmingqi 2026-08-28	19/950	2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11	董八千 2026-08-26	11/550	2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4	winsaint 2026-08-27	9/450	2026-08-28 11:06 by maolC
[基金申请] 出来了 +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
[基金申请] 国合里面能看到了 +7	一怀馨秋 2026-08-26	7/350	2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6	wittyboy 2026-08-26	14/700	2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr