This is a free format group specifying global switches.
SCFTYP together with MPLEVL specifies the wavefunction. You may choose from: = RHF Restricted Hartree Fock calculation (default for EVEN number of electrons) = UHF Unrestricted Hartree Fock calculation (default for ODD number of electrons) = ROHF Restricted open shell Hartree-Fock. (high spin, see GVB for low spin) = GVB Generalized valence bond wavefunction or OCBSE type ROHF. (needs $SCF input) = MCSCF Multiconfigurational SCF wavefunction (this requires $DRT input) = CI Configuration Interaction calculation (this requires $DRT input) MPLEVL = n chooses Moller-Plesset perturbation level. The default is 0 to skip. This is implemented only for n=2, only for RHF, UHF, and ROHF wave functions, and only for ENERGY, TRUDGE,SURFACE, and FFIELD run types. RUNTYP specifies the type of computation, for example at a single geometry point: = ENERGY Single point energy. (default) = GRADIENT Single point energy plus gradient. = HESSIAN Compute energy second derivatives, and perform harmonic vibrational analysis. See the $FORCE group. multiple geometry options: = OPTIMIZE Optimize the molecular geometry using analytic energy gradients. See $STATPT. = TRUDGE Non-gradient total energy minimization. See groups $TRUDGE and $TRURST. = SADPOINT Locate saddle point (transition state). See the $STATPT group. = IRC Follow intrinsic reaction coordinate. See the $IRC group. = GRADEXTR Trace gradient extremal. See the $GRADEX group. = DRC Follow dynamic reaction coordinate. See the $DRC group. = SURFACE Scan linear cross sections of the potential energy surface. See $SURF. single geometry property options: = PROP Properties will be calculated. A $DATA deck and converged $VEC group should be input. Optionally, Boys localization can be done. See $ELPOT, etc. = MOROKUMA Performs dimer energy decomposition. See the $MOROKM group. = TRANSITN Find radiative transition moment. See the $TRANST group. = SPINORBT Find spin-orbit coupling constant. See the $TRANST group. = FFIELD applies finite electric fields, most commonly to extract polarizabilities. See the $FFCALC group. = TDHF analytic computation of time dependent polarizabilities. See the $TDHF group.
EXETYP = RUN Actually do the run. (default) = CHECK Wavefunction and energy will not be evaluated. This lets you speedily check input and memory requirements. See the overview section for details. = DEBUG Massive amounts of output are printed, useful only if you hate trees. = routine Maximum output is generated by the routine named. Check the source for the routines this applies to. MAXIT = Maximum number of SCF iteration cycles. Pertains only to RHF, UHF, ROHF, or GVB runs. See also MAXIT in $MCSCF. (default = 30) ICHARG = Molecular charge. (default=0, neutral) MULT = Multiplicity of the electronic state = 1 singlet (default) = 2,3,... doublet, triplet, and so on.
ECP = effective core potential control. = NONE all electron calculation (default). = READ read the potentials in $ECP group. = SBK use Stevens, Basch, Krauss, Jasien, Cundari potentials for all heavy atoms (Li-Rn are available). = HW use Hay, Wadt potentials for all the heavy atoms (Na-Xe are available).
COORD = choice for molecular geometry in $DATA. = UNIQUE only the symmetry unique atoms will be given, in Cartesian coords (default). = HINT only the symmetry unique atoms will be given, in Hilderbrandt style internals. = CART Cartesian coordinates will be input. = ZMT GAUSSIAN style internals will be input. = ZMTMPC MOPAC style internals will be input.
UNITS = distance units, any angles which are entered in $DATA must be in degrees. = ANGS Angstroms (default) = BOHR Bohr atomic units NZVAR = Coordinate switch. = 0 Use Cartesian coordinates (default). = M If COORD=ZMT or ZMTMPC and a $ZMAT is not given: the internal coordinates will be those defining the molecule in $DATA. In this case, $DATA must not contain any dummy atoms. M is usually 3N-6 (or 3N-5 for linear). = M For other COORD choices, or if $ZMAT is given: the internal coordinates will be those defined in $ZMAT. This allows the use of more sophisticated internal coordinate choices. M is ordinarily 3N-6 (3N-5), unless linear bends are used in the $ZMAT.
LOCAL = controls orbital localization. = NONE Skip localization (default). = BOYS Do Foster-Boys localization. = RUEDNBRG Do Edmiston-Ruedenberg localization. = POP Do Pipek-Mezey population localization. See the $LOCAL group. Localization does not work for SCFTYP's GVB or CI.
MOLPLT = flag that produces an input deck for a molecule drawing program distributed with GAMESS. (default is .FALSE.) PLTORB = flag that produces an input deck for an orbital plotting program distributed with GAMESS. (default is .FALSE.) AIMPAC = flag to create an input deck for Bader's atoms in molecules properties code. (default=.FALSE.) For information about this program, contact Richard F.W. Bader Dept. of Chemistry McMaster University Hamilton, Ontario L8S-4M1 Canada bader@sscvax.cis.mcmaster.ca RPAC = flag to create the input files for Bouman and Hansen's RPAC electronic excitation and NMR shieldings program. RPAC works only with RHF wavefunctions. Contact Prof. Aage Hansen in Copenhagen (nahaeh@vm.uni-c.dk) about this program. (default is .FALSE.) FRIEND = string to prepare input to other quantum programs, choose from = HONDO for HONDO 8.2 = MELDF for MELDF = GAMESSUK for GAMESS (UK Daresbury version) = GAUSSIAN for Gaussian 9x = ALL for all of the above
For the most part, the default is the only sensible value, and unless you are sure of what you are doing, these probably should not be touched.
NPRINT = Print/punch control flag See also EXETYP for debug info. (options -7 to 5 are primarily debug) = -7 Extra printing from Boys localization. = -6 debug for geometry searches = -5 minimal output = -4 print 2e-contribution to gradient. = -3 print 1e-contribution to gradient. = -2 normal printing, no punch file = 1 extra printing for basis,symmetry,ZMAT = 2 extra printing for MO guess routines = 3 print out property and 1e- integrals = 4 print out 2e- integrals = 5 print out SCF data for each cycle. (Fock and density matrices, current MOs = 6 same as 7, but narrow terminal output This option isn't perfect. = 7 normal printing and punching (default) = 8 more printout than 7. The extra output is (AO) Mulliken and overlap population analysis, eigenvalues, Lagrangians, ... = 9 everything in 8 plus Lowdin population analysis, final density matrix.
IREST = restart control options (for OPTIMIZE run restarts, see $STATPT) Note that this option is unreliable! = -1 reuse dictionary file from previous run, useful with GEOM=DAF and/or GUESS=MOSAVED. Otherwise, this option is the same as 0. = 0 normal run (default) = 1 2e restart (1-e integrals and MOs saved) = 2 SCF restart (1-,2-e integrls and MOs saved) = 3 1e gradient restart = 4 2e gradient restart GEOM = select where to obtain molecular geometry = INPUT from $DATA input (default for IREST=0) = DAF read from DICTNRY file (default otherwise)
As noted in the first chapter, binary file restart is not a well tested option!
NOSYM = 0 the symmetry specified in $DATA is used as much as possible in integrals, SCF, gradients, etc. (this is the default) = 1 the symmetry specified in the $DATA group is used to build the molecule, then symmetry is not used again. Some GVB or MCSCF runs (those without a totally symmetric charge density) require you request no symmetry. INTTYP = POPLE use fast Pople routines for sp integral blocks, and HONDO Rys polynomial code for all other integrals. (default) = HONDO use HONDO/Rys integrals for all integrals. This option produces slightly more accurate integrals but is also slower. NORMF = 0 normalize the basis functions (default) = 1 no normalization NORMP = 0 input contraction coefficients refer to normalized Gaussian primitives. (default) = 1 the opposite. ITOL = primitive cutoff factor (default=20) = n products of primitives whose preexponential factor is less than 10**(-n) are skipped. ICUT = n integrals less than 10.0**(-n) are not saved on disk. (default = 9)