计算化学公社

标题: Multiwfn结合CP2K计算弱相互作用时基组选择问题 [打印本页]

作者
Author: renj864878110 时间: 2023-5-21 18:40
标题: Multiwfn结合CP2K计算弱相互作用时基组选择问题
sob老师，我用CP2K计算的体系中有第6周期元素如Bi，对于计算弱作用较省时的选择是6-31G*，6-311G**对于Bi元素不能用，我该选择哪种基组给Bi呢？如果给它选择DZVP-MOLOPT-SR-GTH, 但是出错了。多谢指点！

作者
Author: sobereva 时间: 2023-5-21 21:10
6-31G*计算弱相互作用是开玩笑

DZVP-MOLOPT-SR-GTH对Bi明明有定义，BASIS_MOLOPT文件里都给明了
Bi DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q5
1
2 0 2 5 2 2 1
   1.666928293168 -0.433606016760 -0.100663867810 -0.092396269342 -0.108190299502 -0.195956091275
   1.317941720436  0.723433454454  0.259612627134  0.184381669404  0.210739707604  0.351618767590
   0.222221383065 -0.715303457579 -0.216478402471 -0.507231634301 -0.568218339722 -1.047717462114
   0.074524438950 -0.308379410009 -1.423645113253 -0.421377693158 -0.757828317889 -1.887710553530
   0.066725427306  0.159245921554  2.308274018558  0.106621540317  1.788226073856  0.443921215549

有报错时始终自觉交代清楚报错

对于计算相互作用能目的可以用BASIS_MOLOPT_UZH里的TZ档次的MOLOPT基组

作者
Author: renj864878110 时间: 2023-5-22 18:17

sobereva 发表于 2023-5-21 21:10
6-31G*计算弱相互作用是开玩笑

DZVP-MOLOPT-SR-GTH对Bi明明有定义，BASIS_MOLOPT文件里都给明了

感谢sob老师，我是想试试用Multiwfn+CP2K分析一个周期性结构的IRI, 因带有MOLOPT基组在我的单机上耗时长，想找个耗时少的基组。我知道6-31G*分析弱作用是不对的！可能是表达错误，多谢指教！

作者
Author: sobereva 时间: 2023-5-23 02:02

renj864878110 发表于 2023-5-22 18:17
感谢sob老师，我是想试试用Multiwfn+CP2K分析一个周期性结构的IRI, 因带有MOLOPT基组在我的单机上耗时长 ...

6-31G*用于Multiwfn做IRI分析目的毫无问题，但是不能拿6-31G*算弱相互作用能。计算弱相互作用大家首先会联想到是算相互作用能

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