Keywords used in MOPAC2007

&

Turn next line into keywords

+

Add another line of keywords

0SCF

Read in data, then stop

1ELECTRON

Print final one-electron matrix

1SCF

Do one scf and then stop

AIDER

Read in ab-initio derivatives

AIGIN

Geometry must be in gaussian format

AIGOUT

In arc file, include ab-initio geometry

ALLBONDS Print final bond-order matrix, including bonds to hydrogen

ALLVEC

Print all vectors (keywords vectors also needed)

AUTOSYM

Symmetry to be imposed automatically

AUX Output auxiliary information for use by other programs

AM1

Use the AM1 hamiltonian

BAR=n.nn

reduce bar length by a maximum of n.nn%

BCC (Unique)

Solid is body-centered cubic (used by BZ)

BFGS

Optimize geometries using bfgs procedure

BIGCYCLES=n

Do a maximum of n big steps

BIRADICAL

System has two unpaired electrons

BONDS

Print final bond-order matrix

CAMP

Use Camp-King converger in SCF

CARTAB

Print point-group character table

C.I.=n
C.I.=(n,m)

A multi-electron configuration interaction specified

CHARGE=n

Charge on system = n (e.g. NH4 = +1)

CHARST

Print details of working in CHARST

CIS

C.I. uses 1 electron excitations only

CISD

C.I. uses 1 and electron excitations

CISDT

C.I. uses 1, 2 and 3 electron excitations

COMPFG

Print heat of formation calculated in COMPFG

COSWRT

Write details of the solvent accessible surface to a file

CUTOFP=n.nn

Madelung distance cutoff is n .nn Angstroms

CYCLES=n

Do a maximum of n steps

DCART

Print part of working in DCART

DDMAX=n.nn

See EF code

DDMIN=n.nn

Minimum trust radius in a EF/TS calculation

DEBUG

Debug option turned on

DEBUGPULAY

Print working in PULAY

DENOUT

Density matrix output (channel 10)

DENSITY

Print final density matrix

DERI1

Print part of working in DERI1

DERI2

Print part of working in DERI2

DERITR

Print part of working in DERIT

DERIV

Print part of working in DERIV

DERNVO

Print part of working in DERNVO

DFORCE

Force calculation specified, also print force matrix.

DFP

Use Davidson-Fletcher-Powell method to optimize geometries

DISEX=n.nn

Distance for interactions in fine grid in COSMO

DMAX=n.nn

Maximum stepsize in eigenvector following

DOUBLET

Doublet state required

DRC
DRC=n.nnn

Dynamic reaction coordinate calculation

DUMP=nn.nn

Write restart files every n seconds

ECHO

Data are echoed back before calculation starts

EF

Use ef routine for minimum search

EIGEN Print canonical eigenvectors instead of LMOs in MOZYME calculations

EIGS

Print all eigenvalues in ITER

ENPART

Partition energy into components

EPS=n.nn

Dielectric constant in COSMO calculation

ESP

Electrostatic potential calculation

ESPRST

Restart of electrostatic potential

ESR

Calculate RHF spin density

EXCITED

Optimize first excited singlet state

EXTERNAL=name

Read parameters off disk

FIELD=(n.nn,m.mm,l.ll)

An external electric field is to be used

FILL=n

In RHF open and closed shell, force M.O. n to be filled

FLEPO

Print details of geometry optimization

FMAT

Print details of working in FMAT

FOCK

Print last Fock matrix

FREQCY

Print symmetrized Hessian in a FORCE calculation

FORCE

Force calculation specified

GEO-CHECK

Apply extra safety checks

GEO-OK

Override some safety checks

GNORM=n.nn

Exit when gradient norm drops below n .n

GRADIENTS

Print all gradients

GRAPH

Generate unformatted file for graphics

GRAPHF

Generate formatted file for graphics suitable for  Jmol and MOPETE.

HESS=n

Options for calculating Hessian matrices in EF

H-PRIORITY
H-PRIORITY=n.nn

