GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Eigenvalue   189 is   0.00D+00 should be greater than     0.000001 Eigenvector:
                         A1        A3        A2        A14       A10
  1                    0.34057   0.34057   0.34057  -0.32109  -0.21590
                         R27       R26       A13       R24       D29
  1                    0.19454  -0.16251  -0.15775  -0.14252   0.13054
NTrRot=   140 NTRed=   329 NAtoms=    65 NSkip=   140 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /home/zyang/g09/l103.exe at Thu Mar  5 14:27:58 2015.
Job cpu time:  0 days  1 hours 54 minutes 38.0 seconds.
File lengths (MBytes):  RWF=    408 Int=      0 D2E=      0 Chk=    113 Scr=      1
-------------------------------------------------------------

单纯的QM 这种错误，用opt（cartesian）就可以解决，但是对于

蛋白晶体结构，这种方法就不适用，No microiterations with Cartesian

-----------------------------

这是高斯软件的一个bug，主要是结构中存在线性结构，尝试修改小数点后的数字0位随机数，但是没用。

----------------------------------

经过尝试，发现对于我研究的这个体系，

g09 A 版本有这个问题，用更高版本的gaussian ,可能会避开这个问题。

用g09 D版本可以避开这个问题。

----------------------

增加一些原子，或者减少一些原子也可以避开这个问题，

比如我在高层中增加配体周围一些水分子。