[AMBER] QM/MM Convergence questions

From: 辛志宏 <xzhfood.njau.edu.cn>
Date: Thu, 5 Mar 2020 17:07:44 +0800 (GMT+08:00)

Dear Amber,

An error occur when I run a QM/MM with amber and gausssin09, I do not know how to solve the problems, any suggestion will be much appreciate, thank you for your help in advanced, the error information is as follows:

The QM/MM input file:

Initial min of our structure QMMM
&cntrl
  imin = 1, maxcyc = 5000, ncyc=2000,
  ntb = 0,
  cut = 20.0, ntc=2, ntf=2,
  ntpr=1000,
  restraint_wt=2.0,
  ifqnt = 1
&end
&qmmm
iqmatoms=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,715,716,717,718,719,720,721,722,723,724,725,757,758,759,760,761,762,763,764,765,766,767,768,769,770,771,772,773,774,775,776,777,778,2365,2366,2367,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2382,8628,8629,8630,8679,8680,8681
  qmcharge = 0,
  spin = 2,
  qm_theory = 'EXTERN',
  qmshake=1,
  qm_ewald=1, qm_pme=1,
  diag_routine=0,
  writepdb=1
  qmcut = 20.0
&end
&gau
  method = 'B3LYP',
  basis = '6-31G*',
  num_threads = 10,
  mem = '6GB'
&end


Error form the gau_job.log file,
.....................................
RMSDP=1.93D-04 MaxDP=3.33D-02 DE= 1.64D-04 OVMax= 8.63D-02

 Cycle 128 Pass 1 IDiag 1:
 RMSU= 2.75D-05 CP: 1.00D+00 -5.05D-01 -8.69D-01 2.19D+00
 E= -2424.60356615312 Delta-E= 0.001037414186 Rises=F Damp=F
 DIIS: error= 8.55D-04 at cycle 125 NSaved= 20.
 NSaved=20 IEnMin= 5 EnMin= -2424.61776698937 IErMin=19 ErrMin= 4.94D-04
 ErrMax= 8.55D-04 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.31D-04 BMatP= 1.14D-04
 IDIUse=3 WtCom= 9.97D-01 WtEn= 2.98D-03
 Rare condition: small coef for last iteration: 0.000D+00
 Coeff-Com: -0.702D-01-0.293D-01-0.544D-01-0.305D+00 0.621D+00-0.423D+00
 Coeff-Com: -0.423D-01 0.140D+00 0.716D-01 0.997D-01-0.177D+00-0.897D-01
 Coeff-Com: 0.104D+00-0.256D+00 0.537D+00-0.734D+00 0.809D-01 0.274D+00
 Coeff-Com: 0.127D+01-0.203D-01
 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
 Coeff-En: 0.000D+00 0.100D+01
 Coeff: -0.702D-01-0.293D-01-0.544D-01-0.305D+00 0.621D+00-0.423D+00
 Coeff: -0.423D-01 0.140D+00 0.716D-01 0.997D-01-0.177D+00-0.897D-01
 Coeff: 0.104D+00-0.256D+00 0.537D+00-0.734D+00 0.809D-01 0.274D+00
 Coeff: 0.127D+01-0.203D-01
 Gap= 0.009 Goal= None Shift= 0.000
 Gap= 0.010 Goal= None Shift= 0.000
 RMSDP=6.02D-05 MaxDP=1.08D-02 DE= 1.04D-03 OVMax= 2.84D-02

>>>>>>>>>> Convergence criterion not met.
 SCF Done: E(UB3LYP) = -2424.60356615 A.U. after 129 cycles
            NFock=128 Conv=0.60D-04 -V/T= 2.2388
 <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 1.5874 S= 0.8555
 <L.S>= 0.000000000000E+00
 KE= 1.957254069352D+03 PE=-1.340126621204D+04 EE= 5.052337417625D+03
 Annihilation of the first spin contaminant:
 S**2 before annihilation 1.5874, after 1.2129
 Convergence failure -- run terminated.
 Error termination via Lnk1e in /home/xzhfood/gaussian/g09/l502.exe at Thu Mar 5 03:47:42 2020.
 Job cpu time: 0 days 5 hours 47 minutes 3.7 seconds.
 File lengths (MBytes): RWF= 465 Int= 0 D2E= 0 Chk= 12 Scr= 1

Tianyou Xin
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Mar 05 2020 - 01:30:02 PST
Custom Search