计算化学公社

标题: gaussian对阴离子体系部分加弥散和全部加弥散过渡态能垒差距太大 [打印本页]

作者
Author: zhibisuyusheng 时间: 2025-6-14 14:24
标题: gaussian对阴离子体系部分加弥散和全部加弥散过渡态能垒差距太大
体系中有阴离子需要加弥散函数，但是由于体系部分结构原子数太多用6-31+g*算的很慢，所以我只对TFA阴离子的C O 加了弥散，算完结果发现全部加弥散的话过渡态能垒是17.45，只对TFA阴离子的C O加弥散能垒是4.27，请问这种结构是不能只对阴离子部分加弥散吗

作者
Author: 18217265596 时间: 2025-6-14 15:34
粗粗看了看，为什么用M06?
弥散部分:一般加弥散对单点精度的影响总是好的，只不过会算的慢，但是加弥散会改变势能面，不加弥散的过渡态和加弥散的应该是不同的
另外：一切误差都是计算精度不足造成的，你这结构上dlpno ccsd(t）带一个低级别弥散应该是算的动的
算几个高精度单点就知道低精度下的误差是哪里来的了

作者
Author: zhibisuyusheng 时间: 2025-6-14 15:52

18217265596 发表于 2025-6-14 15:34
粗粗看了看，为什么用M06?
弥散部分:一般加弥散对单点精度的影响总是好的，只不过会算的慢，但是加弥散会 ...

我同个势能面有铜，112原子，所以不能用太高级别，试过用m06 6-31+g* lanl2dz直接算不动了

作者
Author: 18217265596 时间: 2025-6-14 16:54

zhibisuyusheng 发表于 2025-6-14 15:52
我同个势能面有铜，112原子，所以不能用太高级别，试过用m06 6-31+g* lanl2dz直接算不动了

112原子的话双杂化也能算动的，计算资源不足建议去orca中RI加速scf会快很多，用multiwfn生成输入文件基本没有学习成本，你只算个单点是很合适的

作者
Author: sobereva 时间: 2025-6-15 05:20
下文明确说了，优化和单点的级别通常应当明确区分。绝对不要用个对opt freq浪费、对单点又明显不够的基组做opt freq然后直接拿得到的能量算势垒
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）

6-31G*这种区区2-zeta档次的基组，对于算能量问题根本不够，加了弥散依然太次。认真看下文
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈弥散函数和“月份”基组
http://sobereva.com/119

算单点建议用def2-TZVP，带负电荷显著的用ma-def2-TZVP。为节约时间，不重要区域的原子可以降低基组。

倘若你确保运行方式无误（恰当并行了，没有用Windows版Gaussian），对于单点来说还出现“m06 6-31+g* lanl2dz直接算不动了”，说明该用个明显更像样的机子。

没事甭用M06做opt freq，本来耗时就高，还对积分格点要求高、容易出虚频，用了int=superfine更是耗时暴增。结合SMD做opt freq更是纯粹自找不快，又难收敛又容易出虚频，下文都明确说了
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

并且把此文也认真看了
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

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