24小时热门版块排行榜    

查看: 3971  |  回复: 8

jugengfans

金虫 (正式写手)

[求助] 怎样从outcar中看d轨道占据数

直接到结尾处看 total charge吗?

为什么d轨道的总电子数大于5?

发自小木虫IOS客户端

[ Last edited by ljw4010 on 2016-5-8 at 14:29 ]
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

(LORBIT=11 will make VASP write the site and L-decomposed charge density, n^I_d, in the OUTCAR file).

然后怎么看呢?
2楼2016-05-07 05:53:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

今天又想了一下,OUTCAR里面有两组数据 一组是total charge
s          p         d         total
但是这里面的d轨道居然占了5个以上的电子!不科学。。。。

如果以magnetic moment为标准呢?
d轨道的magnetic moment大概在2.8左右,倒是和实际情况挺符合

问题是为什么total charge里面 d轨道占了5个电子呢!百思不得其解。

有没有更直接的看某个原子、某个轨道占据数的方法?
3楼2016-05-07 11:44:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by jugengfans at 2016-05-07 11:44:32
今天又想了一下,OUTCAR里面有两组数据 一组是total charge
s          p         d         total
但是这里面的d轨道居然占了5个以上的电子!不科学。。。。

如果以magnetic moment为标准呢?
d轨道的magnet ...

有啥不科学的,很可能是WS半径划分的时候出现交叠,原子球内的电子数被double count了。不过要把WS半径做的reasonable个人觉得有难度,其他population的方法也各有各的局限性。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2016-05-08 08:04:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 卡开发发 at 2016-05-08 08:04:10
有啥不科学的,很可能是WS半径划分的时候出现交叠,原子球内的电子数被double count了。不过要把WS半径做的reasonable个人觉得有难度,其他population的方法也各有各的局限性。...

total charge

# of ion     s       p       d       tot
----------------------------------------
  1        0.382   0.546   5.063   5.991
  2        0.332   0.508   4.876   5.716
  3        0.355   0.568   4.918   5.841
  4        0.349   0.553   4.867   5.769
  5        0.354   0.567   4.881   5.803
  6        0.349   0.557   4.865   5.772
  7        0.354   0.567   4.881   5.802
  8        0.349   0.553   4.867   5.769
  9        0.355   0.568   4.918   5.841
10        0.332   0.508   4.876   5.716
11        0.382   0.546   5.062   5.989
12        1.581   3.368   0.000   4.949
13        1.581   3.368   0.000   4.949
14        1.574   3.437   0.000   5.011
15        1.574   3.437   0.000   5.011
16        1.571   3.433   0.000   5.004
17        1.571   3.433   0.000   5.004
18        1.570   3.439   0.000   5.009
19        1.570   3.439   0.000   5.009
20        1.570   3.449   0.000   5.019
21        1.570   3.449   0.000   5.019
22        1.569   3.437   0.000   5.006
23        1.569   3.437   0.000   5.006
24        1.570   3.449   0.000   5.020
25        1.570   3.449   0.000   5.020
26        1.570   3.439   0.000   5.009
27        1.570   3.439   0.000   5.009
28        1.571   3.433   0.000   5.004
29        1.571   3.433   0.000   5.004
30        1.574   3.437   0.000   5.011
31        1.574   3.437   0.000   5.011
32        1.581   3.368   0.000   4.949
33        1.581   3.368   0.000   4.949
------------------------------------------------
tot       38.497  81.418  54.075 173.990



magnetization (x)

# of ion     s       p       d       tot
----------------------------------------
  1        0.015   0.007   2.138   2.160
  2       -0.008  -0.012  -2.392  -2.412
  3        0.015   0.019   2.338   2.372
  4       -0.013  -0.018  -2.507  -2.538
  5        0.015   0.020   2.478   2.514
  6       -0.013  -0.019  -2.482  -2.514
  7        0.015   0.020   2.479   2.514
  8       -0.013  -0.018  -2.507  -2.538
  9        0.015   0.019   2.338   2.372
10       -0.008  -0.012  -2.392  -2.412
11        0.015   0.007   2.138   2.160
12        0.001  -0.078   0.000  -0.077
13        0.001  -0.078   0.000  -0.077
14        0.000   0.004   0.000   0.003
15        0.000   0.004   0.000   0.003
16        0.000  -0.020   0.000  -0.020
17        0.000  -0.020   0.000  -0.020
18        0.001   0.058   0.000   0.059
19        0.001   0.058   0.000   0.059
20        0.000  -0.031   0.000  -0.032
21        0.000  -0.031   0.000  -0.032
22        0.001   0.050   0.000   0.051
23        0.001   0.050   0.000   0.051
24        0.000  -0.031   0.000  -0.032
25        0.000  -0.031   0.000  -0.032
26        0.001   0.057   0.000   0.059
27        0.001   0.057   0.000   0.059
28        0.000  -0.020   0.000  -0.020
29        0.000  -0.020   0.000  -0.020
30        0.000   0.004   0.000   0.003
31        0.000   0.004   0.000   0.003
32        0.001  -0.078   0.000  -0.077
33        0.001  -0.078   0.000  -0.077
------------------------------------------------
tot        0.044  -0.158   1.628   1.514

谢谢回复,但是请看我这个结果,d轨道电荷数和d轨道的总磁矩差的有点多。。。是磁性的原子,有3个3d电子。

因为都需要把波函数投影到d轨道算占据数,为什么磁矩和d轨道占据的电子差距那么大呢?
WS元胞半径在算lm decomposed charge和磁矩的时候,应该是一样的吧。
5楼2016-05-08 08:33:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

ljw4010: 屏蔽内容, 重复 2016-05-08 14:29:17
本帖内容被屏蔽

6楼2016-05-08 08:33:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

这是一个AB2型结构,A为阳离子,B为阴离子。  A与B成键的时候失去电子,只剩下3个核外电子。。
magnetic moment 2.138个。这个比较符合3个核外电子的构向
然而,total charge的第一行 3d轨道居然有5.几个电子!

难道说,total charge可能也包含离子键成键后,离域的那部分电子。。。
但是那部分电子应该被阴离子B夺走了,远离A离子才对,为什么还会投影到A的轨道上呢?
7楼2016-05-08 08:42:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
6楼: Originally posted by jugengfans at 2016-05-08 08:33:18
total charge

# of ion     s       p       d       tot
----------------------------------------
  1        0.382   0.546   5.063   5.991
  2        0.332   0.508   4.876   5.716
  3        0.35 ...

磁矩总的是alpha+beta,磁矩是alpha-beta,显然前者会被误差放大。
不一定挂在论坛,计算问题问题欢迎留言。
8楼2016-05-08 09:44:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
7楼: Originally posted by jugengfans at 2016-05-08 08:42:05
这是一个AB2型结构,A为阳离子,B为阴离子。  A与B成键的时候失去电子,只剩下3个核外电子。。
magnetic moment 2.138个。这个比较符合3个核外电子的构向
然而,total charge的第一行 3d轨道居然有5.几个电子!
...

可能啊,甚至有些TmO(Tm=过渡金属)的体系如果用WS划分会出现O上有3d电子的情况。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2016-05-08 09:45:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见