|
|
各位老师大家好！在做一个离子体系（阴离子B -1，阳离子肼+1）的氢转移过渡态的时候，遇到如下问题：
阳离子为肼离子，在M062X/6-311++G**下进行结构优化(还试过了PBE0,B3LYP,MN15等等，基组也试过了很多了，不管能不能用，都试的那种），不管怎么优化，肼离子上的H原子都会到阴离子上，首先说明下初始模型是我从晶胞抠出来的，但是，我用Molclus进行构象搜索（做了200个），得到的结构跟初始模型是一样的。然后我就想着直接做过渡态（不优化反应物，直接做），但是一直振荡，帖子中所有的办法都试过了，所有挨个试的。后来，我采用了溶剂模型，iefpcm跟smd都试过了，在溶剂化模型的作用下，H原子不到阴离子上了，但是这个情况是不管怎么找过渡态都到不了了（产物不能稳定存在）。我同样采用构象搜索的办法用HB（阴离子加H,模拟成氢转移完成的产物）与N2H4进行搜索，HB上的H都回到阳离子上了。现在的情况总结一下就是，不加溶剂模型，反应物不能存在（H会直接到阴离子上），加了，产物就存在不了了。溶剂模型我用的参数是eps=51.7(肼的）。我从晶胞中接着抠出一些分子想做个显式模型，但是这一段时间算下来，结果是一样的。之后，我固定了其中的重原子，这个过渡态马上就找到了。
想请教各位老师，现在这种情况怎么解决，请各位老师给个办法，谢谢各位老师！
|
|
发表于 Post on 2024-5-26 20:46:59
收藏 Add to favorites
楼主 Author