2/2 CASSCF

Roos B O, Taylor P R and Siegbahn P E M 1980 A complete active space SCF method (CASSCF) using a density matrix formulated super-CI approach Chem. Phys. 48 157-73... 

Techniques have been developed within the CASSCF method to characterize the critical points on the excited-state PES. Analytic first and second derivatives mean that minima and saddle points can be located using traditional energy optimization procedures. More importantly, intersections can also be located using constrained minimization [42,43]. Of particular interest for the mechanism of a reaction is the minimum energy path (MEP), defined as the line followed by a classical particle with zero kinetic energy [44-46]. Such paths can be calculated using intrinsic reaction coordinate (IRC) techniques... 

For this reason, there has been much work on empirical potentials suitable for use on a wide range of systems. These take a sensible functional form with parameters fitted to reproduce available data. Many different potentials, known as molecular mechanics (MM) potentials, have been developed for ground-state organic and biochemical systems [58-60], They have the advantages of simplicity, and are transferable between systems, but do suffer firom inaccuracies and rigidity—no reactions are possible. Schemes have been developed to correct for these deficiencies. The empirical valence bond (EVB) method of Warshel [61,62], and the molecular mechanics-valence bond (MMVB) of Bemardi et al. [63,64] try to extend MM to include excited-state effects and reactions. The MMVB Hamiltonian is parameterized against CASSCF calculations, and is thus particularly suited to photochemistry. 

To use direct dynamics for the study of non-adiabatic systems it is necessary to be able to efficiently and accurately calculate electronic wave functions for excited states. In recent years, density functional theory (DFT) has been gaining ground over traditional Hartree-Fock based SCF calculations for the treatment of the ground state of large molecules. Recent advances mean that so-called time-dependent DFT methods are now also being applied to excited states. Even so, at present, the best general methods for the treatment of the photochemistry of polyatomic organic molecules are MCSCF methods, of which the CASSCF method is particularly powerful. 

CASSCF is a version of MCSCF theory in which all possible configurations involving the active orbitals are included. This leads to a number of simplifications, and good convergence properties in the optimization steps. It does, however, lead to an explosion in the number of configurations being included, and calculations are usually limited to 14 elections in 14 active orbitals. 

A simple example would be in a study of a diatomic molecule that in a Hartree-Fock calculation has a bonded cr orbital as the highest occupied MO (HOMO) and a a lowest unoccupied MO (LUMO). A CASSCF calculation would then use the two a electrons and set up four CSFs with single and double excitations from the HOMO into the a orbital. This allows the bond dissociation to be described correctly, with different amounts of the neutral atoms, ion pair, and bonded pair controlled by the Cl coefficients, with the optimal shapes of the orbitals also being found. For more complicated systems... 

The present high cost of full CASSCF direct dynamics means that it is not possible to use such calculations to run large numbers of trajectories. As a result it cannot be used to build up experience of the types of effects to be found in dynamical studies of organic photochemistry, and in their interpretation. This problem can be remedied by performing calculations using the MMVB force field [63,64]. 

Jones et al. [144,214] used direct dynamics with semiempirical electronic wave functions to study electron transfer in cyclic polyene radical cations. Semiempirical methods have the advantage that they are cheap, and so a number of trajectories can be run for up to 50 atoms. Accuracy is of course sacrificed in comparison to CASSCF techniques, but for many organic molecules semiempirical methods are known to perform adequately. 

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

The energies of this Cl and of the other ones calculated in this work are listed in Table III. The calculated CASSCF values of the energies of the two lowest electronically states are 9.0 eV (5i, vertical) and 10.3 eV ( 2, vertical) [99]. They are considerably higher than the expenmental ones, as noted for this method by other workers [65]. In all cases, the computed conical intersections lie at much lower energies than the excited state, and are easily accessible upon excitation to Si. In the case of the H/allyl Cl, the validity confirmation process recovered the CHDN and 1,3-CHDN anchors. An attempt to approach the third anchor [BCE(I)] resulted instead in a biradical, shown in Figure 43. The bhadical may be regarded as a resonance hybrid of two allyl-type biradicals. 

The CASSCF(8,8)/DZV Energies of Some Stable Molecules and Conical Intersections Relevant to 1,4-cyclohexadiene (CHDN) Photochemistry (kcal/mole with respect to CHDN)... 

Another approach is spin-coupled valence bond theory, which divides the electrons into two sets core electrons, which are described by doubly occupied orthogonal orbitals, and active electrons, which occupy singly occupied non-orthogonal orbitals. Both types of orbital are expressed in the usual way as a linear combination of basis functions. The overall wavefunction is completed by two spin fimctions one that describes the coupling of the spins of the core electrons and one that deals with the active electrons. The choice of spin function for these active electrons is a key component of the theory [Gerratt ef al. 1997]. One of the distinctive features of this theory is that a considerable amount of chemically significant electronic correlation is incorporated into the wavefunction, giving an accuracy comparable to CASSCF. An additional benefit is that the orbitals tend to be... 

Fig. 5.35 Geometry predicted by CASSCF ab initio calculations of the two possible transition structure geometries for the Diels-Alder reaction between ethene and butadiene. (Figure adapted from Houk KN, J Gonzalez and Y Li 1995. Pericyclic Reaction Transition States Passions and Punctilios 1935-1995. Accounts of Chemical Research 28 81-90.)...

---