【求助】关于使用VASP优化臭氧分子的问题

wrtgcn

最近想要做臭氧分子在催化剂表面吸附分解过程的研究，优化臭氧分子的时候碰到了一些问题。臭氧的电子结构比较特殊，为自旋极化单重态。我做臭氧分子优化的时候测试了ISPIN=1（闭壳），ISPIN=2（默认磁矩1,1,1），以及ISPIN=2设置磁矩（-1,0,1）的结果，把盒子设置为不对称形状且采用ISYM=0。发现采用默认磁矩得到了最低能量，其次是闭壳，而（-1,0,1）初始磁矩的能量结果反而是最高的。三个结果最终的总磁矩都是0.
想请问大家，采用默认磁矩得到的结果合理吗？

archer

我认为初始给的只是猜测，要看最后收敛的结果。VASP算这种分子可能会出现收敛问题，得到很离谱的结果。

Accelerator

臭氧的基态是闭壳层单重态。

sobereva

Accelerator 发表于 2019-4-23 23:14
臭氧的基态是闭壳层单重态。

结论和泛函有关系。在实验结构下，常用的B3LYP、M06、M06-2X等结论都是基态是自旋极化单重态，某些纯泛函和HF成分低的泛函给出的基态可能是闭壳层单重态。

wrtgcn

archer 发表于 2019-4-23 22:48
我认为初始给的只是猜测，要看最后收敛的结果。VASP算这种分子可能会出现收敛问题，得到很离谱的结果。

我测试的3个计算最终收敛结果都是磁矩处处为0.结构相差很小，能量的差值也不是很大。选取能量最低的结果您觉得合理吗？

wrtgcn

Accelerator 发表于 2019-4-23 23:14
臭氧的基态是闭壳层单重态。

直接当成闭壳单重态计算就可以吗？我用的方法是PAW_PBE。

archer

wrtgcn 发表于 2019-4-24 09:02
我测试的3个计算最终收敛结果都是磁矩处处为0.结构相差很小，能量的差值也不是很大。选取能量最低的结果 ...

磁矩是不是还有01-1这种初始排布，我觉得是可以用能量最低的，你可以试算下吸附结构。