| 查看: 4328 | 回复: 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 | 表界面计算 |
» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有8人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
6楼2019-06-23 15:31:56
2楼2016-03-09 21:57:53
3楼2016-04-17 19:56:10
4楼2016-04-17 19:57:12
