计算化学公社

标题: ISMEAR=-5 轨道占据数为负值，合理吗？ [打印本页]

作者
Author: helpme 时间: 2017-6-5 11:20
标题: ISMEAR=-5 轨道占据数为负值，合理吗？
计算了一个例子，一层Ni原子，
发现ISMEAR=-5时仍然有轨道占据为负值。请问这个正常吗？
这里的负占据怎么理解？

谢谢！

附件是输入文件以及OUTCAR。

作者
Author: 卡开发发 时间: 2017-6-5 20:22
试试
“
ISMEAR = 0 or 1 or 2
SIGMA = 0.1（最好测试下）
MAGMOM = 2（或者再稍大一些）
”
看看是否能够改善结果。
金属体系还是不要用ISMEAR=5吧，K点的充分性不一定能够保证，另外有磁性的体系磁矩按大点猜没问题。

作者
Author: helpme 时间: 2017-6-6 08:53
用ISMEAR=0; SIMGA=0.1试了，没有负数占据了。谢谢！
看来是把 ISMEAR=-5 过于神话了。
体系是一层Ni原子，又是金属又是二维（不适合四面体的构成），所以不适合ISMEAR=-5.

作者
Author: 卡开发发 时间: 2017-6-6 11:43

helpme 发表于 2017-6-6 08:53
用ISMEAR=0; SIMGA=0.1试了，没有负数占据了。谢谢！
看来是把 ISMEAR=-5 过于神话了。
体系是一层Ni原子 ...

四面体的网格应该是k空间的。金属系统（甚至其他某些磁性系统）建议还是采用展宽的方法比较好，四面体方法约束太多了，比如不可约k点大于4等。

作者
Author: lycheeho 时间: 2024-3-27 14:23
本帖最后由 lycheeho 于 2024-3-27 14:24 编辑

卡开发发发表于 2017-6-6 11:43
四面体的网格应该是k空间的。金属系统（甚至其他某些磁性系统）建议还是采用展宽的方法比较好，四面体方 ...

老师您好！我用VASP算孤立分子的单点能，分子基态是开壳层四重态，但我用ISPIN=1，相当于考虑其为闭壳层的情况，结果在EIGENVAL文件中发现占据情况如下，请问这样合理吗？

作者
Author: 卡开发发 时间: 2024-3-27 14:29

lycheeho 发表于 2024-3-27 14:23
老师您好！我用VASP算孤立分子的单点能，分子基态是开壳层四重态，但我用ISPIN=1，相当于考虑其为闭壳层 ...

不好说，这里可能会存在两种情况:
1、一般局域/半局域泛函（如LDA/GGA/MetaGGA）是否会低估多重度的情况，如果低估，这样可能会导致即便正确猜测自旋多重度也会算出来闭壳层的可能，那么就要考虑方法的选择是否合理。
2、局域/半局域泛函也能定性正确，那么这种情况需要对原子的自旋密度进行正确的初猜（如使用MAGMOM或使用NUPDOWN），同时使用ISPIN=2，这样才可能保证计算合理。

作者
Author: lycheeho 时间: 2024-3-27 14:42

卡开发发发表于 2024-3-27 14:29
不好说，这里可能会存在两种情况:
1、一般局域/半局域泛函（如LDA/GGA/MetaGGA）是否会低估多重度的情况 ...

谢谢您！我开ISPIN=2可以正常算出正常的四重态基态，电子填充如下图，我现在就是想人为控制其收敛成闭壳层结构，或者说分子失去3个自由基，所以就直接尝试的不考虑自旋的ISPIN=1，请问这么做合理吗？或者有没有其他更好的办法？谢谢！

作者
Author: 卡开发发 时间: 2024-3-27 15:19

lycheeho 发表于 2024-3-27 14:42
谢谢您！我开ISPIN=2可以正常算出正常的四重态基态，电子填充如下图，我现在就是想人为控制其收敛成闭壳 ...

分数占据的情况你可以试试看把SIGMA设置小一些再用ISPIN=1试试，设置ISPIN=1没有什么问题。

作者
Author: lycheeho 时间: 2024-3-27 16:03

卡开发发发表于 2024-3-27 15:19
分数占据的情况你可以试试看把SIGMA设置小一些再用ISPIN=1试试，设置ISPIN=1没有什么问题。

好的，谢谢您！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)