关于之前那个水解反应机理的贴子

(前情提要)
前几天有个关于化学的贴子里有个老哥想问一下这个水解反应的可能的机理
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-l3hyKeT1kSe5-80.jpg.thumb.jpg[/img]

然后讨论里我看到有3种可能的机理
第一是我猜的
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-k3r4K8T1kSa6-4q.jpg[/img]

第二是
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-jfscK9T1kScp-4m.jpg[/img]


第三是
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-2lmaKcT3cShy-46.jpg[/img]

其中第二和第三的protonation可能和提出者的想法不太一样，但是transition state是差不多

前几天事情比较忙，就没怎么仔细算，今天折腾了一下，可能有一些比较粗糙的成果(甚至能不能定性都要打个问号)。想着既然都算了就发出来看看，结果那个贴子好像没法回复了，只好另开一个新贴。

所有计算的级别都是B3LYP/6-31G*
首先是我做了几个计算确认了第一个质子应该是在protonated在pyrazole的氮上的。
几何优化后的
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-82v8K3T8S2a-34.jpg[/img]
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-2maqK13T3cSsg-mc.jpg.thumb.jpg[/img]
比
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-7px1K3T8S2a-34.jpg[/img]
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-gy7uK14T3cSsg-mc.jpg.thumb.jpg[/img]
能量低了20.8 kcal/mol

比较有趣的是后面一种可能会导致一种之前没想到过的副产物
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-l9owZ1bT3cSxe-is.gif.thumb.jpg[/img]
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-icsvZ1cT3cSxe-is.gif.thumb.jpg[/img]
这个氧拿了质子直接带着ethanol就跑了[s:ac:哭笑]

然后今天在痛苦的找transition state
mech3是我觉得最有意思的，所以就先试了试。它的四元环应该是最不稳定的，试了几次就放弃了，打算看看mech1和mech2的情况，等这两个算完了再试试看。
mech2的ts比较顺利，很快就找到了ts，虚频在-103.99 cm-1,振动模式看起来应该是反应对应的ts。
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-95coXfZ4bT3cSxs-j0.gif.thumb.jpg[/img]

mech1就有点难，目前只找到了一个结构，只有一个虚频，并且振动模式反应比较接近
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-at4lXdZ3pT3cSzo-k2.gif.thumb.jpg[/img]
但是跑几何优化的时候一直跑不出来优化后的transition state，而且我现在在想是不是需要的氯离子帮忙。
比如
[img]https://img.nga.178.com/attachments/mon_202011/14/-7Q5-5ijlXjZ5eT3cSzo-k2.gif.thumb.jpg[/img]
(这个结构并不是最优化的ts，我做了个freq结果看到了一车的虚频[s:a2:大哭])
但是似乎这个结构已经比优化完了的mech2的ts要低5.12 kcal/mol了