计算化学公社

标题: DFT+U计算Fe2O3中O的PDOS曲线出现向左平移情况 [打印本页]

作者
Author: nusiew 时间: 2022-10-1 16:38
标题: DFT+U计算Fe2O3中O的PDOS曲线出现向左平移情况
本帖最后由 nusiew 于 2022-10-1 19:06 编辑

大家好，如下是我采用cp2k计算α-Fe2O3中O的PDOS，分别采用了不同的U值（Multiwfn默认设置+3.9eV，文献中U值一般设为5eV，计算得到的带隙大概在2eV），计算方法总体为smear+UKS+MAGNETIZATION，2*2*2扩胞，Γ点计算，Fe的磁矩分别设为+3和-3（α-Fe2O3磁序为AFM，总磁矩为0）。
（1）U值设为5eV时我计算的带隙总体与文献一致，但是为啥我计算出来的PDOS曲线横坐标在-1eV~1eV之间为0，而文献中基本都是横坐标在0~2eV之间PDOS为0呢？是不是我的inp文件有问题？
（2）当U取3.9eV和5eV时，SCF自洽计算出现警告（ *** WARNING in dft_plus_u.F:1458 :: DFT+U energy contibution is negative *** *** possibly due to unphysical Mulliken charges! ***），U取0.1eV和2eV则没有该警告，请问这种情况的原因以及解决方案。

(, 下载次数 Times of downloads: 2)

相关问题总结：
求助：CP2K中的DFT+U
求助：CP2K用DFT+U出现WARNING DFT+U energy contibution is negative possibl
Fe2O3 计算DOS 对称性问题请教

作者
Author: sobereva 时间: 2022-10-2 04:32
你的展宽设得略宽
你的结果和文献2也差不多
有些文献虽然标的是相对于Ef的能量，但Ef取的是不对的或者有误导性的。文献1把价带顶直接当Ef，这是有问题的。

作者
Author: sigma-jlu 时间: 2023-4-23 20:17
本帖最后由 sigma-jlu 于 2023-4-23 20:20 编辑

请问您“DFT+U energy contibution is negative”报错消失，与跟文献中的pdos接近这个无法兼顾的矛盾后来解决没有？我刚才运行了您给的inp，的确也很耗时。

作者
Author: sigma-jlu 时间: 2023-4-23 20:37
magnetization为什么是-3 ？
因为每个Fe 3+的d轨道被半填充，spin应该是5 uB（中子散射实验结果约4.9uB）。

作者
Author: sobereva 时间: 2023-4-24 03:58

sigma-jlu 发表于 2023-4-23 20:37
magnetization为什么是-3 ？
因为每个Fe 3+的d轨道被半填充，spin应该是5 uB（中子散射实验结果约4.9uB） ...

只有&BS才能控制初猜波函数中原子的净电荷，用magnetization设置初猜的话原子是中性的，照着Fe3+来设没有意义
magnetization的值并不需要特别讲究，只不过会影响SCF收敛性，能收敛到恰当的态即可，不收敛的时候可以再调。

作者
Author: sobereva 时间: 2023-4-24 03:59

sigma-jlu 发表于 2023-4-23 20:17
请问您“DFT+U energy contibution is negative”报错消失，与跟文献中的pdos接近这个无法兼顾的矛盾后来解 ...

前面我说了，别盲目拿文献说事，很多文献里Ef都取得不对。提问前先看别人的回帖

作者
Author: sigma-jlu 时间: 2023-4-24 09:30

sobereva 发表于 2023-4-24 03:59
前面我说了，别盲目拿文献说事，很多文献里Ef都取得不对。提问前先看别人的回帖

多谢提醒。我以后注意，看完所有回帖，仍有疑问或补充的再发帖。

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