高斯计算阴离子脱氢的过渡态失败

wyx123

老师好，我想设计阴阳离子（阴离子-1价 阳离子+2价）逐渐接近成单键，然后成双键同时阴离子脱氢得过程。第一步找到了单键成双键的过渡态,irc后优化得到的结构符合预期（图1）第二步直接脱氢的过渡态计算失败（图2）； 于是设计了图三中脱出来的氢和另一种阴离子（-1价）结合， irc 后得到的两端结构相似（图3），第二个阴离子也跑到另一侧与阳离子结合了，计算失败。关键词分别为opt=(calcfc,ts,noeigen) freq b3lyp/6-31g(d) 和 irc=(calcfc,maxpoints=50,lqa) b3lyp/6-31g(d)

请问老师 我设计的脱氢过程不合理吗，或者阴离子先脱氢再与阳离子结合？

file:///C:/Users/13662/AppData/Local/Temp/msohtmlclip1/01/clip_image002.png图1

file:///C:/Users/13662/AppData/Local/Temp/msohtmlclip1/01/clip_image004.png图2

file:///C:/Users/13662/AppData/Local/Temp/msohtmlclip1/01/clip_image006.png图3

dnlx

你似乎在模拟这样的反应：铜离子 +　碳酸氢根　－－＞　碳酸铜 + 氢离子

先假定你要模拟的产物是碳酸铜离子对，而不是碳酸铜沉淀物。否则模拟方法就错了。

如果确实模拟的是形成碳酸铜离子对：
1. 这个反应可能是两步反应，真正形成离子对的依然是碳酸根和铜离子，或者形成CuHCO3+离子、再电离出氢离子。
2. 这个反应很可能是热力学而不是动力学控制的，需要计算相关物质的热力学量，而不是寻找过渡态。
3. 应该用隐式溶剂模型，而且离子周围要加上足够个数的显式溶剂分子。

wyx123

dnlx 发表于 2023-12-15 12:40
你似乎在模拟这样的反应：铜离子 +　碳酸氢根　－－＞　碳酸铜 + 氢离子

先假定你要模拟的产物是碳酸铜 ...

老师 我要模拟的是沉淀物 该怎样模拟呢

dnlx

wyx123 发表于 2023-12-15 20:30
老师 我要模拟的是沉淀物 该怎样模拟呢

没做过，不知道是不是属于分子动力学模拟。

wyx123

dnlx 发表于 2023-12-15 22:26
没做过，不知道是不是属于分子动力学模拟。

谢谢老师