计算化学公社

标题: 高斯任务激发态优化报错 [打印本页]

作者
Author: HERO 时间: 2016-6-4 15:41
标题: 高斯任务激发态优化报错
在进行激发态优化和频率计算时出现错误：Cycle 128  Pass 1  IDiag  1: RMSU=  9.90D-08 CP:  9.99D-01 -3.00D+00 -3.00D+00 -3.00D+00  2.69D+00
                  CP: -3.00D+00  6.51D-01 -2.54D-01  3.00D+00  3.00D+00
                  CP:  3.00D+00  3.00D+00  3.00D+00  3.00D+00  2.05D+00
                  CP:  1.33D+00  1.95D+00  2.10D+00
E= -2741.22831835809    Delta-E=       0.000026644580 Rises=F Damp=F
DIIS: error= 1.63D-04 at cycle 118 NSaved=  18.
NSaved=18 IEnMin= 2 EnMin= -2741.22961063042    IErMin= 1 ErrMin= 1.55D-04
ErrMax= 1.63D-04 EMaxC= 1.00D-01 BMatC= 1.75D-05 BMatP= 1.47D-05
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com:  0.157D+00-0.329D+00 0.582D+00-0.626D-01-0.109D+01 0.860D+00
Coeff-Com: -0.106D+01 0.366D+00 0.168D+00 0.101D+01-0.790D+00-0.595D+00
Coeff-Com:  0.298D+00-0.562D+00 0.218D+00-0.197D+01 0.308D+01 0.709D+00
Coeff:    0.157D+00-0.329D+00 0.582D+00-0.626D-01-0.109D+01 0.860D+00
Coeff:    -0.106D+01 0.366D+00 0.168D+00 0.101D+01-0.790D+00-0.595D+00
Coeff:    0.298D+00-0.562D+00 0.218D+00-0.197D+01 0.308D+01 0.709D+00
Gap=    0.099 Goal= None Shift= 0.000
RMSDP=5.82D-06 MaxDP=4.66D-04 DE= 2.66D-05 OVMax= 2.06D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RM062X) =  -2741.22831836    A.U. after  129 cycles
         Convg  = 0.5819D-05          -V/T =  2.0037
KE= 2.731175841958D+03 PE=-1.253626991990D+04 EE= 4.045640455029D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in /public/g09/l502.exe at Sat Jun  4 05:16:55 2016.
Job cpu time: 15 days 19 hours  0 minutes 34.4 seconds.
File lengths (MBytes):  RWF= 488 Int=    0 D2E=    0 Chk=    45 Scr=    1

第一激发态优化输入文件如下：
%nprocshared=15
%mem=32GB
%chk=TMP-S1.chk
#p opt freq m062x/TZVP td=(singlets,root=1) integral=ultrafine

TMP-S1

0 1
  Si       0.621711000000    0.347262000000    1.253534000000
  O       -0.460310000000    0.837395000000    2.346464000000
  O       1.944895000000    1.311714000000    1.300825000000
  O       -0.033731000000    0.388775000000    -0.264413000000
  Si       3.063242000000    2.115059000000    0.389530000000
  O       3.841219000000    1.012152000000    -0.572568000000
  Si       3.366686000000    -0.352302000000    -1.358902000000
  O       1.907223000000    -0.059378000000    -2.065634000000
  Si       0.311463000000    -0.350667000000    -1.711488000000
  Si       2.596833000000    -1.962870000000    1.236583000000
  O       3.161722000000    -1.568847000000    -0.269582000000
  O       1.149691000000    -1.197501000000    1.486386000000
  H       -1.390566000000    0.609495000000    2.147801000000
  H       2.377711000000    3.106750000000    -0.464178000000
  H       4.052143000000    2.728393000000    1.287738000000
  H       4.364780000000    -0.727229000000    -2.356268000000
  H       -0.517108000000    0.272824000000    -2.752956000000
  H       0.139762000000    -1.811833000000    -1.599555000000
  H       2.355294000000    -3.413923000000    1.273165000000
  H       3.577834000000    -1.522956000000    2.247437000000
  O       -2.875939000000    0.190746000000    1.309949000000
  P       -3.244065000000    -0.055594000000    -0.111979000000
  O       -4.743733000000    -0.559456000000    -0.309459000000
  O       -3.136244000000    1.202304000000    -1.071525000000
  C       -2.323698000000    2.345622000000    -0.752947000000
  H       -2.913556000000    3.224081000000    -1.001531000000
  H       -1.419709000000    2.311527000000    -1.354892000000
  H       -2.060019000000    2.350660000000    0.303710000000
  O       -2.425763000000    -1.199436000000    -0.860472000000
  C       -2.152230000000    -2.403346000000    -0.120833000000
  H       -1.784319000000    -3.125448000000    -0.843150000000
  H       -3.062619000000    -2.778875000000    0.346589000000
  H       -1.390847000000    -2.207867000000    0.635502000000
  C       -5.786796000000    0.180514000000    0.336428000000
  H       -5.867619000000    1.174978000000    -0.104536000000
  H       -5.589141000000    0.259165000000    1.405347000000
  H       -6.706041000000    -0.371073000000    0.164224000000

