24小时热门版块排行榜    

查看: 1221  |  回复: 1
本帖产生 1 个 1ST强帖 ,点击这里进行查看

lascqy

银虫 (小有名气)

[求助] vasp5.2 计算光学性质时的取点密度问题

用vasp5.2,NEDOS=2000,得到的介电常数
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.190277    0.043709    0.043709    0.014921    0.000000    0.000000    0.000000
    0.380555    0.106859    0.106859    0.032319    0.000000    0.000000    0.000000
    0.570832    0.242305    0.242305    0.056362    0.000000    0.000000    0.000000
    0.761110    0.748739    0.748739    0.097848    0.000000    0.000000    0.000000
    0.951387    2.967057    2.967057    0.202469    0.000000    0.000000    0.000000
画出来的图形很不理想
而在论坛上看见的人家的
frequency dependent IMAGINARY DIELECTRIC FUNCTION (RPA, no local field effects)
     E(ev)      X         Y         Z        XY        YZ        ZX
   --------------------------------------------------------------------------------------
     0.000000    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
     0.092884    0.001782    0.001782    0.001782    0.000000    0.000000    0.000000
     0.185768    0.003575    0.003575    0.003575    0.000000    0.000000    0.000000
     0.278652    0.005388    0.005388    0.005388    0.000000    0.000000    0.000000
明显比我的密集很多,求指教
回复此楼

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

纳米方向

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

后天一

木虫 (小有名气)

我们的自由!

【答案】应助回帖

感谢参与,应助指数 +1
uuv2010(金币+3, 1ST强帖+1): 多谢指导! 2012-02-22 22:17:43
lascqy(金币+10): ★★★★★最佳答案 今天刚看到~~正解啊~~~帮大忙了 2012-02-24 11:27:15
介电函数的取点和两个参数有关,NEDOS和ISMEAR
介电函数在OUTCAR中输出循环范围为I=1,NEDOS,ISTEP
NEDOS控制最大点数,vasp在计算介电函数的时候至少会取1000个点,如果你的NEDOS大于1000,介电函数的取点就等于NEDOS,也就是:
NEDOS=MAX(NEDOS_IN,1000)
而步数(省略数据)和NEDOS和ISMEAR有关,如果你的NEDOS为2000,而ISMEAR>=0的话,vasp会省略一部分数据,一般会除10,那你的取点应该是200个点.
具体代码如下
      ISTEP=MIN(10,NEDOS/40)
      IF (KPOINTS%ISMEAR<0) ISTEP = 1

另外,完整的数据在vasprun.xml文件里面,这里面的数据点不会省略,一共NEDOS个点。你可以提取相关信息。

希望对你有帮助,来个1ST最好了,哈哈
2楼2012-02-22 21:22:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lascqy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见