wzkchem5 发表于 2021-7-27 15:21
如果没有特殊原因，个人建议选Davidson, relaxed reference
扫描曲线出现不平滑，说明二面角不是合适的 ...

老师，我明白了，谢谢您！！！

hhwang 发表于 2021-7-27 07:46
感谢老师您的回复，对于第二点，如图片所示，老师我现在选取的能量是MRCI的值，这个数值应该是没有加上Da ...

如果没有特殊原因，个人建议选Davidson, relaxed reference
扫描曲线出现不平滑，说明二面角不是合适的扫描坐标，这个和是不是molpro的结构没有关系，如果是高斯优化的结构，但是仍然固定同一根键，固定到同样的键长，大概率也有这个问题

hhwang 发表于 2021-7-26 12:57
感谢老师您的回复。还有一些疑问的地方想请教您。
1.老师在非平衡点优化时出现虚频，而且频率不太好消去 ...

1. 还是那句话，非平衡点本来理论上就不保证没有虚频，所以即使虚频不小或者距离不长，也不代表要消掉它。只有当你优化平衡点的时候，才需要让虚频数目和理论值一致（极小值点0个，过渡态1个），非平衡点没法这么做，因为非平衡点的虚频数目没有理论值，几个都有可能。
2. 对，注意根据你的需要，读取加Davidson或Pople校正的能量，而不是不加校正的能量
3. 你做柔性扫描的时候有固定任何的自由度吗？如果没有，那和高斯优化再用高斯做柔性扫描是等价的。如果高斯优化的结构和molpro优化的结构定性相符，那这样做也行，但是写文章的时候不应该说在molpro优化结构的基础上做的柔性扫描，以免误导
4. 参见http://sobereva.com/552

wzkchem5 发表于 2021-7-26 00:20
1. 一般的频率分析得到的虚频的数目只有在平衡点（梯度为0的点）才有意义。扫描（不管是柔性还是刚性）曲线 ...

感谢老师您的回复。还有一些疑问的地方想请教您。
1.老师在非平衡点优化时出现虚频，而且频率不太好消去的情况，我询问了做计算的老师，她的意思是在解离曲线上有虚频是很有可能发生的，在5埃或者10埃这样较长的距离下计算出现有几个虚频，如果虚频很小，可以忽略。请问老师可以这样近似处理吗？
我想用这种长距离的点计算出能量和未解离结构的能量的差和高斯计算的能量之间进行能量校正，类比频率校正的不合理性，这种能量校正的方法合理吗？
2.老师，不加零点能是不是直接用molpro计算的单点能，也就是最后输出的MRCI的能量？
3.您提到需要做hindered rotor处理，或者至少做quasi-RRHO处理，老师，我请教了同学，这边hindered rotor的处理方式是用molpro的优化的结构，用高斯中结构优化相同的方法采用高斯软件对其进行柔性扫描，这样处理合理吗？molpro的结构用高斯扫会不会引起更大的误差？
4.quasi-RRHO处理该怎样实现呢？
问题有点多，麻烦您了，非常感谢您！

1. 一般的频率分析得到的虚频的数目只有在平衡点（梯度为0的点）才有意义。扫描（不管是柔性还是刚性）曲线上的绝大部分点都不是梯度为0的点，所以即便有虚频也不说明任何问题。如果要判断分子在垂直于扫描坐标的方向上是否处于极小值点，应该把扫描坐标投影掉再做频率分析，但一般程序不支持这么做，起码高斯应该不支持，molpro不知道是否支持，估计也不支持
2. 扫描曲线不应该加零点能，因为在梯度非0的点计算出来的热力学量没有物理意义。但是解离能需要报一个不加零点能和一个加的（或者把不加的那个体现在扫描曲线图里，而不在正文里写出来，这样也行）
3. CASPT2没有频率校正因子（因为校正因子和活性空间大小有关，每个人的活性空间大小都不一样，所以无法提前拟合出校正因子给大家使用），所以不用校正，因为没法校正
4. 没办法拟合，因为拟合用的每个分子都需要不同的活性空间，你没法证明你拟合用的活性空间和你实际计算的活性空间对频率的影响是可比的。拟合校正因子一般只能用于单参考态方法
5. 需要。即使不做hindered rotor处理，至少也要做quasi-RRHO处理