Unconverge fragment SCF resulted to an absurd energy in NEDA

Unconverge fragment SCF resulted to an absurd energy in NEDA

Postby adebayo » Mon Aug 13, 2018 1:15 pm

Greetings.

I am trying to run NEDA analysis on four molecules that are very similar. Each were fragmented into 3 but one of them did not have the frament 2 converged. Because the SCF of fragment 2 did not converged, it resulted into an absurd value of interaction energy.

Here is the input:
Code: Select all
 $SYSTEM MWORDS=80 kdiag=0  TIMLIM=100000 $END
 $CONTRL  RUNTYP=ENERGY NOSYM=1  MAXIT=200 $END
 $CONTRL SCFTYP=RHF  DFTTYP=B3LYP $END
 $p2p p2p=.t. dlb=.t. $end
 $INTGRL NOPK=1 NWORD=3000000 $END
 $BASIS GBASIS=N311 NGAUSS=6 $END
 $BASIS NDFUNC=2 NDFUNC=1 NPFUNC=2 $END
 $BASIS DIFFSP=.TRUE. $END
 $TRANS nosym=1 $END
 $CONTRL COORD=UNIQUE UNITS=ANGSTROM $END
 $SCF NCONV=5 DIRSCF=.t. FDIFF=.F. DIIS=.T. SOSCF=.F. $END
 $STATPT  METHOD=GDIIS  NSTEP=100000   $END
 $ELDENS  IEDEN=1  $END
 $ELPOT   IEPOT=1  $END
 $ELFLDG IEFLD=2 $END
 $CUBE    CUBE=.T. mesh=coarse $END
 $CONTRL ICHARG=0 MULT=1 $END
 $NBO $END
 $DEL
NAO
NEDA
(3,4,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36)
(1,2,5,6,37)
(38,39,40,41,42,43,44,45,46,47,48)
END
 $END

 $DATA
NEDA analysis of molecule
C1
O      8.0      1.5470000000    0.7460000000   -1.0030000000
O      8.0      2.2740000000    0.0050000000    0.9820000000
N      7.0     -5.1960000000   -0.8690000000   -0.2070000000
N      7.0     -3.7650000000    1.0210000000    0.4620000000
N      7.0      3.6860000000    0.3680000000   -0.7710000000
H      1.0      3.6530000000    0.6760000000   -1.7160000000
C      6.0     -2.8870000000   -1.2320000000    0.2360000000
C      6.0     -4.1380000000   -0.3900000000    0.6500000000
H      1.0     -4.4250000000   -0.5430000000    1.6860000000
C      6.0     -3.3490000000   -1.9720000000   -1.0470000000
H      1.0     -3.6210000000   -2.9940000000   -0.8050000000
H      1.0     -2.5770000000   -2.0000000000   -1.8060000000
C      6.0     -1.8540000000   -0.1620000000   -0.0450000000
C      6.0     -4.5960000000   -1.2110000000   -1.4810000000
H      1.0     -4.3390000000   -0.3280000000   -2.0710000000
H      1.0     -5.2770000000   -1.8160000000   -2.0680000000
C      6.0     -2.4120000000    1.0990000000    0.1310000000
C      6.0     -2.4370000000   -2.2050000000    1.3290000000
H      1.0     -1.6030000000   -2.8110000000    0.9880000000
H      1.0     -2.1190000000   -1.6760000000    2.2210000000
H      1.0     -3.2490000000   -2.8750000000    1.5960000000
C      6.0     -6.4060000000   -0.0910000000   -0.2710000000
H      1.0     -6.2610000000    0.8980000000   -0.7050000000
H      1.0     -6.8230000000    0.0270000000    0.7240000000
H      1.0     -7.1390000000   -0.6210000000   -0.8700000000
C      6.0     -0.5340000000   -0.3000000000   -0.4100000000
H      1.0     -0.0850000000   -1.2670000000   -0.5490000000
C      6.0     -4.2950000000    1.9960000000    1.3900000000
H      1.0     -4.2490000000    2.9900000000    0.9610000000
H      1.0     -3.7560000000    2.0070000000    2.3390000000
H      1.0     -5.3340000000    1.7800000000    1.5900000000
C      6.0     -1.6430000000    2.2380000000   -0.0350000000
H      1.0     -2.0570000000    3.2200000000    0.1000000000
C      6.0      0.2310000000    0.8430000000   -0.5780000000
C      6.0     -0.3080000000    2.0960000000   -0.3920000000
H      1.0      0.3140000000    2.9610000000   -0.5310000000
C      6.0      2.4870000000    0.3370000000   -0.1390000000
C      6.0      4.9520000000    0.0130000000   -0.2690000000
C      6.0      6.0250000000    0.1060000000   -1.1510000000
H      1.0      5.8580000000    0.4350000000   -2.1630000000
C      6.0      5.1740000000   -0.4120000000    1.0350000000
H      1.0      4.3590000000   -0.4870000000    1.7230000000
C      6.0      7.3020000000   -0.2210000000   -0.7380000000
H      1.0      8.1180000000   -0.1430000000   -1.4350000000
C      6.0      6.4610000000   -0.7370000000    1.4330000000
H      1.0      6.6220000000   -1.0660000000    2.4440000000
C      6.0      7.5310000000   -0.6470000000    0.5600000000
H      1.0      8.5240000000   -0.9030000000    0.8840000000
 $END



