关于LDA+U方法计算氧化物生成能的问题

jiangning198511

对于金属氧化物而言，一般使用LDA+U方法可以在不大幅增加计算量的情况下得到比较合理的电子结构，但对于计算反应能而言，就有个问题，比如计算氧化铁的生成能(E=E(Fe2O3)-2E(Fe)-1.5E(O2)), 这里面Fe2O3的能量计算中需要使用LDA+U，而Fe的bulk计算中是否也应该使用LDA+U？希望大佬可以解答一下@卡开发发

卡开发发

确实有人说不同体系相同元素只有用相同U能够比较，但这也不是严格说法。从（http://hjkgrp.mit.edu/content/in ... -u-variations-dftur）可以看到，U是一个和体系核坐标相关的值，即便构成体系的元素一样，不同的排列、距离其实也就是不同的化学环境，U值会不一样。因此涉及到化学环境剧烈变化的体系，我没办法说是否非得一定DFT+U，这问题具有很强的争议性；反之，化学环境变化很小的情况，固定某个元素U值可以作为一种近似方案。

jiangning198511

卡开发发 发表于 2020-1-2 13:25
确实有人说不同体系相同元素只有用相同U能够比较，但这也不是严格说法。从（http://hjkgrp.mit.edu/content ...

确实，但不用相同的U值得话，能量差就会有U值差异带来的影响，就好像A用了LDA, B用了GGA，A-B的结果没有意义，在量化计算中，能量差在同一个泛函下比较才有意义

卡开发发

jiangning198511 发表于 2020-1-2 13:31
确实，但不用相同的U值得话，能量差就会有U值差异带来的影响，就好像A用了LDA, B用了GGA，A-B的结果没有 ...

是。但化学环境变化剧烈的情况，采用相同U只也不见得合理，这也是前面回复的内容。在Cococcioni他们做自洽线性响应的文章，好像是按照某种取平均的方式处理的吧。