24小时热门版块排行榜    

查看: 2222  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 286求调剂 +8 Faune 2026-03-21 8/400 2026-03-21 17:16 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 08工科 320总分 求调剂 +6 梨花珞晚风 2026-03-17 6/300 2026-03-21 03:40 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
请填处理意见