24小时热门版块排行榜    

查看: 2375  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-16 18:00 by i3KBRD19ERkG
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:48 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:03 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:15 by i3KBRD19ERkG
[找工作] 浙江师范大学是怎么坑我的 +7 recruit123 2026-09-09 7/350 2026-09-16 15:07 by 防腐剂科技部
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:27 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见