24小时热门版块排行榜    

查看: 4028  |  回复: 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变量得到的数值与实验有较大差距。
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 表界面计算

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

对于未来,我一点也不不担心,因为我相信,时光会把我变得更好!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingeraaa

银虫 (小有名气)

用GW主要是带隙校正,这样算出来的能带才比较准,得出来的光学数据才可信

发自小木虫Android客户端
对于未来,我一点也不不担心,因为我相信,时光会把我变得更好!
6楼2019-06-23 15:31:56
已阅   回复此楼   关注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

金虫 (小有名气)

3楼2016-04-17 19:56:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scuedu2

金虫 (小有名气)

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

发自小木虫Android客户端
4楼2016-04-17 19:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085600材料与化工 +5 安全上岸! 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
信息提示
请填处理意见