$SCFMI group          (optional, relevant if SCFTYP=RHF)
 
 
    The Self Consistent Field for Molecular Interactions
(SCF-MI) method is a modification of the usual Roothaan
equations that avoids basis set superposition error (BSSE)
in intermolecular interaction calculations, by expanding
each monomer's orbitals using only its own basis set.
Thus, the resulting orbitals are not orthogonal.  The
presence of a $SCFMI group in the input triggers the use
of this option.
 
    The implementation is limited to ten monomers, treated
at the RHF level.  The energy, gradient, and therefore
semi-numerical hessian are available.  The SCF step may be
run in direct SCF mode, and parallel calculation is also
enabled.  The calculation must use Cartesian Gaussian AOs
only, not spherical harmonics.  The SCF-MI driver differs
from normal RHF calculations, so not all converger methods
are available.  Finally, this option is not compatible with
electron correlation treatments (DFT, MP2, CI, or CC).
 
    The first 3 parameters must be given.  All atoms of a
fragment must appear consecutively in $DATA.
 
NFRAGS    = number of distinct fragments present.  Both
            the supermolecule and its constituent monomers
            must be well described as closed shells by RHF
            wavefunctions.
 
NF        = an array containing the number of doubly
occupied
            MOs for each fragment.
 
MF        = an array containing the number of atomic basis
            functions located on each fragment.
 
ITER      = maximum number of SCF-MI cycles, overriding
            the usual MAXIT value.  (default is 50).
 
DTOL      = SCF-MI density convergence criteria.
            (default is 1.0d-10)
 
ALPHA     = possible level shift parameter.
            (default is 0.0, meaning shifting is not used)
 
DIISON    = a flag to active the DIIS convergence.
            (default is .TRUE.)
 
MXDIIS    = the maximum number of previous effective Fock
and
            overlap matrices to be used in DIIS
(default=10)
 
DIISTL    = the density change value at which DIIS starts.
            (default=0.01)
 
A Huckel guess is localized by the Boys procedure onto each
fragment to provide starting orbitals for each:
 
ITLOC     = maximum number of iteration in the localization
            step (Default is 50)
 
CNVLOC    = convergence parameter for the localization.
            (default is .01).
 
IOPT      =   prints additional debug information.
          = 0 standard outout (default)
          = 1 print for each SCF-MI cycle MOs, overlap
              between the MOs, CPU times.
          = 2 print some extra informations in secular
              systems solution.
 
   "Modification of Roothan Equations to exclude BSSE
       from Molecular Interaction calculations"
    E. Gianinetti, M. Raimondi, E. Tornaghi
    Int. J. Quantum Chem. 60, 157-166 (1996)
 
   "Implementation of Gradient optimization algorithms
     and Force Constant computations in BSSE-free direct
     and conventional SCF approaches"
A. Famulari, E. Gianinetti, M. Raimondi, M. Sironi
    Int. J. Quantum Chem. 69, 151-158 (1997)