计算化学公社

标题: 求助PSI4中怎么引用aug-cc-pvtz-pp基组 [打印本页]

作者
Author: luckylife 时间: 2022-3-3 14:55
标题: 求助PSI4中怎么引用aug-cc-pvtz-pp基组
老师，我想在PSI4中引用aug-cc-pvtz-pp基组，因为我用高斯结构优化以及能量计算时用的是此基组，为了保持一致，所以想从外部引用此基组。以下是我的输入文件和输出文件。但是PSI4计算结束后看输出文件，发现输出文件上面Br原子写的基组是aug-cc-pvtz-pp，继续下滑往下看，发现Br原子写的基组又变成了def2-qzvpp，想知道这是为什么以及怎么修改。

作者
Author: hebrewsnabla 时间: 2022-3-3 14:59
本帖最后由 hebrewsnabla 于 2022-3-3 15:04 编辑

那是jkfit辅助基组。你对Br自定义了基组，但是没有指定辅助基组，有可能psi4默认就是对于Br借用这个辅助基组。

你可以自己指定辅助基组，或者，既然你要和高斯一致的话，应该关掉jkfit。

作者
Author: luckylife 时间: 2022-3-3 15:10

hebrewsnabla 发表于 2022-3-3 14:59
那是jkfit辅助基组。你对Br自定义了基组，但是没有指定辅助基组，有可能psi4默认就是对于Br借用这个辅助基 ...

老师，那请问我应该怎么关掉jkfit呢，需要输入什么关键词呢？

作者
Author: luckylife 时间: 2022-3-3 15:12

luckylife 发表于 2022-3-3 15:10
老师，那请问我应该怎么关掉jkfit呢，需要输入什么关键词呢？

或者说怎么指定辅助基组呢，我看了PSI4的相关说明，实在是有些搞不明白了，不知道该怎么操作。

作者
Author: sobereva 时间: 2022-3-3 22:42
对于SAPT目的还不如直接用def2-QZVP呢，省得你折腾

作者
Author: luckylife 时间: 2022-3-4 09:36

sobereva 发表于 2022-3-3 22:42
对于SAPT目的还不如直接用def2-QZVP呢，省得你折腾

老师，那可以H原子、C原子用aug-cc-pvtz，Br原子采用def2-QZVP吗？还是说它们全用def2-QZVP呢？

作者
Author: 喵星大佬 时间: 2022-3-4 09:47

luckylife 发表于 2022-3-4 09:36
老师，那可以H原子、C原子用aug-cc-pvtz，Br原子采用def2-QZVP吗？还是说它们全用def2-QZVP呢？

不开DF可以

作者
Author: luckylife 时间: 2022-3-4 15:47

喵星大佬发表于 2022-3-4 09:47
不开DF可以

老师，那除了DF这个关键词，还可以选用什么SCF-type模式呢？

作者
Author: sobereva 时间: 2022-3-4 19:01

luckylife 发表于 2022-3-4 09:36
老师，那可以H原子、C原子用aug-cc-pvtz，Br原子采用def2-QZVP吗？还是说它们全用def2-QZVP呢？

全用

作者
Author: luckylife 时间: 2022-3-5 15:46

sobereva 发表于 2022-3-4 19:01
全用

好的，谢谢老师

作者
Author: 星斗如盘 时间: 2022-9-25 18:07

sobereva 发表于 2022-3-3 22:42
对于SAPT目的还不如直接用def2-QZVP呢，省得你折腾

请教一下，结构优化使用的是6-31g(d)，EDA可以使用def2-QZVP吗？

作者
Author: sobereva 时间: 2022-9-26 05:25

星斗如盘发表于 2022-9-25 18:07
请教一下，结构优化使用的是6-31g(d)，EDA可以使用def2-QZVP吗？

如果是弱相互作用体系，优化用6-31G(d)太low，这类体系至少6-311G**或def-TZVP档次优化

作者
Author: 星斗如盘 时间: 2022-9-26 12:53

sobereva 发表于 2022-9-26 05:25
如果是弱相互作用体系，优化用6-31G(d)太low，这类体系至少6-311G**或def-TZVP档次优化

不是弱相互体系，是想研究对过渡态进行EDA分析，请问100个原子级别，想做SAPT能量分解，需要什么级别的方法和基组？
目前使用SAPT0/jun-cc-pVDZ（铜）级别进行计算，发现jun-cc-pVDZ不支持镍，而SAPT2+(3)δMP2/aug-cc-pVTZ（金）级别算不动，结构优化是b3lyp-D3/6-31g(d)(SDD for Ni进行)水平下进行的

作者
Author: sobereva 时间: 2022-9-27 06:00

星斗如盘发表于 2022-9-26 12:53
不是弱相互体系，是想研究对过渡态进行EDA分析，请问100个原子级别，想做SAPT能量分解，需要什么级别的方 ...

SAPT只适合弱相互作用范畴，根本没法用于片段划分涉及截断化学键的情况

作者
Author: zjxitcc 时间: 2022-9-27 12:49

星斗如盘发表于 2022-9-26 12:53
不是弱相互体系，是想研究对过渡态进行EDA分析，请问100个原子级别，想做SAPT能量分解，需要什么级别的方 ...

这种情况适合用GKS-EDA ，见《GKS-EDA计算简介》

作者
Author: PESPES 时间: 2023-4-13 17:19

luckylife 发表于 2022-3-3 15:10
老师，那请问我应该怎么关掉jkfit呢，需要输入什么关键词呢？

同学，你会在psi4里用赝势基组了么？还是用的全电子基组？

作者
Author: PESPES 时间: 2023-4-14 19:57

luckylife 发表于 2022-3-3 15:10
老师，那请问我应该怎么关掉jkfit呢，需要输入什么关键词呢？

scf_type pk 可以试试

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