| 查看: 1606 | 回复: 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 |
» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
等效位置的态密度分析
已经有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