Semiempirical calculations, design

Once a PES has been computed, it is often fitted to an analytic function. This is done because there are many ways to analyze analytic functions that require much less computation time than working directly with ah initio calculations. For example, the reaction can be modeled as a molecular dynamics simulation showing the vibrational motion and reaction trajectories as described in Chapter 19. Another technique is to fit ah initio results to a semiempirical model designed for the purpose of describing PES s. 

Alchemy 2000 (we tested Version 2.05) is a graphic interface for running molecular mechanics and semiempirical calculations. Calculations can be done with the built-in Tripos force field or by calling the MM3 or MOPAC programs, which are included with the package. Alchemy is designed by Tripos and sold by SciVision. 

T vo main streams of computational techniques branch out fiom this point. These are referred to as ab initio and semiempirical calculations. In both ab initio and semiempirical treatments, mathematical formulations of the wave functions which describe hydrogen-like orbitals are used. Examples of wave functions that are commonly used are Slater-type orbitals (abbreviated STO) and Gaussian-type orbitals (GTO). There are additional variations which are designated by additions to the abbreviations. Both ab initio and semiempirical calculations treat the linear combination of orbitals by iterative computations that establish a self-consistent electrical field (SCF) and minimize the energy of the system. The minimum-energy combination is taken to describe the molecule. 

Reasonable correlation between measured and calculated TPA was observed for a series of chromophores with medium size [123, 131, 132, 224, 227, 233, 284, 285, 289, 295-300]. However, problems can arise for larger molecules requiring larger Cl space. Under these circumstances, the calculated quantities A/(0. M(y, M02, oi, and 02 may not reliably describe TPA according to Eq. (18). Despite these problems, semiempirical calculations of Im[ (— o. co)] have been the most applied method for a theoretical design of chromophores with large TP absorptivity [131, 132, 233, 283-285, 287-289, 301]. Hence, ZINDO/S-MRDCI has become an attractive theoretical tool for the calculation of TPA due to the modest calculation time. 

Having established that intramolecular diene cyclozirconation can be carried out under conditions of either kinetic or thermodynamic control, and having shown that semiempirical calculations (ZINDO) can be used to predict the relative stabilities of diastereomeric zirconacycles, we next undertook the computationally based design of a diene such that cyclozirconation would be directed toward a desired diastereomer. 

C. Aleman and M. Orozco,/. Comput.-AidedMol. Design, 6,331 (1992). On the Suitability of Semiempirical Calculations as Sources of Force Field Parameters. 

Ah initio trajectory calculations have now been performed. However, these calculations require such an enormous amount of computer time that they have only been done on the simplest systems. At the present time, these calculations are too expensive to be used for computing rate constants, which require many trajectories to be computed. Semiempirical methods have been designed specifically for dynamics calculations, which have given insight into vibrational motion, but they have not been the methods of choice for computing rate constants since they are generally inferior to analytic potential energy surfaces fitted from ah initio results. 

One very popular technique is an adaptation of the Born model for orbital-based calculations by Cramer and Truhlar, et. al. Their solvation methods (denoted SMI, SM2, and so on) are designed for use with the semiempirical and ah initio methods. Some of the most recent of these methods have a few parameters that can be adjusted by the user in order to customize the method for a specific solvent. Such methods are designed to predict ACsoiv and the geometry in solution. They have been included in a number of popular software packages including the AMSOL program, which is a derivative of AMPAC created by Cramer and Truhlar. 

The conductor-like screening model (COSMO) is a continuum method designed to be fast and robust. This method uses a simpler, more approximate equation for the electrostatic interaction between the solvent and solute. Line the SMx methods, it is based on a solvent accessible surface. Because of this, COSMO calculations require less CPU time than PCM calculations and are less likely to fail to converge. COSMO can be used with a variety of semiempirical, ah initio, and DFT methods. There is also some loss of accuracy as a result of this approximation. 

Most of the semiempirical methods are not designed to correctly predict the electronic excited state. Although excited-state calculations are possible, particularly using a CIS formulation, the energetics are not very accurate. However, the HOMO-LUMO gap is reasonably reproduced by some of the methods. 

Another search trial [54] for the above two molecules was done by using additional information on the charges of the atoms in each molecule, charge in this case means the net charge obtained by a semiempirical molecular orbital calculation. Two modes of charge designation can be used in the system. 

---