 |
|
|
$EOMINP inputs.
Any publication describing the results of CC calculations
obtained using GAMESS should reference the appropriate
papers, which are listed on the output of every run, and in
chapter 4 of this manual.
Analytic gradients are not available, so use CCTYP only for
RUNTYP=ENERGY, TRUDGE, SURFACE, or maybe FFIELD, or request
numerical derivatives.
Generally speaking, the Renormalized energies are obtained
at similar cost to the standard values, while Completely
Renormalized energies cost twice the time. For usage tips
and more information about resources on the various Coupled
Cluster methods, see Section 4, 'Further Information'.
CIMTYP chooses a Cluster-In-Molecule (CIM) calculation.
= NONE skip CIM computation, i.e., perform
a canonical calculation (default).
= SECIM perform a single-environment CIM (SECIM)
computation.
= DECIM perform a dual-environment CIM (DECIM)
computation.
= GSECIM perform a generalized SECIM (GSECIM)
computation. The $CIMFRG must be
included as well.
See also $CIMINP and, optionally, $CIMFRG and $CIMATM.
If CIMTYP is given, SUBMTD in $CIMINP is required. Only
RUNTYP=ENERGY and SCFTYP=RHF or ROHF work when CIMTYP is
given. See SUBMTD in $CIMINP for more details.
* * * * *
RELWFN Selects all-electron scalar relativity treatment.
See the $RELWFN input group for more information,
including nuclear derivative availability.
= NONE use the basic Schrodinger equation (default)
= LUT-IOTC local unitary transformation modification
of IOTC, due to H.Nakai, J.Seino, Y.Nakajima.
This is the fastest and most numerically
reliable scalar relativity method, so it is
preferred over RESC, DK, or IOTC.
= IOTC infinite-order two-component method of
M. Barysz and A.J. Sadlej.
= DK Douglas-Kroll transformation, available at
the 1st, 2nd, or 3rd order.
= RESC relativistic elimination of small component,
the method of T. Nakajima and K. Hirao,
available at 2nd order only.
= NESC normalised elimination of small component,
the method of K. Dyall, 2nd order only.
* * * * *
RUNTYP specifies the type of computation, for
example at a single geometry point:
= ENERGY Molecular energy. (default)
= GRADIENT Molecular energy plus gradient.
= HESSIAN Molecular energy plus gradient plus
second derivatives, including harmonic
harmonic vibrational analysis.
See the $FORCE and $CPHF input groups.
For FMO, use FMOHESS instead of HESSIAN.
= FMOHESS the same as HESSIAN, for FMO runs,
supported only for RHF, R-DFT, UHF,
U-DFT, and ROHF.
= GAMMA Evaluate up to 3rd nuclear derivatives,
by finite differencing of Hessians.
See $GAMMA, and also NFFLVL in $CONTRL.
multiple geometry options:
= OPTIMIZE Optimize the molecular geometry using
analytic energy gradients. See $STATPT.
= TRUDGE Non-gradient total energy minimization.
See $TRUDGE and $TRURST.
= SADPOINT Locate saddle point (transition state).
See $STATPT.
= MEX Locate minimum energy crossing point on
the intersection seam of two potential
energy surfaces. See $MEX.
= CONICAL Locate conical intersection point on
the intersection seam of two potential
energy surfaces. See $CONICL.
= IRC Follow intrinsic reaction coordinate.
See $IRC.
= VSCF anharmonic vibrational corrections.
See $VSCF.
= DRC Follow dynamic reaction coordinate.
See $DRC.
= MD molecular dynamics trajectory, see $MD.
= GLOBOP Monte Carlo-type global optimization.
See $GLOBOP.
= OPTFMO genuine FMO geometry optimization using
nearly analytic gradient. See $OPTFMO.
= GRADEXTR Trace gradient extremal. See $GRADEX.
= SURFACE Scan linear cross sections of the
potential energy surface. See $SURF.
single geometry property options:
= COMP composite thermochemistry calculation,
including G3MP2. See $COMP input.
= G3MP2 evaluate heat of formation using the
G3(MP2,CCSD(T)) methodology. See test
example exam43.inp for more information.
= PROP Molecular properties will be calculated.
Orbital localization can be requested as
well. See $ELPOT, etc.
Converged orbitals must be input in a
$VEC input, which suffice to reproduce
the wavefunction only for simple SCF:
RHF, UHF, ROHF, or DFT counterparts.
GVB also works (CICOEF may be needed).
All other calculations must instead use
RUNTYP=ENERGY to regenerate the density
matrix.
= RAMAN computes Raman intensities, see $RAMAN.
= NACME non-adiabatic coupling matrix element
between two or more state averaged MCSCF
wavefunctions. The calculation has no
specific input group, but must use only
SCFTYP=MCSCF with CISTEP=ALDET or ORMAS.
= NMR NMR shielding tensors for closed shell
molecules by the GIAO method. See $NMR.
= EDA Perform energy decomposition analysis.
Give one of $MOROKM or $LMOEDA inputs.
= QMEFPEA QM/EFP solvent energy analysis,
see $QMEFP.
= TRANSITN Compute radiative transition moment or
spin-orbit coupling. See $TRANST.
= FFIELD applies finite electric fields, most
commonly to extract polarizabilities.
See $FFCALC.
= TDHF analytic computation of time dependent
polarizabilities. See $TDHF.
= TDHFX extended TDHF package, including nuclear
polarizability derivatives, and Raman
and Hyper-Raman spectra. See $TDHFX.
= MAKEFP creates an effective fragment potential,
for SCFTYP=RHF or ROHF only.
See $MAKEFP, $DAMP, $DAMPGS, $STONE, ...
= FMO0 performs the free state FMO calculation.
See $FMO.
* * * * * * * * * * * * * * * * * * * * * * * * * * * * *
Note that RUNTYPs which require the nuclear gradient are
GRADIENT, HESSIAN, OPTIMIZE, SADPOINT,
GLOBOP, IRC, GRADEXTR, DRC, and RAMAN
These are efficient with analytic gradients, which are
available only for certain CI or MP2 calculations, but no
CC calculations, as indicated above. See NUMGRD.
* * * * * * * * * * * * * * * * * * * * * * * * * * * * *
NUMGRD Flag to allow numerical differentiation
of the energy. Each gradient requires
the energy be computed twice (forward
and backward displacements) along each
totally symmetric modes. It is thus
recommended only for systems with just a
few symmetry unique atoms in $DATA.
The default is .FALSE.
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.
Note that you must set PARALL=.TRUE.
in $SYSTEM to test distributed memory
allocations.
= 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.
* * * * * * *
ICHARG = Molecular charge. (default=0, neutral)
MULT = Multiplicity of the electronic state
= 1 singlet (default)
= 2,3,... doublet, triplet, and so on.
ICHARG and MULT are used directly for RHF, UHF, ROHF.
For GVB, these are implicit in the $SCF input, while
for MCSCF or CI, these are implicit in $DRT/$CIDRT or
$DET/$CIDET input. You must still give them correctly.
* * * the next three control molecular geometry * * *
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.
= PRINAXIS Cartesian coordinates will be input,
and transformed to principal axes.
Please read the warning just below!!!
= ZMT GAUSSIAN style internals will be input.
= ZMTMPC MOPAC style internals will be input.
= FRAGONLY means no part of the system is treated
by ab initio means, hence $DATA is not
given. The system is defined by $EFRAG.
Note: the choices PRINAXIS, ZMT, ZMTMPC require input of
all atoms in the molecule. They also orient the molecule,
and then determine which atoms are unique. The
reorientation is likely to change the order of the atoms
from what you input. When the point group contains a 3-
fold or higher rotation axis, the degenerate moments of
inertia often cause problems choosing correct symmetry
unique axes, in which case you must use COORD=UNIQUE rather
than Z-matrices.
Warning: The reorientation into principal axes is done
only for atomic coordinates, and is not applied to the axis
dependent data in the following groups: $VEC, $HESS, $GRAD,
$DIPDR, $VIB, nor Cartesian coords of effective fragments
in $EFRAG. COORD=UNIQUE avoids reorientation, and thus is
the safest way to read these.
Note: the choices PRINAXIS, ZMT, ZMTMPC require the use
of a group named $BASIS to define the basis set. The first
two choices might or might not use $BASIS, as you wish.
UNITS = distance units, any angles must be in degrees.
= ANGS Angstroms (default)
= BOHR Bohr atomic units
NZVAR = 0 Use Cartesian coordinates (default).
= M If COORD=ZMT or ZMTMPC, and $ZMAT is not given:
the internal coordinates will be those defining
the molecule in $DATA. In this case, $DATA may
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 more sophisticated
internal coordinate choices. M is ordinarily
3N-6 (3N-5), unless $ZMAT has linear bends.
NZVAR refers mainly to the coordinates used by OPTIMIZE
or SADPOINT runs, but may also print the internal's
values for other run types. You can use internals to
define the molecule, but Cartesians during optimizations!
* * * * * * *
Pseudopotentials may be of two types: ECP (effective core
potentials) which generate nodeless valence orbitals, and
MCP (model core potentials) producing valence orbitals with
the correct radial nodal structure. At present, ECPs have
analytic nuclear gradients and Hessians, while MCPs have
analytic nuclear gradients.
PP = pseudopotential selection.
= NONE all electron calculation (default).
= READ read ECP potentials in the $ECP input.
= SBKJC use Stevens, Basch, Krauss, Jasien,
Cundari ECP potentials for all heavy
atoms (Li-Rn are available).
= HW use Hay, Wadt ECP potentials for heavy
atoms (Na-Xe are available).
= MCP use Huzinaga's Model Core Potentials.
The correct MCP potential will be chosen
to match the requested MCP valence basis
set (see $BASIS).
* * * * * * *
LOCAL = controls orbital localization.
= NONE Skip localization (default).
= BOYS Do Foster-Boys-like localization.
= RUEDNBRG Do Edmiston-Ruedenberg localization.
= POP Do Pipek-Mezey population localization.
= SVD Do single value decomposition, to project
the molecular orbitals onto atoms. This
is available only for SCFTYP=RHF, ROHF,
and MCSCF (full space or ORMAS). The
ORIENT keyword in $LOCAL is pertinent.
See the related $LOCAL input.
Localization is not available for SCFTYP=GVB.
DFTB only works with LOCAL=POP (and NONE).
* * * * * * *
ISPHER = Spherical Harmonics option
= -1 Use Cartesian basis functions to construct
symmetry-adapted linear combination (SALC)
of basis functions. The SALC space is the
linear variation space used. (default)
= 0 Use spherical harmonic functions to create
SALC functions, which are then expressed
in terms of Cartesian functions. The
contaminants are not dropped, hence this
option has EXACTLY the same variational
space as ISPHER=-1. The only benefit to
obtain from this is a population analysis
in terms of pure s,p,d,f,g functions.
= +1 Same as ISPHER=0, but the function space
is truncated to eliminate all contaminant
Cartesian functions [3S(D), 3P(F), 4S(G),
and 3D(G)] before constructing the SALC
functions. The computation corresponds
to the use of a spherical harmonic basis.
QMTTOL = linear dependence threshhold
Any functions in the SALC variational space whose
eigenvalue of the overlap matrix is below this
tolerence is considered to be linearly dependent.
Such functions are dropped from the variational
space. What is dropped is not individual basis
functions, but rather some linear combination(s)
of the entire basis set that represent the linear
dependent part of the function space. The default
is a reasonable value for most purposes, 1.0E-6.
When many diffuse functions are used, it is common
to see the program drop some combinations. On
occasion, in multi-ring molecules, we have raised
QMTTOL to 3.0E-6 to obtain SCF convergence, at the
cost of some energy.
MAXIT = Maximum number of SCF iteration cycles. This
pertains only to RHF, UHF, ROHF, or GVB runs.
See also MAXIT in $MCSCF. (default = 30)
* * * interfaces to other programs * * *
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, see the URL
http://www.chemistry.mcmaster.ca/aimpac
DGRID = flag to add extra digits in molecular orbitals to
the log file for use by Kohout's DGrid program:
http://www2.cpfs.mpg.de/~kohout/dgrid.html
This is one of the modern alternatives to the old
AIMPAC codes, in the QTAIM/ELF arena.
(default .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
PLTORB, MOLPLT, and AIMPAC decks are written to file
PUNCH at the end of the job. Thus all of these correspond
to the final geometry encountered during jobs such as
OPTIMIZE, SAPDOINT, IRC...
In contrast, selecting FRIEND turns the job into a
CHECK run only, no matter how you set EXETYP. Thus the
geometry is that encountered in $DATA. The input is
added to the PUNCH file, and may require some (usually
minimal) massaging.
PLTORB and MOLPLT are written even for EXETYP=CHECK.
AIMPAC requires at least RUNTYP=PROP.
* * *
NFFLVL used to determine energies and gradients away
from equilibrium structures, at the coordinates
given in $DATA. The method will use a Taylor
expansion of the potential surface around the
stationary point. See $EQGEOM, $HLOWT, $GLOWT.
This may be used with RUNTYP=ENERGY or GRADIENT.
= 2 uses only Hessian information, which gives a
reasonable energy, but not such a good gradient.
= 3 uses Hessian and 3rd nuclear derivatives in the
Taylor expansion, producing more accurate values
for the energy and for the gradient.
* * * computation control switches * * *
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 wider 132 columns 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.
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 input
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.
ETOLLZ = threshold to label molecular orbitals by Lz
values. Small matrices of the Lz operator are
diagonalized for the sets of MOs whose orbital
energies are degenerate to within ETOLLZ. This
option may be used in molecules with distorted
linear symmetry for approximate labelling.
Default: 1.0d-6 for linear, 0 (disable) if not.
INTTYP selects the integral package(s) used, all of which
produce equally accurate results. This is therefore
used only for debugging purposes.
= BEST use the fastest integral code available for
any particular shell quartet (default):
s,p,L or s,p,d,L rotated axis code first.
ERIC s,p,d,f,g precursor transfer equation
code second, up to 5 units total ang. mom.
Rys quadrature for general s,p,d,f,g,L,
or for uncontracted quartets.
= ROTAXIS means don't use ERIC at all, e.g. rotated
axis codes, or else Rys quadrature.
= ERIC means don't use rotated axis codes, e.g.
ERIC code, or else Rys quadrature.
= RYSQUAD means use Rys quadrature for everything.
GRDTYP = BEST use Schlegel routines for spL gradient
blocks, and Rys quadrature for all
other gradient integrals. (default)
= RYSQUAD use Rys quadrature for all gradient
integrals. This option is only slightly
more accurate, but is rather slower.
HSSTYP = BEST use faster routines.
= GENERAL use slower code (default).
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 exponential
factor is less than 10**(-n) are skipped.
ICUT = n integrals less than 10.0**(-n) are not
saved on disk. (default = 9). Direct
SCF will calculate to a cutoff 1.0d-10
or 5.0d-11 depending on FDIFF=.F. or .T.
ISKPRP = 0 proceed as usual
1 skip computation of some properties which
are not well parallelised. This includes
bond orders and virial theorem, and can help
parallel scalability if many CPUs are used.
Note that NPRINT=-5 disables most property
computations as well, so ISKPRP=1 has no
effect in that case. (default: 0)
* * * restart options * * *
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 integrals 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!
==========================================================
==========================================================
generated on 7/7/2017