激发态频率计算

赵云跳槽

我需要计算激发态的频率
为了防止任务过大，我把计算分为了两步，先是优化S1，#p opt pbe1pbe/6-31g(d) td=(nstates=3,root=1)，
计算后完成后再计算频率 #p opt freq td=(nstates=3,root=1) pbe1pbe/6-31g(d)
现在有点疑问：
（1）我明明使用的优化完成后的构型重新生成的gjf文件（没有使用chk文件，防止坐标或文件丢失），为何在第二步计算中还是需要多次优化收敛再计算频率？
（2）第二步经过多次收敛后得到的能量与第一次的不同，较第一步低，怎样理解？
（3）明明已经找到了S1-opt，为何后续计算还要进行能量以及激发计算，而不是直接进行频率计算？

Gaoktp

防止坐标文件丢失你就把第一步的chk文件另外复制一个,
计算频率也不需要写opt，写了opt你要在输出文件中找找有没有 “Optimization completed”有的话说明你第二步的opt结束了，是在计算频率了，频率本来就慢！！参见  http://bbs.keinsci.com/forum.php ... tid=1221&extra=

sobereva

第一步用opt，第二步用freq就完事了。

1 这是很正常的。几何优化时的Hessian是不断更新的，最后一步的Hessian和你重新优化时默认基于价力场生成的初始Hessian是不同的，所以即便之前已经收敛了，重新再优化的时候也并非能立刻收敛。另外，重新生成新输入文件时还有坐标精度的轻微丢失。

所以第二步没必要再写opt。防止把chk搞坏一方面像lao7说的那样自行备份一下，要么就利用g09后期版本加入的%oldchk设定让计算一开始就自动把之前的chk复制为当前任务的chk。

2 比较几何结构，看是否有明显偏差。如果是，可能收敛到了不同极小点。如果只是有轻微的数值差异可以忽视，是数值精度问题造成的。

3 做激发态频率计算自然要做电子激发计算。