Input PDF
Input PDF
2-1
(5 December 2014)
*********************************
*
*
* Section 2 - Input Description *
*
*
*********************************
Input Description
2-2
*
name
function
module:routine
-----------------------Molecule, basis set, wavefunction specification:
$CONTRL chemical control data
INPUTA:START
$SYSTEM computer related options
INPUTA:START
$BASIS
basis set
INPUTB:BASISS
$DATA
molecule, geometry, basis set
INPUTB:MOLE
$ZMAT
internal coordinates
ZMATRX:ZMATIN
$LIBE
linear bend coordinates
ZMATRX:LIBE
$SCF
HF-SCF wavefunction control
SCFLIB:SCFIN
$SCFMI
SCF-MI input control data
SCFMI :MIINP
$DFT
density functional theory
DFT
:DFTINP
$TDDFT
time-dependent DFT
TDDFT :TDDINP
$CIS
singly excited CI
CISGRD:CISINP
$CISVEC vectors for CIS
CISGRD:CISVRD
$MP2
2nd order Moller-Plesset
MP2
:MP2INP
$RIMP2
resolution of the identity MP2
RIMP2 :RIDRVR
$AUXBAS RI-MP2's basis set specification RIMP2 :RIDRVR
$CCINP
coupled cluster input
CCSDT :CCINP
$EOMINP equation of motion CC
EOMCC :EOMINP
$MOPAC
semi-empirical specification
MPCMOL:MOLDAT
$GUESS
initial orbital selection
GUESS :GUESMO
$VEC
orbitals
(formatted) GUESS :READMO
$MOFRZ
freezes MOs during SCF runs
EFPCOV:MFRZIN
$DFTB
DFTB input
DFTBLB:INPUT
$DFTBSK Slater-Koster table input
DFTBSK:SKTAB
Note that MCSCF and CI input is listed below.
Potential energy surface options:
$STATPT
$TRUDGE
$TRURST
$FORCE
$CPHF
$CPMCHF
$MASS
$HESS
$GRAD
$DIPDR
$ALPDR
$VIB
$VIB2
$VSCF
$VIBSCF
$GAMMA
$EQGEOM
STATPT:SETSIG
TRUDGE:TRUINP
TRUDGE:TRUDGX
HESS :HESSX
CPHF :CPINP
MCPCGX:MCPCGX
VIBANL:RAMS
HESS :FCMIN
HESS :EGIN
HESS :DDMIN
RAMAN :ADMIN
HESS :HSSNUM
HESS :HSSFUL
VSCF :VSCFIN
VSCF :VGRID
HESS :GAMMXX
HESS :FFCARX
Input Description
$HLOWT
$GLOWT
$IRC
$DRC
$MEX
$CONICL
$MD
$RDF
$GLOBOP
$GLBFRG
$GRADEX
$SURF
2-3
HESS :FFCARX
HESS :FFCARX
RXNCRD:IRCX
DRC
:DRCDRV
MEXING:MEXINP
MDEFP :MDX
MDEFP :RDFX
GLOBOP:GLOPDR
GLOBOP:GLOPDR
GRADEX:GRXSET
SURF :SRFINP
Interpretation, properties:
$LOCAL
$TRUNCN
$ELMOM
$ELPOT
$ELDENS
$ELFLDG
$POINTS
$GRID
$PDC
$MGC
$RADIAL
$MOLGRF
$STONE
$COMP
$RAMAN
$NMR
$MOROKM
$LMOEDA
$QMEFP
$FFCALC
$TDHF
$TDHFX
LOCAL :LMOINP
EFPCOV:TRNCIN
PRPLIB:INPELM
PRPLIB:INPELP
PRPLIB:INPELD
PRPLIB:INPELF
PRPLIB:INPPGS
PRPLIB:INPPGS
PRPLIB:INPPDC
PRPPOP:RADWFN
PARLEY:PLTMEM
PRPPOP:STNRD
COMP :COMPX
RAMAN :RAMANX
NMR
:NMRX
MOROKM:MOROIN
MOROKM:MMOEDIN
EFINP :QMEFPAX
FFIELD:FFLDX
TDHF :TDHFX
TDX:FINDTDHFX
Solvation models:
$EFRAG
use effective fragment potential
$FRAGNAME specifically named fragment pot.
$FRGRPL inter-fragment repulsion
$EWALD
Ewald sums for EFP electrostatics
$MAKEFP generate effective fragment pot.
$PRTEFP simplified EFP generation
$DAMP
EFP multipole screening fit
$DAMPGS initial guess screening params
$PCM
polarizable continuum model
EFINP :EFINP
EFINP :RDSTFR
EFINP :RDDFRL
EWALD :EWALDX
EFINP :EFPX
EFINP :PREFIN
CHGPEN:CGPINP
CHGPEN:CGPINP
PCM
:PCMINP
Input Description
$PCMGRD
$PCMCAV
$TESCAV
$NEWCAV
$IEFPCM
$PCMITR
$DISBS
$DISREP
$SVP
$SVPIRF
$COSGMS
$SCRF
2-4
PCMCV2:PCMGIN
PCM
:MAKCAV
PCMCV2:TESIN
PCM
:DISREP
PCM
:IEFDAT
PCMIEF:ITIEFIN
PCMDIS:ENLBS
PCMVCH:MORETS
SVPINP:SVPINP
SVPINP:SVPIRF
COSMO :COSMIN
SCRF :ZRFINP
ECPLIB:ECPPAR
MCPINP:MMPRED
INPUTB:RWFINP
PRPLIB:INPEF
INT2A :INTIN
QMFM :QFMMIN
TRANS :TRFIN
FMOIO :FMOMIN
FMOIO :FMOPIN
FMOIO :FMOXYZ
FMOGRD:OPTFMO
FMOIO :FMOLMO
FMOIO :FMOBON
FMOIO :EMINOU
FMOIO :EDIN
FMOGRD:RSTOPT
INPUTA:GDDINP
Polymer model:
$ELG
ELGLIB:ELGINP
DC SCF input
DC correlation method input
subsystem definition for SCF
subsystem definition for MP2/CC
restart data for DC-MP2
restart data for DC-CC
DCLIB
DCLIB
DCLIB
DCLIB
DCMP2
DCCC
:DCINP
:DCCRIN
:DFLCST
:DFLCST
:RDMPDC
:RDCCDC
Input Description
2-5
clusters in molecules
$CIMINP
$CIMATM
$CIMFRG
CIMINF:CIMINP
CIMINF:CIMINP
CIMINF:CIMPRT
GAMESS:WFNCI
ALDECI:DETINP
ALDECI:DETINP
ALGNCI:GCIINP
ALGNCI:GCIINP
ORMAS :FCINPT
CEEIS :CEEISIN
CEEIS :RDCEEIS
ALGNCI:GCIGEN
GMCPT :OSRDDAT
GMCPT :OSMKREF
GMCPT :OSMKREF
GMCPT :OSMKREF
FSODCI:SOCINP
GUGDRT:ORDORB
GUGDRT:ORDORB
MCSCF :MCSCF
MP2
:MRMPIN
DEMRPT:DMRINP
MCQDPT:MQREAD
EXCORR:GETEXC
IVOCAS:IVODRV
IVOCAS:ORBREAD
GUGSRT:GUGSRT
GUGEM :GUGAEM
GUGDGA:GUGADG
GUGDM :GUGADM
GUGDM2:GUG2DM
LAGRAN:CILGRN
TRFDM2:TRF2DM
DIAB:DIABINP
TRNSTN:TRNSTX
Input Description
$CONTRL
2-6
==========================================================
$CONTRL
group
(note:
= UHF
= ROHF
= GVB
= MCSCF
= NONE
Input Description
$CONTRL
2-7
DFTTYP = NONE
= XXXXXX
TDDFT
= NONE
= EXCITE
= SPNFLP
= POL
* * * * *
VBTYP
= NONE
= VB2000
MPLEVL =
= 0
= 2
Input Description
$CONTRL
2-8
* * * * *
CITYP
=
= NONE
= CIS
= SFCIS
= ALDET
= ORMAS
= FSOCI
= GENCI
= GUGA
PMTD1
CCTYP
Input Description
$CONTRL
2-9
= CCSD(TQ)
= CR-CC(Q)
Input Description
= CR-EOM
= CR-EOML
$CONTRL
2-10
ionization processes,
= IP-EOM2
ionized EOMCC with up to 2h1p
excitations (i.e., IP-EOMCCSD)
= IP-EOM3A ionized EOMCC with all 1h and 2h1p,
and active-space 3h2p excitations
(i.e., IP-EOMCCSDt)
= EA-EOM2
electron-attached EOMCC with up to 2p1h
excitations (i.e., EA-EOMCCSD)
= EA-EOM3A electron-attached EOMCC with all 1p and
2p1h, and active-space 3p2h excitations
(i.e., EA-EOMCCSDt).
Labels "p" and "h" in the description of IP and EA EOMCC
methods refer to particles (unoccupied correlated orbitals)
and holes (occupied correlated orbitals). EA and IP runs
produce both ground and excited states of systems obtained
by attaching an electron to or removing an electron from
the underlying CCSD reference ground state, using the EOMCC
formalism. Thus, EA and IP runs read $CCINP as well as
$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'.
Input Description
$CONTRL
2-11
CIMTYP
= FMOHESS
Input Description
= GAMMA
$CONTRL
2-12
= OPTIMIZE
= TRUDGE
= SADPOINT
= MEX
= CONICAL
= IRC
= VSCF
= DRC
= MD
= GLOBOP
= OPTFMO
= GRADEXTR
= SURFACE
= COMP
= G3MP2
= PROP
Input Description
= RAMAN
= NACME
= NMR
= EDA
= QMEFPEA
= TRANSITN
= FFIELD
= TDHF
= TDHFX
= MAKEFP
= FMO0
$CONTRL
2-13
* * * * * * * * * * * * * * * * * * * * * * * * * * * * *
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
Input Description
$CONTRL
2-14
= DEBUG
= routine
ICHARG =
Molecular charge.
MULT
=
= 1
= 2,3,...
(default=0, neutral)
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
Input Description
$CONTRL
= ZMT
= ZMTMPC
= FRAGONLY
2-15
NZVAR
= 0
= M
= M
Input Description
$CONTRL
2-16
=
= NONE
= READ
= SBKJC
= HW
= MCP
pseudopotential selection.
all electron calculation (default).
read ECP potentials in the $ECP input.
use Stevens, Basch, Krauss, Jasien,
Cundari ECP potentials for all heavy
atoms (Li-Rn are available).
use Hay, Wadt ECP potentials for heavy
atoms (Na-Xe are available).
use Huzinaga's Model Core Potentials.
The correct MCP potential will be chosen
to match the requested MCP valence basis
set (see $BASIS).
* * * * * * *
LOCAL
=
=
=
=
=
=
Input Description
= 0
= +1
$CONTRL
2-17
Input Description
$CONTRL
2-18
Input Description
$CONTRL
2-19
NOSYM
=
=
=
=
=
=
=
=
=
=
=
-7
-6
-5
-4
-3
-2
1
2
3
4
5
=
=
7
8
= 0
= 1
ETOLLZ
Input Description
$CONTRL
2-20
NORMF
= 0
= 1
NORMP
= 0
= 1
ITOL
=
= n
ICUT
= n
ISKPRP = 0
proceed as usual
Input Description
1
$CONTRL
2-21
IREST
=
= -1
=
=
=
=
=
GEOM
0
1
2
3
4
=
select where to obtain molecular geometry
= INPUT from $DATA input (default for IREST=0)
= DAF
read from DICTNRY file (default otherwise)
Input Description
$SYSTEM
2-22
==========================================================
$SYSTEM
group
(optional)
Input Description
$SYSTEM
2-23
PARALL =
KDIAG
=
= 0
= 1
= 2
= 3
COREFL =
Input Description
$SYSTEM
2-24
MODIO
Input Description
$SYSTEM
2-25
==========================================================
Input Description
$BASIS
2-26
==========================================================
$BASIS
group
(optional)
Input Description
$BASIS
2-27
Input Description
$BASIS
2-28
= SPK-AnZP = SPKrnZP -
= SPKrAnZP -
Input Description
$BASIS
2-29
Input Description
$BASIS
2-30
Input Description
$BASIS
2-31
http://setani.sci.hokudai.ac.jp/sapporo/Welcome.do
The mapping between the data base names and the keywords
used in GAMESS is (for n=D,T,Q):
data base name
keyword
Sapporo-nZP
SPK-nZP
Sapporo-nZP+diffuse
SPK-AnZP
Sapporo-DK-nZP
SPKrnZP
Sapporo-DK-nZP+diffuse
SPKrAnZP
Sapporo-nZP-2012
SPK-nZC
Sapporo-nZP-2012+diffuse
SPK-nZCD
Sapporo-DK-nZP-2012
SPKrnZC
Sapporo-DK-nZP-2012+diffuse
SPKrnZCD
Input Description
$BASIS
2-32
MCP-DZP
MCP-TZP
MCP-QZP
MCP/NOSeC-V-TZP+diffuse
MCP/NOSeC-V-QZP+diffuse
MCP-ATZP
MCP-AQZP
MCP/NOSeC-CV-DZP
MCP/NOSeC-CV-TZP
MCP/NOSeC-CV-QZP
MCPCDZP
MCPCTZP
MCPCQZP
MCP/NOSeC-CV-DZP+diffuse
MCP/NOSeC-CV-TZP+diffuse
MCP/NOSeC-CV-QZP+diffuse
MCPACDZP
MCPACTZP
MCPACQZP
GBASIS = IMCP-SR1 and IMCP-SR2 valence basis sets to be used with the improved
MCPs with scalar relativistic effects.
Input Description
$BASIS
2-33
MNDO
AM1
PM3
RM1
DFTB
selects
selects
selects
selects
selects
Input Description
$BASIS
2-34
=
=
=
=
Input Description
= DUNNING
= HUZINAGA
= HONDO7
$BASIS
2-35
Input Description
$BASIS
2-36
$end
$ecp
Cr-ecp SBKJC
O-ecp none
...snipped...
O-ecp none
H-ecp none
...snipped...
H-ecp none
$end
* * *
EXTFIL = a flag to read basis sets from an external file,
Input Description
$BASIS
2-37
Input Description
$DATA
2-38
==========================================================
$DATA group
$DATAS group
$DATAL group
(required)
(if NESC chosen, for small component basis)
(if NESC chosen, for large component basis)
TITLE
----------------------------------------------------------2-
GROUP, NAXIS
Input Description
$DATA
2-39
----------------------------------------------------------4-
Input Description
$DATA
2-40
Input Description
if
if
if
if
if
COORD=UNIQUE,
COORD=HINT,
COORD=CART,
COORD=ZMT,
COORD=ZMTMPC,
$DATA
follow
follow
follow
follow
follow
card
card
card
card
card
2-41
sequence
sequence
sequence
sequence
sequence
U
U
C
G
M
NAME, ZNUC, X, Y, Z
NAME
Input Description
$DATA
2-42
Input Description
$DATA
2-43
Input Description
$DATA
2-44
----------------------------------------------------------5G-
ATOM
ATOM
i1 BLENGTH
ATOM
i1 BLENGTH
i2 ALPHA
ATOM
i1 BLENGTH
i2 ALPHA
i3 BETA i4
It can be
Input Description
$DATA
2-45
----------------------------------------------------------9G-
----------------------------------------------------------10G-
STRING VALUE
Input Description
$DATA
2-46
----------------------------------------------------------5M-
ATOM
ATOM BLENGTH
Input Description
$DATA
2-47
Input Description
$ZMAT
2-48
==========================================================
$ZMAT
group
Input Description
$ZMAT
2-49
AUTO
Input Description
$ZMAT
2-50
Input Description
$ZMAT
1.0
0.0
1.0
1.0
0.0
-1.0
2-51
0.0
1.0
0.0
Input Description
$LIBE
2-52
==========================================================
$LIBE
group
Input Description
$SCF
2-53
==========================================================
$SCF
group
Input Description
DAMP
SHIFT
RSTRCT
DEM
=
=
=
=
$SCF
2-54
defaults for
ab initio:
semiempirical:
EXTRAP
T
T
DAMP
F
F
SHIFT RSTRCT
F
F
F
F
DIIS
F/T
F
SOSCF
T/F
F
Input Description
$SCF
2-55
* * * * * * * * * * * * * * * * * * *
more info on the convergence methods,
the 'Further Information' section.
* * * * * * * * * * * * * * * * * * *
Input Description
$SCF
2-56
---- orbital modification options ---The four options UHFNOS, VVOS, MVOQ, and ACAVO are
mutually exclusive. The latter 3 require RUNTYP=ENERGY,
and should not be used with any correlation treatment.
UHFNOS = flag controlling generation of the natural
orbitals of a UHF function. (default=.FALSE.)
VVOS
= 0
= n
ACAVO
PACAVO =
Input Description
$SCF
2-57
NHOMO
NLUMO
MOM
Input Description
$SCF
2-58
----- GVB wavefunction input ----The next parameters define the GVB wavefunction. See
also MULT in the $CONTRL input. The GVB wavefunction
assumes orbitals are in the order core, open, pairs.
NCO
NSETO
NO
NPAIR
The number of geminal pairs in the -GVBfunction. Maximum of 12. The default
corresponds to open shell SCF (default=0).
CICOEF =
ALPHA
BETA
Input Description
$SCF
2-59
0
1
2
VTSCAL =
Input Description
$SCF
2-60
SCALF
MAXVT
VTCONV =
* * * * * * * * * * * * * * * * *
more discussion of GVB/ROHF input
the 'further information' section
* * * * * * * * * * * * * * * * *
==========================================================
Input Description
$SCFMI
2-61
==========================================================
$SCFMI
group
NF
occupied
MF
ITER
DTOL
Input Description
$SCFMI
2-62
ALPHA
DIISON
MXDIIS
and
CNVLOC
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)
Input Description
$DFT
2-63
==========================================================
$DFT
group
Input Description
$DFT
2-64
Many choices are given below, perhaps the most sensible are
local DFT: SVWN
pure DFT GGA: BLYP, PW91, B97-D, PBE/PBEsol
hybrid DFT GGA: B3LYP, X3LYP, PBE0
pure DFT meta-GGA: revTPSS
hybrid DFT meta-GGA: TPSSh, M06
but of course, everyone has their own favorite!
pure exchange functionals:
= SLATER
Slater exchange
= BECKE
Becke 1988 exchange
= GILL
Gill 1996 exchange
= OPTX
Handy-Cohen exchange
= PW91X
Perdew-Wang 1991 exchange
= PBEX
Perdew-Burke-Ernzerhof exchange
These will be used with no correlation functional at all.
pure correlation functionals:
= VWN
Vosko-Wilk-Nusair correlation, using
their electron gas formula 5 (aka VWN5)
= VWN3
Vosko-Wilk-Nusair correlation, using
their electron gas formula 3
= VWN1RPA Vosko-Wilke-Nusair correlation, using
their e- gas formula 1, with RPA params.
= PZ81
Perdew-Zener 1981 correlation
= P86
Perdew 1986 correlation
= LYP
Lee-Yang-Parr correlation
= PW91C
Perdew-Wang 1991 correlation
= PBEC
Perdew-Burke-Ernzerhof correlation
= OP
One-parameter Progressive correlation
These will be used with 100% HF exchange, if chosen.
= SVWN
Input Description
$DFT
2-65
Input Description
$DFT
2-66
Input Description
$DFT
2-67
Input Description
$DFT
2-68
Input Description
$DFT
2-69
= a flag to turn on the computation of the E(3) nonadditive energy term. Pertinent only for DFT-D3,
it forces DC=.TRUE. (default=.FALSE.)
DCSR
DCS6
Input Description
DCS8
$DFT
2-70
The old keywords DCPAR and DCEXP were replaced by DCS6 and
DCSR in 2010. Similarly, DCOLD has morphed into IDCVER.
- - The Local Response Dispersion (LRD) correction includes
atomic pair-wise -C6/R**6, -C8/R**8, and -C10/R**10 terms,
whose coefficients are computed from the molecular system's
electron density and its nuclear gradient. The nuclear
gradient assumes the dispersion coefficients do not vary
with geometry, which causes only a very small error in the
gradient. Optionally, 3 and 4 center terms may be added,
at the 1/R**6 level; in this case, nuclear gradients may
not be computed at all.
Since the three numerical parameters are presently known
only for the long-range exchange corrected BOP functional,
calculations may specify simply DFTTYP=LCBOPLRD. The
"LCBOPLRD" functional will automatically select the
following:
DFTTYP=BOP
LC=.TRUE. MU=0.47
LRDFLG=.TRUE. LAMBDA=0.232 KAPPA=0.600 RZERO=3.22
leaving only the choice for MLTINT up to you.
References for LRD are
T.Sato, H.Nakai J.Chem.Phys. 131, 224104/1-12(2009)
T.Sato, H.Nakai J.Chem.Phys. 133, 194101/1-9(2010)
LRDFLG = flag choosing the Local Response Dispersion (LRD)
C6, C8, and C10 corrections. Default=.FALSE.
MLTINT = flag to add the 3 and 4 center 6th order terms,
the default=.FALSE. Note that nuclear gradients
are not available if these multi-center terms
are requested.
Three numerical parameters may be input. The defaults
shown are optimized for the BOP functional with the LC
correction for long-range exchange.
Input Description
$DFT
2-71
MU
Input Description
$DFT
2-72
CMP2
NLEB
Input Description
$DFT
2-73
JANS
NTHE
NPHI
Input Description
$DFT
2-74
NLEB0
NTHE0
NPHI0
Input Description
$DFT
2-75
technical parameters:
THRESH = threshold for ignoring small contributions to the
Fock matrix. The default is designed to produce
no significant energy loss, even when the grid is
as good as "army grade". If for some reason you
want to turn all threshhold tests off, of course
requiring more CPU, enter 1.0e-15.
default: 1.0e-4/Natoms/NRAD/NTHE/NPHI
GTHRE
Input Description
$DFT
2-76
AUXFUN = AUX0
= AUX3
THREE
Input Description
$TDDFT
2-77
==========================================================
$TDDFT
group
(relevant if TDDFT chosen in $CONTRL)
Input Description
$TDDFT
2-78
Input Description
$TDDFT
2-79
Input Description
$TDDFT
2-80
MULT
TDPRP
TPA
= a flag requesting two-photon absorption crosssections. These are computed for each of the
NSTATE excited states, after first evaluating
Input Description
$TDDFT
2-81
NONEQ
BETA
PFREQ
Input Description
$TDDFT
2-82
NLEB
NTHE
NPHI
SG1
Input Description
$TDDFT
2-83
Input Description
$DFTB
2-84
==========================================================
$DFTB
group
Input Description
NDFTB
$DFTB
2-85
DAMPXH =
DAMPEX =
SRSCC
ITYPMX
Input Description
= 0
= 1
= 2
$DFTB
2-86
Broyden's method.
Interpolation is applied for atomic
(or shell-resolved when SRSCC=.TRUE.)
charges, but not Hamiltonian matrix.
(reserved)
DIIS for charges.
Default: 0.
ETEMP
DISP
DISPPR
HUBDER
Input Description
$DFTB
2-87
SPNCST
Input Description
$DFTB
2-88
Input Description
$DFTSBK
2-89
==========================================================
$DFTBSK
group
(required if GBASIS=DFTB)
Input Description
$CIS
2-90
==========================================================
$CIS
group
Input Description
$CIS
2-91
IROOT
CISPRP =
HAMTYP =
=
=
MULT
Input Description
$CIS
2-92
- - - - - - - - - - - DIAGZN =
=
=
DGAPRX =
NGSVEC =
MXVEC
NDAVIT =
DAVCVG =
CHFSLV =
=
=
RDCISV =
MNMEDG =
Input Description
$CIS
2-93
==========================================================
$CISVEC
group
Input Description
$MP2
2-94
==========================================================
$MP2
group
MP2PRP=
Input Description
$MP2
2-95
OSPT=
CODE
RHF
RHF
energy gradient
SERIAL
DDI
IMS
RIMP2
SERIAL
DDI
IMS
-
UHF
UHF
ROHF
ROHF
ROHF
energy gradient
energy gradient energy
OSPT=ZAPT
ZAPT
RMP
SERIAL SERIAL
SERIAL
SERIAL
DDI
DDI
DDI
DDI
RIMP2
-
The default for serial runs (p=1) is CODE=IMS for RHF, and
CODE=SERIAL for UHF or ROHF (provided PARALL is .FALSE. in
$SYSTEM). When p>1 (or PARALL=.TRUE.), the default becomes
CODE=DDI. However, if FMO is in use, the default for
closed shell parallel runs is CODE=IMS. The "SERIAL" code
for OSPT=RMP will run with modest scalability when p>1.
The many different MP2 programs are written for different
hardware situations. Here N is the number of atomic basis
functions, and O is the number of correlated orbitals in
the run:
The original SERIAL programs use N**3 memory, and have
larger disk files and generally takes longer than CODE=IMS.
The IMS program uses N*O**2 memory, and places most of its
data on local disks (so you must have good disk access),
and will run in parallel...ideal for small clusters. Using
this program on a node where the disks are of poor quality
(SATA-type) and with many cores accessing that single disk
may be very I/O bound. Adding more memory can make this
program run more efficiently. Network traffic is modest
when running in parallel.
Input Description
$MP2
2-96
CUTOFF =
Input Description
$MP2
2-97
for some reason you wish to use the other two closed shell
codes for SCS-MP2 gradients.
The following keywords apply to any CODE=SERIAL MP2 run, or
to parallel ROHF+MP2 runs using OSPT=RMP:
LMOMP2=
CPHFBS =
NWORD =
METHOD= n
AOINTS=
==========================================================
Input Description
$RIMP2
2-98
===========================================================
$RIMP2
group
USEDM
Input Description
$RIMP2
2-99
IVMTD
Input Description
$RIMP2
2-100
Input Description
$AUXBAS
2-101
===========================================================
$AUXBAS
group
Input Description
$CCINP
2-102
==========================================================
$CCINP
group
Input Description
$CCINP
2-103
NFZV
MAXCC
ICONV
Input Description
$CCINP
2-104
IREST
Input Description
$CCINP
2-105
Input Description
$EOMINP
2-106
==========================================================
$EOMINP
group
(optional, for CCTYP=EOM-CCSD, CR-EOM, or CR-EOML)
(optional, for CCTYP=EA-EOMx or IP-EOMx)
(optional for CCSD properties, or CCTYP=CR-CCL)
Input Description
$EOMINP
2-107
Input Description
$EOMINP
2-108
GROUP
NSTATE
Input Description
IROOT
$EOMINP
2-109
Input Description
$EOMINP
2-110
Input Description
$EOMINP
2-111
Input Description
$EOMINP
2-112
MACT
= fine tuning
For MINIT=1
For MINIT=1
For MINIT=2
For MINIT=2
The default
Input Description
$EOMINP
2-113
MICEOM
MICCI
Input Description
$EOMINP
2-114
==========================================================
Input Description
$MOPAC
2-115
==========================================================
$MOPAC
group
Input Description
$GUESS
2-116
==========================================================
$GUESS
group
Input Description
= FMO
= HUCSUB
= DMREAD
$GUESS
2-117
NORB
Input Description
$GUESS
2-118
TOLE
Input Description
$GUESS
2-119
Input Description
2-120
==========================================================
$VEC
group
$DM
group
$MOFRZ
group
IFRZ
==========================================================
Input Description
$STATPT
2-121
==========================================================
$STATPT
group
Pick from
NR
RFO
QA
Input Description
$STATPT
2-122
DXMAX
TRMAX
TRMIN
--- the next three control mode following --IFOLOW = Mode selection switch, for RUNTYP=SADPOINT.
For METHOD=RFO or QA, the mode along which the
energy is maximized, other modes are minimized.
Usually referred to as "eigenvector following".
For METHOD=SCHLEGEL, the mode whose eigenvalue
is (or will be made) negative. All other
curvatures will be made positive.
Input Description
$STATPT
2-123
Input Description
$STATPT
2-124
IHREP
Input Description
$STATPT
2-125
1
0
-1
-2
NPUN
3
2
1
0
-1
-2
Print
Print
Print
Print
INTR
INTR
INTR
INTR
at
at
at
at
all,
all,
all,
0+L,
orbitals
orbitals
orbitals
orbitals
at all
at 0+L (default)
never
never
---- the next parameters control harmonic constraints --Harmonic constraints can be added to the current geometry
by setting ALL the keywords below. For instance, to
harmonically constrain the distance between atom 3 and 12
to a distance of 2.0 Angstrom and a force constant of 500
kcal/mol, the following example can be used:
IHMCON(1)=1,3,12 SHMCON(1)=2.0 FHMCON(1)=500.0
The default is all zeros which means do not do this.
IHMCON = array of coordinates to constrain. The input
is similar to IZMAT in $ZMAT, a code integer,
and the atoms involved in the coordinate.
The code integer may only be 1, for stretches.
SHMCON = equilibrium constraint values for the distances
Input Description
$STATPT
2-126
RMIN
Input Description
RMAX
$STATPT
2-127
RLIM
Input Description
$TRUDGE
2-128
==========================================================
$TRUDGE
group
IEX
Input Description
$TRUDGE
2-129
Input Description
$TRURST
2-130
==========================================================
$TRURST
group
Input Description
$FORCE
2-131
==========================================================
$FORCE
group
Input Description
$FORCE
2-132
(default is .FALSE.)
= 1
= 2
VIBSIZ =
Input Description
$FORCE
2-133
= 1
= 0
NPUN
= 2
= 1
= 0
----- the rest control normal coordinate analysis ---VIBANL = flag to activate vibrational analysis.
(the default is .TRUE. for RUNTYP=HESSIAN, and
otherwise is .FALSE.)
SCLFAC = scale factor for vibrational frequencies, used
in calculating the zero point vibrational energy.
Some workers correct for the usual overestimate
in SCF frequencies by a factor 0.89. ZPE or other
methods might employ other factors, see
J.P.Merrick, D.Moran, L.Radom
J.Phys.Chem.A 111, 11683-11700 (2007).
The output always prints unscaled frequencies, so
this value is used only during the thermochemical
analysis. (Default is 1.0)
TEMP
FREQ
Input Description
$FORCE
2-134
Input Description
$CPHF
2-135
==========================================================
$CPHF
group
Input Description
$CPMCHF
2-136
==========================================================
$CPMCHF
group
(relevant for analytic RUNTYP=HESSIAN,NACME,CONICAL)
= maximum iterations.
(default=300)
CPTOL
PRCCHG
Input Description
$CPMCHF
2-137
NACST
--- the next three choose the other CPMCHF convergers --If none is selected, conjugate gradient (CG) is used.
GCRODR
GCRO
Input Description
$CPMCHF
2-138
(default=.FALSE.)
MICIT
KICIT
Input Description
$CPMCHF
2-139
IPDIR
ITERB
Input Description
$MASS
2-140
==========================================================
$MASS
Input Description
$HESS $GRAD
2-141
==========================================================
$HESS
group
(relevant for RUNTYP=HESSIAN if RDHESS=.TRUE.)
(relevant for RUNTYP=IRC if FREQ,CMODE not given)
(relevant for RUNTYP=OPTIMIZE,SADPOINT if HESS=READ)
$GRAD
group
(relevant for RUNTYP=OPTIMIZE or SADPOINT)
(relevant for RUNTYP=HESSIAN when RDHESS=.TRUE.)
This
Input Description
$DIPDR $ALPDR
2-142
==========================================================
$DIPDR
group
$ALPDR
group
Input Description
$VIB $VIB2
2-143
==========================================================
$VIB
group
$VIB2
group
Input Description
$VSCF
2-144
==========================================================
$VSCF
group
Input Description
$VSCF
2-145
Input Description
$VSCF
2-146
AMP
STPSZ
Input Description
$VSCF
2-147
VCI
IEXC2
IEXC
0
1
2
3
1
1
2
IEXC2
0
0
0
0
1
2
1
Input Description
$VSCF
2-148
VCFCT
MPDIP
Input Description
CMODE
$VSCF
2-149
$VIBSCF
group
Input Description
2-150
==========================================================
$GAMMA
group
required if RUNTYP=GAMMA
$EQGEOM
group
$HLOWT
$GLOWT
group
group
These are the lower triangular parts of the hessian and 3rd
derivative matrices, read in the same format as printed by
an earlier RUNTYP=GAMMA.
==========================================================
Input Description
$IRC
2-151
==========================================================
$IRC
group
GCUT
RCUT
ACUT
MXOPT
IHUPD
Input Description
$IRC
2-152
GA
PACE = AMPC4
Input Description
$IRC
2-153
----- The next two are used by all PACE choices ----STRIDE = Determines how far apart points on the reaction
path will be. STRIDE is used to calculate the
step taken, according to the PACE you choose.
The default is good for the GS2 method, which is
very robust. Other methods should request much
smaller step sizes, such as 0.10 or even 0.05.
(default = 0.30 sqrt(amu)-Bohr)
NPOINT = The number of IRC points to be located in this
run. The default is to find only the next point.
(default = 1)
----- constraint ----Of course, applying a constraint to the saddle point search
and the reaction path means that you are not locating the
true saddle, nor following the true reaction path.
IFREEZ = array of Cartesian coordinates to freeze. The
IRC stepper works in mass-weighted Cartesian
space, making it impossible to freeze internal
coordinates. An input of IFREEZ(1)=4,8 means to
freeze the x coordinate of the 2nd atom and the
y coordinate of the 3rd atom, that is, we count
coordinates x1,y1,z1,x2,y2,z2,x3,y3,z3,...
----- The next two let you choose your output volume ----Let F mean the first IRC point found in this run,
and L mean the final IRC point of this run.
Let INTR mean the internuclear distance matrix.
NPRT
1
0
-1
-2
NPUN
1
0
-1
-2
Print
Print
Print
Print
INTR
INTR
INTR
INTR
at
at
at
at
all,
all,
all,
F+L,
orbitals
orbitals
orbitals
orbitals
Input Description
$IRC
2-154
----- The next two tally the reaction path results. The
defaults are appropriate for starting from a saddle
point, restart values are automatically punched out.
NEXTPT = The number of the next point to be computed.
STOTAL = Total distance along the reaction path to next
IRC point, in mass weighted Cartesian space.
Input Description
$IRC
2-155
* * * * * * * * * * * * * * * *
hints about IRC tracking, see
'further information' section.
* * * * * * * * * * * * * * * *
==========================================================
Input Description
$DRC
2-156
==========================================================
$DRC
group
Input Description
initial
$DRC
2-157
NVEL
Input Description
$DRC
2-158
ENM
Input Description
$DRC
2-159
HESS
HESS2
Input Description
$DRC
2-160
1
0
-1
NPUN
2
1
0
-1
==========================================================
Input Description
$MEX
2-161
===========================================================
$MEX
group
(relevant if RUNTYP=MEX)
STPSZ
NRDMOS =
=
=
=
=
(default = 0.1D+00)
Input Description
$MEX
2-162
NMOS1
NMOS2
NPRT
= Printing orbitals
= 0 No orbital printed out except at the first
geometry (default)
= 1 Orbitals are printed each geometry. If MCSCF
is used, CI expansions are also printed.
TGRMS
===========================================================
Usage notes:
1. Normally $CONTRL will not give SCFTYP or MULT keywords.
SCF1 and SCF2 can be given in any order. The combinations
permitted ensure roughly equal sophistication in the
treatment of electron correlation.
2. After reading $MEX, SCFTYP and MULT will be set to the
more complex of the two choices, which is considered to be
RHF < ROHF < UHF < GVB < MCSCF. This permits the $SCF
input defining a GVB wavefunction to be read and tested for
correctness, in a GVB+ROHF run. Since only one SCFTYP is
stored while reading the input, you might need to provide
some keywords that are normally set by default for the
other (such as ensuring DIIS is selected in $SCF if either
of the states is UHF).
Input Description
$MEX
2-163
Input Description
$CONICL
2-164
===========================================================
$CONICL
group
(relevant if RUNTYP=CONICAL)
Input Description
$CONICL
2-165
SIGMA
Input Description
$CONICL
2-166
Input Description
$MD
2-167
===========================================================
$MD
group
(relevant if RUNTYP=MD)
= MD integrator selection.
= FROG (leapfrog). This is less accurate, and lacks
the special ensemble stepper option NVTNH.
= VVERLET (velocity Verlet) - default.
DT
NVTNH
Input Description
$MD
2-168
velocity Verlet.
= 0 means use NVE Verlet stepping
= 1 means use NVT Verlet stepping
= 2 means use Nose/Hoover chain NVT Verlet stepping
The default is 2 if either NVT option RSTEMP or
RSRAND is chosen, but is 0 otherwise.
NSTEPS = number of MD time steps to be found in this run,
default=10000.
TTOTAL = total time elapsed in the previous part of a MD
trajectory which is being restarted (READ=.TRUE.).
The default means this trajectory is a new one, or
perhaps the start of a production phase of the MD.
(default=0.0 seconds)
* * *
BATHT
default=.FALSE.
Input Description
$MD
2-169
DELR
NPROP
Input Description
CCMS
$MD
2-170
USAMP
IUSTYP
Input Description
$MD
2-171
MREMD
LEVERY
Input Description
$MD
2-172
MBT
MBR
QRAND
READ
Input Description
$RDF
2-173
==========================================================
$RDF
group
water
1
1
water
water
1
2
1
3
2
1
3
1
water
Input Description
STOP
HH
dum
dum
dum
dum
STOP
$end
water
2
2
2
3
3
2
3
3
$RDF
water
2-174
==========================================================
Input Description
$GLOBOP
2-175
==========================================================
$GLOBOP
group
(relevant to RUNTYP=GLOBOP)
JSTRND =
RIORD
=
=
=
Input Description
$GLOBOP
2-176
TEMPF
NTEMPS =
MCTYP
NFRMOV =
NFMORV =
MCMIN
NGEOPT =
NTRAN
NROT
NBLOCK =
Input Description
$GLOBOP
2-177
NBLOCK=NGEOPT/(NTRAN+NROT)
If NBLOCK is input, the number of geometries at
each temperature will be taken as
NGEOPT=NBLOCK*(NTRAN+NROT)
Each block has NTRAN translational steps followed
by NROT rotational steps.
NAIFG
AIMOVE =
SCALE
ALPHA
BOLTWT =
=
=
NSTMIN =
Input Description
$GLOBOP
2-178
DACRAT =
UPDFAC =
SEPTOL =
Input Description
$GLOBOP
2-179
=
=
IFXFRG =
NPRBND =
NSMTP
SMTEMP =
==========================================================
Input Description
$GLBFRG group
$GLOBOP
2-180
(relevant to RUNTYP=GLOBOP)
==========================================================
This group defines the sets of ab initio atoms in $DATA
that are treated as groups in Monte Carlo runs. It also
defines any groups that are to be frozen, not moved during
the Monte Carlo search.
AIFRG
Input Description
$GLOBOP
IBNDS
PRSEP
2-181
==========================================================
Input Description
$GRADEX
2-182
==========================================================
$GRADEX
group
DPRED
STPT
Input Description
$GRADEX
2-183
---- the next parameters apply only to METHOD=SR ---DELCOR = the corrector step should be smaller than this
value before the next predictor step is taken.
(default = 0.001)
MYSTEP = maximum number of micro iteration allowed to
bring the corrector step length below DELCOR.
(default=20)
SNUMH
Input Description
$GRADEX
2-184
* * * * * * * * * * * * * * * * *
the 'further information' section
some help with GRADEXTR runs.
* * * * * * * * * * * * * * * * *
Input Description
$SURF
2-185
==========================================================
$SURF
group
Input Description
$SURF
2-186
IGRP1
ORIG1
DISP1
Input Description
$LOCAL
2-187
==========================================================
$LOCAL
group
Input Description
$LOCAL
2-188
Input Description
$LOCAL
2-189
NOUTA
Input Description
$LOCAL
2-190
* * *
NINA
MOINA
NINB
MOINB
Input Description
$LOCAL
2-191
Input Description
$LOCAL
2-192
IJMO
MOIJ
ZIJ
IPROT
Input Description
$LOCAL
2-193
Input Description
$LOCAL
2-194
atomic units.
For more information, see
I.Adamovic, M.S.Gordon
==========================================================
* *
For
the
the
* *
* * * * * * * * * * * * * * * *
hints about localizations, and
LCD energy decomposition, see
'further information' section.
* * * * * * * * * * * * * * * *
==========================================================
Input Description
$TRUNCN
2-195
==========================================================
$TRUNCN
group
IMOPR
Input Description
$TRUNCN
2-196
completely unaltered.
- - There are now two ways to specify what orbitals are to
be truncated. The most common usage is for preparation of
a buffer zone for QM/MM computations, with an Effective
Fragment Potential representing the non-quantum part of
the system. This input is NATAB, NATBF, ICAPFR, ICAPBF,
in which case the $DATA input must be sorted into three
zones. The first group of atoms are meant to be treated
in later runs by full quantum mechanics, the second
group by frozen localized orbitals as a 'buffer', and the
third group is to be substituted later by an effective
fragment potential (multipoles, polarizabilities, ...).
Note that in the DOPROJ=.TRUE. run, all atoms are still
quantum atoms.
NATAB
NATBF
Input Description
$TRUNCN
2-197
shorter:
IZAT
ILAT
==========================================================
Input Description
$ELMOM
2-198
==========================================================
$ELMOM
group
(not required)
= 0
1
2
3
WHERE
= COMASS
NUCLEI
POINTS
CUM
Input Description
$ELMOM
2-199
==========================================================
Input Description
$ELPOT
2-200
==========================================================
$ELPOT
group
(not required)
= COMASS
NUCLEI
POINTS
GRID
PDC
center of mass
at each nucleus (default)
at points given in $POINTS
at grid given in $GRID
at points controlled by $PDC.
Input Description
$ELDENS
2-201
==========================================================
$ELDENS
group
(not required)
MORB
WHERE
= COMASS
NUCLEI
POINTS
GRID
center of mass
at each nucleus (default)
at points given in $POINTS
at grid given in $GRID
Input Description
$ELFLDG
2-202
==========================================================
$ELFLDG
group
(not required)
WHERE
= COMASS
NUCLEI
POINTS
- center of mass
- at each nucleus (default)
- at points given in $POINTS
* *
For
see
* *
* * * * * * * * * * * * * * * * * *
conversion factors, and references
the 'further information' section.
* * * * * * * * * * * * * * * * * *
Input Description
$POINTS $GRID
2-203
==========================================================
$POINTS
group
(not required)
which properties
group must
by an integer
The next NPOINT
the X, Y, and Z
==========================================================
$GRID
group
(not required)
Input Description
$PDC
2-204
==========================================================
$PDC
group
=
determines the points to be used, choose
GEODESIC to use a set of points on several fused
sphere van der Waals surfaces, with points
selected using an algorithm due to Mark
Spackman. The results are similar to those
from the Kollman/Singh method, but are
less rotation dependent. (default)
CONNOLLY to use a set of points on several fused
sphere van der Waals surfaces, with points
selected using an algorithm due to Michael
Connolly. This is identical to the method
used by Kollman & Singh (see below)
CHELPG to use a modified version of the CHELPG
algorithm, which produces a symmetric
grid of points for a symmetric molecule.
CONSTR = NONE
Input Description
$PDC
2-205
DELR
Input Description
$PDC
2-206
QPOLE
Input Description
$MGC
2-207
==========================================================
$MGC
group
NPOINT
CP
Input Description
$RADIAL
2-208
==========================================================
$RADIAL
group
RMAX
Input Description
$RADIAL
2-209
Input Description
$MOLGRF
2-210
==========================================================
$MOLGRF
group
MESH
BOUND
CHECK
Input Description
$STONE
2-211
==========================================================
$STONE
group
(optional)
NRAD <nrad>
NANG <nang>
SMOOTH <nbecke>
SMRAD <nbckrd>
Input Description
i
name
$STONE
2-212
---------------------------------------------------------ATOMS
---------------------------------------------------------BONDS
---------------------------------------------------------CMASS
Input Description
$STONE
2-213
nmo, where
QMMMBUF
NMO
==========================================================
The second and third moments on the printout can be
converted to Buckingham's tensors by formula 9 of
A.D.Buckingham, Quart.Rev. 13, 183-214 (1959)
These can in turn be converted to spherical tensors
by the formulae in the appendix of
S.L.Price, et al. Mol.Phys. 52, 987-1001 (1984)
Input Description
$RAMAN
2-214
==========================================================
$RAMAN
group
Input Description
$COMP
2-215
==========================================================
$COMP
group
= NUMBER OF OPTIMIZATIONS
= INDEX OF OPTIMIZATION USED FOR HESSIAN CALCULATION
= NUMBER OF SINGLE POINT ENERGY CALCULATIONS
TOTAL # ENERGY CALCULATIONS = IOP + IENE
=
=
=
=
=
=
=
=
Input Description
$COMP
2-216
=
=
=
=
=
=
=
=
ISPHER
GBASIS
NGAUSS
NDFUNC
NFFUNC
NPFUNC
DIFFSP
DIFFS
SPHERICAL HARMONICS
INTERNAL BASIS SETS
VALUES
(D POLARIZATION FUNCTION) VALUES
(F POLARIZATION FUNCTION) VALUES
(P POLARIZATION FUNCTION) VALUES
(SP DIFFUSE FUNCTION) OPTION
(S DIFFUSE FUNCTION) OPTION
Input Description
$COMP
2-217
a) select RUNTYP=G3MP2
b) select RUNTYP=COMP and $COMP's RMETHD=G32CCSD
c) select RUNTYP=COMP and explicitly define it:
$COMP
iop=2, ihs=1, iene=3, icalc=2, ienref=5,
cenhlc=-0.009170, czpscl=0.8929, cheats=g3mp2,
icmplv(2)=2,
icmplv(4)=2,
icmplv(5)=2,
ccctyp(3)=ccsd(t),
icisph(4)=1,
cgbass(1)=n31,n31,n31,g3l,n31,
icngss(1)=6,6,6,0,6,
icndfc(1)=1,1,1,0,1,
icnacr(2)=0,
icnbcr(2)=0,
icnder(2)=1,
icenca(1)=3,4,
icencb(1)=5,5
$END
In case of failures when the program shifts basis sets, try
running without point group symmetry: NOSYM=1 in $CONTRL.
===========================================================
Input Description
$NMR
2-218
==========================================================
$NMR
group
PPARA
PEVEC
PITER
PRMAT
Input Description
$NMR
2-219
PTEINT
TEDBG
==========================================================
Input Description
$MOROKM
2-220
==========================================================
$MOROKM group
(relevant if RUNTYP=EDA)
Input Description
$MOROKM
2-221
EQUM
RDENG
==========================================================
The present implementation has some quirks:
1. The initial guess of the monomer orbitals is not
controlled by $GUESS. The program first looks for a $VEC1,
$VEC2, ... group for each monomer. The orbitals must be
Input Description
$MOROKM
2-222
Input Description
$MOROKM
2-223
Input Description
$LMOEDA
2-224
==========================================================
$LMOEDA
group
(relevant if RUNTYP=EDA)
Input Description
$LMOEDA
2-225
Input Description
$QMEFP
2-226
===========================================================
$QMEFP
group
STEP2
EMULT
Input Description
$QMEFP
2-227
Those QM methods which are not based on fully selfconsistent solutions of the QM/EFP interaction Hamiltonian
(namely TDDFT, CIS, MP2, CCSD, EOM-CCSD) provide results
which include the EFP's perturbation by the correlated
density, and/or a particular excited state's density. This
approach is termed "Method 2" in the following references:
1. P.Arora, L.V.Slipchenko, S.P.Webb, A.DeFusco, M.S.Gordon
J.Phys.Chem.A 114, 6742-6750(2010)
2. A.DeFusco, J.Ivanic, M.W.Schmidt, M.S.Gordon
J.Phys.Chem.A 115, 4574-4582(2011)
===========================================================
Input Description
$FFCALC
2-228
==========================================================
$FFCALC
group
Input Description
$FFCALC
2-229
Input Description
$FFCALC
2-230
Input Description
$TDHF
2-231
==========================================================
$TDHF
group
FREQ
Input Description
$TDHF
2-232
ATOL
BTOL
* * *
BSHG
BEOPE
BOR
INIG
GTHG
GOKE
Input Description
$TDHF
2-233
Input Description
$TDHF
2-234
Input Description
$TDHFX
2-235
==========================================================
$TDHFX
group
USE_Q
Input Description
$TDHFX
2-236
Input Description
$TDHFX
2-237
ALPHA w =
compute the dynamic polarizability:
alpha(-w;w).
BETA w1 w2 / BETA_NI w1 w2 =
compute the dynamic first hyperpolarizability:
beta(-w1-w2;w1,w2).
GAMMA w1 w2 w3 / GAMMA_NI w1 w2 w3 =
compute the dynamic second hyperpolarizability:
gamma(-w1-w2-w3;w1,w2,w3).
POCKELS w / POCKELS_NI w =
compute electro-optic Pockels effect: beta(-w;w,0).
OR w / OR_NI w =
optical rectification: beta(0;w,-w).
SHG w / SHG_NI w =
second harmonic generation: beta(-2w;w,w).
KERR w / KERR_NI w =
DC Kerr effect: gamma(-w;w,0,0).
ESHG w / ESHG_NI w =
electric field induced 2nd harm gen: gamma(-2w;w,w,0).
THG w / THG_NI w =
third harmonic generation: gamma(-3w;w,w,w).
DFWM w / DFWM_NI w =
degenerate four wave mixing gamma(-w;w,-w,w).
See the review
D.P.Shelton, J.E.Rice Chem.Rev. 94, 3-29(1994)
for more information on the quantities just above. The
next options are nuclear derivatives of some of the above.
DMDX_NI =
compute the dipole derivative matrix,
the geometrical first derivative of MU.
DADX w / DADX_NI w =
compute the polarizability derivative matrix, the
Input Description
$TDHFX
2-238
Input Description
$TDHFX
2-239
HRAMAN 0.02
FREE 0.02
FREE 0.02 0.02
HRAMAN 0.03
$END
References:
"Time Dependent Hartree-Fock schemes for analytic
evaluation of the Raman intensities"
O.Quinet, B.Champagne J.Chem.Phys. 115, 6293-6299(2001).
"Analytical TDHF second derivatives of dynamic electronic
polarizability with respect to nuclear coordinates.
Application to the dynamic ZPVA correction."
O.Quinet, B.Champagne, B.Kirtman
J.Comput.Chem. 22, 1920-1932(2001).
"Analytical time-dependent Hartree-Fock schemes for the
evaluation of the hyper-Raman intensities"
O.Quinet, B.Champagne J.Chem.Phys. 117, 2481-2488(2002).
errata: JCP 118, 5692(2003)
"Analytical time-dependent Hartree-Fock evaluation of the
dynamically zero-point averaged (ZPVA) first
hyperpolarizability"
O.Quinet, B.Kirtman, B.Champagne
J.Chem.Phys. 118, 505-513(2003).
Computer quirks:
1. This package uses file numbers 201, 202, ... but some
compilers (chiefly g77) may not support unit numbers above
99. The remedy is to use a different computer or compiler.
2. If you experience trouble running this package under
AIX, degrade the optimization of subroutine JDDFCK in
hss2b.src, by placing this line
@PROCESS OPT(2)
immediately before JDDFCK, recompile hss2b, and relink.
==========================================================
Input Description
$EFRAG
2-240
==========================================================
$EFRAG
group
(optional)
Input Description
$EFRAG
2-241
=CART
=INT
POLMETHD=SCF
=FRGSCF
Input Description
$EFRAG
2-242
dipolar interaction.
POSITION=OPTIMIZE Allows full optimization within the
ab initio part, and optimization of
the rotational and translational
motions of each fragment. (default)
=FIXED
Allows full optimization of the
ab initio system, but freezes the
position of the fragments. This
makes sense only with two or more
fragments, as what is frozen is the
fragments' relative orientation.
FIXED may be used with RUNTYP being
OPTIMIZE, SADPOINT, HESSIAN and IRC.
=EFOPT
the same as OPTIMIZE, but if the
fragment gradient is large, up to
5 geometry steps in which only the
fragments move may occur, before
the geometry of the ab initio piece
is relaxed. This may save time by
reusing the two electron integrals
for the ab initio system.
NBUFFMO = n
Input Description
$EFRAG
SWR1, SWR2
2-243
RCUT
=
=
=
=
=
switch
switch
switch
switch
switch
off
off
off
off
off
Input Description
$EFRAG
2-244
ISCRPOL
=
fragment-fragment polarization screening.
= 0 damping is controlled by POLSCR sections in
the $FRAGNAME inputs. If not found, there
will be no screening. If POLSCR is found,
you must also use ISCRELEC=0 and SCREEN3.
= 1 damping will use a Tang-Toennis style
Gaussian formula,
(1-exp(aR**2)(1+aR**2)
where the default value of a=0.6. In order
to change the 'a' parameter, give
POLAB <a's value>
STOP
in the $FRAGNAME input. A smaller value
may be useful for ionic EFPs. (default)
ISCRDISP =
fragment-fragment dispersion screening
= 0 Use Tang-Toennies damping, with a fixed
parameter a=1.5.
= 1 use an overlap based damping factor,
1-S**2(1-2ln|S|+2ln**2|S|)
instead. There is no parameterization, so
there's no other input. (default)
Input Description
$EFRAG
2-245
The following keywords are for use with the EFP2-AI (a.k.a.
EFP2-QM) dispersion calculation, that is, the calculation
of the dispersion energy in a mixed system containing one
or more EFP2 fragment(s) and a molecule modeled with a
fully ab initio method (e.g. Hartree-Fock).
QMDISP
= 0
= 1
ISCRQMDS
=-1
= 0
= 1
NODSGRD
Input Description
$EFRAG
2-246
FRAGNAME=XXX
NAME, X, Y, Z
NAME, I, DISTANCE, J, BEND, K, TORSION
(COORD=CART)
(COORD=INT)
Input Description
$EFRAG
2-247
Input Description
$FRAGNAME
2-248
==========================================================
$FRAGNAME
group
(required for each FRAGNAME given in $EFRAG)
input
input
input
input
input
input
input
input
input
input
input
monopoles,
follow input sequence
dipoles,
follow input sequence
quadrupoles,
follow input sequence
octopoles,
follow input sequence
electrostatic screening,
follow input seq.
polarizable points,
follow input sequence
polarizability screening, follow input seq.
fitted "repulsion",
follow input sequence
Pauli exchange,
follow input sequence
dispersion,
follow input sequence
charge transfer,
follow input sequence
-EM-ED-EQ-EO-ES-P-PS-R-PE-D-CT-
Input Description
$FRAGNAME
2-249
COORDINATES
NAME
is a unique string identifying the point.
X, Y, Z are the Cartesian coordinates of the point, and
must be in Bohr units.
WEIGHT, ZNUC are the atomic mass and nuclear charge, and
should be given as zero only for points which are
not nuclei.
In EFP1 potentials, the true nuclei will appear twice, once
for defining the positive nuclear charge and its screening,
and a second time for defining the electronic distributed
multipoles.
Repeat line -3- for each expansion point, and terminate
the list with a "STOP".
---------------------------------------------------------Note: the multipole expansion produced by RUNTYP=MAKEFP
comes from Stone's distributed multipole analysis (DMA).
An alternative expansion, from a density based multipole
expansion (DBME) performed on an adaptive grid is placed in
the job's PUNCH file. This alternative multipole expansion
may be preferable if large basis sets are in use (the DMA
expansion is basis set sensitive). The DBME values can be
inserted in place of the DMA values, for -EM-, -ED, -EQ-,
and -EO- sections, if you wish. Experience suggests that
DBME multipoles are about as accurate as those obtained
using DMA.
-EM1-
MONOPOLES
Input Description
$FRAGNAME
2-250
Input Description
-EO2-
$FRAGNAME
2-251
SCREEN
NAME, A, B
Input Description
$FRAGNAME
2-252
NAME, A, B
POLARIZABLE POINTS
NAME, X, Y, Z
Input Description
$FRAGNAME
2-253
POLSCR
NAME, B
Input Description
$FRAGNAME
2-254
REPULSIVE POTENTIAL
See also the $FRGRPL input group, which defines the fit for
the EFP1-EFP1 repulsion term.
REPULSIVE POTENTIAL signals the start of the subgroup
containing the fitted exchange repulsion potential, for the
interaction between the fragment and the ab initio part of
the system. This term also accounts, in part, for other
effects, since it is a fit to a remainder. The fitted
potential has the form
N
sum
i
-R2-
C * exp[-D * r**2]
i
i
NAME, X, Y, Z, N
Input Description
-R3-
$FRAGNAME
2-255
C, D
Input Description
-D2-
$FRAGNAME
2-256
NAME, X, Y, Z
Input Description
$FRAGNAME
2-257
Input Description
$FRGRPL
2-258
==========================================================
$FRGRPL
group
PAIR=FRAG1 FRAG2
NAME1 NAME2 A B
*or*
NAME1 NAME2 'EQ' NAME3 NAME4
Input Description
$FRGRPL
2-259
Input Description
$EWALD
2-260
==========================================================
$EWALD
group
KMAX
==========================================================
Input Description
$MAKEFP
2-261
==========================================================
$MAKEFP
group
(relevant if RUNTYP=MAKEFP)
FRAG
CHTR
CTVVO
DISP
Input Description
$MAKEFP
2-262
Input Description
$PRTEFP
2-263
==========================================================
$PRTEFP
group
(optional)
Input Description
$PRTEFP
2-264
Input Description
$DAMP
2-265
==========================================================
$DAMP
group
Input Description
$DAMP
2-266
Input Description
$DAMP
2-267
MAXIT
THRSH
==========================================================
Input Description
$DAMPGS
2-268
==========================================================
$DAMPGS group
Input Description
$PCM
2-269
==========================================================
$PCM
group
(optional)
Input Description
$PCM
2-270
Input Description
$PCM
2-271
THF
DMSO (or DMETSOX)
SMD has many additional solvents, see below.
The default solvent name is "INPUT" which means you must
give the numerical values defining some other solvent, as
described below.
* * * non-SMD calculations * * *
The next set of parameters controls the computation:
parameterization of the solvents, ICOMP which has an impact
on the PCM electrostatics, and other keywords related to
cavitation, dispersion, and repulsion corrections: ICAV,
IDISP, IREP/IDP.
------ICOMP
Input Description
$PCM
2-272
IDP
WB
ETA2
Input Description
$PCM
of the solvent.
2-273
(default=1.75)
Input Description
$PCM
2-274
Input Description
$PCM
2-275
=
option for spheres forming the cavity:
= 0 centers spheres on each nucleus in the quantum
solute, and every atom in EFP. (default)
= N use N initial sphere, whose centers XE, YE, ZE
and radii RIN must be specified in $PCMCAV.
MXTS
Input Description
$PCM
2-276
Input Description
$PCMCAV
2-277
==========================================================
$PCMCAV
group
(optional)
Input Description
$PCMCAV
2-278
The source for the van der Waals radii is "The Elements",
2nd Ed., John Emsley, Clarendon Press, Oxford, 1991, except
for C,N,O where the Pisa group's experience with the best
radii for PCM treatment of singly bonded C,N,O atoms is
taken. The radii for a few transition metals are given by
A.Bondi, J.Phys.Chem. 68, 441-451(1964).
ALPHA
Input Description
$TESCAV
2-279
==========================================================
$TESCAV
group
(optional)
METHOD =
Input Description
$TESCAV
2-280
MTHALL =
1, 2, 3, or 4 (default = 4)
All values in the array METHOD are set to MTHALL
MTHAUT =
0 or 1 (default = 0)
If RUNTYP=OPTIMIZE and frozen atoms are defined
by IFCART, MTHAUT=1 will select METHOD=1 for
frozen atoms. See also AUTFRE and NTSFRZ.
BONDRY =
NTSFRZ =
==========================================================
Input Description
$NEWCAV
2-281
==========================================================
$NEWCAV
group
(optional)
PM
NEVAL
Input Description
$PCMGRD
2-282
==========================================================
$PCMGRD
group
(optional)
It
==========================================================
Input Description
$IECPCM
2-283
==========================================================
$IEFPCM
group
(optional)
Input Description
$PCMITR
2-284
==========================================================
$PCMITR
group
MXITR1 =
MXITR2 =
THRSLS =
DENSLS =
IDIRCT =
Input Description
$PCMITR
2-285
=
=
=
=
=
RCUT1
mid-range
Default=15.0 A
molecule, the
off.
RCUT2
RASC
RABI
Input Description
$PCMITR
2-286
==========================================================
Input Description
$DISBS
2-287
==========================================================
$DISBS
group
(optional)
XYZE
NADD=7 NKTYP(1)= 0 0 1 1
XYZE(1)=2.9281086
0.0
2.9281086
0.0
2.9281086
0.0
2.9281086
0.0
2.9281086
0.0
2.9281086
0.0
2.9281086
0.0
2 2 3
.0001726
.0001726
.0001726
.0001726
.0001726
.0001726
.0001726
0.2
0.05
0.2
0.05
0.75
0.2
0.2 $END
==========================================================
Input Description
$DISREP
2-288
==========================================================
$DISREP
group
(optional)
Input Description
$DISREP
2-289
DKA
RWA
Other elements have DKA and RWA values of 0.0 and so must
be given in the input deck, or the dispersion/repulsion
energy will be 0. For EFP/PCM calculations, only QM atoms
need DKA and RWA values to calculate the DIS-REP energy.
---- the following data must given for ILJ=1:
RHO
EPSI
Input Description
$SVP
2-290
==========================================================
$SVP
group
(optional)
Input Description
$SVP
2-291
Input Description
$SVP
2-292
Input Description
$SVP
2-293
Input Description
$SVP
2-294
RHOSOLV
DISDMP
DISLIN
EXCLIN
NGSLGR
Input Description
$SVP
2-295
= flag to turn on an alternative (to DEFESR) semiempirical correction for local electrostatic
effects based on the electric field's normals to
the surface cavity. This also adds cavitation/
dispersion/solvent structure (CDS) effects drawn
from the SMD model, see SMD in $PCM.
(Default=.FALSE.)
* * *
Input Description
$SVP
2-296
IPNRF
==========================================================
Input Description
$SVPIRF
2-297
==========================================================
$SVPIRF
group
Input Description
$COSGMS
2-298
==========================================================
$COSGMS
group
(optional)
DISEX
OUTCHG
Input Description
$COSGMS
2-299
Input Description
$SCRF
2-300
==========================================================
$SCRF
group
(optional)
==========================================================
Additional information on the SCRF model can be
found in the Further Information chapter.
Input Description
$ECP
2-301
==========================================================
$ECP
group
Input Description
$ECP
2-302
*** Card sets -2- and -3- are repeated LMAX+1 times
***
NGPOT
CLP,NLP,ZLP
SBKJC
NONE
SBKJC
NONE
Input Description
$ECP
2-303
Input Description
$MCP
2-304
==========================================================
$MCP
ANAT
is a 8 character name for the MCP atom.
It must match the name given for that atom
in the $DATA input group.
ZEFF, MCPFMT
FORMAT(F10.2, A8)
ZEFF
(ACOEF(L), L=1,NOAN)
FORMAT(MCPFMT)
(AEXPN(L), L=1,NOAN)
FORMAT(MCPFMT)
Input Description
$MCP
2-305
(NINT(L), L=1,NOAN)
NINT(L)
FORMAT(10I3)
(C(I), I=1,NG(IS))
FORMAT(MCPFMT)
expansion coefficients of the core orbital
COORD=HINT
PP=MCP $END
0.0
0.0 - O K
<<<< this is an MCP atom
<<<< (311/311/1) basis
-0.16370140
0.22673090E-01
-0.26304451
0.19109693
0.58040872
0.50918856
1.0000000
1.0000000
1.0000000
1.0000000
Input Description
$MCP
O 8.0
LC
1.2031
MCP READ
L
3
1 44.242510
2 6.2272700
3 1.4361751
L
1
1 0.40211473
L
1
1 0.12688798
D 1; 1 1.154 1.0
H 1.0
TZV
P 1; 1
PCC
1.1012
1.100
2-306
0.0
0.0 - O K
<<<< this is an MCP atom
<<<< (311/311/1) basis
-0.13535836
0.17372951E-01
-0.30476423
0.16466813
0.43955753
0.46721611
1.0000000
1.0000000
1.0000000
1.0000000
121.875 0.0 + O K I
<<<< not an MCP atom, TZV+pol basis
1.0
$END
$MCP
NR (2S/2P)
1 16
6.00(4D15.8)
.31002267
.27178756E-01
S(2)P(4)
<<<< MCP for the atom O
Input Description
25.973731
0 0
41.361784
57480.749
829.06758
31.542773
1.6077662
.85822477D-04
.76877460D-02
.22782741D+00
.55937988D-02
$MCP
2-307
13.843290
17270.167
346.04791
14.815300
.77153240
.18173691D-03
.20823429D-01
.33492260D+00
.16121923D-02
5766.9282
151.12147
7.0298236
.37052330
.84803428D-03
.52424753D-01
.28833079D+00
.10915544D-04
2107.0076
68.233250
3.3561489
.17799002
.25439914D-02
.11864010D+00
.93046197D-01
.21431633D-03
$END
==========================================================
Input Description
$RELWFN
2-308
==========================================================
$RELWFN
group
(optional)
Input Description
$RELWFN
2-309
There are two possible basis set choices for IOTC, DK,
or RESC calculations. The Sapporo relativistic segmented
contractions are available for elements down to Xenon (see
GBASIS=SPKrnZP or SPKrAnZP). DK3 basis sets for H-Lr
obtained at U. of Tokyo exist in the form of general
contractions,
http://www.riken.jp/qcl/
publications/dk3bs/periodic_table.html
which gives the EPAPS data published by
T.Tsuchiya, M.Abe, T.Nakajima, K.Hirao
J.Chem.Phys. 115,4463-4472(2001)
A program to extract this web page into GAMESS's format
is provided with GAMESS, see file ~/gamess/tools/dk3.f.
Light to medium atom main group (H-Kr) DK2 bases exist,
look for the names cc-pVnZ_DK on
http://www.emsl.pnl.gov:2080/forms/basisform.html
For NESC, you must provide three basis sets, for the
large and small components and an averaged one, which are
given in $DATAL, $DATAS, $DATA, respectively. The only
possible choice for these basis sets is due to Dyall, and
these are available from
http://www.emsl.pnl.gov:2080/forms/basisform.html Their
names are similar to cc-pVnZ(pt/sf/lc), pt=point or
fi=finite nucleus, sf for spin-free and the final field is
lc=large component ($DATAL), sc=small component ($DATAS),
and wf is a typo for Foldy-Wouthuysen 2e- basis ($DATA).
In GAMESS you can only use point nucleus approximation. The
need to input three basis sets means that you cannot use
$BASIS input, and you must use COORD=UNIQUE style input in
the various $DATA's. The three $DATA input groups must
contain identical information except for the primitive
expansion coefficients, as the three basis sets must have
the same exponents. In case the option to treat only some
atoms relativistically is chosen, all non-relativistic
atoms must have identical basis input in all three groups.
The finite size of nuclei is not taken into account, so
do not use any basis set obtained including this effect.
For NESC, the one electron part of the spin-orbit
operator can be corrected, while for RESC, one can compute
spin-orbit coupling with relativistic corrections to both
one and two electron SOC integrals, unless internal
uncontraction is requested (in this case only 1 electron
Input Description
$RELWFN
2-310
Input Description
$RELWFN
2-311
Input Description
$RELWFN
2-312
Input Description
$RELWFN
2-313
==========================================================
Input Description
$EFIELD
2-314
==========================================================
$EFIELD
This
external
general,
energies
group
(not required)
EVEC
SYM
==========================================================
Restrictions: analytic hessians are not available, but
numerical hessians are. Because an external field causes a
molecule with a dipole to experience a torque, geometry
optimizations must be done in Cartesian coordinates only.
Internal coordinates eliminate the rotational degrees of
freedom, which are no longer free.
A nuclear hessian calculation will have two rotational
modes with non-zero "frequency", caused by the torque. A
gas phase molecule will rotate so that the dipole moment is
anti-parallel to the applied field. To carry out this
rotation during geometry optimization will take many steps,
and you can help save much time by inputting a field
opposite the molecular dipole. There is also a stationary
point at higher energy with the dipole parallel to the
field, which will have two imaginary frequencies in the
hessian. These will appear as the first two modes in a
hessian run, but will not have the i for imaginary included
on the printout since they are rotational modes.
sign conventions:
Dipole vectors are considered to point from the negative
end of the molecule to the positive end. Thus HCl at the
MP2/aug-cc-pVDZ level's geometry of R=1.2831714 has a
positive dipole, if we place Cl at the origin and H along
Input Description
$EFIELD
2-315
Input Description
$INTGRL
2-316
==========================================================
$INTGRL
group
(optional)
Input Description
$INTGRL
2-317
NORDER = 0 (default)
= 1 Sort integrals into canonical order. There
is little point in selecting this option, as
no part of GAMESS requires ordered integrals.
See also NSQUAR through NOMEM.
NSQUAR = 0 Sorted integrals will be in triangular
canonical order (default)
= 1 instead sort to square canonical order.
NDAR
= Number of direct access logical records to be
used for the integral sort (default=2000)
LDAR
= Length of direct access records (site dependent)
NBOXMX = 200
Maximum number of bins.
NWORD =
0
Memory to be used (default=all of it).
NOMEM =
0
If non-zero, force external sort.
The following parameters control integral restarts.
IST=JST=KST=LST=1
NREC=1
INTLOC=1
Values shown are defaults, and mean not restarting.
==========================================================
Input Description
$FMM
2-318
==========================================================
$FMM
group
NS
IWS
IDPGD
IEPS
==========================================================
Input Description
$TRANS
2-319
==========================================================
$TRANS
group
Input Description
$TRANS
2-320
Zero
Input Description
$FMO
2-321
==========================================================
$FMO
group
IEFMO
MODEFM
Input Description
$FMO
2-322
Input Description
$FMO
2-323
TDTYP
CITYP
NFRAG
LAYER
Input Description
$FMO
2-324
Input Description
$FMO
2-325
Input Description
NACUT
$FMO
2-326
Input Description
$FMO
2-327
Input Description
$FMO
2-328
Input Description
$FMO
2-329
Input Description
$FMO
2-330
Input Description
$FMOPRP
2-331
==========================================================
$FMOPRP
group
CONV
Input Description
$FMOPRP
2-332
unconverged SCF.
IJVEC
Input Description
$FMOPRP
2-333
Input Description
$FMOPRP
2-334
Input Description
$FMOPRP
2-335
FMO documentation.
Default: 9999,9999 (which means do not use).
NAODIR = a parameter to decide whether to enforce DIRSCF.
Useful for incore integral runs in parallel.
NAODIR is the number of AO orbitals that is
expected to produce 100,000,000 non-zero
integrals. Using this and assuming NAO**3.5
dependence, the program will then guess how many
integrals will each n-mer have and whether they
will fit into the available memory. If they are
determined not to fit, DIRSCF will be set true.
This option overwrites MCONV but not MCONFG.
If set to 0, then the default in-core integral
strategy is used. (default=0)
VDWRAD = array of van der Waals radii in Angstrom, one for
each atom in the periodic table. Reasonable values
are set only for a few light atoms and otherwise a
value of 2.5 is used. VDWRAD values are used only
to compute distance between fragments and thus
somewhat affect all distance-based approximations.
II. Parameters defining parallel execution
MODPAR = parallel options (additive options)
(default: 13, which is 1+4+8)
1 turns on/off heavy job first strategy (reduces
waiting on remaining jobs at barrier points)
(see also 8)
4 broadcast all fragments done by a group at once
rather than fragment by fragment.
8 alters the behavior of fragment initialixation:
if set, fragments are always done in the reverse
order (nfg, nfg-1, ...1) because distance
calculation costs decrease in the same order and
they usually prevail over making Huckel orbitals
or running free monomer SCF. Note that during
SCC (monomer SCF) iterations the order in which
monomers are done is determined by MODPAR=1.
16 if set, hybrid orbital projectors will not be
parallelized (may be useful on slow networks)
32 reserved
64 Broadcast F40 for FMO restarts. F40 should only
Input Description
$FMOPRP
2-336
Input Description
$FMOPRP
2-337
Input Description
$FMOPRP
2-338
Input Description
$FMOPRP
2-339
Input Description
$FMOPRP
2-340
Input Description
$FMOPRP
2-341
R0BDA
E0BDA
EFMO0
Input Description
$FMOPRP
2-342
Input Description
$FMOXYZ
2-343
==========================================================
$FMOXYZ
group
Input Description
$FMOXYZ
2-344
H 1
O 8
$end
II. advanced.
This allows you to mix basis sets, have multiple layers or
a non-standard without involving EXTFIL.
1. Do not define $BASIS.
2. Put each atom type in $DATA, followed by basis set,
either explicit or built in.
The names of atoms in $DATA have the following format,
where brackets indicate optional parameters:
S[.N][-L]
N and L may be omitted (taking the default value of 1),
S is the atom name (discarded upon reading),
N is the basis set ordinal number,
L is the layer.
S[.N][-L] may not exceed 8 characters.
Example: 2-layer water dimer. In the first layer, you want
to use STO-3G for the first molecule and your own basis set
for the second. In the second layer, you want to use 6-31G
and 6-31G* for the first and second molecules,
respectively.
$DATA
water dimer (H2O)2
C1
H-1 1
! explanation: layer 1, basis 1 (STO-3G) for Hydr.
sto 3
O-1 8
sto 3
Input Description
$FMOXYZ
2-345
1 5.0 1 1
H-2 1
n31 6
O-2 8
n31 6
Input Description
$OPTFMO
2-346
==========================================================
$OPTFMO
group
(relevant if RUNTYP=OPTFMO)
Input Description
$OPTFMO
2-347
Input Description
$FMOLMO
2-348
==========================================================
$FMOHYB
group
(optional, for FMO runs)
(this group was previously known as $FMOLMO)
Input Description
$FMOBND
2-349
==========================================================
$FMOBND
group
Input Description
$FMOEFP
2-350
==========================================================
$FMOEFP group
Input Description
$FMOEFP
2-351
(default=2)
Input Description
2-352
==========================================================
$FMOENM
group
$FMOEND
group
$OPTRST
group
Input Description
$GDDI
2-353
==========================================================
$GDDI
group
Input Description
runs.
$GDDI
2-354
(default is DLB)
Input Description
$ELG
2-355
===========================================================
$ELG
group
NCT
IPRI
= PRINT LEVEL
0 minimum printing (default)
3 debugging printing
LDOS
Input Description
$ELG
2-356
ATOB
Input Description
$DANDC
2-357
==========================================================
$DANDC
group
Input Description
$DANDC
2-358
SUBTYP =
= ATOM
Input Description
$DANDC
2-359
Input Description
$DANDC
2-360
=
Fermi function cutoff factor (default=15.0).
= p The value of Fermi function less than 10**(-p)
is considered as 0. The value greater than
[1 - 10**(-p)] is considered as 1.
NDCPRT =
= +1
= +2
= +4
= +8
IORBD
Input Description
$DANDC
2-361
DIIQTR = .TRUE.
Input Description
$DANDC
2-362
==========================================================
Input Description
$DCCORR
2-363
==========================================================
$DCCORR
group
(optional)
relevant for MPLEVL=2
relevant for CCTYP=LCCD, CCD, CCSD, CCSD(T), R-CC
Input Description
$DCCORR
2-364
WOCC
ONLYOC =
ITPART =
ISTCOR =
restart option for DC-MP2/CC.
= 0 does DC-MP2/CC calculation from the beginning
(default).
= n reads subsystem correlation energies
corresponding to subsystem 1-(n-1) from input
Input Description
$DCCORR
2-365
==========================================================
Input Description
$SUBSCF $SUBCOR
2-366
==========================================================
$SUBSCF group
$SUBCOR group
1,3,4,5
2,6,7,8
2,6,7,8,9,11
1,3,4,10,12,13,14
==========================================================
Input Description
$MP2RES $CCRES
2-367
==========================================================
$MP2RES group
$CCRES group
0.000E+00
0.000E+00
0.000E+00
-0.133110332082E+00
-0.130740147906E+00
-0.130483660838E+00
-0.132440119981E+00
-0.131009546477E+00
-0.130832334600E+00
-0.131377855474E+00
Input Description
$CIMINP
2-368
==========================================================
$CIMINP
group
NONE)
Input Description
$CIMINP
2-369
ZETA1
ZETA2
= 1
Input Description
= 2
$CIMINP
2-370
ETA
FCORE
MOFIX
Input Description
$CIMINP
2-371
= HIGHER
Input Description
$CIMINP
= LMO
2-372
orbitals (default).
use the localized molecular orbitals
(SUBMTD=CCSD, NONE, or HF only). IF
SUBMTD=NONE or HF, the only method that
can be accepted in $CIMFRG or $CIMATM is
CCSD.
= 1
Input Description
= 2
$CIMINP
2-373
==========================================================
* * * * * * * * * * * * * * * * * * * * *
For information on Cluster-in-Molecules,
see the 'further information' section
* * * * * * * * * * * * * * * * * * * * *
Input Description
$CIMATM $CIMFRG
2-374
==========================================================
$CIMATM
group
(optional for CIMTYP=SECIM in $CONTRL;
relevant to multi-level SECIM calculations)
$CIMFRG
group
Input Description
$CIMATM $CIMFRG
2-375
Input Description
$FFDATA $FFDB
2-376
==========================================================
$FFDATA group
$FFPDB group
JOUT
KOUT
LOUT
Input Description
ITSTAT =
=
=
=
=
$FFDATA $FFDB
2-377
PRES0
=
=
=
=
=
maximum number
maximum number
maximum number
maximum number
maximum number
(i.e. improper
MXWAGG = maximum number
of MM atoms.
of bonds.
of bond angles.
of dihedral rotation angles.
of dihedral bending angles.
torsion in CHARMM).
of wagging angles.
Input Description
$FFDATA $FFDB
2-378
specifies how
0
No adding
1
Add water
2
Add water
ITYPWAT=
=
=
=
IADDNA1=
IADDK1 =
IADDCA2=
IADDMG2=
IADDCL1=
number
number
number
number
number
of
of
of
of
of
CENTER = X, Y, Z
Na+
K+
Ca2+
Mg2+
Cl-
ions
ions
ions
ions
ions
randomly
randomly
randomly
randomly
randomly
added.
added.
added.
added.
added.
Default=0.
Default=0.
Default=0.
Default=0.
Default=0.
Input Description
$FFDATA $FFDB
2-379
Input Description
$FFDATA $FFDB
2-380
Input Description
$FFDATA $FFDB
2-381
Input Description
$FFDATA $FFDB
2-382
NRMSD
NGYRA
Input Description
$FFDATA $FFDB
2-383
NFOLD
RFOLD
IDOCHG =
IDOPOL =
IDOLJ =
IDOCMAP=
Input Description
$FFDATA $FFDB
2-384
Input Description
$FFDATA $FFDB
2-385
COORDINATES
NAME, NUC, X, Y,
NAME
NUC
X,Y,Z
STOP
Z
= The name of the atom.
= nuclear charge of the atom
= Cartesian coordinates in angstrom
PARAMETERS
NAME, MASS, Q, POL, SIGMA, EPSILON, SIGMA2, EPSILON2
NAME
= The name of the atom.
Q
= Force field charge (e) on the atom.
POL
= Polarizability of the atom, in A**3
SIGMA
= Lennard-Jones parameter in angstrom.
CHARMM uses RMIN/2. QuanPol will
automatically check whether it is
SIGMA or RMIN/2.
EPSILON = Lennard-Jones parameter in kcal/mol
SIGMA2, EPSILON2
= the LJ parameters for select 1-4 cases
in CHARMM. Give zeros if not these.
Give the same number of lines as specified in
the COORDINATES section
STOP
QMMMREP
N, C1, Z1, C2, Z2, C3, Z3, C4, Z4
N
= number of Gaussian type potentials.
currently only 4 is allowed.
C1 - C4 = strength factor of the potential
Z1 - Z4 = radial factor of the potential
give the same number of lines as specified in
the COORDINATES section
STOP
BOND
Input Description
$FFDATA $FFDB
2-386
Input Description
$FFDATA $FFDB
ATOM2
ATOM3
ATOM4
DBFC
=
=
=
=
DIHB0
2-387
STOP
CMAP
SERIAL#, ATOM1, ATOM2, ATOM3, ATOM4, ATOM5, ITYPE
SERIAL# = serial number of the CHARMM correction
map phi,psi couples. this is only
for notation purpose.
ATOM1
= serial number in COORDINATES section
for the first atom in the phi angle
(the carbonyl carbon of an amino acid
residue) of the peptide backbone in a
phi,psi couple.
ATOM2
= serial number in COORDINATES section
for the second atom in the phi angle
(peptide N atom) of a phi,psi couple.
this is the first atom of the psi angle.
ATOM3
= serial number in COORDINATES section
for the third atom in the phi angle
(the alpha C atom) of a phi,psi couple.
this is the 2nd atom of the psi angle.
ATOM4
= serial number in COORDINATES section
for the fourth atom in the phi angle
(a carbonyl carbon) of a phi,psi couple.
this is the third atom of the psi angle.
ATOM5
= serial number in COORDINATES section
for the fourth atom in the psi angle
(the next peptide N atom).
ITYPE
= specifies the two amino acid
residues a phi,psi couple belongs to.
=1 alanine-alanine
(most cases)
=2 alanine-proline
=3 glycine-glycine/glycine-proline
=4 proline-alanine
=5 proline-proline
(here alanine stands for any amino acid
residue other than glycine and proline)
Input Description
$FFDATA $FFDB
2-388
STOP
When $FFPDB is given, CMAP is automatically generated for
restart jobs, via $FFDATA input.
WAGGING
SERIAL#, ATOM2, ATOM3, ATOM4, ATOM1, WFC
SERIAL# = serial number of the wagging angle.
this is only for notation purpose.
ATOM2
= serial number in COORDINATES section
for the second atom in the wagging
angle.
ATOM3
= same as ATOM2, but for the third atom.
ATOM4
= same as ATOM2, but for the fourth atom.
ATOM1
= same as ATOM2, but for the first atom.
WFC
= wagging force constant kcal/mol/rad**2
STOP
Use a $END line to end $FFDATA.
---------------------------------------------------------**** the following inputs pertain to $FFPDB ****
-2- Simply paste a PDB text file into $FFPDB will work.
(1). H atoms must be added beforehand and appear at the
correct places. Currently QuanPol cannot add H atoms
or any other missing atoms.
(2). PDB format is enforced. Sequential numbers are not
used by QuanPol. Chemical symbols are used.
(3). Rename the atom at the very end of each chain as 'OXT'
or 'HXT', and delete the 'TER' lines. Multiple chains
are allowed.
(4). SSBOND lines are required to define S-S bonds.
Use a $END line to end $FFPDB.
==========================================================
Input Description
$CIINP
2-389
* * * *
and CI
section
* * * *
==========================================================
$CIINP
group
Input Description
$CIINP
2-390
Input Description
2-391
==========================================================
Input Description
2-392
NACT
NELS
SZ
STSYM =
Input Description
2-393
Input Description
2-394
Input Description
2-395
IROOT
Input Description
2-396
Input Description
2-397
Input Description
2-398
Input Description
$ORMAS
2-399
==========================================================
$ORMAS
group
MSTART
MAXE
Input Description
$ORMAS
2-400
GVAL
Input Description
$ORMAS
2-401
QCORR
FDIRCT
Input Description
$CEEIS
2-402
==========================================================
$CEEIS
group
Input Description
$CEEIS
2-403
Input Description
$CEEIS
2-404
$CEDATA
group
Input Description
$GCILST
2-405
==========================================================
$GCILST
group
Input Description
$GCILST
2-406
Input Description
$GMCPT
2-407
==========================================================
$GMCPT
group
Input Description
$GMCPT
2-408
SZ:
STSYM:
Input Description
$GMCPT
2-409
NEXCIT
NSPACE
MSTART
MINE
MAXE
Input Description
$GMCPT
2-410
IROOT
ISPINA
KNOSYM
KNOSPN
MXITER
THRCON
THRENE
MAXBAS
Input Description
MDI
$GMCPT
2-411
KXGMC
IWGT
KFORB
KROT
Input Description
$GMCPT
2-412
KSZDOE
THRGEN
THRHDE
Input Description
$GMCPT
2-413
EDSHFT
5. miscellaneous data
CEXCEN = string defining the units for the excitation
energy. Choose from these 4 strings (any case):
eV (default), cm-1, Kcal/mol, KJ/mol
DDTFPT = a flag requesting the distributed data integral
transformation be used, if the run is parallel.
This option requires MEMDDI in $SYSTEM. If there
is not enough memory to allow this, turn this
option off to use an alternate parallel
transformation (DEFAULT=.TRUE.).
Note: There are additional technical parameters for $GMCPT,
documented only in the source code file gmcpt.src.
---In case it is desirable for the GMC-QDPT program to
reproduce results obtained by the DETMRPT/MCQDPT programs:
a) use a CAS-SCF reference in the MCSCF step
b) select REFTYP=ORMAS here, and enter NSPACE=1, giving
only one value for MSTART, MINE, MAXE
c) retain the entire CAS reference in the internal
determinant's perturbation space, THRWGT=-1.0
d) select the original external determinant space's
perturbation treatment, KROT=.FALSE.
e) use equal alpha/beta orbital energies, KSZDOE=.FALSE.
f) in multi-state mode, select KXGMC off, to reproduce
those program's 0-th order reference states
g) ensure ISA is turned off, THRDE= -1.0
Input Description
$GMCPT
2-414
Input Description
$PDET/$ADDDET/$REMDET
2-415
==========================================================
$PDET
group
$ADDDET
$REMDET
group
group
Input Description
$SODET
2-416
==========================================================
$SODET
group
(required if CITYP=FSOCI)
NSOST
==========================================================
Input Description
$DRT $CIDRT
2-417
==========================================================
$DRT group
$CIDRT group
FOCI
SOCI
Input Description
$DRT $CIDRT
2-418
Input Description
$DRT $CIDRT
2-419
(default=10000)
==========================================================
Input Description
$MCSCF
2-420
==========================================================
$MCSCF
group
(for SCFTYP=MCSCF)
SOSCF
Input Description
$MCSCF
2-421
MICIT
NWORD
FORS
Input Description
$MCSCF
2-422
FINCI
Input Description
$MCSCF
2-423
MCFMO
--- the next three refer to FOCAS optimizations --CASDII = threshold to start DIIS (default=0.05)
CASHFT = level shift value (default=1.0)
NRMCAS = renormalization flag, 1 means do Fock matrix
Input Description
$MCSCF
2-424
Input Description
MOFRZ
$MCSCF
2-425
--- the last few FULLNR options are seldom used --NORB
NOROT
DROPC
==========================================================
Input Description
$MRMP
2-426
==========================================================
$MRMP
group
Input Description
$MRMP
2-427
Input Description
$MRMP
2-428
Input Description
$DETPT
2-429
==========================================================
$DETPT group
NEXT
NOS
Default=.TRUE.
Input Description
$DETPT
2-430
IPTST
WPTST
==========================================================
Input Description
$MCQDPT
2-431
==========================================================
$MCQDPT
group
Controls 2nd order MCQDPT (multiconfiguration quasidegenerate perturbation theory) runs, if requested by
MPLEVL=2 in $CONTRL. MCQDPT2 is implemented only for FORS
(aka CASSCF) wavefunctions. The MCQDPT method is a
multistate, as well as multireference perturbation theory.
The implementation is a separate program, interfaced to
GAMESS, with its own procedures for determination of the
canonical MOs, CSF generation, integral transformation, CI
in the reference CAS, etc. Therefore some of the input in
this group repeats data given elsewhere, particularly for
$DET/$DRT.
Analytic gradients are not available. Spin-orbit
coupling may be treated as a perturbation, included at the
same time as the energy perturbation. If spin-orbit
calculations are performed, the input groups for each
multiplicity are named $MCQD1, $MCQD2, ... rather than
$MCQDPT. Parallel calculation is enabled.
When applied to only one state, the theory is known as
multi-reference Moller-Plesset (MRMP), but the term MCQDPT
is used when this theory is used in its multi-state form.
Please note that this perturbation theory is not the same
thing as the CASPT2 theory, and should -NEVER- be called
that. A more complete discussion may be found in the
'Further Information' chapter.
Most values will inherit sensible defaults for the state
symmetry and the orbital space counts from the $DET or $DRT
input defining the MCSCF: however for multi-state runs, the
user probably has to supply the desired state and weighting
information.
In case of diabatic state generation at the MCQDPT
level, the settings for state selection and weights will be
inherited from the $DIABAT input, to be the same as used
for the Diabatic MO generation. Thus diabatization runs
will probably not give any input here, although they might
override NMOFZC/NMODOC defaults.
*** MCSCF reference wavefunction ***
Input Description
$MCQDPT
NEL
MULT
NMOACT =
NMOFZC =
NMODOC =
NMOFZV =
2-432
If the input file does not provide a $DET or $DRT, the user
must give NMOFZC, NMODOC, and NMOACT correctly here.
STSYM = The symmetry of the target electronic state(s).
See $DET for possible values: use AP/APP in Cs, not
primes. This must be given, and need not match the
state symmetry used in optimizing the orbitals by
$DET or $DRT, although it often does.
Default is the totally symmetric representation.
NOSYM
Input Description
$MCQDPT
2-433
XZERO
EDSHFT =
(default=0.0,0.0)
Input Description
$MCQDPT
2-434
= 0 skip canonicalization
(default when DIABAT=.TRUE.).
= 1 determine the canonical Fock orbitals.
Input Description
$MCQDPT
2-435
Input Description
$MCQDPT
2-436
THRGEN =
THRENE =
THRCON =
MDI
MXBASE =
NSOLUT =
NSTOP
LPOUT
Input Description
$MCQDPT
timing information.
2-437
(default=0)
PARAIO =
DELSCR =
Input Description
$EXCORR
2-438
===========================================================
$EXCORR
group
NFRZC
NINACT
NCORR
NFRZV
PUNTOL
Input Description
$EXCORR
2-439
RUNR12
SINGLS
F12EXP
NTHRDS
===========================================================
Input Description
$CASCI
2-440
==========================================================
$CASCI
group
This group carries out the Improved Virtual Orbital Complete Active Space CI method of Freed, Chaudhuri, and
co-workers. IVO-CASCI starts with a RHF reference, and
then generates IVOs, which are used in a CI computation
within an active space chosen by the user. The input
consists of this group, a $MCQDPT, and perhaps a $IVOORB
input, along with SCFTYP=RHF and MPLEVL=2. MULT in $CONTRL
applies to the SCF reference, while MULT in $MCQDPT selects
the spin of the IVO-CASCI state(s). Doublets are treated
by using a cation RHF reference.
IVOCAS = a flag to turn on IVO-CASCI computation. This is
usually the only input required (default=.FALSE.)
MOLIST = a flag to request complete control over the
active space specification. The default uses the
parameters in $MCQDPT to select from the IVOs
with the lowest energy. (default=.FALSE.)
DEGENR = a flag to indicate the HOMO is degenerate.
The program should set this for you.
PRINT
Input Description
$CASCI
2-441
0.0
0.0
.0000000000
.9813851814
.0289123030
.4758735367
$IVOORB
group
Input Description
$CISORT $GUGEM
2-442
$CISORT
group
LDAR
NOMEM
==========================================================
$GUGEM
group
= not used.
==========================================================
Input Description
$GUGDIA
2-443
==========================================================
$GUGDIA
group
Input Description
$GUGDIA
2-444
EREF
==========================================================
Input Description
$GUGDM
2-445
==========================================================
$GUGDM
group
==========================================================
Input Description
$GUGDM2
2-446
==========================================================
$GUGDM2
group
NOMEM
NDAR
LDAR
Input Description
$LAGRAN
$LAGRAN
group
2-447
NWORD =
NDAR
= 4000
LDAR
==========================================================
Input Description
$TRFDM2
2-448
==========================================================
$TRFDM2
group
NWORD =
= 2000
LDAR
NBOXMX= 200
==========================================================
Usually neither $LAGRAN nor $TRFDM2 are given. Since these
groups are normally used only for CI gradient runs, we list
here the restrictions on GUGA CI gradients:
a) SCFTYP=RHF, only
b) no FZV orbitals in $CIDRT, all MOs must be used.
c) the derivative integrals are computed in the 2nd
derivative code, which is limited to spd basis sets.
d) the code does not run in parallel.
e) Use WSTATE in $GUGDM2 to specify the state whose
gradient is to be found. Use IROOT in $GUGDM to
specify the state whose other properties will be
found. These must be the same state!
f) excited states often have different symmetry than the
ground state, so think about GROUP in $CIDRT.
g) the gradient can probably be found for any CI for
which you have sufficient disk to do the CI itself.
Time is probably about 2/3 additional.
See also $CISGRD for CI singles gradient
Input Description
$DIABAT
2-449
==========================================================
$DIABAT group
Input Description
$DIABAT
2-450
Input Description
NGRST
$DIABAT
2-451
Input Description
$DIABAT
2-452
NMLAP
Input Description
$DIABAT
2-453
THVIR
MAXIT
Input Description
$DIABAT
2-454
* * *
Diabatization runs can read additional input groups which
are not well described here. The $REFCSF and $DFMVEC input
groups are often given. They are usually prepared by a
three-fold DIABAT=.TRUE. run at some reference geometry,
using both REFGRP=.FALSE. and REFMOS=.FALSE.
$REFCSF: list of dominant configurations at the reference
geometry, which is created by a REFGRP=.FALSE. run. Note
that the SLCTTH threshold assists in deciding how many CSFs
are placed in the reference groups. The $REFCSF data is
then read by all REFGRP=.TRUE. runs at the various other
geometries, in the same format as it is generated.
$DFMVEC: a set of temporary DMOs for ordering the DMOs at
the current geometry. Usually these are DMOs prepared by a
3-fold way calculation at one reference geometry, but could
be from a 4-fold way calculation at a geometry very close
to the current one. This contains only active MOS, namely
NDOC+NALP+NAOS+NBOS+NVAL from $DRT.
$DPSVEC contains NDLAP resolution DMOs.
$DIAVEC contains NMLAP reference DMOs.
Note: $DFMVEC, $DPSVEC, $DFMVEC are typical $VEC type
inputs. Each is read ignoring the MO index, so you might
prepare the order reference MOs from converged natural
orbitals by simply deleting all doubly occupied orbitals,
and keeping all active orbitals. Similarly, reference or
resolution orbitals may be plucked from any desired orbital
set: natural, canonical, localized...
At present, there is no way for GAMESS to generate the
orientation information, although this can be read in.
Therefore their contents are not well described:
$LCLDC and $LCLVL are orientation data for DPSVEC and
DIAVEC, respectively, required if ORIENT=.TRUE.
==========================================================
Input Description
$TRANST
2-455
===========================================================
$TRANST
group
Input Description
$TRANST
2-456
Input Description
$TRANST
2-457
Input Description
$TRANST
2-458
-1
0
1
-2
Input Description
$TRANST
states is recommended.
Note that $GUESS is not read by this RUNTYP!
be in $VEC1 and possibly $VEC2 input groups.
2-459
Orbitals must
NFZC
NOCC
Input Description
$TRANST
2-460
Input Description
$TRANST
2-461
Input Description
$TRANST
2-462
Input Description
$TRANST
2-463
Input Description
$TRANST
2-464
Input Description
$TRANST
2-465
Input Description
$TRANST
2-466
Input Description
$TRANST
2-467
TOLE
==========================================================
Input Description
$TRANST
* * * * * * * * * * * * * * * * * * *
For information on RUNTYP=TRANSITN,
see the 'further information' section
* * * * * * * * * * * * * * * * * * *
2-468
Input Description
2-469
DRC,
157
DRT,
418
ECP,
302
EFIELD,
315
EFRAG,
241
ELDENS,
202
ELFLDG,
203
ELG,
356
ELMOM,
199
ELPOT,
201
EOMINP,
106
EQGEOM,
151
EWALD,
261
EXCORR,
439
FFCALC,
229
FFDATA,
377
FFPDB,
377
FMM,
319
FMO,
322
FMOBND,
350
FMOEND,
353
FMOENM,
353
FMOHYB,
349
FMOPRP,
332
FMOXYZ,
344
FORCE,
131
FRAGNAME,
249
FRGRPL,
259
GAMMA,
151
GCILST,
406
GDDI,
354
GEN,
392
GLBFRG,
181
GLOBOP,
176
GLOWT,
151
GMCPT,
408
GRAD,
142
GRADEX,
183
GRID,
204
GUESS,
116
GUGDIA,
444
GUGDM,
446
GUGDM2,
447
GUGEM,
443
HESS,
142
HLOWT,
151
IEFPCM,
284
INTGRL,
317
IRC,
152
IVOORB,
442
LAGRAN,
448
LIBE,
52
LMOEDA,
225
LOCAL,
188
MAKEFP,
262
MASS,
141
MCP,
305
MCQDPT,
432
MCSCF,
421
MD,
168
MEX,
162
MGC,
208
MOFRZ,
120
MOLGRF,
211
MOPAC,
115
MOROKM,
221
MP2,
94
MP2RES,
368
MRMP,
427
NEWCAV,
282
NMR,
219
OPTFMO,
347
OPTRST,
353
ORMAS,
400
PCM,
270
PCMCAV,
278
PCMGRD,
283
PCMITR,
285
PDC,
205
PDET,
416
POINTS,
204
PRTEFP,
264
QMEFP,
227
RADIAL,
209
RAMAN,
215
RDF,
174
RELWFN,
309
REMDET,
416
RIMP2,
98
SCF,
53
SCFMI,
61
SCRF,
301
SODET,
417
STATPT,
121
STONE,
212
SUBCOR,
367
SUBSCF,
367
SURF,
186
SVP,
291
SVPIRF,
298
SYSTEM,
22
TDDFT,
77
TDHF,
232
TDHFX,
236
TESCAV,
280
TRANS,
320
TRANST,
456
TRFDM2,
449
TRUDGE,
128
TRUNCN,
196
TRURST,
130
VEC,
120
VEC1,
120
VEC2,
120
VIB,
144
VIB2,
144
VIBSCF,
150
VSCF,
145
ZMAT,
48