应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 怎样从outcar中看d轨道占据数
91/1返回列表
查看: 3719  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jugengfans

金虫 (正式写手)

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

直接到结尾处看 total charge吗？

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

发自小木虫IOS客户端

[ Last edited by ljw4010 on 2016-5-8 at 14:29 ]
回复此楼

» 猜你喜欢
1楼 2016-05-06 21:43:00
已阅   关注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的方法也各有各的局限性。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
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 的主题更新
91/1返回列表
信息提示
关闭
请填处理意见
关闭