计算化学公社

标题: 请问计算声子谱时GGA+U算法中的U值需要加上吗？ [打印本页]

作者
Author: 碧雪霜花剑 时间: 2020-11-13 11:27
标题: 请问计算声子谱时GGA+U算法中的U值需要加上吗？
计算声子谱时，结构优化得到的CONTCAR是通过GGA+U算法得到的，那么在计算声子谱时，INCAR中参数U值需要加上吗？

作者
Author: rabbitkiller 时间: 2020-11-13 12:59
(, 下载次数 Times of downloads: 39)
本人新手，我猜不用，算scf的时候，U的效应已经在波函数和密度里了，求声子谱解的不是KS方程（看上面的图片），只要有密度信息就行了。
（上面的图片来自QE的讲义，我猜VASP算声子谱应该解的是差不多的方程吧）

作者
Author: 碧雪霜花剑 时间: 2020-11-13 15:30

rabbitkiller 发表于 2020-11-13 12:59
本人新手，我猜不用，算scf的时候，U的效应已经在波函数和密度里了，求声子谱解的不是KS方程（看上面的图 ...

我也是初学者，虽然还不是很理解，还是非常谢谢你

作者
Author: 卡开发发 时间: 2020-11-15 06:41
本帖最后由卡开发发于 2020-11-15 06:46 编辑

构型优化做了DFT+U算声子频率显然要加，GGA+U的Hamiltonian和GGA是不一样的，这样两者的势能面也很可能不相同。建议计算声子频率（振动）所使用的理论方法和精度都应与构型优化的完全相同。上面的图只能说明用到了密度导数的信息而已，并不能够说明与密度之外的量无关，光看总能的偏导数能够得知是否+U结果会可能有所不同。

作者
Author: 碧雪霜花剑 时间: 2020-11-15 15:41
本帖最后由碧雪霜花剑于 2020-11-15 15:43 编辑

卡开发发发表于 2020-11-15 06:41
构型优化做了DFT+U算声子频率显然要加，GGA+U的Hamiltonian和GGA是不一样的，这样两者的势能面也很可能不相 ...

感谢感谢，除了在做分子谱的计算以外，我还在进行分子动力学的计算，想请教一下计算分子动力学时加U是必要的吗？

作者
Author: 卡开发发 时间: 2020-11-15 16:15

碧雪霜花剑发表于 2020-11-15 15:41
感谢感谢，除了在做分子谱的计算以外，我还在进行分子动力学的计算，想请教一下计算分子动力学时加U是必 ...

如果你希望优化结构和动力学两个过程用的是相同的势能面，那就得和构型优化过程的计算参数相同。

作者
Author: 碧雪霜花剑 时间: 2020-11-15 17:04

卡开发发发表于 2020-11-15 16:15
如果你希望优化结构和动力学两个过程用的是相同的势能面，那就得和构型优化过程的计算参数相同。

感谢

作者
Author: rabbitkiller 时间: 2020-11-16 19:23

卡开发发发表于 2020-11-15 06:41
构型优化做了DFT+U算声子频率显然要加，GGA+U的Hamiltonian和GGA是不一样的，这样两者的势能面也很可能不相 ...

明白啦，谢谢指教！

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