在CP2K中对Ueff值进行拟合

user_fu

各位老师好，我在CP2K中计算硅酸亚铁表面二价铁原子的U值时，先通过查找文献发现Vasp中对于硅酸亚铁中的Fe原子Ueff是3.6ev，同时我查找资料发现硅酸亚铁中的铁原子在0k时是反铁磁磁性，因此我在硅酸亚铁晶胞中设置了反铁磁性，并在CP2K中设置总体磁性为反铁磁，相关铁原子的磁矩大小相等,方向相反,总磁矩为0,自旋多重度为1，打开UKS，但是开始计算后出现如下报错。通过在群里提问发现，各个计算软件的Ueff值不能通用，只能通过线性方法计算或者通过晶格常数和能带拟合，于是我找到了硅酸亚铁的实验晶格参数，想通过不同Ueff值进行拟合，我初步尝试了Ueff=2ev，发现开始晶胞优化后没有出现warning警告，但是当Ueff设置为3.0ev时，还是会出现Ueff，于是我想可能是合理的Ueff值是在0-3ev之间，但是我不知道这个0-3ev之间的间距该如何确定才能比较合理而不引起审稿人质疑，我也查看了一些文献，发现它们都是Ueff0-5ev进行模拟，然后在其中取更精确的Ueff值进行判定，似乎没有出现大的Ueff发出警告这种情况，因此我想请教各位老师，我是否应该也像文章中这样先在一个大的范围（包括vasp计算出的精确值）中先进行一个粗糙的筛选，然后在小的范围内精确筛选，以及如果在CP2K中计算vasp计算出的精确值附近以及超过时出现warning该如何解决，请各位老师给出宝贵建议。