12.01.2015 Views

Chem3D Users Manual - CambridgeSoft

Chem3D Users Manual - CambridgeSoft

Chem3D Users Manual - CambridgeSoft

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Administrator<br />

The procedures assume you have a basic<br />

understanding of the computational concepts and<br />

terminology of semi-empirical methods, and the<br />

concepts involved in geometry optimization<br />

(minimization) and single-point computations. For<br />

more information see “Computation Concepts” on<br />

page 129.<br />

For help with MOPAC, see the online MOPAC<br />

manual at:<br />

http://www.cachesoftware.com/mopac/Mopac2002<br />

manual/<br />

MOPAC Semiempirical<br />

Methods<br />

The method descriptions that follow represent a<br />

very simplified view of the semi-empirical methods<br />

available in <strong>Chem3D</strong> and CS MOPAC. For more<br />

information see the online MOPAC manual.<br />

Extended Hückel Method<br />

Developed from the qualitative Hückel MO<br />

method, the Extended Hückel Method (EH)<br />

represents the earliest one-electron semi-empirical<br />

method to incorporate both σ and p valence<br />

systems. It is still widely used, owing to its versatility<br />

and success in analyzing and interpreting groundstate<br />

properties of organic, organometallic, and<br />

inorganic compounds of biological interest. Built<br />

into <strong>Chem3D</strong>, EH is the default semi-empirical<br />

method used to calculate data required for<br />

displaying molecular surfaces.<br />

The EH method uses a one-electron Hamiltonian<br />

with matrix elements defined as follows:<br />

H µµ = – I µ<br />

H µν = 0.5K( H µµ + H νν )S µν µ ≠ ν<br />

where I µ is the valence state ionization energy<br />

(VSIE) of orbital µ as deduced from spectroscopic<br />

data, and K is the Wolfsberg-Helmholtz constant<br />

(usually taken as 1.75). The Hamiltonian neglects<br />

electron repulsion matrix elements but retains the<br />

overlap integrals calculated using Slater-type basis<br />

orbitals. Because the approximated Hamiltonian<br />

(H) does not depend on the MO expansion<br />

coefficient C νi , the matrix form of the EH<br />

equations:<br />

can be solved without the iterative SCF procedure.<br />

RHF<br />

The default Hartree-Fock method assumes that the<br />

molecule is a closed shell and imposes spin<br />

restrictions. The spin restrictions allow the Fock<br />

matrix to be simplified. Since alpha (spin up) and<br />

beta (spin down) electrons are always paired, the<br />

basic RHF method is restricted to even electron<br />

closed shell systems.<br />

Further approximations are made to the RHF<br />

method when an open shell system is presented.<br />

This approximation has been termed the 1/2<br />

electron approximation by Dewar. In this method,<br />

unpaired electrons are treated as two 1/2 electrons<br />

of equal charge and opposite spin. This allows the<br />

computation to be performed as a closed shell. A CI<br />

calculation is automatically invoked to correct<br />

errors in energy values inherent to the 1/2 electron<br />

approximation. For more information see<br />

“Configuration Interaction” on page 167.<br />

With the addition of the 1/2 electron<br />

approximation, RHF methods can be run on any<br />

starting configuration.<br />

UHF<br />

HC=SCE<br />

The UHF method treats alpha (spin up) and beta<br />

(spin down) electrons separately, allowing them to<br />

occupy different molecular orbitals and thus have<br />

different orbital energies. For many open and<br />

closed shell systems, this treatment of electrons<br />

166•MOPAC Computations<br />

<strong>CambridgeSoft</strong><br />

MOPAC Semi-empirical Methods

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

Saved successfully!

Ooh no, something went wrong!