应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助pdos的问题
111/2返回列表12下一页
查看: 2252  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小豆sumy

金虫 (小有名气)

[求助] 求助pdos的问题

我计算态密度的结果与文献不相符合，首先是总态密度，文献在4eV的地方就收敛了，而我的在4eV后还有图形。
其次是PDOS，在导带部分应该是p轨道比s轨道贡献大，而我的刚好相反。


初入论坛，金币不多，但请多指教！非常感谢！

我的图



文献图
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
木虫一枚，多多指教
1楼 2011-11-21 16:19:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

likai0106

金虫 (小有名气)

【答案】应助回帖

★ ★
mazuju028(金币+2): 谢谢交流 2011-11-21 18:09:59
mazuju028: 回帖置顶 对题。 2011-11-21 18:10:29
有部分分波程序对于DOSCAR的非自旋处理不准确，特别是p和d的态密度。其中脚本中并没有将pz,px,py加和，而是直接输出的pz。想检验你的程度正确与否，直观的看脚本，如果看不出来就加上自旋算一个，看看结果。
至于4eV上面的部分，建议你看看分波里面这个贡献来自于哪里，之后就能知道大致其问题的所在了。
赞一下(1人) 回复此楼
3楼2011-11-21 16:36:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

小豆sumy

金虫 (小有名气)
就是价带部分很符合，导带部分差别太大
赞一下 回复此楼
木虫一枚，多多指教
2楼2011-11-21 16:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小豆sumy

金虫 (小有名气)
引用回帖:
3楼: Originally posted by likai0106 at 2011-11-21 16:36:48:
有部分分波程序对于DOSCAR的非自旋处理不准确，特别是p和d的态密度。其中脚本中并没有将pz,px,py加和，而是直接输出的pz。想检验你的程度正确与否，直观的看脚本，如果看不出来就加上自旋算一个，看看结果。
至于 ...

谢谢你！我刚刚试过，加自旋后总态密度确实有所改善，在4eV左右就收敛了。但是PDOS还是存在同样的问题。
就是这个图中，在费米能级附近，两边都应该是P轨道贡献大，而在0eV右边导带部分却S比P峰值大，请问这是什么原因呢???　
赞一下 回复此楼
木虫一枚，多多指教
4楼2011-11-21 21:18:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★
mazuju028(金币+1): 谢谢交流 2011-11-22 09:04:29
尝试增加NBANDS
回复此楼

» 本帖已获得的红花（最新10朵）

小豆sumy
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
5楼2011-11-22 00:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-11-23 10:35:20
引用回帖:
4楼: Originally posted by 小豆sumy at 2011-11-21 21:18:18:
谢谢你！我刚刚试过，加自旋后总态密度确实有所改善，在4eV左右就收敛了。但是PDOS还是存在同样的问题。
就是这个图中，在费米能级附近，两边都应该是P轨道贡献大，而在0eV右边导带部分却S比P峰值大，请问这是什 ...

这个是一个原子的还是总的情况啊？首先你先找到对应每种原子的各个轨道态密度对应好，有时不同原子的贡献会不同。比如如果你做ZnO，价带和导带部分都不是一个原子提供的，如果都是一个原子提供的也体现不出来成键的效果啊。
赞一下(1人) 回复此楼
6楼2011-11-22 08:01:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小豆sumy

金虫 (小有名气)
送鲜花一朵
引用回帖:
5楼: Originally posted by iamikaruk at 2011-11-22 00:23:25:
尝试增加NBANDS

还是没有解决，我再看看吧~谢谢~
赞一下 回复此楼
木虫一枚，多多指教
7楼2011-11-22 13:09:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小豆sumy

金虫 (小有名气)
引用回帖:
6楼: Originally posted by likai0106 at 2011-11-22 08:01:13:
这个是一个原子的还是总的情况啊？首先你先找到对应每种原子的各个轨道态密度对应好，有时不同原子的贡献会不同。比如如果你做ZnO，价带和导带部分都不是一个原子提供的，如果都是一个原子提供的也体现不出来成键 ...

嗯，这是一个Central Si 原子的DOS，我是参照文献做的。文献画出了最中心和最边缘硅原子的DOS，但与我画的有出入。
赞一下 回复此楼
木虫一枚，多多指教
8楼2011-11-22 13:11:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小豆sumy

金虫 (小有名气)
引用回帖:
3楼: Originally posted by likai0106 at 2011-11-21 16:36:48:
有部分分波程序对于DOSCAR的非自旋处理不准确，特别是p和d的态密度。其中脚本中并没有将pz,px,py加和，而是直接输出的pz。想检验你的程度正确与否，直观的看脚本，如果看不出来就加上自旋算一个，看看结果。
至于 ...

你好，如果对p轨道存在脚本没有将pz,px,py加和，只输出pz，那要怎么解决呢？
赞一下 回复此楼
木虫一枚，多多指教
9楼2011-11-23 09:56:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

【答案】应助回帖

引用回帖:
9楼: Originally posted by 小豆sumy at 2011-11-23 09:56:51:
你好，如果对p轨道存在脚本没有将pz,px,py加和，只输出pz，那要怎么解决呢？

修改脚本就行了，你打开脚本会看到一个awk的输出其中有一个写的是自旋不为1的情况，如果学过c语言你根据你的知识会看到一个print的语句后面有输出情况其后面有一列是$2,$3,$4,你对应你的DOSCAR就会知道这种输出是有问题的。所以改成DOSCAR中对应的列就行了，比如你的计算如果没有自旋的话应该是输出1+1+3+5+7=17列的，所以改一下就行了。
赞一下 回复此楼
10楼2011-11-23 11:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小豆sumy 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 7/350 2026-08-28 07:06 by hejicker
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[教师之家] 售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
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