| 查看: 975 | 回复: 8 | ||
| 当前主题已经存档。 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[交流]
【求助】VASP程序如何计算体系的电荷转移情况
|
||
|
我计算的是碳体系，用的是超软赝势，每个碳原子应该有2s2p电子，可是我计算出来的碳的电子数大约是1S2p电子，请问各位这对计算体系的能量有影响吗？而且我需要计算管子的电荷转移，所以一定要计算准确，请问应该如何调试？ [ Last edited by wuchenwf on 2009-6-21 at 20:57 ] |
回复此楼
4楼2008-09-08 17:18:32
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
5楼2008-09-09 09:19:17
