计算化学公社

标题: orca计算单点能时使用CCSD(T）-F12，如何给cc-pvdz-F12基组添加aug弥散函数？ [打印本页]

作者
Author: Luffy-D-M 时间: 2022-2-10 18:12
标题: orca计算单点能时使用CCSD(T）-F12，如何给cc-pvdz-F12基组添加aug弥散函数？
各位老师好！
我使用orca计算含C、H、O14个原子的体系的单点能，体系中有一个阴离子，使用Multiwfn自动生成orca的输入文件，采用了CCSD(T)-F12/cc-pvdz-F12的条件，我优化的时候用的m062/aug-cc-pvtz,担心单点能没有采用aug而被审稿人指出来，请问如何给CCSD(T)-F12/cc-pvdz-F12添加aug？我尝试加了aug，但是计算报错说是基组和方法不匹配。我查了一下orca手册，CCSD(T)-F12方法没有关于aug的基组，还有阴离子体系如果不加aug计算单点能的话结果是不是就很不合理？望各位老师指点。

作者
Author: wzkchem5 时间: 2022-2-10 18:35
F12系列基组本身就比较弥散，一般不需要加aug。你可以把cc-pVDZ-F12的基组定义打出来，看看它最弥散的s函数和p函数的exponent是不是和aug-cc-pVDZ的exponent差不多，如果前者不明显比后者大，就证明cc-pVDZ-F12至少和aug-cc-pVDZ差不多弥散

作者
Author: Luffy-D-M 时间: 2022-2-10 18:51

wzkchem5 发表于 2022-2-10 18:35
F12系列基组本身就比较弥散，一般不需要加aug。你可以把cc-pVDZ-F12的基组定义打出来，看看它最弥散的s函数 ...

老师，这个基组的定义怎么打出来，是需要在orca的输入文件里面写入特定的关键词吗？还是直接在正常结束的orca输出文件*.out里面可以找到？谢谢！！

作者
Author: wzkchem5 时间: 2022-2-10 19:50

Luffy-D-M 发表于 2022-2-10 11:51
老师，这个基组的定义怎么打出来，是需要在orca的输入文件里面写入特定的关键词吗？还是直接在正常结束的 ...

需要写 %output print[p_basis] 2 end。但是不需要把整个任务跑完，跑几秒钟就打印出来了，然后把任务杀掉就好

作者
Author: Luffy-D-M 时间: 2022-2-10 19:53

wzkchem5 发表于 2022-2-10 19:50
需要写 %output print 2 end。但是不需要把整个任务跑完，跑几秒钟就打印出来了，然后把任务杀掉就好

好的老师，那就是CCSD(T)-F12/cc-pvdz-F12和CCSD(T)/aug-cc-pvdz分别都跑同一个任务是吗？

作者
Author: Luffy-D-M 时间: 2022-2-10 20:12

wzkchem5 发表于 2022-2-10 19:50
需要写 %output print 2 end。但是不需要把整个任务跑完，跑几秒钟就打印出来了，然后把任务杀掉就好

老师，我算完了，但是我不知道该怎么看，您能帮忙看一下吗？这是输出文件 (, 下载次数 Times of downloads: 15) ，

作者
Author: wzkchem5 时间: 2022-2-10 23:10

Luffy-D-M 发表于 2022-2-10 13:12
老师，我算完了，但是我不知道该怎么看，您能帮忙看一下吗？这是输出文件，

看所有的s轨道里面primitive exponent最小的那个。如果不知道哪列是exponent，看orca说明书"Assigning or Adding Basis Functions to an Element"章节

作者
Author: Luffy-D-M 时间: 2022-2-10 23:35

wzkchem5 发表于 2022-2-10 23:10
看所有的s轨道里面primitive exponent最小的那个。如果不知道哪列是exponent，看orca说明书"Assigning or ...

谢谢老师，我看了orca手册，我看F12的s,p的指数比aug的要大，但是好像大的不是很多，我不知道这个相差多少才大致一样的标准，还得麻烦您指点一下！我附上O原子的基组定义的截图,是比较这两组数据吗。 (, 下载次数 Times of downloads: 13) (, 下载次数 Times of downloads: 15)

作者
Author: wzkchem5 时间: 2022-2-11 02:10

Luffy-D-M 发表于 2022-2-10 16:35
谢谢老师，我看了orca手册，我看F12的s,p的指数比aug的要大，但是好像大的不是很多，我不知道这个相差多 ...

是看最小的，不是看最大的。所以是比较0.07376和0.07896。p函数也可以类似地比较一下最小的exponent。exponent越小越弥散。
所以可以看到对于C、O，两个基组都差不多弥散，只有对于H是aug-cc-pVDZ更弥散一点。所以如果你的氢原子没有明显负电荷，cc-pVDZ-F12可以不用加弥散

作者
Author: Luffy-D-M 时间: 2022-2-11 22:34

wzkchem5 发表于 2022-2-11 02:10
是看最小的，不是看最大的。所以是比较0.07376和0.07896。p函数也可以类似地比较一下最小的exponent。exp ...

谢谢老师，我明白了

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