计算化学公社

标题: z-matrix如何冻结片段TS优化 [打印本页]

作者
Author:
chenxin199261    时间: 2014-11-13 12:03
标题: z-matrix如何冻结片段TS优化
我在用z-matrix优化时候。
只有两个variables,其余都是constants

但在优化过程中,我发现我发现constant的变量,仍然在计算力,导致体系始终不收敛。
请问,这个情况怎么办?


我正在rigid扫描两个变量,实在不行就从二维势能面上找吧

another question:
我用cartesian优化,冻结原子,是否会造成中间的原子收到外围原子的挤压而难以收敛?

作者
Author:
sobereva    时间: 2014-11-14 07:11
虽然constant的变量会计算受力,但是这不是不收敛的原因,因为在判断是否达到收敛限的时候constant变量是不纳入统计的。

这是个简单的例子,从输出信息中会看到最后一步时A1受力明显不为0,但判断受力时只考虑了B1和B2,它们很接近0,所以正常收敛了。

  1. #p b3lyp/6-31G** opt=z-matrix

  2. b3lyp/6-31G** opted

  3. 0 1
  4. O
  5. H 1 B1
  6. H 1 B2 2 A1

  7. B1=0.76533395
  8. B2=0.76533395

  9. A1=123.66503176
复制代码

几何优化不收敛本质上是数值问题层面的原因,巧合因素居多,和是否设定冻结没有必然关系,尽管在设定时越胡来越可能导致收敛的难度增加。

作者
Author:
panger    时间: 2014-11-15 13:37
冻结坐标会不会引入额外的虚频
作者
Author:
sobereva    时间: 2014-11-15 13:58
panger 发表于 2014-11-15 13:37
冻结坐标会不会引入额外的虚频

如果坐标冻结的数值和它实际优化后的值差异较大,则肯定会引入虚频




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3