24小时热门版块排行榜    

查看: 2074  |  回复: 6

nkleof

银虫 (正式写手)

[交流] 【求助】如何选取合适的ISMEAR 已有4人参与

vasp说明文档里提到ISMEAR的取值有几种,常用的有-5,-1,0,N,分别对应tetrahedron method with Blöchl corrections,Fermi-smearing,Gaussian smearing,method of Methfessel-Paxton order N,在不同的体系里应该采取不同的ISMEAR值,比如Avoid to use  ISMEAR>0 for semiconductors and insulators, since it might cause problems.
那么对于一个未知的体系,如何设置ISMEAR值比较好?需要取不同的ISMEAR值都算一下吗?比如计算掺杂或者磁序改变对体系DOS的影响,事先不可能知道体系是金属性还是半导体性或者绝缘体性。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

nkleof(金币+1):谢谢讨论 2010-04-26 11:32
一般金属取1, 非金属取0。
很简单,最快的方法就是用ms算一下,精度取中就行。
真理是一点点接近的!
2楼2010-04-26 11:02:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by gavinliu7390 at 2010-04-26 11:02:10:
一般金属取1, 非金属取0。
很简单,最快的方法就是用ms算一下,精度取中就行。

那么对于未知体系是要用各种ISMEAR值都算一遍吗?
另外非金属情况取-5更好吧?
3楼2010-04-26 11:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


gavinliu7390(金币+1):谢谢交流! 2010-04-26 15:12
nkleof(金币+1):谢谢讨论! 2010-04-26 16:36
先用几个离子步(NSW=10)来跑一下看看.对于不熟悉的体系,我先用ISMEAR=-5跑十步,看看输出的提示和消耗情况,然后来调整参数.自然ISMEAR,SIGMA,POTIM都是调节的对象.SIGMA=0(ncg~=900)比-5(ncg~=4000)跑的快的多.
但是我不明白为什么.
4楼2010-04-26 12:45:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by cenwanglai at 2010-04-26 12:45:36:
先用几个离子步(NSW=10)来跑一下看看.对于不熟悉的体系,我先用ISMEAR=-5跑十步,看看输出的提示和消耗情况,然后来调整参数.自然ISMEAR,SIGMA,POTIM都是调节的对象.SIGMA=0(ncg~=900)比-5(ncg~=4000)跑的快的多.
但 ...

谢谢,看来只能是尝试了。以收敛速度作为ISMEAR值合不合适的评价标准吗?
5楼2010-04-26 16:39:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengxiangm

木虫 (小有名气)

nkleof(金币+1):谢谢讨论! 2010-04-26 17:31
一般用ISMEAR=0,SIGMA=0.05都能得到一个比较合理的值
6楼2010-04-26 17:01:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by zengxiangm at 2010-04-26 17:01:54:
一般用ISMEAR=0,SIGMA=0.05都能得到一个比较合理的值

vasp说明文档里有这样一句  For semiconductors or insulators use the tetrahedron method ( ISMEAR=-5), if the cell is too large (or if you use only a single or two k-points) use  ISMEAR=0 in combination with a small  SIGMA=0.05.
请再详细一点讨论一下“一个比较合理的值”,这个值是指所得到的基态性质吗?也就是说一般非金属都用ISMEAR=0,SIGMA=0.05比较好吧?
7楼2010-04-26 17:38:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nkleof 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 v7Aev8ORQjin 2026-08-22 3/150 2026-08-23 20:58 by 0p9y1isJp2jM
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:22 by lXgU0RvNTDC9
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 19:10 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见