| 查看: 2803 | 回复: 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计算的结果会有多大影响？ |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关遗传算法的问题求助
已经有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人回复
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 ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
