24小时热门版块排行榜    

查看: 1546  |  回复: 8
当前主题已经存档。

tiantian2921

新虫 (小有名气)

[交流] 【求助】分析dos,作图问题

VASP计算完DOS进行分析,分析的是DOSCAR文件,用自编的小程序,但是图像出来非常不平滑,都是尖锐的峰,谁能帮我解决一下吗?

不知道是我编的程序的问题,还是需要对图像平滑,但是平滑后的效果很差,图像不好看
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1,VIP+0):谢谢提示 11-20 20:28
自己搜索以下本版面很多相关的数据处理程序和辅助软件
2楼2009-11-20 20:22:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyong.gan

木虫 (正式写手)

★
aylayl08(金币+1,VIP+0):谢谢提示 11-21 20:51
引用回帖:
Originally posted by tiantian2921 at 2009-11-20 17:28:
VASP计算完DOS进行分析,分析的是DOSCAR文件,用自编的小程序,但是图像出来非常不平滑,都是尖锐的峰,谁能帮我解决一下吗?

不知道是我编的程序的问题,还是需要对图像平滑,但是平滑后的效果很差,图像不好看

可以直接复制到origin里面,得到的图也会有很多毛刺,不过可以在origin里面对图进行平滑处理
3楼2009-11-21 14:44:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1,VIP+0):谢谢提示 11-21 20:51
如果你用的ISMEAR=0的话,可以试着增大一下sigma的值试试
4楼2009-11-21 18:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):讨论奖励! 11-25 11:42
我也遇这样的问题,还没有解决,我试着在origin里面插值进行光滑处理,但是把态密度展到带隙中去了,本来是绝缘体结果成了半金属,我本来是采用ISMEAR=-5,后来试着用=0,增加SIGMA,但是图跟开始相差很大,不知道是什么原因?
5楼2009-11-21 21:54:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian2921

新虫 (小有名气)

引用回帖:
Originally posted by 含草 at 2009-11-21 21:54:
我也遇这样的问题,还没有解决,我试着在origin里面插值进行光滑处理,但是把态密度展到带隙中去了,本来是绝缘体结果成了半金属,我本来是采用ISMEAR=-5,后来试着用=0,增加SIGMA,但是图跟开始相差很大,不知道 ...

以下是VASP手册的原话
6楼2009-11-23 16:43:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian2921

新虫 (小有名气)

★ ★
zxzj05(金币+2,VIP+0):讨论奖励! 11-25 11:42
引用回帖:
Originally posted by 含草 at 2009-11-21 21:54:
我也遇这样的问题,还没有解决,我试着在origin里面插值进行光滑处理,但是把态密度展到带隙中去了,本来是绝缘体结果成了半金属,我本来是采用ISMEAR=-5,后来试着用=0,增加SIGMA,但是图跟开始相差很大,不知道 ...

对不起,误操作,写了一半,就发出去了,重新发,很抱歉

以下是VASP手册的原话,即ISMEAR取值的建议中的两点

1.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.
2.For the calculations of the DOS and very accurate total energy calculations (no relaxation in metals) use the tetrahedron method (ISMEAR=-5)

我理解是这样的,对于计算DOS无论对于金属,半导体或者绝缘体,ISMEAR=-5
而不能取别的值,不然会导致DOS的计算的不正确。不知道我的理解是否准确。

建议你也可以试试,算你的那个体系,ISMEAR=0把SIGMA取小到0.05,看看结果如何。Sample Text
7楼2009-11-23 16:51:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)

★
zxzj05(金币+1,VIP+0):讨论奖励! 11-25 11:43
可以借鉴一下MS中dos图的画法
8楼2009-11-23 23:30:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjhwoods

铜虫 (正式写手)

学习学习了
9楼2010-02-05 05:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tiantian2921 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:25 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 10:52 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 TZqyNTBa07NP 2026-11-13 4/200 2026-11-15 10:08 by xSrQNSSjcgoX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:46 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:29 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:45 by xSrQNSSjcgoX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 07:01 by xSrQNSSjcgoX
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:20 by Equinoxhua
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:16 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
信息提示
请填处理意见