Hamiltonian semiempirical

The cornerstone of semiempirical and ab initio molecular orbital methods is the Harhee equation and its extensions and variants, the Harhee-Fock and Roothaan-Hall equations. We have seen that the Hamiltonian for the hydrogen atom. 

Semiempirical calculations are set up with the same general structure as a HF calculation in that they have a Hamiltonian and a wave function. Within this framework, certain pieces of information are approximated or completely omitted. Usually, the core electrons are not included in the calculation and only a minimal basis set is used. Also, some of the two-electron integrals are omitted. In order to correct for the errors introduced by omitting part of the calculation, the method is parameterized. Parameters to estimate the omitted values are obtained by fitting the results to experimental data or ah initio calculations. Often, these parameters replace some of the integrals that are excluded. 

Unlike semiempirical methods that are formulated to completely neglect the core electrons, ah initio methods must represent all the electrons in some manner. However, for heavy atoms it is desirable to reduce the amount of computation necessary. This is done by replacing the core electrons and their basis functions in the wave function by a potential term in the Hamiltonian. These are called core potentials, elfective core potentials (ECP), or relativistic effective core potentials (RECP). Core potentials must be used along with a valence basis set that was created to accompany them. As well as reducing the computation time, core potentials can include the effects of the relativistic mass defect and spin coupling terms that are significant near the nuclei of heavy atoms. This is often the method of choice for heavy atoms, Rb and up. 

The types of algorithms described above can be used with any ah initio or semiempirical Hamiltonian. Generally, the ah initio methods give better results than semiempirical calculations. HE and DFT calculations using a single deter-... 

The most commonly used semiempirical for describing PES s is the diatomics-in-molecules (DIM) method. This method uses a Hamiltonian with parameters for describing atomic and diatomic fragments within a molecule. The functional form, which is covered in detail by Tully, allows it to be parameterized from either ah initio calculations or spectroscopic results. The parameters must be fitted carefully in order for the method to give a reasonable description of the entire PES. Most cases where DIM yielded completely unreasonable results can be attributed to a poor fitting of parameters. Other semiempirical methods for describing the PES, which are discussed in the reviews below, are LEPS, hyperbolic map functions, the method of Agmon and Levine, and the mole-cules-in-molecules (MIM) method. 

The ah initio methods available are RHF, UHF, ROHE, GVB, MCSCF along with MP2 and Cl corrections to those wave functions. The MNDO, AMI, and PM3 semiempirical Hamiltonians are also available. Several methods for creating localized orbitals are available. 

The final application considered in this chapter is chosen to illustrate the application of a QM-MM study of an enzyme reaction that employs an ab initio Hamiltonian in the quantum region [67]. Because of the computational intensity of such calculations there are currently very few examples in the literahire of QM-MM shidies that use a quanhim mechanical technique that is more sopliisticated than a semiempirical method. MuUiolland et al. [67] recently reported a study of part of the reaction catalyzed by citrate synthase (CS) in wliich the quanhim region is treated by Hartree-Fock and MP2 methods [10,51],... 

The systems discussed in this chapter give some examples using different theoretical models for the interpretation of, primarily, UPS valence band data, both for pristine and doped systems as well as for the initial stages of interface formation between metals and conjugated systems. Among the various methods used in the examples are the following semiempirical Hartree-Fock methods such as the Modified Neglect of Diatomic Overlap (MNDO) [31, 32) and Austin Model 1 (AMI) [33] the non-empirical Valence Effective Hamiltonian (VEH) pseudopotential method [3, 34J and ab initio Hartree-Fock techniques. 

CAChe 5.0, available in 2002, includes a new, more powerful, semiempirical method that uses the PM5 Hamiltonian, a MOPAC 2002 offering, modeling of molecules with up to 20,000 atoms, the inclusion of all main group elements in one semiempirical method, and using MOPAC AMl-d, supports the transition metals Pt, Fe, Cu, Ag, Mo, V, and Pd. Researchers can now import and display, in 3D, proteins from the Protein Data Bank (PDB), optimize proteins, dock ligands, and model reactions on protein molecules. 

Hqm is the Hamiltonian of the QM region and might be based on semiempirical, ah initio molecular orbital or density functional theory (DFT) methods. Hqm/mm represents the interactions between the QM and MM regions, Hmm is the Hamiltonian of the purely MM region, and Hboundary is Hamiltonian for the boundary of the system, if this contribution is included. The corresponding total energy ( tot) of the QM/MM system is ... 

Ab initio molecular orbital methodology or density functional theory [158-160] would be suited for this combined QM/MM approach. However, in order to be able to compute the QM energies along the Monte Carlo simulation, nowadays a semiempirical Hamiltonian, like AMI [161], is a much more computationally efficient method. Before using AMI, the goodness of the semiempirical results in gas phase in comparison with the ab initio ones has to be tested. For systems in which the semiempirical results are poor, the relation... 

In semiempirical calculations variational principle itself is manipulated, or to give another interpretation, the molecular Hamiltonian is modified in order to save computation time. Two semiempirical methods differing in the concept and the degree of sophistication have been applied frequently to ion-molecule complexes and therefore will be mentioned briefly. 

The second kind of semiempirical procedure mentioned here is even cruder. In the extended Hueckel method (EHM 64<65)) the electronic structure of the molecule is simulated by an effective Hamiltonian. The total energy of the molecule is represented by a sum of one electron energies and even the nuclear repulsion terms are not taken into account explicitly. This type of approach can be shown to give an approximate idea of electronic structures and relative energies of unpolar molecules like hydrocarbons, but it fails inevitably when applied to structures with appreciable polarity 66>. Therefore any application of EHM calculations to interactions between polar molecules or ions should be regarded with a good deal of scepticism. 

Small monosaccharides have molecular sizes at the upper limit of the range that is currently treatable with initio methods. An exaiqple of the application of initio calculations to carbohydrates is given in the paper by Garrett and Serianni in this volume. Semiempirical quantum mechanical calculations, which use simplified molecular Hamiltonians with parameters taken from experiment, extend quantum mechanical calculations to larger molecules. However, the reliability is reduced compared to the best ab initio results. 

The first two-step calculations of the P,T-odd spin-rotational Hamiltonian parameters were performed for the PbF radical about 20 years ago [26, 27], with a semiempirical accounting for the spin-orbit interaction. Before, only nonrelativistic SCF calculation of the TIF molecule using the relativistic scaling was carried out [86, 87] here the P,T-odd values were underestimated by almost a factor of three as compared to the later relativistic Dirac-Fock calculations. The latter were first performed only in 1997 by Laerdahl et al. [88] and by Parpia [89]. The next two-step calculation, for PbF and HgF molecules [90], was carried out with the spin-orbit RECP part taken into account using the method suggested in [91]. 

---