Heat of formation takes priority in DRC

HYPERFINE

Hyperfine coupling constants to be calculated

INT

Make all coordinates internal coordinates

IRC
IRC=n

Intrinsic reaction coordinate calculation

ISOTOPE

Force matrix written to disk (channel 9 )

ITER

Print details of working in ITER

ITRY=nn

Set limit of number of SCF iterations to n

IUPD=n

Mode of Hessian update in eigenvector following

KINETIC=n.nnn

Excess kinetic energy added to DRC calculation

KING

Use Camp-King converger for SCF

LARGE

Print expanded output

LBFGS Use the low-memory version of the BFGS optimizer

LET

Override certain safety checks

LEWIS Print the Lewis structure

LINMIN

Print details of line minimization

LOCALIZE

Print localized orbitals

LOG Generate a log file

MECI

Print details of MECI calculation

MEP=1

Print Orozco-Luque electrostatic potentials for cross-section

MEP=2

Print Orozco-Luque electrostatic potentials for Connolly surface

MERS=(n1,n2,n3)

Keyword for BZ

MICROS=n

Use specific microstates in the C.I.

MINMEP

Minimize MEP minima in the plane defined

MMCCR1

Multiplier for Molecular Mechanics correction for 5 membered rings

MMCCR2

Exponent for Molecular Mechanics correction for 5 membered rings

MMCCROK

Apply Molecular Mechanics corrections to 5 membered rings

MMOK

Use molecular mechanics correction to CONH bonds

MNDO

Use the MNDO hamiltonian

MNDOD

Use the MNDO-d hamiltonian

MODE=n

In EF, follow Hessian mode no. n

MOLDAT

Print details of working in MOLDAT

MOLSYM

Print details of working in MOLSYM

MOPAC

Use old MOPAC definition for 2nd and 3rd atoms

MS=n

In MECI, magnetic component of spin

MULLIK

Print the Mulliken population analysis

N**2

In excited state COSMO calculations, set the value of N**2

NLLSQ

Minimize gradients using NLLSQ

NOANCI

Do not use analytical C.I. derivatives

NOLOG

Suppress log file trail, where possible

NOMM

Do not use molecular mechanics correction to CONH bonds

NONET

NONET state required

NONR

Do not use Newton-Raphson method in EF

NOREOR

In symmetry work, use supplied orientation

NOSYM

Point-group symmetry set to C1

NOTHIEL

Do not use Thiel's FSTMIN technique

NOXYZ

Do not print Cartesian coordinates

NSPA=n

Sets number of geometric segments in COSMO

NSURF

Number of surfaces in an ESP calculation

OCTET

Octet state required

OLDCAV

In COSMO, use the old Solvent Accessible Surface calculation

OLDENS

Read initial density matrix off disk

OLDFPC

Use the old fundamental physical constants

OLDGEO

Previous geometry to be used

OMIN=n.nn

In TS, minimum allowed overlap of eigenvectors

OPEN(n1,n2)

Open-shell RHF calculation requested

OPT-C

Optimize all carbon atom coordinates

OPT-H

Optimize all hydrogen atom coordinates

OPT-N

Optimize all nitrogen atom coordinates

OPT-O

Optimize all oxygen atom coordinates

OPT-S

Optimize all sulfur atom coordinates

OPTB  <not implemented>
OPTR  <not implemented>

P=n.nn

An applied pressure of n.nn Newtons/m2 to be used

PDB

Input geometry is in protein data bank format

PDB=(text)

User defined chemical symbols in protein data base

PDBOUT

Output geometry in pdb format

PECI

C.I. involves paired excitations only

PI

Resolve density matrix into s, p, and d components

pKa Print the pKa for ionizable hydrogen atoms attached to oxygen atoms

PL

Monitor convergence of density matrix in ITER

PM3

Use the MNDO-PM3 Hamiltonian

PM6

Use the PM6 Hamiltonian

PMEP

Complete semiempirical MEP calculation

PMEPR

Complete semiempirical MEP in a plane to be defined

