24小时热门版块排行榜    

查看: 1789  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ip0371

木虫 (正式写手)

[求助] 画出每一层原子的PDOS图以及获取各原子的磁矩 已有1人参与

各位同仁:

如下图附件所示,研究体系是一个由左、中、右三个区域拼接成的一个晶体(即磁性隧道结):
中间的粉红色代表Al原子;
灰色表示Cr原子;
黄色表示Si原子;
蓝色表示Co原子。

我用VASP计算它并得到DOSCAR文件后,使用split_dos脚本将其分割成
DOS0、DOS1、DOS2...

这里需要请教大家的是:
1.如何获得中间区域三层Al原子(粉红色表示的原子)每一层原子的PDOS图?
2.怎样获得界面处每一层原子的磁矩大小?当然,计算前会在INCAR中设置ISPIN=2。貌似MAGMOM这个参数是设置初始磁矩的,但是不清楚如何设?

不懂这些技巧和使用方法,还请大家多多指教!谢谢了

画出每一层原子的PDOS图以及获取各原子的磁矩
pic01.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ip0371

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 郭芳8888 at 2014-06-13 11:17:21
分裂开后dos0是总的,而dos1 dos2....就是单个原子的,导出后可以画出能态密度图,也就是dos0是整个体系的能态密度图,而dos1 dos2 ...是每个原子的;
在OUTCAR中 输入/magnetization 可以查看每个原子的磁矩,顺序 ...

谢谢你的回答

追问:

按你说的DOS1、DOS2、DOS3...就是指的单个原子的分波态密度(还是局域态密度啊?),找到处于同一层的所有原子,用sum脚本将这几个原子态密度叠加后就是这一层原子的总的态密度吗?不好意思分波态密度和局域态密度想请你简单区别解释一下可以吗?呵呵呵,谢谢啦
3楼2014-06-13 11:41:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

郭芳8888

金虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-06-24 11:43:31
分裂开后dos0是总的,而dos1 dos2....就是单个原子的,导出后可以画出能态密度图,也就是dos0是整个体系的能态密度图,而dos1 dos2 ...是每个原子的;
在OUTCAR中 输入/magnetization 可以查看每个原子的磁矩,顺序是对应POSCAR的原子坐标次序。MAGMOM设置的是初始磁矩,一般尽量设置的接近理论值,当然你设置为其他的,只是计算时间要长一点。
希望对你有用。
2楼2014-06-13 11:17:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ip0371

木虫 (正式写手)

知道的朋友请帮忙啊~~~~~~~~~~

DOSCAR里的数据怎么有10列之多啊??(在INCAR里设置的ISPIN=2)

4楼2014-06-13 15:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见