| 查看: 3908 | 回复: 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 | 表界面计算 |
» 猜你喜欢
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
2楼2016-03-09 21:57:53
6楼2019-06-23 15:31:56
7楼2019-06-23 15:33:18