计算化学公社

标题: 求助：DFT+U加spin后，同结构同初始磁矩输入却得到不同的磁矩与能量 [打印本页]

作者
Author: ykr 时间: 2022-5-12 04:07
标题: 求助：DFT+U加spin后，同结构同初始磁矩输入却得到不同的磁矩与能量
本帖最后由 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？这样也太麻烦了。

作者
Author: 文申的量子化学 时间: 2022-8-7 18:42
同学，你解决类没

作者
Author: yuexiatx 时间: 2022-8-18 15:18
计算每一种mag下的能量，最低者为基态
应当合理利用NUPDOWN

作者
Author: Tang-Tommy 时间: 2023-10-19 10:05
您好，请问您解决了吗？

作者
Author: Tang-Tommy 时间: 2023-10-19 10:07

yuexiatx 发表于 2022-8-18 15:18
计算每一种mag下的能量，最低者为基态
应当合理利用NUPDOWN

舍友，你也遇到过相应的问题了吗

作者
Author: 我静静等待 时间: 2023-10-19 20:37
可以看一下我帖子下面的回复 http://bbs.keinsci.com/thread-34510-1-1.html

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