Chemistry 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Chemistry
标  题: 优化稳定几何构型出现小的虚频zz from厦大
发信站: BBS 哈工大紫丁香站 (Wed Jun 30 10:14:15 2004)

from http://210.34.15.126/cgi-bin/topic.cgi?forum=3&topic=864

>Q1:                                                                          
                                          
 我在做irc计算时，在优化irc所指向的产物时，虽然收敛但在进行频率计算时，却出现

－60的虚颦，不知应算minimal,还是transition state. 恳请赐教。

A:
1.记得帖子上说过：虚频小于100可不考虑。如是这样的话，你这个就应对应极小值。记

得SUPI曾说过消除虚频是将虚频所对应简正模式的坐标加到输入的相应原子的坐标中。那

将相应原子的坐标减去虚频所对应简正模式的坐标值，是不是会是虚频扩大？
2.这句话是在zixia上请教的大侠对我说的，版权不在于我，呵呵！
其实，对于虚频我现在仍然很头痛，尤其是低激发态的分子的几何构型优化，总是会有
虚频存在，动用了以前的一些经验全然不顶用，555

3. gaussian计算有个问题就是可能有些对称结构的频率可能出现小虚频
这个时候就是稍微降低一下它的对称性
4. 这个情况我也遇以过,后面降低对称性消频成功了.

5.很小的虚频可以用Grid命令通过增大积分的网格数消除。
 int(grid=ultrafine)
一般仅对小虚频有用。－－我算的吸附体系，虚频在20几，只能稍微减小一点，但是去
不掉。

>Q2 比如一个平面结构，Cs对称性， 要是降低了就成C1了。                         
                                          
那写文章时说是Cs还是C1？ or quasi－Cs ？

A Cs对称性，C1只是gaussian计算的bug而已

>Q3 请教一下，什么是“虚频所对应简正模式的坐标值”，在输出文件的什么地方去找? 
                                         
A 在输出文件中查找“ Frequencies” 就会看到。

>Q4  感谢诸位大侠的讨论，另我受益匪浅。我在优化过程中，没有加任何对称性限制   
                                          
得到的结构是C1对称的。是一个平面结构，-60的虚颦对应两个分子之间垂直于平面的微

弱摆动。不知道，从一个过渡态是否可能指向一个van der waals 过渡态。如何保证体
系较低的对称性。

A 这种情况可以尝试一下把这个虚频消除掉,把构型朝虚频振动的方向改变一点再作优化.


>Q4  如何降低对称性?                                                          
                                          
比如一个xy平面结构,有z方向的虚频, 但是在z方向,分子中原子的受力还是为0啊.
我试过用nosymm, 结果还是平面结果,对称性不会降低啊.---当然虚频也就没有除掉了.
A  初始构型就把那个沿虚频振动的原子向z方向偏一下。比如说有对称性时二面角为0度

，现在初始值给成10度。

--
╔═══════════════════╗
║★★★★★友谊第一  比赛第二★★★★★║
╚═══════════════════╝


※ 来源:·哈工大紫丁香 http://bbs.hit.edu.cn·[FROM: 202.118.229.*]