| 查看: 2414 | 回复: 7 | ||||
corbin银虫 (小有名气)
|
[交流]
VASP的GW计算很耗时间吗？已有4人参与
|
| VASP的GW计算很耗时间吗？我提交了一个任务，ZnO 4个原子的原胞，居然可以用到22G内存，几十个小时了，还没算完，请问有人用VASP做过GW计算吗？ |
回复此楼» 收录本帖的淘贴专辑推荐
 收藏 |
» 猜你喜欢
今天能不能出来名单
已经有10人回复
-
博后网上办公系统维护
已经有8人回复
-
博后面上今天下午会公布吗？有无消息？
已经有38人回复
-
对象没有，还非常想读博，难以抉择
已经有47人回复
-
联培博士文章第一单位署名问题交流~
已经有19人回复
-
过柱子，产品，杂质在是 pe:ea=100:1 也一起出来？
已经有5人回复
-
酚羟基甲基化
已经有5人回复
-
刚刚收到科研之友邮件
已经有43人回复
-
【45岁以上博士】柔性人才引进项目，有补贴
已经有17人回复
-
怎么带研究生？
已经有24人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算过渡态问题
已经有13人回复
- vasp计算能带，刚开始就服务器就卡了
已经有20人回复
- 我们用VASP计算材料的力学性能，是绝对零度的。有意义吗？正解者有红花相送哦
已经有8人回复
- vasp计算的半导体结构，bands带隙比dos里面的大，何故？
已经有33人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- vasp 中能用 GW算能带吗？
已经有7人回复
- 【求助】VASP中的GW方法能计算能带吗？
已经有5人回复
- 【求助】VASP计算光学性质时的GW方法（RPA, no local field effects)
已经有5人回复
- 【求助】band structures using hybrid functionals 的Kpoints问题
已经有18人回复
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有51人回复
- 【求助】vasp5.2的GW怎么算？
已经有14人回复
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3717.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
2楼2011-05-27 20:28:09
tider
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1950.4
- 散金: 40
- 红花: 5
- 帖子: 313
- 在线: 87.2小时
- 虫号: 1190614
- 注册: 2011-01-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-05-28 00:06:27
corbin
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 457.1
- 散金: 4
- 帖子: 113
- 在线: 93.5小时
- 虫号: 1045382
- 注册: 2010-06-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2011-05-28 19:13:52
liuhuimeiwoo
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1469.1
- 散金: 10
- 红花: 3
- 帖子: 235
- 在线: 40.2小时
- 虫号: 911880
- 注册: 2009-11-24
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2011-05-31 10:35:26
corbin
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 457.1
- 散金: 4
- 帖子: 113
- 在线: 93.5小时
- 虫号: 1045382
- 注册: 2010-06-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我发现改了KPOINT 从8x8x8变为6x6x6，而且将GW计算的参数ENCUTGW改为150eV之后，计算速度还是可以接受的，但是现在不太清除到底是KPOINT对计算的影响大，还是ENCUTGW对计算的影响大 另外，对于GW的计算步骤还是有点不理解 先做DFT计算，vasp说明书上给的Si的例子 DFT-INCAR System = Si NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 ! small sigma is required to avoid partial occupancies LOPTICS = .TRUE. 然后再做GW0计算 GW0-INCAR System = Si NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 LOPTICS = .TRUE. ALGO = GW0 ; NOMEGA = 50 NELM = 4 说明书上的意思是，由于GW计算需要DFT计算得到的WAVECAR，并且在计算过程中需要文件WAVEDER 我的问题是，从第二步GW0计算的INCAR设置来看，由于目录中存在WAVECAR，所以ISTART=1，那么ICHARG=0(?)，也就是说GW0计算中不需要读取DFT计算的CHARGCAR？ 问这个问题的原因是在网上看到做GW0计算的Bandstructure的INCAR设置， 在这份设置中，ISTART=1，ICHARG=11，然后讲KPOINT改为倒空间中的路径，相对于DFT计算的K点个数，大大减少，这样不影响GW计算结果吗？ |
6楼2011-05-31 15:31:14
|
本帖内容被屏蔽 |
7楼2011-05-31 15:40:51
corbin
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 457.1
- 散金: 4
- 帖子: 113
- 在线: 93.5小时
- 虫号: 1045382
- 注册: 2010-06-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我在网上找了Si的例子，如下： 首先是自述文件： ======================================================================================= Here is two VASP-GW examples. - BN : Graphene-like BoronNitride 2d structure. - Bulk silicon. In each folder, input files are available (INCAR1,INCAR2,POSCAR,POTCAR,KPOINTS) and the procedure is: <1> cp INCAR1 INCAR <2> $vasp <3> mkdir DFT <4> cp * DFT <5> cp INCAR2 INCAR <6> $vasp INCAR1 is for initial DFT calculation and INCAR2 is for GW0 part. Lets examine the results in BN folder. - KPOINTS should be "Automatic mesh". You can examine the kpoints determined by vasp (BZ.png,BZ.eps,draw_bz.m) - After step <2>, WAVEDER, WAVECAR files are generated by INCAR1 and they are used in step <6> by INCAR2. - GW0 OUTCAR contains lots of information about GW... - In DFT folder, there is a WEIRD band structure (bands.eps) which is plotted along "Automatic mesh" kpoints. - In DFT_band folder, there is a NICE band structure (bands.eps) which is plotted with "Line-mode". You can examine also the kpoints determined by vasp (BZ.png,BZ.eps). - GW_band folder is the tricky part. Here I select the kpoints from "Automatic mesh" along G-K-M direction. You can examine .eps and .m files there. As in GW_bands.eps, band gap increases from 4.58 to 6.23 with GW0. I do not know whether this procedure is CORRECT or there exists easier methods to ger GW band structure along special paths in BZ. The procedure is exactly the same for Silicon folder. However, GW_band part is more difficult since we need to deal with 3-dimensional Kpoints mesh. So I just plot the e-values along "Automatic mesh". The band gap is ill defined for this case. However we can calculate the shift of QP energies (scissor operation ??). The "band gap" in DFT folder is 0.58 eV and this increases to 1.28 eV as in GW_band folder. This increase is close to the value in PHYSICAL REVIEW B 75, 235102 (2007). ======================================================================================= INCAR1 ======================================================================================= System = GW ISMEAR = 0 SIGMA = 0.05 LOPTICS = .TRUE. ALGO = Fast NBANDS = 160 ======================================================================================= INCAR2 ======================================================================================= System = GW ISMEAR = 0 SIGMA = 0.05 LOPTICS = .TRUE. ALGO = GW0 NBANDS = 160 NOMEGA = 64 NELM = 4 ======================================================================================= POSCAR ======================================================================================= scaled 0.5306060285750001 -5.0992500000000000 0.0000000000000000 5.0992500000000000 0.0000000000000000 5.0992500000000000 5.0992500000000000 -5.0992500000000000 5.0992500000000000 0.0000000000000000 Si 2 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 ======================================================================================= KPOINTS ======================================================================================= Automatic mesh 0 G 6 6 6 ======================================================================================= POTCAR就不贴了，用的是PAW 按照自述文件的步骤，可以成功完成GW0计算 |
8楼2011-05-31 16:58:03
5