24小时热门版块排行榜

返回列表

当前主题已经存档。

loovfnd

至尊木虫 (著名写手)

[资源] 【转帖】网上看到，关于优化时对称性

关于对称性
1。gaussian中输入什么对称性，一般优化的结果仍然还是那个对称性，比如CO2，如果初始两个CO键长输入不是完全相等（比如一个1.214,一个1.215），那么程序就会判断为C∞v 对称，那么优化结果虽然键长几乎相等，但仍然认为是C∞v ，这个从振动频率或者分子轨道对称性上可以看出来。--我们知道，CO2实际上是直线的两边对称的构型，其对称性应该是D∞v 。因此，为了得到高的对称性，必须输入的时候，精确地输入数值，比如sqrt(2)，就要保留很多的小数点，180.0角，就不能写成179.9。

2。有时计算过程中对称性会变化，比如做过渡态的时候，这时需要用 IOP(2/16=3)，否则计算会出错退出。

3。比如用直角坐标输入一个正三角形构型，其对称性应该是D3h，但是如果输入的小数点后面的数字不够多，那么常常得到的是C2v或其它。为了消除输入文件中坐标的有效位数的影响，得到较高的对称性，可以降低对称性判断的严格性。一般可以用symm=loose，这等价于 IOP (2/17=4, 2/18=3)。还可以减小这4和3这两个数值，使得更加loose，但不能过小，否则会出错。symm=loose只是在第一步判断输入构型的对称性时用到。
此外，也可以用gaussview来调整设置初始构型的对称性。

4。如果要降低对称性，那么可以用symm(PG=C3v)等等来做。使判断出来的对称性为C3v的一个子群。即由PG来限制最高对称性。

附用到的IOP的详细解释。

IOp(2/16)
action taken if the point group changes during an optimization. IOP(2/16=0)Abort the job. IOP(2/16=1)Keep going. IOP(2/16=2)Keep going and leave symmetry on, using the old symmetry. IOP(2/16=3)Keep going and leave symmetry on, using the new symmetry

IOp(17)
Tolerance for distance comparisons in symmetry determination.
0 Default (determined in the symmetry package, currently 1.d-8).
N>0 10-N.
N<0 10N, use the same tolerance for orientation.

IOp(18)
Tolerance for non-distance comparisons in symmetry determination.
0 Default (determined in the symmetry package, currently 1.d-7).
N>0 10 -N.
N<0 10 N, use the same tolerance for orientation.

有时候分子对称性太高经常无法优化收敛，这时可以加一个nosymm，有时挺管用的。

回复此楼