|
|
本帖最后由 ykr 于 2022-5-12 10:30 编辑
我的体系含Cr、O、Ce、Zr、C、H等元素，结构主要是ZrO2（无磁性）。设置了Cr的U为3.5，Ce的U为4.5。其他元素没有设U。Ce、Cr设置初始磁矩都为3。其他所有原子，包括O的磁矩都是0.
计算中发现，同结构同初始磁矩输入文件，做结构优化，第一个离子步得到总磁矩mag=2，有可能计算后磁矩也在2左右。为方便区别，我这里记这样的结果为“首次计算”。
但是，如果用如200步优化后的结构作为输入结构，重新开始计算，我记为“第二次计算”，那么第一个离子步得到的磁矩变了，mag=6或者-1。而且“第二次计算”的能量和“首次计算”的能量差距较大，我体系120个原子左右，能量最多可以差几个eV。这样的话，我做催化计算的话就有奇怪的结果发生，如大于3 eV的能垒，反应能大于活化能。
在这种情况下，该取哪个结果？还是说所有的结果都不可信？
实际上我感觉如果对于其他元素，也有可能会有这种问题。我感觉这种问题来源可能出在U值和spin上。但文献中说，为了避免结果能量的差异，都设置为统一的初始磁矩，这样结果有可比性。但我也有测试，改变初始磁矩，仍然是“首次计算”和“第二次计算”有差异。难道要对每个原子分别设置磁矩，比如我做的表面的模型，表面原子的磁矩可能高，体相原子的磁矩低甚至0？这样也太麻烦了。
|
|
发表于 Post on 2022-5-12 04:07:01
收藏 Add to favorites