Serge Gorelsky || AOMix user registration | Publications | FAQ | AOMix | AOMix-CDA | ALP-Vibro | RedS | SWizard | SIG-JK | U2WIN

AOMix
Last updated: The most current version of AOMix:
January 19, 2010 6.46 | software updates and changes | software download | feedback

Events: March 20, 2009: the presentation about the AOMix software at Q-Chem workshop at Pittsburgh Supercomputing Center

AOMix is a user-friendly software for the molecular orbital (MO) analysis. It calculates the MO compositions in terms of the constituent chemical fragments (you can specify them as atoms, groups of atoms, atomic orbitals, fragment molecular orbitals, groups of atomic orbitals, etc.) in the molecule or atom. AOMix automatically processes output files of multiple quantum-chemical packages (see the list below). AOMix also allows you to analyze chemical structure (bonding/antibonding nature of molecular orbitals) using overlap populations (total and per molecular orbital), valence indices, 2-center (Mayer, Lowdin, Wiberg, and bond-order symmetry components) and 3-, 4-, 5- and 6-center bond orders, charge decomposition analysis, CDA (total, per molecular orbital, per symmetry type), and condensed Fukui functions.

Partial and Total Density of States Plot Overlap Population Density of States Plot, also known as Crystal Orbital Overlap Population (COOP) diagram

AOMix generates total, partial, and overlap population density-of-states (DOS) plots (these are available in line-plot and continuous-plot representations, see Figures 1, 2, and 3). Overlap population DOS plots are also refered in the literature as Crystal Orbital Overlap Population (COOP) diagrams.

A very important part of quantum chemical research is to interpret the results in terms of qualitative concepts. This software was developed to help in the analysis of the nature of the chemical bonding in molecular systems and to monitor changes in the electron density distribution upon the electron excitation. Let's say, there is a band in an absorption spectrum of a molecule or an ion at 400 nm which is assigned to a HOMO -> LUMO+2 electron excitation. What does it tell us about properties of this molecule/ion? What do we know about the nature of the corresponding excited state? What will happen with this molecule / ion after the photoexcitation? AOMix helps to answer these questions.

With AOMix, you can


AOMix

  1. calculates the Wiberg bond order indices between molecular fragments (calculated in the canonical orbital basis: (P*S) x (P*S) and in the orthonormal Lowdin orbital basis);

  2. calculates the Mayer bond orders (including the bond order contributions from alpha- and beta-spin occupied orbitals for open-shell systems) between molecular fragments;

  3. if molecular symmetry is present, calculates symmetry components (for each irreducible representation) of Mayer bond orders (including the bond order contributions from alpha- and beta-spin MOs for open-shell systems) between fragments. Symmetry decomposition of bond orders helps in picturing electronic structure and contributions to chemical bonding;

  4. calculates total and free valence of atoms / fragments;

  5. calculates the eigenvalues and eigenvectors of the AO overlap matrix;

  6. performs Lowdin population analysis (LPA) and other PA methods which are based on the Sa x P x S1-a formula (where P is the density matrix and S is the overlap matrix);

  7. calculates the 3-, 4-, 5- and 6-center bond order indices between molecular fragments to help in finding possible multi-center orbital interations.

    AOMix output files also contain:

    • Lowdin and Mulliken alpha-spin electron populations,
    • Lowdin and Mulliken beta-spin electron populations,
    • Lowdin and Mulliken gross electron populations,
    • Lowdin and Mulliken spin densities,
    • (for spin-unrestricted calculations) the overlap matrix between beta-spin molecular orbitals and alpha-spin molecular orbitals (the so-called mutual overlap matrix).

    Here is one of the AOMix output files. In this example, fragments are defined by automatically by the program (each atom is a fragment).

    In addition, if each atom is a fragment, AOMix writes two UCSF Chimera attribute files. The first file (with the atomic attributes) allows easy visualization of the calculated atom-based properties:

    • atomic charges,
    • MPA, NPA and LPA spin densities,
    • total and free valence indices of atoms.
    Atomic valence indices in the anthracene molecule

    Atomic valence indices in the anthracene molecule (C14H10)

    The second file (with the PseudoBond attributes) allows easy visualization of the calculated two-center properties:

    • Lowdin and Wiberg bond orders,
    • Mayer bond orders and its components.
Mayer bond orders in the anthracene molecule

Mayer bond orders for the C-C and C-H bonds in the anthracene molecule (C14H10)

