| 查看: 2830 | 回复: 16 | |||
mazuju028木虫 (正式写手)
|
[求助]
GW计算流程问题
|
|
GW计算在vasp网站上有简单说明，还是有些不明确的地方想和大家交流一下。 先简单罗列一下步骤， 1，较少空带的自洽计算。 System = Si groundstate occupied orbitals ISMEAR = 0 ; SIGMA = 0.05 ! small sigma is required to avoid partial occupancies EDIFF = 1E-8 ! required tight tolerance for groundstate orbitals 2，较多能带的自洽计算。 System = Si unoccupied orbitals ALGO = Exact ! use exact diagonalization of the Hamiltonian NELM = 1 ! since we are already converged stop after one step NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 ! small sigma is required to avoid partial occupancies LOPTICS = .TRUE. 3，G0W0计算。 System = Si NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 LOPTICS = .TRUE. ALGO = GW0 ; NOMEGA = 50 ---疑问如下--------------------------------------------------------------------------- 1，这三步都应该是在同一目录下进行，所以第二步应该是默认读取第一步的波函数，第三步默认读取第二步的波函数？ 2，有关K点设置，是否可以从第一步自洽计算开始就用特殊K点？ 3，第二步中NELM=1如果不设置，对最后GW计算的结果会有多大影响？ |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关遗传算法的问题求助
已经有3人回复
- GW计算的错误提示【问题已解决】
已经有4人回复
- 关于冷却问题的计算
已经有10人回复
- abinit计算晶体生成焓的问题讨论
已经有7人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 求助签证问题以及流程注意事项等
已经有9人回复
- 【求助】aspen流程图中选择单元操作的一个小问题
已经有3人回复
- 【求助】VASP计算光学性质时的GW方法（RPA, no local field effects)
已经有5人回复
- 【求助】GW计算中k点的选择问题
已经有16人回复
- 【求助】GW精度问题
已经有13人回复
ApolloYang
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2209.6
- 散金: 114
- 红花: 4
- 帖子: 390
- 在线: 236.9小时
- 虫号: 582657
- 注册: 2008-07-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-25 15:29:34
mazuju028: 金币+4, ★★★很有帮助, 谢谢。 2012-05-25 17:03:55
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-25 15:29:34
mazuju028: 金币+4, ★★★很有帮助, 谢谢。 2012-05-25 17:03:55
|
NELM = 1 ! since we are already converged stop after one step 看注释，是因为step1进行了自洽运算，故这里设置NELM=1是为了减少计算量 2，有关K点设置，是否可以从第一步自洽计算开始就用特殊K点？ 楼主是为了算能带吗？ 我曾经试过，特殊的k点一开始可以使用的，但是最后GW的计算仍然得不到能带图。 |
2楼2012-05-25 10:06:12
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
3楼2012-05-25 15:57:26
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2770.4
- 散金: 528
- 红花: 10
- 帖子: 943
- 在线: 419.5小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-05-25 17:03:36
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2770.4
- 散金: 528
- 红花: 10
- 帖子: 943
- 在线: 419.5小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2012-05-25 17:04:40
ApolloYang
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2209.6
- 散金: 114
- 红花: 4
- 帖子: 390
- 在线: 236.9小时
- 虫号: 582657
- 注册: 2008-07-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2012-05-26 08:04:27
ApolloYang
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2209.6
- 散金: 114
- 红花: 4
- 帖子: 390
- 在线: 236.9小时
- 虫号: 582657
- 注册: 2008-07-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2012-05-26 08:05:51
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10506.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
8楼2012-05-26 11:55:50
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
9楼2012-05-26 16:16:13
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10506.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
10楼2012-05-26 16:40:33
