|
|
本帖最后由 NicoleDH 于 2025-3-30 15:45 编辑
大家好，
我正在用Gaussian09做一个轴手性化合物的旋转能垒计算，该化合物有两个芳香平面通过手性轴连接而成。我想计算从A构象到B构象的旋转能量和半衰期。
通过对A的优化构象进行二面角扫描后，我选取了能量最高的点作为初猜过渡态结构，柔性扫描的关键词如下：
# opt=modredundant B3LYP/6-31G(d,p) scf=tight int=ultrafine
...
D 5 6 14 27 S 36 10.0
将初猜过渡态结构（TS_guess）做优化，关键词如下。虽然正常计算完毕，也收敛了，只有一个虚频，但是虚频的震动方向并不是手性轴的转动方向，而是一个基团上氢原子的转动。
检查优化后的初猜过渡态（TS_guess_opt），发现TS_guess_opt与TS_guess的二面角差别比较大，而与A却很接近。
# opt=(ts,calcfc,noeigen) freq B3LYP/6-31G(d,p) scf=tight int=ultrafine
请教大家：
1. 二面角扫描的曲线为什么不够平滑？扫描这一步就出问题了吗？
2. 虚频的震动方向不是手性轴的转动方向这个问题该如何解决？我看社长的帖子里说，可能是初猜TS结构不合理或者计算级别不合适。我选的TS_guess中两个芳环平面基本上是平行的，两个围绕轴转动肯定会有这个中间态的。
3. 我还拿能量第二高的点作过初猜过渡态结构，做一样的优化，但出现了一样的问题。
看了社长的帖子，应该暂时没必要做QST2寻找过渡态。
非常感谢！
|
|
发表于 Post on 2025-3-30 15:45:30
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites