24小时热门版块排行榜    

查看: 2627  |  回复: 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的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

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

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的回帖

卡开发发

专家顾问 (著名写手)

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然放榜时间 +4 布布和一二 2026-08-08 4/200 2026-08-09 22:24 by otj2008
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见