计算化学公社

标题: 求助，pp.x计算的问题 [打印本页]

作者
Author: TJU_RYC 时间: 2024-5-11 16:30
标题: 求助，pp.x计算的问题
本帖最后由 TJU_RYC 于 2024-5-13 15:37 编辑

  求助各位老师同学，pp.x计算输入文件 (图a)；计算后输出里只读取了两个赝势文件，还提示找不到原子 (图b)，但是文件夹里三个赝势文件都是有的 (图c)
  最后能计算出cube文件用VESTA打开发现原子不对：原本S的原子变成Co，原本Co的原子变成Li (图d) ，图e为计算结果
  想问问大家；我是哪一步计算出错了？这样计算出的结果对吗？在VESTA里把原子都改正确的原子可以吗？

问题终结：属于是版本问题（qe7.3，已知两处错误：1.生成的pslibrary赝势文件不能用；2.pp.x无法正确读取赝势文件）

感谢卡开发发老师

作者
Author: 卡开发发 时间: 2024-5-11 20:18
问题不一定出在pp.x，要不然先看看pw.x的输入，看看在赝势定义部分是否有不正常的地方？

作者
Author: TJU_RYC 时间: 2024-5-11 22:24

卡开发发发表于 2024-5-11 20:18
问题不一定出在pp.x，要不然先看看pw.x的输入，看看在赝势定义部分是否有不正常的地方？

老师，您看我上传了pw.x输入的赝势文件部分的定义，而且在截图左边也是通过scf计算出的文件，所以很纳闷是不是qe7.3的pp.x有问题

作者
Author: 卡开发发 时间: 2024-5-11 22:46
本帖最后由卡开发发于 2024-5-11 22:47 编辑

TJU_RYC 发表于 2024-5-11 22:24
老师，您看我上传了pw.x输入的赝势文件部分的定义，而且在截图左边也是通过scf计算出的文件，所以很纳闷 ...

输入文件里面一些东西肉眼一下子很难看出问题来。但是确实不好说，QE的IO部分（甚至包括别的）不同版本改的都比较激进，例如7.x版本我印象当中之前某些赝势可能就因为增设了什么字段没办法使用，确实我也很难说存在什么bug。不然换个版本试试看？

作者
Author: TJU_RYC 时间: 2024-5-11 23:14

卡开发发发表于 2024-5-11 22:46
输入文件里面一些东西肉眼一下子很难看出问题来。但是确实不好说，QE的IO部分（甚至包括别的）不同版本改 ...

qe编译太难了，我刚编译7.2想用7.2的pp.x试一下可惜编译失败😂，老师您有qe吗要不分享一个pp.x模块给我看看行不行，在qe根目录，pp/src路径里面，不是qe/bin那个，那个pp.x好像是快捷方式

作者
Author: 卡开发发 时间: 2024-5-12 08:28

TJU_RYC 发表于 2024-5-11 23:14
qe编译太难了，我刚编译7.2想用7.2的pp.x试一下可惜编译失败😂，老师您有qe吗要不分享一个pp.x模 ...

这个编译环境不一样我分享给你可能作用也不大，按道理7.2应该也不是很难编译吧。

作者
Author: TJU_RYC 时间: 2024-5-12 10:07

卡开发发发表于 2024-5-12 08:28
这个编译环境不一样我分享给你可能作用也不大，按道理7.2应该也不是很难编译吧。

想试试行不行，主要是报错了，makeall报错了，改装的库明明都装好了的，就很奇怪

作者
Author: TJU_RYC 时间: 2024-5-12 10:09

卡开发发发表于 2024-5-12 08:28
这个编译环境不一样我分享给你可能作用也不大，按道理7.2应该也不是很难编译吧。

还有就是，老师，如果赝势没有读取正确，是不是pp.x计算了电荷密度不能用呀

作者
Author: 卡开发发 时间: 2024-5-12 14:53

TJU_RYC 发表于 2024-5-12 10:07
想试试行不行，主要是报错了，makeall报错了，改装的库明明都装好了的，就很奇怪

这个问题不好说，要不然有条件我远程看看，或者input发过来我这边测试下。
如果波函数/电荷密度是对的，应该也能用，但这个需要做一些检查才能确定。

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