| 查看: 4314 | 回复: 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 | 表界面计算 |
» 猜你喜欢
科研孤儿太难了
已经有9人回复
-
时间戳又变了
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
8楼2020-08-21 10:37:39
2楼2016-03-09 21:57:53
3楼2016-04-17 19:56:10
4楼2016-04-17 19:57:12
40