AOMix also prints the S0.5 x P x S0.5 matrix for alpha- and beta-spin orbitals.

AOMix can be also used to recover the initial guess from the computed wave function (see Appendix I of the AOMix manual).


AOMix provides:

  1. the MO composition analysis in terms of
    • ALL atomic orbitals in one run (number of fragments = number of molecular orbitals);
    • ALL atoms in one run (number of fragments = number of atoms), these include the automatic s,p,d,f-orbital contribution breakdown (see the example);
    • or any large number of molecular fragments (defined as a list of atoms or/and atomic orbitals);
  2. total and partial density-of-states (TDOS and PDOS) plot data (both continuos and line plots are available);
  3. overlap populations density-of-states (OPDOS) plot data;
  4. MS Windows and UNIX scripts for Gaussian CUBE file generation;
  5. in addition to the regular AOMix output file (for atoms as individual fragments and for atomic orbitals as individual fragments), the HTML file (AOMix-MPA-orb.html) is generated. It contains the information about MO compositions in terms of MOST IMPORTANT atomic orbitals (6 largest contributions with values greater than 0.5%)

Molecular orbitals in the line and continuous DOS plot representations

You can refer to the AOMix manual and these publications to learn about the AOMix capabilities.


The following procedures for the electron population analysis are available in the AOMix package:

  1. Mulliken population analysis (MPA)
  2. modified Mulliken population analysis (MMPA)
  3. C2 population analysis (SCPA)
  4. Lowdin population analysis (LPA) and other methods based on the Sa x P x S1-a formula
  5. charge decomposition analysis (CDA) (using AOMix-CDA)

Supported operating systems  AOMix
MS Windows 2000 +
MS Windows XP +
MS Windows Vista + (*)
MS Windows 7 + (*)

* For processing of Gaussian and Turbomole output files only.




How to refer to the AOMix program in publications:

Compositions of molecular orbitals, overlap populations between molecular fragments, bond orders and density-of-states spectra were calculated using the AOMix program [1,2].

  1. S. I. Gorelsky, AOMix: Program for Molecular Orbital Analysis, http://www.sg-chem.net/, University of Ottawa, 2009.
  2. S. I. Gorelsky, A. B. P. Lever, J. Organomet. Chem. 2001, 635, 187-196.

AOMix can process output files from the following quantum-chemistry packages
(other packages can be also included)

ADF (Scientific Computing & Modelling NV).
AOMix has been tested with the following versions: ADF 2000.01, ADF 2002.03, ADF 2003.01, ADF 2004.01, ADF 2005.01 and ADF 2009.01; please note that ADF calculations with core functions (do not confuse core functions with core orbitals!) cannot be processed by the current version of AOMix.
CNDO/INDO CNDO/INDO (by J.R. Reimers, University of Sydney, Australia)

The General Atomic and Molecular Electronic Structure System GAMESS from Iowa State University and PC GAMESS from Moscow State University.
AOMix has been tested with the following versions: GAMESS-US 6, GAMESS VERSION = 12.
Gaussian 98, Gaussian 03 and Gaussian 09 (Gaussian, Inc.)
HyperChem 4.x-7.x (HyperCube, Inc.)
Jaguar (Schrodinger, Inc.).
AOMix has been tested with the following versions: Jaguar 3.5, Jaguar 4.1, Jaguar 4.2, Jaguar 5.0, Jaguar 5.5, and Jaguar 6.0
MOPAC2009 (James J. P. Stewart, Stewart Computational Chemistry)
AOMix has been tested with the following versions: MOPAC2009, Version 8.311W.
Q-Chem (Q-Chem, Inc.).
AOMix has been tested with the following versions: Q-Chem 3.2
Spartan (Wavefunction, Inc.).
AOMix has been tested with the following versions: Spartan '02, Spartan '03, and Spartan '06.
Turbomole Turbomole (COSMOlogic).
AOMix has been tested with the following versions: TURBOMOLE 6.0
ZINDO ZINDO (by M.C. Zerner, Quantum Theory Project, University of Florida; ZINDO is available in Cerius2 (Accelrys Inc.) and CAChe (CAChe Group, Fujitsu));
AOMix has been tested with the following versions: ZINDO 98.1 and ZINDO 99.1


How to obtain a copy of AOMix:

If you would like to obtain the AOMix program for yourself, your research group or company, please follow these instructions.

WWW.SG-CHEM.NET
All contents copyright © 1997-2010 S. I. Gorelsky
Last updated: January 2010