24小时热门版块排行榜    

查看: 2307  |  回复: 10

Rosk

银虫 (正式写手)

[求助] 电子组态的计算(价电子跑哪里去呢?)

最近用VASP计算了一个团簇投了一篇论文,最近收到审稿意见,编辑给了修改,意见都还可以,但遇到了一个难题,求助大家!

Ti原子电子基态是3d2 4s2,P原子是3s2 3p3.
我们用VASP-GGA(PAW势)计算了TiP两个原子的组合,得到Ti的3d-, 4s-轨道 是2.511 和 0.502,P原子的3s-, 3p轨道 是1.107和 1.385,因此对于TiP来说 2.511 + 0.502 + 1.107 + 1.385 = 5.5;但是从TiP的电子基态Ti3d,Ti4s,P3s,P3p知道电子总数是2+2+2+3=9。 9-5.5=3.5
审稿人用高斯计算了一下TiP,3d-, 4s-orbitals of Ti and 3s-, 3p-orbitals of P is 2.4 + 0.9 + 1.9 + 3.5 = 8.7 (close to 9)。

因此,审稿人问另外那3.5 电荷到哪里去呢?本人对这个原理了解有限,请求大家帮助。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

★ ★
ben_ladeng(金币+2): 谢谢提醒 2011-08-21 09:00:59
vasp里计算电荷要调整好RWIGS,用默认的值基本都偏小,丢掉的电荷跑到RWIGS之外了

» 本帖已获得的红花(最新10朵)

2楼2011-08-20 18:43:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rosk

银虫 (正式写手)

送鲜花一朵
引用回帖:
2楼: Originally posted by hakuna at 2011-08-20 18:43:33:
vasp里计算电荷要调整好RWIGS,用默认的值基本都偏小,丢掉的电荷跑到RWIGS之外了

非常谢谢,不知道是否有合适的设置经验推荐

由于你没点应助,导致我不能评分
3楼2011-08-20 19:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

你在INCAR里设一下RWIGS,在POTCAR里RWIGS值得基础上增加10%试试看,说穿了就是凑一下,通过合理调整RWIGS,让计算出来的总电子数与理论值一致就行。之后或积分DOS或Bader......

» 本帖已获得的红花(最新10朵)

4楼2011-08-20 19:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rosk

银虫 (正式写手)

送鲜花一朵
引用回帖:
4楼: Originally posted by hakuna at 2011-08-20 19:13:05:
你在INCAR里设一下RWIGS,在POTCAR里RWIGS值得基础上增加10%试试看,说穿了就是凑一下,通过合理调整RWIGS,让计算出来的总电子数与理论值一致就行。之后或积分DOS或Bader......

谢谢,可以设置rwigs之后,然后设置LORBIT =11,直接读出这个电荷, 不需要用积分dos可以吗
5楼2011-08-20 19:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

★ ★
ben_ladeng(金币+2): 多多参与讨论哦 2011-08-21 09:02:17
你可以应用bader电荷分析给出合理的电荷的数值,应用vasp可以做到。细节请参见:
http://theory.cm.utexas.edu/vtsttools/bader/

» 本帖已获得的红花(最新10朵)

6楼2011-08-21 00:19:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rosk

银虫 (正式写手)

送鲜花一朵
引用回帖:
6楼: Originally posted by uuv2010 at 2011-08-21 00:19:30:
你可以应用bader电荷分析给出合理的电荷的数值,应用vasp可以做到。细节请参见:
http://theory.cm.utexas.edu/vtsttools/bader/

请问 这个方法可以准确给出每个轨道的电荷数目吗?
7楼2011-08-21 09:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
8楼2011-08-21 09:53:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★ ★
Rosk(金币+5): 2011-08-22 08:25:31
youzhizhe(金币+2): 谢谢交流。 2011-08-22 20:37:50
vasp计算出的分立态密度的电荷值本身就不太准。没有太大意义,只能定性的分析。
可以适当的给它归一化。
例如castep就是把它归一化了。
真理是一点点接近的!
9楼2011-08-21 22:45:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖


Rosk(金币+15): 3 2011-08-22 08:25:36
youzhizhe(金币+1): 谢谢交流。 2011-08-22 20:38:04
因为审稿人用的是高斯基组, 在电荷守恒方面肯定要比vasp好。
真理是一点点接近的!
10楼2011-08-21 22:46:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Rosk 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 16:47 by skyjudy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见