计算化学公社

标题: 高斯混合基组问题：如何给特定原子指定基组 [打印本页]

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

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

提前感谢！

作者
Author: sobereva 时间: 2019-6-16 06:18
非常简单。在gview里把要用高级别基组的区域选成黄色，进atom selection界面把序号拷出来放到gjf定义基组部分就完了。反选之后其它部分定义成小基组

作者
Author: lnxj99 时间: 2019-6-17 07:10

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

已经解决，谢谢您。

作者
Author: houxudong2020 时间: 2020-11-5 20:31

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

作者
Author: sobereva 时间: 2020-11-5 21:38

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

307-308 前头甭空行

作者
Author: houxudong2020 时间: 2020-11-6 10:49

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

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

作者
Author: sobereva 时间: 2020-11-8 19:21

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

上传gjf文件

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

作者
Author: ZHQ 时间: 2022-6-14 15:58
金属是不是还得考虑雁势

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