计算化学公社

标题: 对周期性Pt(111)晶面计算Hirshfeld-I原子电荷时提示缺少Ptir.rad文件而无法计算 [打印本页]

作者
Author: electrochemyjy 时间: 2025-11-25 12:39
标题: 对周期性Pt(111)晶面计算Hirshfeld-I原子电荷时提示缺少Ptir.rad文件而无法计算
各位老师好，我用CP2K对三层Pt(111)晶面结构优化后(PBE/DZVP-MOLOPT-SR-GTH)，计算单点能，产生molden文件(下面提供单点计算输入、输出文件及molden文件截图)，然后计划用multiwfn计算这个体系的Hirshfeld-I原子电荷分布。我首先对molden文件添加了盒子与有效核电荷信息，然后导入multiwfn进行原子电荷计算(已将multiwfn的examples目录下的atmrad复制到当前运行路径下)，然而计算时出现了以下提示：

Error: atmrad/Ptir.rad is needed but was not prepared!
Current charge of atom 1(Pt): 59.98382555
Note: This error implies that this atom has unusual charge. You should manually provide the corresponding .rad file in "atmrad" folder prior to the calculation. See Section 3.9.13 of Multiwfn manual for detail.

我检查atmrad中Pt元素相关的rad文件，发现已有Pt_0.rad，Pt+1.rad，Pt+2.rad，Pt+3.rad，Pt-1.rad这几种电荷状态下的rad文件，我的计算体系电荷设置为0，请问Ptir.rad是什么文件呢？另外，在此之前，我用同样的molden文件在multiwfn计算了Hirshfeld原子电荷，想看下CP2K计算的Hirshfeld原子电荷是否合理，然而发现用multiwfn计算产生的结果很奇怪，电荷输出结果是在60左右数值(附件提供输出结果)。在博客《使用Multiwfn对周期性体系计算Hirshfeld(-I)、CM5和MBIS原子电荷》中，我看到正常的电荷输出应该是有正负值，表示得失电子的状态？我这个结果实在是奇怪，想问下各位老师我在哪个步骤可能出了问题，感谢指点！

(, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 0)

(, 下载次数 Times of downloads: 0)

(, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 2)

作者
Author: sobereva 时间: 2025-11-25 14:37
实测你的inp文件极难收敛，这是我改后的
(, 下载次数 Times of downloads: 2)

这是计算得到的molden文件编辑后的开头的截图
(, 下载次数 Times of downloads: 2)

Multiwfn目前最新版本计算，我这里没有任何问题（为了节约时间，作为测试用的偏大的0.4 Bohr格点间距）

(, 下载次数 Times of downloads: 2)

作者
Author: electrochemyjy 时间: 2025-11-25 15:26

sobereva 发表于 2025-11-25 14:37
实测你的inp文件极难收敛，这是我改后的

谢谢老师！

作者
Author: electrochemyjy 时间: 2025-11-25 18:45

electrochemyjy 发表于 2025-11-25 15:26
谢谢老师！

问题已解决，是因为我把molden文件中的关键词写错了，把[Nval]写成了[Naval]，因此没有正确读入有效电荷数，从而引起结果异常。

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