| 查看: 2019 | 回复: 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 |
» 猜你喜欢
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有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
20