| 查看: 1682 | 回复: 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 |
» 猜你喜欢
为什么nbs上溴 没有产物点出现呢
已经有7人回复
-
拟解决的关键科学问题还要不要写
已经有9人回复
-
最失望的一年
已经有17人回复
-
求推荐博导
已经有4人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 等效位置的态密度分析
已经有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