求助：共轭体系以热力学焓值计算解离能及电离势等任务的算法

wzkchem5

Ace2333 发表于 2022-3-3 16:27
我懂了，只有做结构优化时不加弥散，但是计算单点能除了氢以外都要使用6-311+G(d,p)

结构优化最好也加弥散，因为我看你会涉及到酚氧负离子，这个负离子的电荷还是挺集中的。如果是一个特别弥散的负离子，比如石墨烯自由基阴离子这种，还可以考虑不加弥散，但是这种相对比较局域的氧负离子一般还是应该加弥散
另外注意不是“除了氢以外都要使用6-311+G(d,p)”，6-311+G(d,p)已经包含了“只对非氢原子加弥散，氢原子不加弥散”的意思了，两个加号才是对所有原子都加弥散，所以如果要用6-311+G(d,p)，对所有原子都用6-311+G(d,p)就行了。当然，只对除了氢以外的原子用6-311+G(d,p)，结果也是对的，只不过白费了一两分钟写混合基组的输入文件而已，而且如果文章计算部分这么写也容易让审稿人觉得你不懂计算。
此外我的本意是你结构优化用6-311+G(d,p)，单点能用更大的基组，比如jul-cc-pVTZ。单点能用6-311+G(d,p)有点太小了，最多只是勉强可以接受，如果没有特殊原因尽量选大一些

Ace2333

wzkchem5 发表于 2022-3-4 00:06
结构优化最好也加弥散，因为我看你会涉及到酚氧负离子，这个负离子的电荷还是挺集中的。如果是一个特别弥 ...

谢谢老师的建议与指导，我会跟我导师商量后酌情考虑结构优化的方式，非常感谢老师在百忙之中对我的帮助！