24小时热门版块排行榜    

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

卡开发发

专家顾问 (著名写手)

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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +8 archer.. 2026-04-09 9/450 2026-04-11 10:35 by zhq0425
[考研] 0831生医工第一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 13/650 2026-04-11 10:04 by maddjdld
[考研] 电子信息279求调剂,有书读就行 +5 wwwooden 2026-04-08 8/400 2026-04-11 09:30 by zhq0425
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +4 弱水听文 2026-04-09 4/200 2026-04-10 22:01 by doctff
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 085600材料与化工专硕329 求调剂 +24 额cc 2026-04-06 25/1250 2026-04-09 16:01 by wp06
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 材料工程专业日语生求调剂 +9 111623 2026-04-07 9/450 2026-04-07 23:31 by 一只好果子?
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
信息提示
请填处理意见