计算化学公社

标题: 求助，弱氢键主导的质子解离过渡态TS搜寻问题 [打印本页]

作者
Author: j5888xm 时间: 2020-10-19 17:12
标题: 求助，弱氢键主导的质子解离过渡态TS搜寻问题
本帖最后由 j5888xm 于 2020-10-19 17:15 编辑

求助各位老师，我计算一个两分子体系的过渡态，为 α-H在溶剂下的离解过程，体系内有明显的pi-pi堆积作用，使用B3LYP(D3)/6-31G(d,p)级别下计算，多次尝试无果，均找不到正常的过渡态，过渡态初猜均为半经验柔性扫描的突跃点。

根据sob老师JCC文章的方法，计算所得 α-H与溶剂间的氢键约为6~7kJ/mol。

询问各位老师几个问题：

一、对于弱相互作用明显的体系，只用加了D3矫正的B3LYP是不是还不够，需要考虑换M06-2X(D3)吗？还是说简化分子，然后使用CCSD(T)/cc-pvtz找过渡态。

二、若是弱氢键主导质子离解，那么这个强度的弱氢键和偶极作用有无明显差别。

谢谢大家。

作者
Author: wzkchem5 时间: 2020-10-19 19:55
B3LYP-D3一般够了，如果不够那也是因为BSSE，通过加大基组或counterpoise的方法减小BSSE的情况下，B3LYP-D3一般不致有定性问题。
什么叫柔性扫描的突跃点？如果指的是扫描过程中能量增加然后突然断崖式下跌，那不叫过渡态，说明你选的扫描坐标不合理。建议想办法得到一个更好的过渡态结构初猜

作者
Author: 喵星大佬 时间: 2020-10-19 20:14

wzkchem5 发表于 2020-10-19 19:55
B3LYP-D3一般够了，如果不够那也是因为BSSE，通过加大基组或counterpoise的方法减小BSSE的情况下，B3LYP-D3 ...

他应该是溶剂模型的问题，我看了是个涉及显式溶剂的反应，但是既没有显式溶剂甚至隐式溶剂都没加，显然找不到，气相下整一个单质子的能量还是很高的

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