24小时热门版块排行榜    

查看: 2357  |  回复: 15
【悬赏金币】回答本帖问题,作者yangbs2009将赠送您 40 个金币

yangbs2009

木虫 (小有名气)

[求助] 态密度求助已有1人参与

最近在用vasp计算单层过渡族金属硫化物的磁性和电子结构方面的性质,对得到的态密度图有些疑问,DOS如图所示
1、物质没有磁性的原因是spinup与spindown在相同能量状态占据的电子数相同;
      存在磁性时则有一部分的spindown电子转移到了spinup的态上,spinup的占据态上电子多于spindown, 而spindown的一些轨道由于没有电子占据,因而会产生空的轨道,如果积分导带和价带的所有spinup与spindown的电子,那么spinup与spindown的电子数目应该是相同的
      但是我计算的结果显示up&down电子数目并不相同,Fe的dz2轨道尤为明显
2、根据分子轨道理论,配位化合物的中心原子与周围的配位原子成键分成了bond state与antibonding state,当我在计算一些过渡族金属掺杂时,有的很明显显示了bonding state与antibonding state;但是有一些的DOS仅在一个能量范围内有电子填充,并未形成bonding&antibonding state

我感觉可能是我态密度图的做法有问题,INCAR的参数如下
general:
System =substitute-Fe
ISTART = 1
ICHARG=11
ENCUT  = 500
ISMEAR =-5
SIGMA = 0.2
#IBRION=2
#ISIF= 2
#NSW = 1000
ISPIN =2
POTIM=0.5
PRCE=High
EDIFFG=-0.01
EDIFF=1E-4
LORBIT=11
LREAL=Auto
MAGMOM = 50*0 24*0 1
NEDOS=600
RWIGS=0.953 1.625 1.302
而且对INCAR的设置做过一些更改,如增大NEDOS的值,或者ismear=0以及减小sigma的值,但结果都没有什么太大的影响

以上是我个人的一些理解,得到的dos不准确是因为incar设置不正确还是什么其他的原因我也不是很清楚,希望得到高人的指点,小弟先在此谢过!

附件是替代Fe之后的态密度图,第一个是总的dos,中间的两个是Fe的五个d轨道的pdos,最后的一个是Fe周围一个S的pdos

态密度求助
Fe.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

pangrui1985

铜虫 (小有名气)

不明白,Fe不该是有磁性的吗?
6楼2014-04-01 10:09:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yangbs2009

木虫 (小有名气)

有人可以给一些建议么,着急用啊
2楼2014-03-30 20:52:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

各位大神,拜托帮忙解决一下啊,感激不尽
3楼2014-03-31 21:40:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
建议尝试降低EDIFF进行SCF计算。

“DOS仅在一个能量范围内有电子填充”是什么意思?是说d轨道没有分裂吗?
4楼2014-04-01 00:34:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
4楼: Originally posted by KalaShayminS at 2014-04-01 00:34:30
建议尝试降低EDIFF进行SCF计算。

“DOS仅在一个能量范围内有电子填充”是什么意思?是说d轨道没有分裂吗?

现在的ediff设置的能量精度已经是0.01mev,我感觉这个精度已经足够了,再降低的话计算效率会很慢
不过还是谢谢你

DOS仅在一个能量范围内有电子填充:
分子轨道理论指明配体和中心原子之间是通过共价键作用的,因此中心原子(Fe)为了成共价键键应该从孤立原子的局域态变为扩展态,所以Fe的pdos应该是比较分散的,但是从上面的图中可以看出,大多数Fe-d的电子都在费米面附近,因此很局域,与理论上的东西不一样啊
5楼2014-04-01 07:42:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
6楼: Originally posted by pangrui1985 at 2014-04-01 10:09:53
不明白,Fe不该是有磁性的吗?

我算的是Fe替代原来体系中的金属原子
是有磁性的啊
7楼2014-04-01 12:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by yangbs2009 at 2014-04-01 07:42:49
现在的ediff设置的能量精度已经是0.01mev,我感觉这个精度已经足够了,再降低的话计算效率会很慢
不过还是谢谢你

