计算化学公社

标题: 单点计算是否总是需要带着弥散？ [打印本页]

作者
Author: WaterMirror 时间: 2022-3-2 02:49
标题: 单点计算是否总是需要带着弥散？
尊敬的各位老师:
计算稍微大一点体系的过渡态与结构（150个原子），总是会有弱相互作用，请问各位老师，我在泛函充分考虑弱相互作用的情况下（M06-2X-D3），仅仅是计算体系的单点能获得准确的电子能量，需要一定要带着弥散函数吗？比如使用def2-TZVP是否就足够了？还是必须要用ma-def2-TZVP呢？因为给def2加弥散在高斯中有点麻烦，如果必须要弥散可能就会考虑用6-311+G**或者6-311++G**了。

是偶然看到Sob老师给别人解答的一个问题，没有明白Sob老师的具体意思，有点懵了。请问各位老师Sob老师这句话的意思是只要计算带有弱相互作用体系的单点能，就必须加弥散(def2-TZVP, 6-311G*等不带有明显弥散的基组不能用)，还是只有专门计算弱相互作用时（Couterpoise或者大基组直接用E(AB)-E(A)-E(B))，才必须使用带有明显弥散的基组呢？

谢谢各位老师！

作者
Author: sobereva 时间: 2022-3-2 02:59
3-zeta基组算弱相互作用必须带弥散函数，要么用4-zeta档次基组。

根本就不麻烦，仔细看
给def2以ma-方式加弥散函数的Gaussian格式的基组定义文件（含所有def2支持的元素）
http://sobereva.com/509（http://bbs.keinsci.com/thread-14370-1-1.html）
撑死了写输入文件时候多花你10秒钟改个关键词和写个路径

6-311+G**耗时和ma-def2-TZVP有天壤之别，根本就不是能相提并论的东西

搞清楚你算单点能到底是为了什么。倘若你要对比两个构象能量差，即便目的不是专门算弱相互作用能，但两个构象中弱相互作用有极大的差异（比如有个构象有pi-pi堆积而另一个构象没有），显然算电子能量也得带弥散函数。

作者
Author: WaterMirror 时间: 2022-3-2 03:09

sobereva 发表于 2022-3-2 02:59
3-zeta基组算弱相互作用必须带弥散函数，要么用4-zeta档次基组。

根本就不麻烦，仔细看

哦哦谢谢Sob老师！
老师我还想问一下：我在小基组下优化的结构，得到的fchk1。使用大基组基于小基组得到的几何计算单点，得到fchk2，是否fchk2里面波函数的“质量”要高于fchk1，更有利于后续进一步的分析呢？
如果使用带有弥散的大基组得到的fchk2_1，是否“质量”也比不带弥散的大基组得到的fchk2_2更好，后面分析的时候更加清晰呢？
谢谢Sob老师！

作者
Author: sobereva 时间: 2022-3-2 04:17

WaterMirror 发表于 2022-3-2 03:09
哦哦谢谢Sob老师！
老师我还想问一下：我在小基组下优化的结构，得到的fchk1。使用大基组基于小基组得 ...

是。但是有些波函数分析方法如Mayer键级、Mulliken电荷、CDA分析等不兼容弥散函数（例如原子电荷受基组的影响在《原子电荷计算方法的对比》http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml的图5里有清晰的体现、CDA分析不适合用弥散函数在我的GCDA原文《广义化的电荷分解分析(GCDA)方法》http://www.hanspub.org/journal/P ... .aspx?paperID=16370里专门说了），即便不带弥散函数的时候，对于它们也不是基组越大一定越好

作者
Author: 喵星大佬 时间: 2022-3-2 14:43
may-cc-pVTZ就好了，虽说不是性价比很高的选择，但也不至于说算不动

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