| 查看: 1067 | 回复: 8 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】VASP程序如何计算体系的电荷转移情况
|
|||
|
我计算的是碳体系，用的是超软赝势，每个碳原子应该有2s2p电子，可是我计算出来的碳的电子数大约是1S2p电子，请问各位这对计算体系的能量有影响吗？而且我需要计算管子的电荷转移，所以一定要计算准确，请问应该如何调试？ [ Last edited by wuchenwf on 2009-6-21 at 20:57 ] |
回复此楼» 猜你喜欢
面上提前没消息，有中的吗
已经有17人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
2楼2008-09-07 10:02:24
|
我计算的结果是这样的， 第一次，我设定RWIGS=0.863，（POTCAR）中碳的默认值，在计算得到的OUTCAR中，我发现计算的charge为： # of ion s p d tot ---------------------------------------- 1 0.878 1.915 0.074 2.867 2 0.878 1.915 0.074 2.867 3 0.878 1.915 0.074 2.867 4 0.878 1.914 0.074 2.866 5 0.878 1.915 0.074 2.867 6 0.878 1.915 0.074 2.867 7 0.878 1.915 0.074 2.867 8 0.878 1.915 0.074 2.867 9 0.878 1.915 0.074 2.867 10 0.878 1.915 0.074 2.867 11 0.878 1.914 0.074 2.866 12 0.878 1.914 0.074 2.866 13 0.878 1.915 0.074 2.867 14 0.878 1.915 0.074 2.867 15 0.878 1.915 0.074 2.867 16 0.878 1.915 0.074 2.867 17 0.878 1.915 0.074 2.867 18 0.878 1.914 0.074 2.866 19 0.878 1.915 0.074 2.867 20 0.878 1.915 0.074 2.867 21 0.878 1.915 0.074 2.866 22 0.878 1.914 0.074 2.866 23 0.878 1.914 0.074 2.866 24 0.878 1.914 0.074 2.866 25 0.878 1.915 0.074 2.867 26 0.878 1.915 0.074 2.867 27 0.878 1.915 0.074 2.867 28 0.878 1.915 0.074 2.867 29 0.878 1.915 0.074 2.867 30 0.878 1.915 0.074 2.867 31 0.878 1.914 0.074 2.866 32 0.878 1.914 0.074 2.866 ------------------------------------------------ tot 28.10 61.26 2.36 91.73 也许是因为设定的RWIGS的值太小所以s和p轨道上的电子数太小，第二次我设定RWIGS=0.913,得到的结果为： total charge # of ion s p d tot ---------------------------------------- 1 0.975 2.133 0.100 3.208 2 0.975 2.133 0.100 3.208 3 0.975 2.133 0.100 3.208 4 0.975 2.133 0.100 3.208 5 0.975 2.133 0.100 3.208 6 0.975 2.133 0.100 3.208 7 0.975 2.133 0.100 3.208 8 0.975 2.133 0.100 3.208 9 0.975 2.133 0.100 3.208 10 0.975 2.133 0.100 3.208 11 0.975 2.133 0.100 3.208 12 0.975 2.133 0.100 3.208 13 0.975 2.133 0.100 3.208 14 0.975 2.133 0.100 3.208 15 0.975 2.133 0.100 3.208 16 0.975 2.133 0.100 3.208 17 0.975 2.133 0.100 3.208 18 0.975 2.133 0.100 3.208 19 0.975 2.133 0.100 3.209 20 0.975 2.133 0.100 3.208 21 0.975 2.133 0.100 3.208 22 0.975 2.133 0.100 3.208 23 0.975 2.133 0.100 3.208 24 0.975 2.133 0.100 3.208 25 0.975 2.133 0.100 3.208 26 0.975 2.133 0.100 3.208 27 0.975 2.133 0.100 3.208 28 0.975 2.133 0.100 3.208 29 0.975 2.133 0.100 3.208 30 0.975 2.133 0.100 3.208 31 0.975 2.133 0.100 3.208 32 0.975 2.133 0.100 3.208 ------------------------------------------------ tot 31.21 68.25 3.20 102.67 s轨道上的电子还是很小，而且p轨道上电子又已经超过2了，所以不是调节RWIGS的问题。那么丢失的那些s电子哪里去了，因要进一步计算体系的电荷转移，所以恳请大家帮忙解决一下这个问题。 |
3楼2008-09-08 16:58:52
4楼2008-09-08 17:18:32
5楼2008-09-09 09:19:17
6楼2008-09-10 00:48:48
7楼2008-09-10 00:49:34
8楼2008-09-10 08:21:15
spur
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.21
- 金币: 4994.7
- 散金: 559
- 红花: 5
- 帖子: 866
- 在线: 374.9小时
- 虫号: 215683
- 注册: 2006-03-11
- 性别: GG
- 专业: 第一性原理计算
9楼2008-09-10 08:53:52
50