| 查看: 2026 | 回复: 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 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 等效位置的态密度分析
已经有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
