Note
This website is presently under development. The primary website for MOPAC is still http://openmopac.net until the majority of its contents have been migrated here.
AUX
¶
AUX[([n][,][COMP][,][PRECISION=m][,][MOS=n][,][XP][,][XS][,][XW])]
When
AUX
is used then auxiliary information is output to a file
called <file>.aux. When AUX
(n) is used, output for the
changing geometry is re-directed to output channel n, n is limited
to the range 0 - 99. Thus AUX
(6) would re-direct the changing
geometry to standard output unit 6. Useful values are n = 0 and n
= 6. Avoid using units 2 to 4, and 9 to 30, as these are used by
MOPAC. The <file>.aux files are intended for use by Graphical User
Interfaces (GUI’s), but the files can be converted using a simple
Windows utility, Read_AUX, into ASCII
tables for use by programs such as EXCEL.AUX is supported by Avogadro, Gabedit, and Molden.
When
COMP
is present, all large blocks of data are printed in
compressed form.To increase precision by m digits (m limited to the range 1 to 9),
use
PRECISION=m
.By default, only the 10 highest occupied and 10 lowest unoccupied
M.O.s are printed, unless
LARGE
is present, in which case all
M.O.s are printed To change the default, use MOS=n
, where n is
the number of occupied and virtual M.O.s to be printed. If n is
greater than the number of M.O.s, it will be re-set to all the M.O.s,
so an alternative way to print all the M.O.s is to use MOS=99999
.
To prevent printing of M.O.’s, use MOS=-99999
.The commas are optional.
If the re-direction unit “n” is used together with any of COMP,
PRECISION, or MOS, then “n” must come first, followed by the other
options.
For large systems, the AUX file might still be very large, even if
COMP is used. To reduce the size of the AUX file, use one or more of
the following options:
If XP is present, the density matrix will not be printed.
If XS is present, the overlap matrix will not be printed.
If XW is present, the M.O.s will not be printed. This is an
alternative to
MOS=-99999
.Only data considered to be important to GUI programs will be included in the .aux file. A consequence of this is that some data might be missing. If such data can be identified, and a message sent to MrMOPAC@ATT.net, then the data will be added to the .aux file.
The auxiliary information is a very compact form of the results of a
calculation in a primitive mark-up language. This file is intended for
use by other programs. An example of the auxiliary file for water is
given below. No explanation is necessary: the file is intended to be
self-descriptive.
START OF MOPAC FILE
####################################
# #
# Start of Input data #
# #
####################################
MOPAC_VERSION=MOPAC2007.7.150W
DATE=”Mon Jun 4 12:15:37 2007”
METHOD=PM6
TITLE=” Water”
KEYWORDS=” SYMMETRY”
ATOM_EL[0003]=
O H H
ATOM_X:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.6018 0.7647 0.0000
0.6018 -0.7647 0.0000
AO_ATOMINDEX[0006]=
1 1 1 1 2 3
ATOM_SYMTYPE[0006]=
S PX PY PZ S S
AO_ZETA[0006]=
5.4218 2.2710 2.2710 2.2710 1.2686 1.2686
ATOM_PQN[0006]=
2 2 2 2 1 1
NUM_ELECTRONS=0008
EMPIRICAL_FORMULA=”H2 O”
####################################
# #
# Geometry optimization #
# #
####################################
HEAT_OF_FORM_UPDATED:KCAL/MOL=-0.540108D+02
GRADIENT_UPDATED:KCAL/MOL/ANG=+0.210682D+02
ATOM_X_UPDATED:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.6018 0.7647 0.0000
0.6018 -0.7647 0.0000
HEAT_OF_FORM_UPDATED:KCAL/MOL=-0.542600D+02
GRADIENT_UPDATED:KCAL/MOL/ANG=+0.893385D+01
ATOM_X_UPDATED:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.5693 0.7558 0.0000
0.5693 -0.7558 0.0000
HEAT_OF_FORM_UPDATED:KCAL/MOL=-0.543026D+02
GRADIENT_UPDATED:KCAL/MOL/ANG=+0.225094D+01
ATOM_X_UPDATED:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.5641 0.7630 0.0000
0.5641 -0.7630 0.0000
HEAT_OF_FORM_UPDATED:KCAL/MOL=-0.543063D+02
GRADIENT_UPDATED:KCAL/MOL/ANG=+0.144932D+00
ATOM_X_UPDATED:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.5610 0.7655 0.0000
0.5610 -0.7655 0.0000
####################################
# #
# Final SCF results #
# #
####################################
HEAT_OF_FORMATION:KCAL/MOL=-0.543063D+02
ENERGY_ELECTRONIC:EV=-0.458386D+03
ENERGY_NUCLEAR:EV=+0.139308D+03
POINT_GROUP=C2v
DIPOLE:DEBYE=+0.206870D+01
AREA:SQUARE ANGSTROMS=+0.424512D+02
VOLUME:CUBIC ANGSTROMS=+0.251444D+02
ATOM_X_OPT:ANGSTROMS[0009]=
0.0000 0.0000 0.0000
0.5610 0.7655 0.0000
0.5610 -0.7655 0.0000
ATOM_CHARGES[0003]=
-0.61867 +0.30933 +0.30933
OVERLAP_MATRIX[000021]=
# Lower half triangle only
1.0000 0.0000 1.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 1.0000
0.1618 0.2330 0.3180 0.0000 1.0000 0.1618 0.2330
-0.3180 0.0000 0.2333
1.0000
EIGENVECTORS[000036]=
-0.8288 -0.2451 0.0000 0.0000 -0.3556 -0.3556 0.0000
0.0000 0.8129 0.0000
0.4118 -0.4118 0.4394 -0.8417 0.0000 0.0000 -0.2220
-0.2220 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.3463 0.4812 0.0000
0.0000 -0.5694 -0.5694
0.0000 0.0000 0.5824 0.0000 -0.5748 0.5748
TOTAL_DENSITY_MATRIX[000021]=
# Lower half triangle only
1.7601 -0.3333 1.5370 0.0000 0.0000 1.3216 0.0000
0.0000 0.0000 2.0000
0.3944 0.5480 0.6695 0.0000 0.6907 0.3944 0.5480
-0.6695 0.0000 0.0123
0.6907
M.O.SYMMETRY_LABELS[0006]=
1a1 1b2 2a1 1b1 3a1 2b2
EIGENVALUES[0006]=
-30.402 -18.983 -14.111 -11.905 4.067 6.038
MOLECULAR_ORBITAL_OCCUPANCIES[00006]=
2.0000 2.0000 2.0000 2.0000 0.0000 0.0000
When
COMP
is present, the overlap and density matrices, and
eigenvectors, are printed in compressed form.Compare the previous uncompressed data with the following compressed
data to see how the compression is done.
OVERLAP_INDICES[013]=
1 3 6 10 15 21 18 13 20 17 12 16 11
OVERLAP_COEFFICIENTS[013]=
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 -0.3180 0.3180 0.2336 0.2332
0.2332 0.1619 0.1618
SET_OF_MOS= 1 6
MO_INDICES[004]=
1 6 5 2
MO_COEFFICIENTS[004]=
0.8287 0.3558 0.3557 0.2452
MO_INDICES[003]=
3 5 6
MO_COEFFICIENTS[003]=
0.8129 0.4118 -0.4118
MO_INDICES[004]=
2 1 5 6
MO_COEFFICIENTS[004]=
0.8416 -0.4396 0.2219 0.2219
MO_INDICES[001]=
4
MO_COEFFICIENTS[001]=
1.0000
MO_INDICES[004]=
5 6 2 1
MO_COEFFICIENTS[004]=
0.5694 0.5693 -0.4812 -0.3465
MO_INDICES[003]=
3 6 5
MO_COEFFICIENTS[003]=
0.5824 0.5749 -0.5748
DENSITY_MATRIX_INDICES[014]=
10 1 3 6 21 15 13 18 12 17 16 11 2 20
DENSITY_MATRIX_COEFFICIENTS[014]=
2.0000 1.7599 1.5369 1.3217 0.6908 0.6908 0.6695 -0.6695 0.5480 0.5480
0.3945 0.3945 -0.3334 0.0124
When COMP and MOS are present, a user-defined set of M.O.s is printed. The default is the 10 highest M.O.s and 10 lowest virtual M.O.s.
If LARGE is present, all the M.O.s are printed.