| 查看: 1737 | 回复: 11 | ||||
[交流]
【求助】split_dos 已有2人参与
|
|
用split_dos分出来的DOS，O原子有d轨道，而且值还不小，最大有大约0.3。 而且每个原子都有6列（自旋极化），画出来的s p d感觉很接近，差别很小。 更重要的是用的是稀土Ho-3赝势，但DOS也有6列，而且后面值都很大，看OUTCAR里，没有f电子。 total charge # of ion s p d tot ---------------------------------------- 1 1.562 3.459 0.000 5.021 2 1.561 3.465 0.000 5.026 3 1.561 3.465 0.000 5.026 4 1.562 3.459 0.000 5.021 5 1.562 3.459 0.000 5.021 6 1.561 3.465 0.000 5.026 7 1.561 3.465 0.000 5.026 8 1.562 3.459 0.000 5.021 9 1.564 3.489 0.000 5.053 10 1.564 3.489 0.000 5.053 11 1.564 3.489 0.000 5.053 12 1.564 3.489 0.000 5.053 13 0.313 0.471 4.883 5.667 14 0.314 0.473 4.884 5.670 15 0.341 0.476 6.004 6.821 16 0.341 0.476 6.002 6.819 17 0.330 6.822 1.257 8.409 18 0.330 6.823 1.257 8.410 19 0.330 6.822 1.257 8.409 20 0.330 6.823 1.257 8.410 其中17-20是稀土，不看charge貌似是f电子写到p上面了，同时d对应的应该是s电子。 但是看磁矩magnetization (x) # of ion s p d tot ---------------------------------------- 1 -0.002 -0.070 0.000 -0.073 2 -0.005 -0.044 0.000 -0.050 3 -0.005 -0.044 0.000 -0.049 4 -0.002 -0.070 0.000 -0.073 5 -0.002 -0.070 0.000 -0.073 6 -0.005 -0.044 0.000 -0.050 7 -0.005 -0.044 0.000 -0.049 8 -0.002 -0.070 0.000 -0.073 9 -0.005 -0.058 0.000 -0.063 10 -0.005 -0.059 0.000 -0.063 11 -0.005 -0.058 0.000 -0.063 12 -0.005 -0.059 0.000 -0.063 13 0.011 0.013 3.005 3.029 14 0.011 0.013 2.993 3.018 15 -0.012 -0.015 -3.634 -3.661 16 -0.012 -0.015 -3.635 -3.662 17 0.002 0.002 0.007 0.010 18 0.002 0.002 0.007 0.010 19 0.002 0.002 0.007 0.010 20 0.002 0.002 0.007 0.010 稀土基本没有磁矩，我是设置了稀土磁矩的，不知道这个是怎么回事 赝势里显示Ho-3是Xe4，应该有4f电子的。 请问牛人，这个是什么问题。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
材料专硕306英一数二
已经有6人回复
-
伙伴们，祝我生日快乐吧
已经有26人回复
-
机械专硕325，寻找调剂院校
已经有4人回复
-
333求调剂
已经有7人回复
-
考研化学学硕调剂，一志愿985
已经有4人回复
-
0854控制工程 359求调剂 可跨专业
已经有9人回复
-
梁成伟老师课题组欢迎你的加入
已经有9人回复
-
化学调剂0703
已经有8人回复
-
环境工程调剂
已经有6人回复
-
326求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 等效位置的态密度分析
已经有9人回复
- DOSCAR画态密度
已经有16人回复
- VASP计算出来的PDOS怎么画图啊？求指教
已经有5人回复
- 【求助成功】split_dos分割DOSCAR的问题
已经有14人回复
- 【求助】有没有能处理LORBIT=12的DOSCAR的软件？
已经有6人回复
- 【分享】处理DOSCAR文件的脚本和程序
已经有19人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
2楼2010-05-14 09:51:28
|
SYSTEM = HoMnFeO3 ISTART = 1 ICHARG = 11 ENCUT = 500 GGA = PE #NELMIN = 4 #NSW = 60 #IBRION = 2 PREC = Accurate ALGO = Fast #ISIF = 3 ISMEAR = -5 LORBIT = 12 LREAL = Auto EDIFF = 1E-5 EDIFFG =-1E-2 ISPIN = 2 MAGMOM = 12*0 1*4 1*4 1*-5 1*-5 1*10 1*-10 1*10 1*-10 LMAXMIX = 6 AMIX = 0.2 BMIX = 0.0001 AMIX_MAG = 0.8 BMIX_MAG = 0.0001 |
3楼2010-05-14 10:16:55
4楼2010-05-14 10:19:12
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
5楼2010-05-14 13:41:30
6楼2010-05-14 13:52:51
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
zhang668(金币+2):非常感谢！ 2010-05-15 11:03:24
|
是这样的，split_dos脚本是提取DOSCAR里的前四列（no-spin），或前7列（spin）。你的体系考虑了自旋，就是7列。 如果你设置LORBIT=12，DOSCAR里会有1+1+3+5=10列（no-spin)，或1+2+6+10=19列（spin）。你的体系是19列。 这时，你如果还用split_dos脚本来分割，得到的DOSi （i为原子序号）文件中的7列数据分别对应 : 能量，s(up), s(down), py(up), py(down), pz(up), pz(down). 所以你认为的O的d 轨道实际上是它的pz 轨道，态密度的值不小也就不奇怪了。 如果你不关心角动量态密度的话，可以设置LORBIT=10 |
7楼2010-05-14 14:49:18
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
8楼2010-05-14 14:50:23
9楼2010-05-15 11:05:04
10楼2010-05-15 11:18:24
