17.11.2014 Views

ChemOffice.Com - CambridgeSoft

ChemOffice.Com - CambridgeSoft

ChemOffice.Com - CambridgeSoft

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

esults in better estimates of the energy in systems<br />

where energy levels are closely spaced, and where<br />

bond breaking is occurring.<br />

UHF can be run on both open and closed shell<br />

systems. The major caveat to this method is the<br />

time involved. Since alpha and beta electrons are<br />

treated separately, twice as many integrals need to<br />

be solved. As your models get large, the time for the<br />

computation may make it a less satisfactory<br />

method.<br />

Configuration Interaction<br />

The effects of electron-electron repulsion are<br />

underestimated by SCF-RHF methods, which<br />

results in the overestimation of energies.<br />

SCF-RHF calculations use a single determinant that<br />

includes only the electron configuration that<br />

describes the occupied orbitals for most molecules<br />

in their ground state. Further, each electron is<br />

assumed to exist in the average field created by all<br />

other electrons in the system, which tends to<br />

overestimate the repulsion between electrons.<br />

Repulsive interactions can be minimized by<br />

allowing the electrons to exist in more places (i.e.<br />

more orbitals, specifically termed virtual orbitals).<br />

The multi-electron configuration interaction<br />

(MECI) method in MOPAC addresses this<br />

problem by allowing multiple sets of electron<br />

assignments (i.e., configurations) to be used in<br />

constructing the molecular wave functions.<br />

Molecular wave functions representing different<br />

configurations are combined in a manner analogous<br />

to the LCAO approach.<br />

For a particular molecule, configuration interaction<br />

uses these occupied orbitals as a reference electron<br />

configuration and then promotes the electrons to<br />

unoccupied (virtual) orbitals. These new states,<br />

Slater determinants or microstates in MOPAC, are<br />

then linearly combined with the ground state<br />

configuration. The linear combination of<br />

microstates yields an improved electronic<br />

configuration and hence a better representation of<br />

the molecule.<br />

Approximate Hamiltonians<br />

in MOPAC<br />

There are five approximation methods available in<br />

MOPAC:<br />

• AM1<br />

• MNDO<br />

• MNDO-d<br />

• MINDO/3<br />

• PM3<br />

The potential energy functions modify the HF<br />

equations by approximating and parameterizing<br />

aspects of the Fock matrix. The approximations in<br />

semi-empirical MOPAC methods play a role in the<br />

following areas of the Fock operator:<br />

• The basis set used in constructing the 1-<br />

electron atom orbitals is a minimum basis set<br />

of only the s and p Slater Type Orbitals (STOs)<br />

for valence electrons.<br />

• The core electrons are not explicitly treated.<br />

Instead they are added to the nucleus. The<br />

nuclear charge is termed N effective .<br />

For example, Carbon as a nuclear charge of<br />

+6-2 core electrons for a effective nuclear<br />

charge of +4.<br />

• Many of the 2-electron Coulomb and<br />

Exchange integrals are parameterized based on<br />

element.<br />

Choosing a Hamiltonian<br />

Overall, these potential energy functions may be<br />

viewed as a chronological progression of<br />

improvements from the oldest method, MINDO/3<br />

to the newest method, PM5. However, although the<br />

<strong>ChemOffice</strong> 2005/Chem3D MOPAC <strong>Com</strong>putations • 163<br />

MOPAC Semi-empirical Methods

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!