Optimiser une géométrie de molécule

Pour une optimisation devez entrer des commandes qui concernent le groupe de Controle général (CONTRL) et pour le groupe des Points Stationnaires (STATPT):

Pour optimiser vers un minimum met RUNTYP=OPTIMIZE au lieu de RUNTYP=ENERGY. Vers un état de transition on met RUNTYP=SADPOINT. On défini aussi le seuil d'optimization (OPTTOL) qui par defaut est de 0.0001 et le nombre de pas maximum (NSTEP) qui par défaut vaut 20.

Il est une bonne pratique d'ajouter un calcul de frequence en fin d'optimisation HSSEND=.TRUE. et le groupe qui défini le calcul du hessien $FORCE METHOD=analytic $END

  • $CONTRL RUNTYP=OPTIMIZE $END
  • $STATPT OPTTOL=0.0005 NSTEP=999 HSSEND=.TRUE. $END
  • $FORCE METHOD=analytic $END

! GAMESS input file for geometry optimization
 $CONTRL
  SCFTYP=RHF MAXIT=200
  RUNTYP=OPTIMIZE
  MULT=1   iCHARG=0
 $END

 $STATPT OPTTOL=0.0005 NSTEP=999 HSSEND=.TRUE. $END

 $BASIS GBASIS=N31 NGAUSS=6
  NDFUNC=1 NPFUNC=0 $END
 
 $DATA
 
  C1
  H 1.0 -0.00000 0.00000 -0.35419
  H 1.0 0.00000 -0.00000 0.35419
  $END
 ! Modalites de calcul (ne pas toucher a priori)
 $SYSTEM TIMLIM=2879 MWORDS=250 $END

 $SCF DIRSCF=.T. $END
 $GUESS GUESS=HUCKEL $END