计算化学公社

标题: 关于自旋极化单重态过渡态溶剂化后，能垒变高问题 [打印本页]

作者
Author: btm 时间: 2023-5-31 17:23
标题: 关于自旋极化单重态过渡态溶剂化后，能垒变高问题
   最近发现涉及自旋极化单重态的过渡态（如ISET反应&自由捕获）进行溶剂化后，再去算反应的能垒，能垒就会变高了。由于我比较了其他反应，例如SN2反应，溶剂化前ISET为优势反应路径（ΔΔG=5kcal/mol）,溶剂化后SN2反应变为了优势路径（ΔΔG=2.7kcal/mol）,由于实验自由基钟证明了以自由基路径为主，因此ISET必然应该能垒低才对，使用马库斯公式求得OSET的能垒21也是比SN2能垒低，一篇文献说ISET一定比OSET能垒低，我尝试了更改泛函目前尝试了wb97xd/m06/m062x,基组(C H O N Cl)6-311G(d,p) (Cu)SDD但是并没有用，因此贴出了自己的计算文件的关键部分，希望老师能指正一下。
   （ Ⅰ ）先找稳定波函数：
   %chk=sol-WTS2-1.chk
# stable=opt guess=mix uwb97xd/genecp guess=read SCRF=(read,SMD)

Title Card Required

0 1
H                -1.83065800 -1.81185000 1.30680300

C H O N Cl 0
6-311g(d,p)
****
Cu 0
SDD
****

Cu  0
SDD

eps=x
epsinf=x
（ Ⅱ ）读波函数：
%oldchk=sol-TS2-1.chk
%chk=sol-TS2-1.chk
# uwb97xd/genecp guess=read SCRF=(read,SMD)

Title Card Required

0 1
H                0.39516100 4.14087900 -3.55310700

C H O N Cl 0
6-311g(d,p)
****
Cu 0
SDD
****

Cu 0
SDD

eps=x
epsinf=x

作者
Author: wzkchem5 时间: 2023-5-31 19:43
算SN2的时候，做了显式溶剂化了吗？

作者
Author: btm 时间: 2023-5-31 19:52

wzkchem5 发表于 2023-5-31 19:43
算SN2的时候，做了显式溶剂化了吗？

没有老师，SN2应该注意什么呀

作者
Author: wzkchem5 时间: 2023-5-31 20:31
本帖最后由 wzkchem5 于 2023-5-31 14:13 编辑

btm 发表于 2023-5-31 12:52
没有老师，SN2应该注意什么呀

如果体系里有和溶剂作用非常强的原子，例如有电荷比较集中的离子（Na+, Cl-, [SO4]2-这类），需要把这些离子第一配位层的溶剂分子画出来，因为溶剂模型不适合用来描述和溶质作用特别强的溶剂。否则计算结果不准

作者
Author: btm 时间: 2023-5-31 20:46

wzkchem5 发表于 2023-5-31 20:31
如果体系里有和溶剂作用非常强的原子，例如有电荷比较集中的离子（Na+, Cl-, [SO4]2-这类），需要把这些 ...

好的老师

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