计算化学公社

标题: 高斯计算错误 Problem with coordinate system [打印本页]

作者
Author: 小和尚 时间: 2015-5-22 21:17
标题: 高斯计算错误 Problem with coordinate system
我要优化水合氢离子（附件），用MP2可以得到结果，但是用DFT方法，B3LYP和MO6-2X优化，都出错误。请问怎么可以解决这个问题，非常感谢

输入文件：

%chk=4H2O+.chk
%mem=1200MB
%nprocshared=4
#p m062x/6-311++g(d,p) 6d opt freq

4H2O+

1 1
O                0.00164000 -0.00322000 -0.34664000
H                0.93615900 -0.32524700 -0.13102900
H                -0.18777600 0.96772000 -0.13178300
H                -0.74603100 -0.65054500 -0.13038600
O                -1.88372400 -1.64106100 0.09326700
H                -2.35818200 -2.08355300 -0.61800300
O                2.36930900 -0.80104400 0.09412500
H                2.98469500 -0.99018800 -0.62118900
H                2.81483800 -0.99065700 0.92510700
H                -2.26062800 -1.93412500 0.92837600
O                -0.48765000 2.44620100 0.09346800
H                -0.63032600 3.07505300 -0.62067200
H                -0.54933800 2.92453300 0.92582400

输出文件：
(Enter /home/g09root/g09/l701.exe)
Compute integral first derivatives.
... and contract with generalized density number  0.
Leave Link  701 at Fri May 22 16:39:53 2015, MaxMem=  157286400 cpu:    0.3
(Enter /home/qinsong/g09root/g09/l702.exe)
L702 exits ... SP integral derivatives will be done elsewhere.
Leave Link  702 at Fri May 22 16:39:53 2015, MaxMem=  157286400 cpu:    0.0
(Enter /home/g09root/g09/l703.exe)
Compute integral first derivatives, UseDBF=F ICtDFT=          0.
Integral derivatives from FoFDir, PRISM(SPDF).
Calling FoFJK, ICntrl=    2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
FoFCou: FMM=F IPFlag=          0 FMFlag=    100000 FMFlg1=       800
      NFxFlg=          0 DoJE=F BraDBF=F KetDBF=F FulRan=T
      Omega=  0.000000  0.000000  1.000000  0.000000  0.000000 ICntrl= 2127 IOpCl=  0
      NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0
      I1Cent=          0 NGrid=          0.
Symmetry not used in FoFCou.
Leave Link  703 at Fri May 22 16:39:54 2015, MaxMem=  157286400 cpu:    5.5
(Enter /home/g09root/g09/l716.exe)
Dipole       = 1.17107945D-03-7.01849135D-03 4.18137571D-01
***** Axes restored to original set *****
-------------------------------------------------------------------
Center    Atomic                Forces (Hartrees/Bohr)
Number    Number             X             Y             Z
-------------------------------------------------------------------
   1       8       -0.000345023 0.000024465 0.000151972
   2       1       -0.000124061 0.000235504 0.000077355
   3       1          0.000156227 -0.000310664 -0.000021876
   4       1          0.000301152 0.000115369 -0.000205702
   5       8          0.000204652 0.000078523 -0.000002895
   6       1       -0.000031721 0.000034844 0.000265423
   7       8       -0.000298557 0.000116560 -0.000205894
   8       1          0.000183406 -0.000107224 0.000078080
   9       1          0.000177490 -0.000097347 0.000070819
   10       1       -0.000236277 -0.000170573 -0.000152666
   11       8          0.000125873 -0.000343366 0.000198128
   12       1       -0.000088503 0.000322194 -0.000105048
   13       1       -0.000024658 0.000101717 -0.000147695
-------------------------------------------------------------------
Cartesian Forces:  Max    0.000345023 RMS    0.000181661
Leave Link  716 at Fri May 22 16:39:54 2015, MaxMem=  157286400 cpu:    0.0
(Enter /home/g09root/g09/l103.exe)

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
NTrRot= -1 NTRed= 28 NAtoms= 13 NSkip= -5 IsLin=F
Problem with coordinate system.
Error termination via Lnk1e in /home/g09root/g09/l103.exe at Fri May 22 16:39:54 2015.
Job cpu time:  0 days  0 hours  0 minutes 42.8 seconds.
File lengths (MBytes):  RWF=    14 Int=    0 D2E=    0 Chk=    1 Scr=    1

作者
Author: sobereva 时间: 2015-5-22 21:18
6d是多余的，没任何好处，还会带来数值不稳定性。
写opt=cartesian再试

作者
Author: 小和尚 时间: 2015-5-23 12:56

sobereva 发表于 2015-5-22 21:18
6d是多余的，没任何好处，还会带来数值不稳定性。
写opt=cartesian再试

用这个方法解决了，非常感谢！能请教一下原理吗？谢谢啊

作者
Author: sobereva 时间: 2015-5-23 21:55

小和尚发表于 2015-5-23 12:56
用这个方法解决了，非常感谢！能请教一下原理吗？谢谢啊

默认在冗余内坐标下进行，优化有时容易出问题，特别是优化过程中组成二面角的四个原子恰好排成直线，或者组成键角的三个原子组成直线时，优化肯定报错。分子体系还好，团簇优化时这个问题出现得相对频繁。

作者
Author: jiangxia0909 时间: 2017-3-16 22:17

sobereva 发表于 2015-5-22 21:18
6d是多余的，没任何好处，还会带来数值不稳定性。
写opt=cartesian再试

您好！我优化的是锕系元素的羟基配合物，我也遇到同样的问题，我加了opt=cartesian，输入是
#p b3lyp/genecp 6D 10F SCF=(MaxCyc=1024) Opt=(MaxCyc=1024) Freq=NoRama
n opt=cartesian geom=check guess=read
----------------------------------------------------------------------
但是立马报错说
Illegal IType or MSType generated by parse.
Error termination via Lnk1e in /usr/local/g09/l1.exe
这是怎么回事呢？

作者
Author: sobereva 时间: 2017-3-16 23:29

jiangxia0909 发表于 2017-3-16 22:17
您好！我优化的是锕系元素的羟基配合物，我也遇到同样的问题，我加了opt=cartesian，输入是
#p b3lyp/g ...

先去掉垃圾关键词再说
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

作者
Author: jiangxia0909 时间: 2017-3-17 09:55

sobereva 发表于 2017-3-16 23:29
先去掉垃圾关键词再说
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

sob老师，我将多余的关键词去掉了，还是出现了 Error in internal coordinate system.的问题，加了opt=cartesian，输入是 #p b3lyp/genecp Opt Freq opt=cartesian
--------------------------------------
还是提醒Illegal IType or MSType generated by parse.这是怎么回事呢？

作者
Author: sobereva 时间: 2017-3-17 23:16

jiangxia0909 发表于 2017-3-17 09:55
sob老师，我将多余的关键词去掉了，还是出现了 Error in internal coordinate system.的问题，加了opt=ca ...

上传输入输出文件
opt不要写两遍

欢迎光临计算化学公社 (http://bbs.keinsci.com/)