请教过渡金属配合物催化中自由能计算方法

lllnankai

请教各位老师专家，对含Cu的过渡金属配合物催化的反应，计算各中间体和过渡态的自由能，采用下列方法和步骤:
1. opt+freq用 PBE0-D3(BJ)/def2-SVP （Cu用SDD），结合 scrf=solvent=dichloromethane,得到优化的结构和自由能校正量G(correction)；
2. 高精度气相单点能 用PBE0-D3(BJ)/def2-TZVP ,得到 E(gas)；
3. M052X/6-31g* (Cu用SDD) ，再算气相单点能,得到 E(gas,m052x)  //考虑到Cu是2配位，没有被有机部分充分包围，因此对Cu单独使用SDD赝势基组
4. M052X/6-31g* (Cu用SDD)，结合SMD溶剂模型 scrf=(smd,solvent=dichloromethane)，再算溶液中单点能,得到E(sol,m052x)
最后用 G(correction) + E(gas) + E(sol,m052x) - E(gas,m052x) + 1.89 kcal/mol 得到各物质的自由能，来绘制反应路径，
请教各位老师专家，上述方法是否合理？谢谢！

sobereva

没必要那么折腾，直接参考下文的方式计算溶剂下的自由能就完了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

对于过渡金属配合物，并没有文章证明M05-2X/6-31G*是结合SMD计算溶解自由能的首选，因此没必要刻意专门在这个级别计算溶解自由能。

lllnankai

sobereva 发表于 2026-7-14 12:55
没必要那么折腾，直接参考下文的方式计算溶剂下的自由能就完了
使用Shermo结合量子化学程序方便地计算分子 ...

感谢社长，去学习一下