| 查看: 7082 | 回复: 34 | ||||||||||
zzy870720z荣誉版主 (文坛精英)
|
[交流]
【求助】VASP计算电子态密度已有14人参与
|
|||||||||
|
最近做点测试计算 可是DOS做出来效果很差 如下图  同时附上自己的INCAR 希望高手帮忙，如何做出漂亮的DOS图 并指出那些参数对DOS图有影响 先谢谢了 ENCUT = 300 eV ISTART = 0 ISPIN = 2 ISYM = 0 GGA = 91 VOSKOWN = 1 ISIF = 0 NELM = 100 ALGO = FAST EDIFF = 1E-4 EDIFFG = 0.001 LREAL = .FALSE. NSW = 0 IBRION = -1 ISMEAR = -5 SIGMA = 0.1 NWRITE = 2 NKPR = 1 |
回复此楼» 收录本帖的淘贴专辑推荐
 VASP | 第一性原理 | vasp 学习贴 | 模拟，理论 |
| VASP | 待办事宜 | VASP报错建模问题集 |
» 猜你喜欢
34岁读博士晚吗
已经有22人回复
-
江西双非一本和四川双一流高校如何选择？
已经有19人回复
-
杭电、天津科技、青农和宁波工程学院如何选？
已经有7人回复
-
博士后基金需要结题吗？
已经有8人回复
-
工材E口JQ有消息了吗
已经有4人回复
-
投稿时忘记修改一作
已经有8人回复
-
硕博巨婴，也许才刚刚开始
已经有56人回复
-
食品博士导师
已经有8人回复
-
竟然失眠了
已经有9人回复
-
饶议：什么制度能保障大学普通教师不用为领导拎包，不用看领导脸色
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP静态计算中ismear的选择
已经有15人回复
- vasp dos曲线平滑
已经有3人回复
- siesta计算分波态密度及态密度，能带
已经有6人回复
- VASP计算态密度时K点选择
已经有6人回复
- VASP结果分析求助
已经有6人回复
- vasp静态计算，非静态计算哪个用来做态密度，哪个用来做能带？
已经有7人回复
- vasp计算含有磷原子结构的态密度时磷原子为什么还有d态
已经有14人回复
- vasp计算，DOS图的分析，分波态密度的数据分析~·
已经有4人回复
- VASP 计算DOS，态密度中费米能的作用
已经有3人回复
- vasp测试SIGMA
已经有12人回复
- 用VASP计算Mg的PDOS的时候出现d电子是什么原因？
已经有6人回复
- 怎么根据态密度图积分算区间的电子数
已经有8人回复
- Guassian和VASP计算过渡态的原理和算法上差别在哪里？
已经有9人回复
- VASP: 计算态密度时，设置哪些参数可以使态密度的曲线更平滑，谢谢
已经有13人回复
- vasp计算在进行静态计算时，出错。
已经有3人回复
- vasp计算结果处理态密度要用到的程序
已经有18人回复
- VASP计算态密度出现问题
已经有12人回复
- 【求助】VASP计算的KPOINTS 设置，求助
已经有10人回复
- 【求助】用vasp计算电子态密度出现错误
已经有10人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【求助】VASP怎样精确计算出费米能级处的态密度值？
已经有22人回复
- 【求助】关于vasp态密度计算
已经有10人回复
- 【求助】电子态密度计算dos图中feimi能级的确定
已经有7人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25647.7
- 散金: 602
- 红花: 166
- 帖子: 3771
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★
zzy870720z(金币+2):你的意思是采用gaussian smearing？但是vasp手册上建议态密度用四面体方法取-5？ 2010-06-03 07:54:04
hometownlove(金币+2):谢谢交流 2010-06-09 13:40:08
zzy870720z(金币+2):你的意思是采用gaussian smearing？但是vasp手册上建议态密度用四面体方法取-5？ 2010-06-03 07:54:04
hometownlove(金币+2):谢谢交流 2010-06-09 13:40:08
|
1)增加k点网格的大小； 或者 2）采用ISMEAR=0, SIGMA = 0.1~0.2。 |
» 本帖已获得的红花（最新10朵）
xiaobo13142楼2010-06-03 07:35:35
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25647.7
- 散金: 602
- 红花: 166
- 帖子: 3771
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
gavinliu7390(金币+3):非常感谢详细指导！ 2010-06-03 14:57:20
zzy870720z(金币+5):谢谢侯老师指点，我先试试，有问题再请教 2010-06-03 16:44:47
hometownlove(金币+1):谢谢交流 2010-06-09 13:44:57
gavinliu7390(金币+3):非常感谢详细指导！ 2010-06-03 14:57:20
zzy870720z(金币+5):谢谢侯老师指点，我先试试，有问题再请教 2010-06-03 16:44:47
hometownlove(金币+1):谢谢交流 2010-06-09 13:44:57
|
下面的资料是从The FHIMD Toolkit User's Manual上拷贝下来的： http://www.fhi-berlin.mpg.de/th/fhi98md/  见里面的公式3.2 1、先用vasp计算态密度，比如采用ISMEAR=-5，同时将NEDOS设置的大一些，默认值为301，将它调大，比如设置为1000左右； 2、得到了DOSCAR文件之后，取出total DOS. 可以简单的认为，在采用ISMEAR=-5时计算态密度，上面公式3.1中的F就是delta函数。 3、按上面公式3.1， 3.2，对计算出来的态密度做一个展宽，这个要自己编程： 公式3.1中的权重都为1.0了，不对k进行循环，只是对能量点的个数， 公式3.1里面|...|内含的量已经由vasp计算出来了的。 [ Last edited by valenhou001 on 2010-6-3 at 13:07 ] |
13楼2010-06-03 13:06:09
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7439.4
- 散金: 9654
- 红花: 117
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
★
zzy870720z(金币+2):感谢交流 2010-06-03 16:45:06
hometownlove(金币+1):谢谢交流 2010-06-09 13:45:52
zzy870720z(金币+2):感谢交流 2010-06-03 16:45:06
hometownlove(金币+1):谢谢交流 2010-06-09 13:45:52
|
展宽不就是SIGMA吗?在INCAR中不是已经给出来了?为什么会有""在处理数据的脚本中来约定"这一说呢? 搜索到下面这一段话:"当体系的电子结构在Fermi 能级附近出现多个简并能级时，体系的自洽过程会相当慢，甚至不能收敛。这是由于能级简并，或能级之间相差不大，导致电子“不知道”应该在哪里分布，从而不断地在几个能级之间“跳来跳去”，但无论电子分布在哪个前线能级上整个体系都不能自洽。 所谓拖尾，是指在前线轨道上，电子可以以分数分布，而不再要求一定要以整数分布，从而使电子可以一分为多地分布在能量相近的各个前线轨道，而不需要在各轨道之间“跳来跳去”，使体系更容易自洽。使用该办法时要对拖尾效应设定一个能量值，使电子在能量相差在该值之内的轨道进行分数分布。设定值越小，计算结果越精确，但亦更难自洽。（也叫展宽）" |
14楼2010-06-03 15:08:12
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1446.3
- 散金: 134
- 红花: 24
- 帖子: 685
- 在线: 316.1小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2010-06-03 08:06:08
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25647.7
- 散金: 602
- 红花: 166
- 帖子: 3771
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
6楼2010-06-03 08:40:00
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25647.7
- 散金: 602
- 红花: 166
- 帖子: 3771
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
10楼2010-06-03 10:32:04
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
4楼2010-06-03 08:25:18
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
5楼2010-06-03 08:26:07
7楼2010-06-03 08:51:56
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
8楼2010-06-03 10:02:56
xiaojie7783
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 957.4
- 散金: 76
- 帖子: 234
- 在线: 32.6小时
- 虫号: 702761
- 注册: 2009-02-17
- 性别: GG
- 专业: 金属材料的力学行为
9楼2010-06-03 10:19:07
5