计算化学公社

标题: 求助：只有一个虚频但是并不是手性轴旋转的方向 [打印本页]

作者
Author: NicoleDH 时间: 2025-3-30 15:45
标题: 求助：只有一个虚频但是并不是手性轴旋转的方向
本帖最后由 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寻找过渡态。

非常感谢！

作者
Author: naonao5205 时间: 2025-3-30 18:06
建议考虑是否收敛到正确的电子态，过于扭曲时候可能是双自由基。
TS guess可以先做一个freq看下虚频是否大概合理，然后直接读hessian即可

作者
Author: Solitude198 时间: 2025-3-30 19:09
经典问题

http://bbs.keinsci.com/forum.php ... mp;page=1#pid329722

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