线性响应计算U值问题求助

鱼丸box

卡开发发 发表于 2023-3-18 08:46
简单说，就是把哪些原子作为Hubbard位点，然后考虑何种模型（+U还是+U+V）取决于进行hp.x前pw.x的计算参 ...

谢谢卡老师回复，但是现在的我依然有两个疑问：
（1）关于老师回答中的"响应得到的参数是什么你再按照格式填回pw.x当中对应的参数进行后续计算即可"这句话，按照QE中的例10（LiCoO2）的输入文件在Co的3d轨道上+U，以及Co-3d对O-2p轨道作用上+V，为什么在LiCoO2.Hubbard_parameters.dat文件中会得出O原子的U值呢？也就是说，我需要在后续计算中也对O原子+U吗？
（2）如果只针对+U的话，我们进行线性响应计算U值，这个U值指的是Ueff吗？如果是，我们又该在哪里寻找J值呢？

卡开发发

鱼丸box 发表于 2023-3-18 14:42
谢谢卡老师回复，但是现在的我依然有两个疑问：
（1）关于老师回答中的"响应得到的参数是什么你再按照格 ...

1、如果你要进行+U+V的计算，U应该参考LiCoO2.Hubbard_parameters.dat，V应该参考Hubbard.dat的数据。因为V当中设置了O-2p，那么U也会被算出来。
2、如果你按照DFT+U（也就是旋转不变形式）进行设置，此时U就是Ueff=U-J，但DFT+U+J目前看hp.x程序描述似乎没办法得到单独的J，有可能得通过pw.x在实空间下手动做，具体我不那么确定。

鱼丸box

卡开发发 发表于 2023-3-18 21:11
1、如果你要进行+U+V的计算，U应该参考LiCoO2.Hubbard_parameters.dat，V应该参考Hubbard.dat的数据。因 ...

谢谢卡老师的回答，对于问题（1）我会先进行一些测试，对于（2），我的理解是pwscf的DFT+U相当于VASP中的LDAUTYPE=2，并设置了J=0。

卡开发发

鱼丸box 发表于 2023-3-19 15:27
谢谢卡老师的回答，对于问题（1）我会先进行一些测试，对于（2），我的理解是pwscf的DFT+U相当于VASP中的 ...

确实是这样。

Tang-Tommy

卡开发发 发表于 2023-3-18 08:46
简单说，就是把哪些原子作为Hubbard位点，然后考虑何种模型（+U还是+U+V）取决于进行hp.x前pw.x的计算参 ...

卡老师，请问一下，我使用Vasp对一个尖晶石体系进行计算，Zn、Fe、O元素。用DFT+U进行结构优化，然而结构优化提交任务后，每次得到的能量都不一样，这个是正常的吗？

卡开发发

Tang-Tommy 发表于 2023-10-22 18:10
卡老师，请问一下，我使用Vasp对一个尖晶石体系进行计算，Zn、Fe、O元素。用DFT+U进行结构优化，然而结构 ...

两次计算的OUTCAR和OSZICAR或STDOUT（也就是屏幕输出内容）发来看看，暂时不好说。特定的U值可能确实会导致计算不稳定。

Tang-Tommy

卡开发发 发表于 2023-10-23 00:04
两次计算的OUTCAR和OSZICAR或STDOUT（也就是屏幕输出内容）发来看看，暂时不好说。特定的U值可能确实会导 ...

请问一下老师，怎么把这个文件发给您过目呢？

鱼丸box

博主您好，请问最后您解决了这个问题吗？我使用了多位点计算，但是结果依然是5.6-5.9之间