DOS仅在一个能量范围内有电子填充:
分子轨道理论指明配体和中心原子之间是通过共价键作用的, ...

顶楼图中d轨道明显是分裂的,d轨道从-4到+4都有,费米面附近有两段空白,只是被dz的高峰掩盖掉了,因此我觉得这不是一个问题。你可以看一下纯Fe的d轨道有多宽进行对比。

顶楼的ediff是1e-4,这个不是0.01meV吧?0.01meV一般都是够的。
还有这个INCAR看起来写了很多不必要的默认参数。
8楼2014-04-01 21:58:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
8楼: Originally posted by KalaShayminS at 2014-04-01 21:58:01
顶楼图中d轨道明显是分裂的,d轨道从-4到+4都有,费米面附近有两段空白,只是被dz的高峰掩盖掉了,因此我觉得这不是一个问题。你可以看一下纯Fe的d轨道有多宽进行对比。

顶楼的ediff是1e-4,这个不是0.01meV吧? ...

我的主要问题是,dz2的自旋向上和自旋向下的电子总数明显不同
把能量范围缩小以后同样是,自旋向下的电子要多很多。
还有为什么要和纯净的Fe的d轨道进行比较呢,不同体系当中,Fe的DOS肯定是不同的啊
9楼2014-04-02 15:07:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
8楼: Originally posted by KalaShayminS at 2014-04-01 21:58:01
顶楼图中d轨道明显是分裂的,d轨道从-4到+4都有,费米面附近有两段空白,只是被dz的高峰掩盖掉了,因此我觉得这不是一个问题。你可以看一下纯Fe的d轨道有多宽进行对比。

顶楼的ediff是1e-4,这个不是0.01meV吧? ...

ediff这个是0.1mev,kbT大概是25mev
与kbT相比,我感觉这个精度足够了
10楼2014-04-02 15:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangbs2009 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 13/650 2024-11-15 07:15 by zju2000
[论文投稿] 找不到审稿人 +5 吉祥如意141913 2024-11-13 5/250 2024-11-15 07:09 by 半生梦君
[教师之家] 大学老师 +9 考研一路顺风 2024-11-13 11/550 2024-11-14 21:47 by EndNoted
[基金申请] 博后基金分组排名 +7 攻城2025 2024-11-14 7/350 2024-11-14 21:19 by 实验小白ha
[基金申请] 去年七月底入站的还能申请下一批吗? +3 brightwo 2024-11-14 3/150 2024-11-14 21:17 by 实验小白ha
[论文投稿] angew三个小修,返修过去已经20多天了没动静,有朋友遇到这种情况吗? +3 迟迟未到场 2024-11-14 3/150 2024-11-14 17:53 by TopEdit
[精细化工] 同一个反应相同的反应条件,是不是反应结果相差不大? +4 青霉素 2024-11-11 4/200 2024-11-14 17:33 by zinedine
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[基金申请] 广东省自然科学基金-面上项目消息有了么? +4 pachang 2024-11-14 5/250 2024-11-14 16:41 by owenyaa
[论文投稿] 职称论文 +8 xiaoqi_66 2024-11-09 10/500 2024-11-14 15:14 by nono2009
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[找工作] 咨询一下江西的高校待遇,人文氛围怎么样? +5 akslis2024 2024-11-09 5/250 2024-11-14 13:53 by 啄木鸟、
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[教师之家] 评正教授需要两个国家级项目,有人用子课题糊弄,结果在评审前资格公示时被举报拿下了 +20 瞬息宇宙 2024-11-12 29/1450 2024-11-13 08:49 by Quakerbird
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[有机交流] 求助NMR +5 苯巴比妥! 2024-11-09 5/250 2024-11-12 11:20 by 88817753
[硕博家园] 同步辐射 +4 小张要加油努力 2024-11-10 7/350 2024-11-11 13:47 by mpdfwxgui
[论文投稿] renewable energy 添加作者 15+3 Hebauwww12 2024-11-09 4/200 2024-11-11 08:42 by 北京莱茵润色
信息提示
请填处理意见