用6-31G*算C O H N S，但是显示“”Atomic number out of range for 6-31G basis set

wangzh

老师们好，我进行结构优化时，用6-31G*算C O H N S，但是显示“”Atomic number out of range for 6-31G basis set.”老师们能不能帮我看一下问题在哪里吗？感谢老师们

乐平

用6-31G*算C O H N S

你不诚实哦…… 体系里明明有镧系的 Ce，提问里却只写 C O H N S ……

你可以查看 Gaussian 官网的 Basis sets（链接如下）， 6-31G 系列的基组只含有 H -- Kr 这些元素，镧系没有包含在此基组内，当然会报错 “”Atomic number out of range for 6-31G basis set.”
https://gaussian.com/basissets/


对于镧系锕系这类的元素，你需要用赝势基组。请阅读下面的几个帖子

谈谈量子化学中基组的选择  http://sobereva.com/336
谈谈赝势基组的选用 http://sobereva.com/373
在线基组和赝势数据库一览 http://sobereva.com/309
详解Gaussian中混合基组、自定义基组和赝势基组的输入  http://sobereva.com/60

sobereva

除了看楼上说的博文外，另见
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

你的压缩包里的gau.gjf和里面的out文件根本不对应

wangzh

乐平 发表于 2023-5-10 10:52
你不诚实哦…… 体系里明明有镧系的 Ce，提问里却只写 C O H N S ……

你可以查看 Gaussian 官网的 Ba ...

老师您好，感谢您的指导，这里Ce用的SDD,没有问题，所以我只写了6-31G算的，不好意思，应该写全的，我觉得我的输入是没有问题的，但还是报错了，所以想请教一下错在哪里了？

wangzh

sobereva 发表于 2023-5-10 12:23
除了看楼上说的博文外，另见
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（htt ...

老师您好，问题就在这里，我用molclus调用高斯进行优化的，输入文件是这个，但是输出文件也是我上传的这个，确实不对应

sobereva

wangzh 发表于 2023-5-10 15:08
老师您好，问题就在这里，我用molclus调用高斯进行优化的，输入文件是这个，但是输出文件也是我上传的这 ...

模板文件没写对

搞清楚molclus目录下的template.gjf是干嘛的，下文明确说了
使用molclus程序做团簇构型搜索和分子构象搜索
http://bbs.keinsci.com/thread-577-1-1.html