24小时热门版块排行榜    

查看: 3895  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qingeraaa

银虫 (小有名气)

[求助] GW+BSE计算介电常数 已有1人参与

问题是:ALGO=BSE 计算很快结束,vasprun.xml 中没有介电常数之类的输出。
输入文件:
step1:标准的优化结构;
step2ptic计算
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变量得到的数值与实验有较大差距。
回复此楼
对于未来,我一点也不不担心,因为我相信,时光会把我变得更好!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scuedu2

金虫 (小有名气)

3楼2016-04-17 19:56:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

qingeraaa

银虫 (小有名气)

出现的错误是:MKL ERROR : Parameter 13 was incorrect on entry to ZGEMM
求解啊
对于未来,我一点也不不担心,因为我相信,时光会把我变得更好!
2楼2016-03-09 21:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scuedu2

金虫 (小有名气)

请问一下介电常数各项异性怎么处理呢,就是三个方向差很多,不能取平均

发自小木虫Android客户端
4楼2016-04-17 19:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jx在路上

金虫 (著名写手)

【答案】应助回帖

步骤确定正确?第二步难道不是要自洽?
5楼2016-07-24 19:07:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见