24小时热门版块排行榜    

查看: 2040  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhang668

木虫 (正式写手)

[交流] 【求助】split_dos 已有2人参与

用split_dos分出来的DOS,O原子有d轨道,而且值还不小,最大有大约0.3。
而且每个原子都有6列(自旋极化),画出来的s p d感觉很接近,差别很小。
更重要的是用的是稀土Ho-3赝势,但DOS也有6列,而且后面值都很大,看OUTCAR里,没有f电子。
total charge     

# of ion     s       p       d       tot
----------------------------------------
  1        1.562   3.459   0.000   5.021
  2        1.561   3.465   0.000   5.026
  3        1.561   3.465   0.000   5.026
  4        1.562   3.459   0.000   5.021
  5        1.562   3.459   0.000   5.021
  6        1.561   3.465   0.000   5.026
  7        1.561   3.465   0.000   5.026
  8        1.562   3.459   0.000   5.021
  9        1.564   3.489   0.000   5.053
10        1.564   3.489   0.000   5.053
11        1.564   3.489   0.000   5.053
12        1.564   3.489   0.000   5.053
13        0.313   0.471   4.883   5.667
14        0.314   0.473   4.884   5.670
15        0.341   0.476   6.004   6.821
16        0.341   0.476   6.002   6.819
17        0.330   6.822   1.257   8.409
18        0.330   6.823   1.257   8.410
19        0.330   6.822   1.257   8.409
20        0.330   6.823   1.257   8.410
其中17-20是稀土,不看charge貌似是f电子写到p上面了,同时d对应的应该是s电子。
但是看磁矩magnetization (x)

# of ion     s       p       d       tot
----------------------------------------
  1       -0.002  -0.070   0.000  -0.073
  2       -0.005  -0.044   0.000  -0.050
  3       -0.005  -0.044   0.000  -0.049
  4       -0.002  -0.070   0.000  -0.073
  5       -0.002  -0.070   0.000  -0.073
  6       -0.005  -0.044   0.000  -0.050
  7       -0.005  -0.044   0.000  -0.049
  8       -0.002  -0.070   0.000  -0.073
  9       -0.005  -0.058   0.000  -0.063
10       -0.005  -0.059   0.000  -0.063
11       -0.005  -0.058   0.000  -0.063
12       -0.005  -0.059   0.000  -0.063
13        0.011   0.013   3.005   3.029
14        0.011   0.013   2.993   3.018
15       -0.012  -0.015  -3.634  -3.661
16       -0.012  -0.015  -3.635  -3.662
17        0.002   0.002   0.007   0.010
18        0.002   0.002   0.007   0.010
19        0.002   0.002   0.007   0.010
20        0.002   0.002   0.007   0.010
稀土基本没有磁矩,我是设置了稀土磁矩的,不知道这个是怎么回事
赝势里显示Ho-3是Xe4,应该有4f电子的。
请问牛人,这个是什么问题。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

split-dos不是就分割出分波态密度吗?
LORBIT=应该为几?
6楼2010-05-14 13:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

zhang668

木虫 (正式写手)

SYSTEM = HoMnFeO3
ISTART = 1
ICHARG = 11
ENCUT = 500
GGA = PE
#NELMIN = 4
#NSW    = 60
#IBRION = 2
PREC   = Accurate
ALGO   = Fast
#ISIF   = 3
ISMEAR = -5
LORBIT = 12
LREAL = Auto
EDIFF = 1E-5
EDIFFG =-1E-2
ISPIN = 2
MAGMOM = 12*0 1*4 1*4 1*-5 1*-5 1*10 1*-10 1*10 1*-10
LMAXMIX = 6
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001
3楼2010-05-14 10:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

稀土没有f电子的原因是赝势的原因,Ho-3确实f电子是芯态。
现在的主要问题就是分割的DOS,O原子d轨道上值比较大。而且每个原子的s p d都很接近,不知道是不是哪里出了问题,当然也有轨道杂化的可能。
4楼2010-05-14 10:19:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)

zhang668(金币+2): 2010-05-14 13:51:15
问题在于你在INCAR里设置了LORBIT = 12,算出来的是角动量态密度。这时候就不能再用split_dos脚本来分割DOSCAR.
5楼2010-05-14 13:41:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:48 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:34 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:53 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 18:07 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 17:34 by hZiFeudZyoGR
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见