MCSCF analytical derivatives

The various response tensors are identified as terms in these series and are calculated using numerical derivatives of the energy. This method is easily implemented at any level of theory. Analytic derivative methods have been implemented using self-consistent-field (SCF) methods for a, ft and y, using multiconfiguration SCF (MCSCF) methods for ft and using second-order perturbation theory (MP2) for y". The response properties can also be determined in terms of sum-over-states formulation, which is derived from a perturbation theory treatment of the field operator — [iE, which in the static limit is equivalent to the results obtained by SCF finite field or analytic derivative methods. 

The ground state force field, vibrational normal modes and frequencies have been obtained with MCSCF analytic gradient and hessian calculations [176]. Frequencies computed with the DZ basis set are compared with experimental ones in Table 16. The T - So transition moments were obtained using distorted benzene geometries with atomic displacements along the normal modes, and with the derivatives in Eq. 97 obtained by numerical differentiation. The normal modes active for phosphorescence in benzene are depicted in Fig. 12. The final formula for the radiative lifetime of the k spin sublevel produced by radiation in all (i/f) bands is (ZFS representation x,y,z is used [49]) ... 

NACMEs based on state-averaged MCSCF wave functions and analytic derivative methods. This should provide NACME s of similar overall accuracy to that obtained for the adiabatic potentials. 

Importantly for direct dynamics calculations, analytic gradients for MCSCF methods [124-126] are available in many standard quantum chemistiy packages. This is a big advantage as numerical gradients require many evaluations of the wave function. The evaluation of the non-Hellmann-Feynman forces is the major effort, and requires the solution of what are termed the coupled-perturbed MCSCF (CP-MCSCF) equations. The large memory requirements of these equations can be bypassed if a direct method is used [233]. Modem computer architectures and codes then make the evaluation of first and second derivatives relatively straightforward in this theoretical framework. 

The selection rules for the QM harmonic oscillator pennit transitions only for An = 1 (see Section 14.5). As Eq. (9.47) indicates diat the energy separation between any two adjacent levels is always hm, the predicted frequency for die = 0 to n = 1 absorption (or indeed any allowed absorption) is simply v = o). So, in order to predict die stretching frequency within the harmonic oscillator equation, all diat is needed is the second derivative of the energy with respect to bond stretching computed at die equilibrium geometry, i.e., k. The importance of k has led to considerable effort to derive analytical expressions for second derivatives, and they are now available for HF, MP2, DFT, QCISD, CCSD, MCSCF and select other levels of theory, although they can be quite expensive at some of the more highly correlated levels of theoiy. 

Lengsfield III, B.H. and Yarkony, D.R. (1986). On the evaluation of nonadiabatic coupling matrix elements using SA-MCSCF/CI wavefunctions and analytic gradient methods. Ill Second derivative terms, J. Chem. Phys. 84, 348-353. 

The analytical expressions for the property derivatives may be obtained by substituting (n) for H(n) in the appropriate expressions given in Sections III— VI, since these expressions were derived with no assumptions about the internal structure of the Hamiltonian. For example, second-order MCSCF properties are given by [Eq. (68)]... 

One of the advantages of this formulation is that the fjfim factors have an analytical form for regular cavities, hence it is relatively easy to get energy derivatives at fixed cavity (Rinaldi et al., 1992). The Nancy s SCRF procedure (available as QCPE program, Rinaldi and Pappalardo, 1992) has been recently supplemented by many additional options, i.e. calculation of Gdis and Gcav, calculation of Gei at the MP2, GVB and MCSCF levels of the quantum theory (Chipot et al., 1992 Rinaldi and Rivail, 1995). 

Calculations of analytic excited state properties for correlated methods have been reported by several groups [107-118]. Excited state dynamic properties from cubic response theory were first obtained by Norman et al. at the SCF level [55] and by Jonsson et al. at the MCSCF [56] level, and in a subsequent study a polarizable continuum model was applied to account for solvation effects [119]. Hattlg et al. presented a general theory for excited state response functions at the CC level using a quasi-energy formulation [120] which was subsequently implemented and applied at the CCSD level [121, 122]. The first ID DFT calculation of dynamic excited state polarizabilities, which we will shortly review here, was presented in [58] for pyrimidine and -tetrazine utilizing the double residue of the cubic response function derived in Section 2.7.3. 

Analytical second derivatives for closed-shell (or unrestricted Hartree-Fock (UHF)) SCF wavefunctions are used routinely now. The extension to the MCSCF case is relatively new, however. In contrast to the first derivatives, the coupled perturbed SCF equations have to be solved in order to calculate the second and third energy derivatives. The closed-shell case is relatively straightforward, and will be discussed. The multiconfigurational formalism is... 

Hence, if we are to make bE = 0, we can avoid these terms. But to do so we have to have E optimum with respect to the location of the atomic basis functions, t (R) the MO coefficients, c(R) and the Cl coefficients, C(R). The first cannot be satisfied unless the atomic orbital basis set is floated off the atomic centers to an optimum location [105], while the second requires optimum MO coefficients, and the third optimum Cl coefficients. In practice, we will introduce atomic orbital derivatives explicitly, so the AOs can follow their atoms. Now focusing only on the MO and Cl coefficients, in SCF we have optimum MOs and no Cl term. In MCSCF, both terms would vanish, whUe in Cl, the MO derivatives would remain, but the Cl coefficients contribution would vanish. In the non-variational coupled-cluster theory, neither will vanish and this means that CC theory forces us into some new considerations for analytical forces. 

The derivative coupling can be calculated for Cl or for MCSCF wave-functions using analytic gradients from the expression ... 

---