24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2303  |  回复: 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

木虫 (正式写手)

【答案】应助回帖


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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考博] 申博/考博 +3 啃面包的小书虫 2026-04-17 4/200 2026-04-17 23:54 by 阳阳阳^_^
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 药学求调剂 +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
信息提示
请填处理意见