POINT=n

Number of points in reaction path

POINT1=n

Number of points in first direction in grid calculation

POINT2=n

Number of points in second direction in grid calculation

POLAR

Calculate first, second and third order polarizabilities

POTWRT

In ESP, write out electrostatic potential to unit 21

POWSQ

Print details of working in POWSQ

PRECISE

Criteria to be increased by 100 times

PRESSURE Apply pressure or tension to a solid or polymer

PRNT=n

Print details of geometry optimization in EF

PRTINT

Print interatomic distances

PRTMEP

MEP contour data output to <filename>.mep

PULAY

Use Pulay's converger to obtain a SCF

QPMEP

Charges derived from Wang-Ford type AM1 MEP

QUARTET

Quartet state required

QUINTET

Quintet state required

RAPID In MOZYME geometry optimizations, only use atoms being optimized in the SCF

RECALC=n

In EF, recalculate Hessian every n steps

RELSCF

Default SCF criterion multiplied by n

RESEQ Re-arrange the atoms to match the PDB convention
RESIDUES Label each atom in a polypeptide with the amino acid residue

RESTART

Calculation restarted

RHF Use Restricted Hartree-Fock methods
RM1 Use the RM1 Hamiltonian

RMAX=n.nn

In TS, maximum allowed ratio for energy change

RMIN=n.nn

In TS, minimum allowed ratio for energy change

ROOT=n

Root n to be optimized in a C.I. calculation

RSCAL

In EF, scale p-RFO to trust radius

RSOLV=n.nn

Effective radius of solvent in COSMO

SADDLE

Optimize transition state

SCALE

Scaling factor for van der waals distance in ESP

SCFCRT=n.nn

Default SCF criterion replaced by the value supplied

SCINCR=n.nn

Increment between layers in ESP

SEPTET

Septet state required

SETUP

Extra keywords to be read from setup file

SEXTET

Sextet state required

SHIFT=n.nn

a damping factor of n defined to start SCF

SIGMA

Minimize gradients using SIGMA

SINGLET

Singlet state required

SLOG=n.nn

In BFGS optimization, use fixed step of length n .nn

SLOPE

Multiplier used to scale MNDO charges

SNAP

Increase precision of symmetry angles

SPIN

Print final UHF spin matrix

STATIC

Calculate Polarizability using electric fields

STEP

Step size in path

STEP1=n.nnn

Step size n for first coordinate in grid calculation

STEP2=n.nnn

Step size n for second coordinate in grid calculation

STO3G

Deorthogonalize orbitals in STO-3G basis

SUPER

Print superdelocalizabilities

SYBYL

Output a file for use by Tripos's SYBYL program

SYMAVG

Average symmetry equivalent ESP charges

SYMOIR

Print characters of eigenvectors and print number of I.R.s

SYMTRZ

Print details of working in subroutine SYMTRZ.

SYMMETRY

Impose symmetry conditions

T=n[M,H,D]

A time of n seconds requested

THERMO
THERMO(nnn)
THERMO(nnn,mmm)
THERMO(nnn,mmm,lll)

Perform a thermodynamics calculation

TIMES

Print times of various stages

T-PRIORITY
T-PRIORITY=n.nn

Time takes priority in DRC

TRANS
TRANS=n

The system is a transition state (used in thermodynamics calculation)

TRIPLET

Triplet state required

TS

Using EF routine for TS search

UHF

Use the Unrestricted Hartree-Fock method

UNSAFE

Use reduced memory (may affect SCF convergence)

VDW(text)

Van der waals radius for atoms in COSMO, ESP, and PMEP defined by user

VECTORS

Print final eigenvectors

VELOCITY

Supply the initial velocity vector in a DRC calculation

WILLIAMS

Use Williams surface

X-PRIORITY=n.nn

Geometry changes take priority in DRC

XENO Allow non-standard residues in proteins to be labeled.

XYZ

Do all geometric operations in Cartesian coordinates

Z=n Number of mers in a cluster