计算化学公社

标题: 请教QE包含动能项的赝势 [打印本页]

作者
Author: ghifi37 时间: 2024-2-14 20:06
标题: 请教QE包含动能项的赝势
想利用QE进行meta-GGA的计算，因其只支持模守恒的赝势，于是就尝试了dojo和sg15，计算时都出现了"kinetic-energy density file not found"。
请问这个动能密度是像VASP一样直接包含在赝势文件中还是一个单独的文件？
找了dojo和sg15都没提及动能项问题，请问有什么解决办法吗？

作者
Author: 卡开发发 时间: 2024-2-15 10:38
按道理不应该吧，这个跳错字段是在read_scf这个soubroutine里面，如果从头开始进行自洽计算应该不会有相关信息，正常输出的样式应该是PW/Example09里面的样子。

这里的kinetic energy density应该不是对应VASP的POTCAR的内容，理论上说应该PAW做metaGGA才必须要core kineic-energy denisty的信息，目前UPF格式根本就没这个字段。

QE应该有人专门拟合了SCAN的模守恒赝势，暂时不能确定一般GGA或者其他metaGGA的赝势直接使用有什么后果，建议仔细测试

。

作者
Author: ghifi37 时间: 2024-2-15 14:55

卡开发发发表于 2024-2-15 10:38
按道理不应该吧，这个跳错字段是在read_scf这个soubroutine里面，如果从头开始进行自洽计算应该不会有相关 ...

哈哈，老哥所言极是！我是为了能加速下计算速度，给了个其它泛函收敛的波函数作初猜。80来个原子，331的k点，给了32核，结果一晚上的输出还是 starting charge xx, renormalised to xx，都没更新。我杀了试试看从头算会不会快点。
主要是SCAN的赝势太少，还一个过渡金属的都没有，不行我就混合着用！谢谢老哥！

作者
Author: ghifi37 时间: 2024-2-15 20:25
从头算还是半天没更新，投降了！
换VASP几分钟就出一个电子步，还是有它流行的道理的！

作者
Author: 卡开发发 时间: 2024-2-15 22:51

ghifi37 发表于 2024-2-15 20:25
从头算还是半天没更新，投降了！
换VASP几分钟就出一个电子步，还是有它流行的道理的！

那个我不确定是跳错结束还是warning，有可能只是warning，不过metaGGA的稳定性本身很不好应该要小心一些。除了mGGA的问题，其实赝势各方面都还得小心，程序要提升体验是不容易的，短的用户流程和算法的鲁棒性里面有不少门道。

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