U参数如何确定？

lyhnh

卡开发发 发表于 2020-6-3 23:12
没那么简单，印象中好像是利用1/X0~Uin这样的关系来加速迭代，如果每个α≠0点的响应只要bare的计算而不 ...

谢谢

卡开发发

lyhnh 发表于 2020-6-4 00:07
谢谢

我没记错的话，方式也变了，老方法是外推，新方法似乎是用某种迭代关系。

zhukangg

这个U参数的测定适用于f电子吗？我测量Pr元素的时候，发现占据数随α的变化并不是线性的。（可能是我对占据数的定义不同，我是将f电子的态密度积分到费米能级得到的占据数），请老师指教

卡开发发

zhukangg 发表于 2021-1-22 16:50
这个U参数的测定适用于f电子吗？我测量Pr元素的时候，发现占据数随α的变化并不是线性的。（可能是我对占据 ...

适合。指的是能带在U的局域投影态上的占据数，通常进行DFT+U计算应该会print出来。PDOS中的可能会受到WS半径、展宽等各方面影响。

zhukangg

卡开发发 发表于 2021-1-23 07:14
适合。指的是能带在U的局域投影态上的占据数，通常进行DFT+U计算应该会print出来。PDOS中的可能会受到WS ...

那请问这个占据数要在OUTCAR和DOSCAR中grep哪个关键词呢，我没有找到。谢谢老师

卡开发发

zhukangg 发表于 2021-1-23 18:31
那请问这个占据数要在OUTCAR和DOSCAR中grep哪个关键词呢，我没有找到。谢谢老师

增加LDAUPRINT=1或2看看，默认那个占据矩阵好像是不打印出来的（也许LDAUTYPE=3的时候不见得能打印，这点我不太确定）。按道理其实用LORBIT=11得到的电荷也行，前面几页给的链接是采用这种方法算的。如果线性程度不太好，那么可以减小alpha的取值范围试试。

另外，QE其实也使用了更好的方法，直接使用DFPT方法进行扰动，也就不必采用超胞方法了。具体可以参考hp.x的使用，不过我还没有测试过，有时间测试可以交流下。

文申的量子化学

卡老师，您好，问一下关于U值的设置是单个原子加U还是，整个体系中的稀土金属都加U了