计算化学公社

标题: 求助赝势基组计算与全电子标量相对论计算的W@Au12团簇Mayer键级差异明显问题 [打印本页]

作者
Author: BLL130 时间: 2022-10-5 20:17
标题: 求助赝势基组计算与全电子标量相对论计算的W@Au12团簇Mayer键级差异明显问题
本帖最后由 BLL130 于 2022-10-5 20:19 编辑

本人在近期采用G16 C.02计算闭壳层Ih点群的W@Au12团簇时，发现赝势基组计算与全电子标量相对论计算的Mayer键级结果相差极大，以下是计算结果（泛函均为revTPSS，以W-Au的Mayer键级为例）：
SDD赝势基组：键级为0.05；
DKH2结合SARC基组（DKH2版本）：键级为1.08；
DKH4结合x2c-TZVPall基组：键级为-2.90（这个明显是有问题的）。
其它Mayer键级（例如Au-Au键级）也存在明显差别。

团簇结构采取revTPSS/SDD优化，已振动分析检验无虚频。SDD基组采取Gaussian内置的，其余基组是从BSE上拷贝而来的。三个计算均采用相应级别下的单点计算任务生成的chk文件，使用formchk转化为fch文件后，采用Multiwfn进行Mayer键级的计算。三个计算均经过波函数检测，为稳定波函数。

附三个单点计算的gjf输入文件如下。更多相关文件请见（链接: https://pan.baidu.com/s/1hZtYG_walsAMLUWO-bzI3w 提取码: 14so），内含每个计算的gjf输入文件、log输出文件、fch文件、bndmat键级矩阵文件。

请问为什么会造成这种差别呢？计算的时候采取哪种结果更加符合实际呢？感谢您的解答！

作者
Author: sobereva 时间: 2022-10-5 20:43
建议改用Multiwfn算Wiberg键级，三个波函数文件结果就很一致了，都是0.7几
当前问题属于Mayer键级自身的问题，在极个别情况下可能不合理，之前我研究聚炔烃的时候也遇到Mayer键级存在的明显问题。

作者
Author: BLL130 时间: 2022-10-5 20:48

sobereva 发表于 2022-10-5 20:43
建议改用Multiwfn算Wiberg键级，三个波函数文件结果就很一致了，都是0.7几
当前问题属于Mayer键级自身的问 ...

好的，感谢sob老师的回复！

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