|
|
本帖最后由 Shana 于 2023-7-30 23:05 编辑
各位老师好：
小弟最近看到一篇文献（https://doi.org/10.1021/acs.chemmater.2c03051），使用quantum espresso对Eu2+的d电子和f电子同时加U。小弟之前一直使用vasp，而在vasp中只能对一个轨道加U，因此对这种方法十分感兴趣。论文中的计算方法如参考文献图片所示：
但在实际操作中，小弟却遇到了两个问题：
1. 文献中所使用的Eu的UPF文件，是基于PAW。小弟从引文48中找到对应的赝势文件，在QE中使用时，却会提示如下的错误：Error in routine offset_atom_wfc (1): The pseudo for Gd does not contain labels for atomic orbitals! Please add them by hand in the pseudo.参见报错信息1图片。
小弟的输入文件设置参见输入文件图片，参照的是qe的官方文件。
Eu的PAW文件参见Eu-PAW.UPF图片，是使用atompaw所生成的。我比较了一下超软赝势的文件（参见Eu-SSSP.UPF图片），发现在开头部分，超软赝势给出了电子轨道信息，而PAW没有给出这一信息，这是否是计算报错的原因？
小弟应该如何修改PAW的赝势文件呢？
2. 小弟将Eu的UPF文件换成SSSP的赝势文件后，由于没有f电子，小弟便模仿文献，尝试对未占据的5d轨道加U，小弟一开始设置Hubbard_occ(3,1) = 5，得到如下的报错信息Hubbard manifold with zero occupations is not allowed。小弟将参数改为0以及注释掉这一行之后，都会得到相同的报错信息。参见报错信息2图片。
小弟疑惑的是，文献中提到能够对unoccpuied 5d orbital 加U，那么是不是小弟少了相关的参数设置？或者应该修改UPF文件？
请各位老师指教。
|
|
发表于 Post on 2023-7-30 22:55:14
收藏 Add to favorites
楼主 Author