24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2050  |  回复: 13

klwang

木虫 (著名写手)

[求助] DOSCAR 中Mg的PDOS有d轨道态密度已有1人参与

用VASP计算含Mg和Nd的化合物,考虑自旋。
Mg的膺势不含d轨道,为什么画出来的Mg的PDOS有d轨道成分(见图)?请大侠指点。

我理解的DOSCAR的数据顺序 ,并且根据此进行的数据归总和分析。希望这个顺序没有错误。
Energy(eV)        s_1_up        s_1_dn        p_1_up        p_1_dn        p_2_up        p_2_dn        p_3_up        p_3_dn        d_1_up        d_1_dn        d_2_up        d_2_dn        d_3_up        d_3_dn        d_4_up        d_4_dn        d_5_up        d_5_dn        f_1_up        f_1_dn        f_2_up        f_2_dn        f_3_up        f_3_dn        f_4_up        f_4_dn        f_5_up        f_5_dn        f_6_up        f_6_dn        f_7_up        f_7_dn

另外,我没有用VASP的tool kits。因为有明显的bug,我自己写的PERL程序生成了原子轨道PDOS,保留了原有数据格式,然后用MATLAB把每个原子轨道PDOS按元素加起来。
PERL程序是没有问题的,我和原始数据比对过了,没有出入。对于MATLAB,我是新手,就是把文件读入矩阵,然后矩阵相加。结果就出现了很奇怪的d轨道成分在Mg的PDOS里。做法如有错误,请大侠指点。或请大侠分享分析方法和程序。

谢谢!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Mg-PDOS.jpg
  • 2015-08-13 12:21:44, 1.29 M

» 猜你喜欢

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
klwang(ljw4010代发): 金币+13, 代发奖励 2015-08-15 08:26:56
有种可能性,在做WS半径的时候Nd有一部分被划入Mg的球体中。

» 本帖已获得的红花(最新10朵)

不一定挂在论坛,计算问题问题欢迎留言。
4楼2015-08-14 03:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖


klwang: 金币+1, 有帮助 2015-08-14 22:59:27
你得看d轨道的成分多少。。和p轨道相比是不是算少。。。POTCAR里面只是原子的价电子,并不代表计算完之后不会有d轨道的成分,因为在和其他原子杂化的时候出现d的成分,我认为很正常。
凝聚态理论与计算物理
8楼2015-08-14 22:31:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
klwang: 金币+5, 有帮助, 谢谢大侠回复 2015-08-13 23:26:17
没有d,一般来说DOSCAR也会有d,但是是否有f得看你的元素。你应该把INCAR、DOSCAR的一部分贴出来看看,才好分析
凝聚态理论与计算物理
2楼2015-08-13 12:51:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klwang

木虫 (著名写手)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-08-12 16:51:15
没有d,一般来说DOSCAR也会有d,但是是否有f得看你的元素。你应该把INCAR、DOSCAR的一部分贴出来看看,才好分析

谢谢回复。
我的材料里面有Nd,含有d和f轨道,但是Mg是没有的,不知道为什么结果这么奇怪。
另外,我注意到OUTCAR里,费米面的电子数不为整数,请问大侠遇到过这种情况吗?

System =  Ba24NdO72Ta12Mg5Zr5
ISTART = 0        
ISMEAR = 0        
SIGMA = 0.05      
PREC  = Accurate  
ENCUT = 615 eV   
LREAL  = Auto     
EDIFF = 1.0E-5   
EDIFFG = -1.0E-2   
ISPIN = 2         
NBANDS = 800      
MAXMIX = 40
NPAR   = 4        
ALGO   = Normal
IBRION = 2        
NSW    = 1        
ISIF   = 3        

LORBIT = 11      
EMIN = -65        
EMAX =  80        
NEDOS = 800      

LDAU = .TRUE.  
LDAUTYPE = 2   
LDAUL = 0 3 0 0 0 0
LDAUU = 0 6 0 0 0 0
LDAUJ = 0 0 0 0 0 0
LDAUPRINT = 2     

MAGMOM = 24*0 2*3 72*0 12*0 5*0 5*0         

LCHARG  = .True.  
LWAVE  = F

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : mg_doscar.dat
  • 2015-08-13 23:23:09, 1.81 M
3楼2015-08-13 23:24:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖


klwang: 金币+1, 有帮助, 抱歉,您的回复不能接受。Mg的d轨道没有在POTCAR里,是不应该出现这种情况的。 2015-08-14 22:21:37
引用回帖:
3楼: Originally posted by klwang at 2015-08-13 23:24:42
谢谢回复。
我的材料里面有Nd,含有d和f轨道,但是Mg是没有的,不知道为什么结果这么奇怪。
另外,我注意到OUTCAR里,费米面的电子数不为整数,请问大侠遇到过这种情况吗?

System =  Ba24NdO72Ta12Mg5Zr5
IS ...

有d正常吧。。毕竟不是晶体,不是原子,有d的成分正常。
凝聚态理论与计算物理
5楼2015-08-14 09:09:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

4樓正解
Simpler, stronger, more tolerant
6楼2015-08-14 14:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klwang

木虫 (著名写手)

送红花一朵
引用回帖:
4楼: Originally posted by 卡开发发 at 2015-08-13 07:25:13
有种可能性,在做WS半径的时候Nd有一部分被划入Mg的球体中。

谢谢回复,因为您没有用应助回复,我没有办法给您金币。送上红花一朵。
7楼2015-08-14 22:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

而且VASP用的是平面波的基,最后投影出d的成分,很正常吧。
凝聚态理论与计算物理
9楼2015-08-14 22:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klwang

木虫 (著名写手)

引用回帖:
8楼: Originally posted by wmy8802217 at 2015-08-14 02:31:18
你得看d轨道的成分多少。。和p轨道相比是不是算少。。。POTCAR里面只是原子的价电子,并不代表计算完之后不会有d轨道的成分,因为在和其他原子杂化的时候出现d的成分,我认为很正常。

问题是我的材料是块材料,不是cluster或molecule,实验上表明是半导体或绝缘体,离子键,Mg没有d轨道参与计算,结果出现d轨道的电子态,您不觉得这种杂化不是很奇怪吗?如果Mg有d轨道参与计算,我不否认您的说法。
10楼2015-08-14 22:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 klwang 的主题更新
信息提示
请填处理意见