夙雪夜归 发表于 2021-7-3 16:33
您看是不是这个意思：
1.对于闭壳层体系，如果不放心ISPIN=1或者2两种情况是否会导致不同的结果，可以进 ...

1、确定是闭壳层的情况，没有特殊原因一般不用开。
2、你的理解是正确的。只开ISPIN=2仅仅对应于MAGMOM=NATOMS*1.0的特殊情况，但并不一定这样就能找到能量最低的情况。

卡开发发 发表于 2021-7-3 16:15
上述测试指的是针对你闭壳层结构自旋密度异常这件事提的。当然，如果确定是否用ispin=1还是2确实可以看下 ...

您看是不是这个意思：
1.对于闭壳层体系，如果不放心ISPIN=1或者2两种情况是否会导致不同的结果，可以进行测试比较一下。但是应该不会有区别或者区别微弱到可以忽略（我猜应该是这样）；
2.对于开壳层体系，包括自由基和带有磁性的体系，都应该在ISPIN=2的情况下进行计算；磁性体系可以在文献中找别人用的磁矩设置方法，或者自行测试得出一合理的结果；而涉及自由基、不是磁性体系的计算也应当猜测不同magmom来看是否对结果会产生影响。

夙雪夜归 发表于 2021-7-3 15:54
1.您说的是对的，vesta里显示的Max和Min确实都是很小的数值。这里的这个数值大小我想应该可以认为是浮点 ...

上述测试指的是针对你闭壳层结构自旋密度异常这件事提的。当然，如果确定是否用ispin=1还是2确实可以看下是否开启，但不要忽略对magmom的猜测。

夙雪夜归 发表于 2021-7-3 15:26
感谢您的回复，下图是我添加了  MAGMOM = 2*0  后重新提交任务得到的结果，结果出现了变化，但感觉还是不 ...

1、其实也正常，因为直接设置ispin=2所有原子全按照1猜的磁矩，其实即便全按照0来猜测，达到设定的阈值之后自旋密度也未必能全部都是严格的0，毕竟浮点数的精度是有限的，你现在可以检查一下vesta中isosurface的值，应该是在一个可忽略的量级。

2、如果有自由基，那么不应该使用ispin=1，不同的磁矩情况下的势能面的形状未必是一样的。如果你担心上述问题，你可以再次测试，计算合理的情况自由基这个单电子对自旋密度的贡献应该远比上述因数值误差引起的大得多。

MAGMOM所有原子全按照0猜试试。本来就没磁性的体系按照ISPIN=1做为好。