| 查看: 4325 | 回复: 7 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
GW+BSE计算介电常数 已有1人参与
|
||||
|
问题是：ALGO=BSE 计算很快结束，vasprun.xml 中没有介电常数之类的输出。 输入文件： step1:标准的优化结构； step2  ptic计算INCAR2 system = fcc Si ISMEAR = 0; SIGMA = 0.05; POTIM=0.1 LWAVE=T LCHARG=T NSW=0 ISIF=2 ENCUT=500 ALGO=exact NBANDS=120 LOPTICS=.TRUE. NEDOS=2000 EDIFF=1e-5 EDIFFG=-0.01 step3: GW INCAR3 system = fcc Si ISMEAR = 0; SIGMA = 0.05; POTIM=0.1 LWAVE=T LCHARG=T NSW=0 ISIF=2 IBRION=2 ENCUT=500 NBANDS=120 LOPTICS=.TRUE. NEDOS=2000 EDIFF=1e-5 EDIFFG=-0.01 ALGO= GW0 ; NOMEGA = 72 NELM=1 LSPECTRAL= .TRUE. STEP4: BSE method system = fcc Si ISMEAR = 0; SIGMA = 0.05; POTIM=0.1 LWAVE=T LCHARG=T NSW=0 ISIF=2 IBRION=2 ENCUT=500 NBANDS=120 LOPTICS=.TRUE. NEDOS=2000 EDIFF=1e-5 EDIFFG=-0.01 ALGO= BSE ; NOMEGA = 72 LSPECTRAL= .TRUE. #NPAR= 20 #NELM= 5 NBANDSO= 12 NBANDSV= 12 #OMEGAMAX= 6 其实我想算的是想石墨烯一样的二维体系，但是一样的错误。OUTCAR的最后是： WARNING: aliasing errors must be expected set NGX to 26 to avoid them WARNING: aliasing errors must be expected set NGY to 26 to avoid them WARNING: aliasing errors must be expected set NGZ to 26 to avoid them exchange correlation table for LEXCH = 8 RHO(1)= 0.500 N(1) = 2000 RHO(2)= 100.500 N(2) = 4000 NQ= 1 0.0000 0.0000 0.0000, parallelization over k-points Bands included in local field effects VB(min)= 6 VB(max)= 4 CB(min)= 5 CB(max)= 3 parameters for HF 0.0000000 0.0000000 求解： 是不是只有GW+BSE这样的方法才可以得到比较精确的介电常数？ 直接用DFT，加Optic变量得到的数值与实验有较大差距。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | 表界面计算 |
» 猜你喜欢
放榜前的不淡定
已经有13人回复
-
今天放榜没戏了吧
已经有6人回复
-
让我中一个面上吧！
已经有11人回复
-
时间戳今天，20号变了
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有24人回复
-
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
3楼2016-04-17 19:56:10
2楼2016-03-09 21:57:53
4楼2016-04-17 19:57:12
jx在路上
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 3373.7
- 散金: 932
- 红花: 5
- 帖子: 1321
- 在线: 823.4小时
- 虫号: 1710564
- 注册: 2012-03-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2016-07-24 19:07:57
40
