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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 怎样从outcar中看d轨道占据数
91/1返回列表
查看: 3938  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +5 布布和一二 2026-08-10 6/300 2026-08-11 10:15 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +15 Lanmanbaby 2026-08-09 23/1150 2026-08-11 09:51 by demonstreets
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 为什么网上很多人说本周 12号出结果 +5 瞬息宇宙 2026-08-10 5/250 2026-08-11 09:29 by 铁帽子农民
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 11/550 2026-08-11 09:14 by xiaruohan
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