zjxitcc 发表于 2025-9-29 09:04
（1）我不知道是我眼神不好，还是你上传错了，你上传的输入文件看起来是SA-CASSCF做限制性优化，没有任何CA ...

抱歉老师是写错了，确实是SA-CASSCF做限制性优化，我的问题是第二个，感谢老师的指导，我去研究一下。

（1）我不知道是我眼神不好，还是你上传错了，你上传的输入文件看起来是SA-CASSCF做限制性优化，没有任何CASPT2相关的东西。问题与输入文件不对应。如果你将来使用CASPT2，处理势能面交叉情形注意要使用XMS-CASPT2。
（2）势能面交叉点有无数个，因为两个高维曲面的重叠部分是一个次（次）高维曲面，往往不是一个点，而是千千万万个点，这其中能量最低的点称为极小能量圆锥交叉点（MECI，仅针对>2个原子的S0/S1情况）或极小能量交叉点（MECP，广义情况）。因此计算类型需要严格区分：势能面扫描 与 寻找MECI点 是两类不同的计算类型/计算任务，不可以混为一谈。

你可以猜想S0/S1 MECI点 很像某一个二面角旋转势能曲线上的某一个结构，但不能保证它就在那条旋转势能曲线上，顶多是接近而已。即使这条柔性扫描曲线上有一个结构它满足S0/S1能量差为零，只能说这个结构是势能面交叉点，但不一定是势能面极小能量交叉点（能量差为零，且能量局域极小，两个条件）。

所以你的问题要分成两种情况：
（i）你实际上想找的是MECI点，并不想做柔性扫描，或者说柔性扫描是次要任务，那么你应当首先在OpenMolcas手册里搜索minimum energy conical intersection或者MECI，有相应优化MECI结构的输入文件和关键词。
（ii）你实际上不想找MECI点，只想完成这一系列柔性扫描，那么在完成后你不能说优化出的结构是MECI点，因为你没有去优化真正MECI点的结构。问题转化为：如何让当前的限制性结构优化收敛。我个人的建议是算一次精确数值Hessian，然后读入Hessian，往往有助于（限制性）结构优化收敛。当然它没有Gaussian那种opt=calcfc/rcfc方便，需要你自己琢磨。