24小时热门版块排行榜    

查看: 3767  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 07化学280分求调剂 +4 722865 2026-03-23 4/200 2026-03-24 00:01 by chixmc
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 材料与化工考研调剂 +4 孅華 2026-03-22 4/200 2026-03-23 16:13 by 一休哥FU
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-17 10/500 2026-03-23 13:05 by 醉在风里
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
信息提示
请填处理意见