计算化学公社

 找回密码 Forget password
 注册 Register
Views: 5048|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[PSI4] 求助PSI4中怎么引用aug-cc-pvtz-pp基组

[复制链接 Copy URL]

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

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



202203031454153593..png (68.13 KB, 下载次数 Times of downloads: 67)

202203031454153593..png

202203031453333316..png (225.91 KB, 下载次数 Times of downloads: 69)

202203031453333316..png

test.inp

4.2 KB, 下载次数 Times of downloads: 16

test.out

72.24 KB, 下载次数 Times of downloads: 4

910

帖子

1

威望

7887

eV
积分
8817

Level 6 (一方通行)

2#
发表于 Post on 2022-3-3 14:59:18 | 只看该作者 Only view this author
本帖最后由 hebrewsnabla 于 2022-3-3 15:04 编辑

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

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

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-3-3 15:10:06 | 只看该作者 Only view this author
hebrewsnabla 发表于 2022-3-3 14:59
那是jkfit辅助基组。你对Br自定义了基组,但是没有指定辅助基组,有可能psi4默认就是对于Br借用这个辅助基 ...

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

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-3-3 15:12:25 | 只看该作者 Only view this author
luckylife 发表于 2022-3-3 15:10
老师,那请问我应该怎么关掉jkfit呢,需要输入什么关键词呢?

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

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

5#
发表于 Post on 2022-3-3 22:42:34 | 只看该作者 Only view this author
对于SAPT目的还不如直接用def2-QZVP呢,省得你折腾
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-3-4 09:36:15 | 只看该作者 Only view this author
sobereva 发表于 2022-3-3 22:42
对于SAPT目的还不如直接用def2-QZVP呢,省得你折腾

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

1657

帖子

5

威望

4562

eV
积分
6319

Level 6 (一方通行)

喵星人

7#
发表于 Post on 2022-3-4 09:47:14 | 只看该作者 Only view this author
luckylife 发表于 2022-3-4 09:36
老师,那可以H原子、C原子用aug-cc-pvtz,Br原子采用def2-QZVP吗?还是说它们全用def2-QZVP呢?

不开DF可以

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
muuu2333 + 4

查看全部评分 View all ratings

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-3-4 15:47:28 | 只看该作者 Only view this author

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

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

9#
发表于 Post on 2022-3-4 19:01:19 | 只看该作者 Only view this author
luckylife 发表于 2022-3-4 09:36
老师,那可以H原子、C原子用aug-cc-pvtz,Br原子采用def2-QZVP吗?还是说它们全用def2-QZVP呢?

全用
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2022-3-5 15:46:30 | 只看该作者 Only view this author

好的,谢谢老师

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

11#
发表于 Post on 2022-9-25 18:07:04 | 只看该作者 Only view this author
sobereva 发表于 2022-3-3 22:42
对于SAPT目的还不如直接用def2-QZVP呢,省得你折腾

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

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

12#
发表于 Post on 2022-9-26 05:25:22 | 只看该作者 Only view this author
星斗如盘 发表于 2022-9-25 18:07
请教一下,结构优化使用的是6-31g(d),EDA可以使用def2-QZVP吗?

如果是弱相互作用体系,优化用6-31G(d)太low,这类体系至少6-311G**或def-TZVP档次优化
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

13#
发表于 Post on 2022-9-26 12:53:33 | 只看该作者 Only view this author
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进行)水平下进行的

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

14#
发表于 Post on 2022-9-27 06:00:54 | 只看该作者 Only view this author
星斗如盘 发表于 2022-9-26 12:53
不是弱相互体系,是想研究对过渡态进行EDA分析,请问100个原子级别,想做SAPT能量分解,需要什么级别的方 ...

SAPT只适合弱相互作用范畴,根本没法用于片段划分涉及截断化学键的情况
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4114

帖子

4

威望

8875

eV
积分
13069

Level 6 (一方通行)

MOKIT开发者

15#
发表于 Post on 2022-9-27 12:49:57 | 只看该作者 Only view this author
星斗如盘 发表于 2022-9-26 12:53
不是弱相互体系,是想研究对过渡态进行EDA分析,请问100个原子级别,想做SAPT能量分解,需要什么级别的方 ...

这种情况适合用GKS-EDA ,见《GKS-EDA计算简介
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 03:02 , Processed in 0.682150 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list