Here is part of the output of NEDA analysis
Code: Select all
--------------
  Fragment  2:
 --------------
NONZERO     BLOCKS
 ITER EX DEM     TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR      INTEGRALS    SKIPPED
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0  0     -242.0568830148  -242.05688301482216.770834457   0.159014489    10879861458  498197238
...
199198  0     -243.7646092885    -0.0000125487   0.547739012   0.000957450    11083207920  488146261
 200199  0     -243.7647175519    -0.0001082634   1.615013602   0.000957450    11078864774  488287964

 SCF IS UNCONVERGED, TOO MANY ITERATIONS
     TIME TO FORM FOCK OPERATORS=   94524.7 SECONDS (     472.6 SEC/ITER)
     FOCK TIME ON FIRST ITERATION=     454.4, LAST ITERATION=     469.7
     TIME TO SOLVE SCF EQUATIONS=     375.9 SECONDS (       1.9 SEC/ITER)

 FINAL R-B3LYP ENERGY IS        0.0000000000 AFTER 200 ITERATIONS

...
--------------
  Fragment  3:
 --------------

 NONZERO     BLOCKS
 ITER EX DEM     TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR      INTEGRALS    SKIPPED
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0  0     -230.8083678281  -230.8083678281   2.205980575   0.319478089    10127693385  509190885
  14 13  0     -231.1243401628    -0.0000000001   0.000002771   0.000000392     9911682596  510992004

          -----------------
          DENSITY CONVERGED
          -----------------
     TIME TO FORM FOCK OPERATORS=    5880.2 SECONDS (     420.0 SEC/ITER)
     FOCK TIME ON FIRST ITERATION=     429.3, LAST ITERATION=     418.7
     TIME TO SOLVE SCF EQUATIONS=      26.1 SECONDS (       1.9 SEC/ITER)

 FINAL R-B3LYP ENERGY IS     -231.1243401628 AFTER  14 ITERATIONS

...

Natural Energy Decomposition Analysis (Summary):

                                                            Component
                    Energy(wfn)        Energy(wfn)          (kcal/mol)
 ------------------------------------------------------------------------------
 C20H23N3O2    -1090.3425815(scf) -1088.4885249(loc)      CT =-1163.438
                                                          ES = -657.352
                                                         POL = -157.716
                                                          XC =  -83.237
 1. C13H17N2(+  -614.1917321(def)  -614.4987244(cp)  DEF(SE) =  192.641( 48.382)
 2. CHNO2(-2)   -242.0568830(def)     0.0000000(cp)  DEF(SE) =*********( 53.759)
 3. C6H5(+)     -230.8083678(def)  -231.1243402(cp)  DEF(SE) =  198.276( 47.169)
                                                              ---------
                                                           E =*********


 Electrical (ES+POL+SE) :   -665.758
   Charge Transfer (CT) :  -1163.438
       Core (XC+DEF-SE) :-151734.585
                        ------------
  Total Interaction (E) :-153563.780



Kindly help.

Thanks.
adebayo
 
Posts: 5
Joined: Fri Feb 09, 2018 2:41 am

Re: Unconverge fragment SCF resulted to an absurd energy in

Postby ericg » Sun Aug 26, 2018 8:53 pm

GAMESS can often have convergence problems. I’d recommend trying DAMP=.T. in $SCF along with DIIS=.T. This often works for me. SOSCF generally has poor convergence behavior, so disabling it, as you have, is usually helpful.

Eric
ericg
 
Posts: 298
Joined: Sat Dec 29, 2012 9:31 am


Return to General NBO Discussion

Who is online

Users browsing this forum: No registered users and 1 guest

cron