在基态结构基础上，计算激发能如下：
Excited State 1:    Singlet-A    8.2543 eV  150.20 nm  f=0.0033  <S**2>=0.000
   102 ->109       0.11476
   103 ->108       0.16158
   103 ->109       0.23071
   103 ->110       -0.12505
   103 ->111       0.11225
   104 ->105       0.45591
   104 ->106       -0.21187
   104 ->107       -0.21019
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2741.16449160
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    8.3178 eV  149.06 nm  f=0.0002  <S**2>=0.000
   102 ->106       -0.12781
   102 ->109       -0.10823
   103 ->105       0.35015
   103 ->106       -0.25592
   103 ->107       -0.20316
   104 ->106       0.11633
   104 ->108       0.21598
   104 ->109       0.27633
   104 ->110       -0.16748

Excited State 3:    Singlet-A    8.3579 eV  148.34 nm  f=0.0001  <S**2>=0.000
   102 ->105       0.14549
   102 ->108       -0.14898
   102 ->111       -0.11014
   103 ->105       -0.19278
   103 ->108       -0.18751
   103 ->111       0.10051
   104 ->105       0.28963
   104 ->106       0.42157
   104 ->110       -0.11039

Excited State 4:    Singlet-A    8.4534 eV  146.67 nm  f=0.0038  <S**2>=0.000
   102 ->105       0.49364
   102 ->106       0.10120
   102 ->111       -0.23610
   103 ->105       -0.10275
   103 ->108       0.15328
   104 ->106       -0.16053
   104 ->109       0.15605
   104 ->111       -0.16076

Excited State 5:    Singlet-A    8.5046 eV  145.78 nm  f=0.0023  <S**2>=0.000
   102 ->105       0.11831
   102 ->106       0.11282
   103 ->105       0.27752
   103 ->106       0.35216
   103 ->109       -0.10927
   103 ->110       -0.13444
   103 ->111       0.18113
   104 ->106       0.10511
   104 ->108       -0.28768
   104 ->111       0.18629

Excited State 6:    Singlet-A    8.7086 eV  142.37 nm  f=0.0028  <S**2>=0.000
   102 ->106       0.43804
   102 ->107       -0.15157
   102 ->110       0.27393
   102 ->113       -0.11435
   102 ->117       -0.11816
   103 ->106       -0.14354
   103 ->110       -0.10372
   104 ->106       0.13551
   104 ->110       0.16743

请问下问题出现在哪里？  怎样给处理。谢谢！

作者
Author: 小范范1989 时间: 2016-6-4 16:46
我的感觉哈
#p opt freq m062x/TZVP td=(singlets,root=1) integral=ultrafine
优化和freq不能同时进行吧？
可能这里的错误是不收敛，但是opt td和freq能同时计算吗？

作者
Author: HERO 时间: 2016-6-4 16:51

小范范1989 发表于 2016-6-4 16:46
我的感觉哈
#p opt freq m062x/TZVP td=(singlets,root=1) integral=ultrafine
优化和freq不能同时进行吧 ...

在做别的体系的时候有这样设置。我试下吧，先优化下，在频率计算。谢谢哈

作者
Author: 小范范1989 时间: 2016-6-4 16:59

HERO 发表于 2016-6-4 16:51
在做别的体系的时候有这样设置。我试下吧，先优化下，在频率计算。谢谢哈

我很水哈，我只是之前opt td和freq同时计算会出错。
基态的 opt freq没问题。
我说的不一定对哈

作者
Author: HERO 时间: 2016-6-4 17:09

小范范1989 发表于 2016-6-4 16:59
我很水哈，我只是之前opt td和freq同时计算会出错。
基态的 opt freq没问题。
我说的不一定对哈

呵呵我以前也是这样的基态opt fre一起，激发态分开，后来发现有的体系激发态可以一齐算。试试吧

作者
Author: liyuanhe211 时间: 2016-6-4 20:59
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: drdavis 时间: 2016-6-4 22:50

小范范1989 发表于 2016-6-4 16:46
我的感觉哈
#p opt freq m062x/TZVP td=(singlets,root=1) integral=ultrafine
优化和freq不能同时进行吧 ...

好像是gaussian09 A版本不能一起算，貌似从B开始修复了这个bug

作者
Author: sobereva 时间: 2016-6-5 00:14
凡是耗时较多的任务不要搞在一起算
当前体系，TDDFT用有限差分做freq是极为耗时的，显然应当opt和freq分成两个任务来做，

但当前出错这个这没关系，凡是看到类似报错，一定要看6L提供的链接，没有其它任何有用的解决方法。

激发态优化和振动分析完全没必要用TZVP级别的基组，6-31G*足矣，基组小了收敛通常更容易。

作者
Author: HERO 时间: 2016-6-5 08:47
本帖最后由 HERO 于 2016-6-5 08:56 编辑

sobereva 发表于 2016-6-5 00:14
凡是耗时较多的任务不要搞在一起算
当前体系，TDDFT用有限差分做freq是极为耗时的，显然应当opt和freq分成 ...

谢谢sob，请问下6L提供的链接是啥意思？试下分开做，把基组调整小下。

作者
Author: HERO 时间: 2016-6-5 10:24

liyuanhe211 发表于 2016-6-4 20:59
解决SCF不收敛问题的方法
http://sobereva.com/61

谢谢哈学习中

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