| 查看: 1307 | 回复: 11 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 阿勇 的 4 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
阿勇铜虫 (小有名气)
|
[交流]
【求助】VASP中 原子能级问题
|
||
|
想用VASP计算得到H，O，Zn各个原子轨道能级，想用他们的轨道能级来分析他们的成键. 首先，我分别在20*20*20 的大的supecell里 放入了单个原子，k点选G点，计算得到他们的轨道能量分别是: H: k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -6.3540 1.00000 (H: 1s) 2 -0.0494 0.00000 3 0.2821 0.00000 4 0.2918 0.00000 O: k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -23.7275 2.00000 (O:2s) 2 -9.1859 1.33333 (O:2p) 3 -9.1859 1.33333 (O:2p) 4 -9.1859 1.33333 O:2p) 5 -0.1450 0.00000 6 0.1190 0.00000 7 0.3438 0.00000 8 0.3443 0.00000 Zn: k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -10.4277 2.00000 (Zn:3d) 2 -10.4277 2.00000 (Zn:3d) 3 -10.4277 2.00000 (Zn:3d) 4 -10.4277 2.00000 (Zn:3d) 5 -10.4277 2.00000 (Zn:3d) 6 -6.1993 2.00000 (Zn:4s) 7 -1.2219 0.00000 8 -1.2219 0.00000 9 -1.2219 0.00000 我想问的是把H，O,Zn 的各个能级放在一起比较可不可行，也就是说，由VASP计算得到的各个原子对应的能量他们的参考零点是否相同？ 我还把O,Zn原子放在同一个supercell 里计算，O Zn之间的距离大概有17A，认为这样得到的轨道能级能反映出他们各自的原子轨道能级，而其他们轨道能量的参考标准应该相同。但计算完后的结果有些不可思议： O, Zn : k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -21.6514 2.00000 (O:2s) 2 -11.7000 2.00000 (Zn:3d) 3 -11.7000 2.00000 (Zn:3d) 4 -11.6998 2.00000 (Zn:3d) 5 -11.6998 2.00000 (Zn:3d) 6 -11.6998 2.00000 (Zn:3d) 7 -7.2443 1.829140 (Zn:4s) 8 -7.1836 1.39029 (O:2p) 9 -7.1836 1.39029 (O:2p) 10 -7.1836 1.39029 (O:2p) 11 -2.0005 0.00000 12 -2.0005 0.00000 13 -2.0005 0.00000 .... 结果发现 Zn 4s 与O 2p之间有电荷转移, 两个原子相隔这么远还能有电荷转移？ 后来我把胞扩大到30*30*30 后发现还是只不过转移量小了一些而已 |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有18人回复
-
垃圾破二本职称评审标准
已经有12人回复
-
职称评审没过，求安慰
已经有17人回复
-
EST投稿状态问题
已经有7人回复
-
谈谈两天一夜的“延安行”
已经有15人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
happyjwx
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3094.5
- 散金: 70
- 帖子: 168
- 在线: 174.5小时
- 虫号: 934744
- 注册: 2009-12-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★ ★ ★
阿勇(金币+2,VIP+0):先谢谢了 1-11 13:15
aylayl08(金币+5,VIP+0):谢谢详细解释 1-11 19:22
阿勇(金币+2,VIP+0):先谢谢了 1-11 13:15
aylayl08(金币+5,VIP+0):谢谢详细解释 1-11 19:22
|
我认为你的计算是有问题的！ 首先，VASP是采用平面波基组而非LCAO，所以根本不要去考虑什么mulliken布居分布，根本不存在。电荷转移只能表征在空间上，而不能说是在哪个原子上面。空间上的转移可以通过计算电荷密度差来研究，原子布居分部目前只能用bader势来看。 然后，你看到的所谓电荷转移，其实是出现了轨道分数占据态。能带中不应该出现分数占据，出现分数坐标说明你的计算中使用了smearing方法来加速SCF收敛，比如说ISMEAR=2，这个时候虽然收敛的很好很速度，但是得到的结果根本不是能量最低。结构已然不对怎么能算对电子结构呢？要算能带结构，ISMEAR要设置为-5才行！ 最后，你的原子距离调整到很大，体系的总轨道就是两原子轨道的简单叠加，可以分别算一个各个原子的能级，分开来看比较好。 [ Last edited by happyjwx on 2010-1-11 at 12:30 ] |
3楼2010-01-11 12:27:58
westmonster
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5637.3
- 帖子: 401
- 在线: 236.4小时
- 虫号: 456224
- 注册: 2007-11-11
- 专业: 凝聚态物性 II ：电子结构
2楼2010-01-11 11:38:57
阿勇
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184.9
- 散金: 61
- 帖子: 86
- 在线: 16.5小时
- 虫号: 33779
- 注册: 2004-01-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-01-11 13:18:10
happyjwx
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3094.5
- 散金: 70
- 帖子: 168
- 在线: 174.5小时
- 虫号: 934744
- 注册: 2009-12-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2010-01-11 13:21:22