高斯混合基组问题：如何给特定原子指定基组

lnxj99

大家好，我最近遇到一个问题：在做一个体系，里面有上百个碳原子，其中只有个别和金属配位，我想减小计算量，只对个别几个用TZVP，其余的用SVP。这个情况该如何去定义呢？
读了卢天老师这个帖子（http://sobereva.com/60）之后，发现G09里面基组是叠加而非覆盖。
从而目前想到这样两个方法：
1)用数字标记挨个给，但是帖子中例子是一个CH4，总共只有5个原子，对于这种有上百个C原子的，如果全部用数字给，会非常长，不清楚高斯读入是否会出现问题；就算不出现问题，也非常麻烦。
2)采用叠加的方法，先给所有碳原子定义成SVP，然后再对特定那几个碳原子增加基函数，使其从SVP变成TZVP。但是我不清楚是否可以这样做，以及具体怎么做。

就此向大家求助，请问是否有更简便的办法？

提前感谢！

ZHQ

金属是不是还得考虑雁势

influenza

卢老师您好，您说的这个atom selection工具，是不是没有办法进行按照原子的种类进行批量选择？比如我想对金属相连的氢使用极化函数，但不想对分子内其它的氢加极化函数，那怎么选？

sobereva

houxudong2020 发表于 2020-11-6 10:49
这样，我也试了，还是不行，它报错在 1-37这里就语法错误了。

上传gjf文件

houxudong2020

sobereva 发表于 2020-11-5 21:38
307-308 前头甭空行

这样，我也试了，还是不行，它报错在 1-37这里就语法错误了。

sobereva

houxudong2020 发表于 2020-11-5 20:31
老师，您好，我也看了您之前的帖子，但是
H                                                   0   -4 ...

307-308 前头甭空行

houxudong2020

sobereva 发表于 2019-6-16 06:18
非常简单。在gview里把要用高级别基组的区域选成黄色，进atom selection界面把序号拷出来放到gjf定义基组部 ...

老师，您好，我也看了您之前的帖子，但是
H                                                   0   -4.80985996   14.60527220   -9.77630283 L
H                                                   0   -6.09031454   15.51020803   -9.77630283 L

1-37
6-31g*
****

307-308
6-31g*
****
我用上面这种方式定义时（我只复制了部分），还是一直出错，我想对体系中的非核心原子用小基组核心原子用大基组，问题出在哪里呢？

lnxj99

sobereva 发表于 2019-6-16 06:18
非常简单。在gview里把要用高级别基组的区域选成黄色，进atom selection界面把序号拷出来放到gjf定义基组部 ...

已经解决，谢谢您。

sobereva

非常简单。在gview里把要用高级别基组的区域选成黄色，进atom selection界面把序号拷出来放到gjf定义基组部分就完了。反选之后其它部分定义成小基组