ORCA 过渡态与虚频实操

下载完整输入、实际结果、原生振动动画和绘图脚本。解压进入资料包目录，先提取已经计算好的结果：

bash analyze.sh

图片在 reference/figures，路径CSV、能垒表和完整数据在 reference/analysis。在新解压目录重新计算整套流程：

bash submit.sh

脚本按两档基组提交两个数组任务。每个任务依次准备两个端点和直接过渡态、做NEB与IRC、重新优化IRC两端，最多同时运行3个独立单进程ORCA计算，再自动分析；新数据在 work，依赖几何和Hessian从新算出的上一步传递。已有 reference 保留。

使用 orca/6.1.1 厂商预编译CPU程序，依赖专用OpenMPI4.1.8（GNU13.3环境）；单次计算 nprocs=1，不在外层加 mpirun。本例在P100节点使用CPU执行，不是ORCA GPU加速。附件模板每个数组任务申请3个CPU核心和1张P100，GPU未参与ORCA计算，资源按相应分区规则计费。不要在登录节点运行这些计算；已有结果可以直接后处理。

1 有一个看起来平面的结构，还不能直接称为过渡态#

本例氨翻转的直接猜测是N和三个H位于同一平面。使用 OptTS 让程序优化鞍点，Freq 检查最终Hessian；还要用IRC确认它连接的两个极小点。

完整实际输入：

! PBE0 def2-SVP def2/J RIJCOSX DefGrid3 VeryTightSCF TightOpt OptTS Freq
%pal nprocs 1 end
%maxcore 500
%geom
 MaxIter 100
 Calc_Hess true
 Recalc_Hess 5
end
* xyz 0 1
N 0 0 0.0000000000
H 0.945 0 0
H -0.4725 0.8183940066 0
H -0.4725 -0.8183940066 0
*

Calc_Hess true 从实际Hessian出发寻找负曲率方向，Recalc_Hess 5 定期重算；初始猜测仍需接近目标鞍点。这里是一个小分子的对称翻转，不能保证任意化学反应都能从随手摆出的平面结构收敛。

原生运行方式：

module load orca/6.1.1
"$ORCA_ROOT/orca" case.inp > case.out
grep 'HAS CONVERGED' case.out
grep -n 'VIBRATIONAL FREQUENCIES' case.out

2 看最后一段频率，区分真实振动和刚体运动#

优化过程中可能先计算初始Hessian，输出不止一段频率。判断最终结构时，看最后一段 VIBRATIONAL FREQUENCIES，不要把初始猜测的频率当成最终结果。

NH₃有4个原子、12个笛卡尔自由度。孤立非线性分子移除3个平移和3个转动后，有6个内部振动；前6个接近0的刚体模式不用于判断反应鞍点。

这里每个最终TS都有一个显著负频率，两个端点都没有显著虚频。非常小的负值应结合优化、网格和数值精度检查，不能把任意负小数都当作化学反应方向。

附件读取原生Hessian，独立做质量加权、去除刚体运动和对角化，与原生频率比较；图中显示的是实际结果。振动模式动画用于看位移方向，不是有限温度MD轨迹，也不代表真实热振幅。

3 用原生工具生成虚频模式动画#

在过渡态结果目录执行：

"$ORCA_ROOT/orca_pltvib" case.hess 6

得到 case.hess.v006.xyz，本例有20帧，可用支持多帧XYZ的软件打开。本例原生输出的0–5为刚体模式，6对应虚频；其他体系要先查看实际模式编号，再决定提取哪个模式。

这段动画中，N相对H平面来回穿过，符合NH₃伞形翻转。附件还核对了动画位移与独立求得的负曲率本征向量。

注意：orca_pltvib 的6是本例输出的笛卡尔模式编号；IRC中默认 hessMode 0 选择其内部最低模式，不能把两个编号体系直接互换。

4 虚频检查之后做什么#

有一个虚频只说明局部曲率符合一阶鞍点，仍要确认这是不是目标反应。下一步将最终 case.xyz 和 case.hess 一起传给IRC，方法、基组、电荷、自旋与原子顺序保持一致。

如果虚频主要是无关的侧链摆动、分子间转动或其他结构变化，就不能据此报告目标反应能垒。先检查位移方向，再回到初始路径或TS猜测调整。

下一篇给出完整IRC输入、两端再优化和三种能垒。基础优化和普通频率见ORCA优化与频率。