Post Hartree-Fock Methods (Lecture 2)

1. Post Hartree-Fock Methods(Lecture 2) NSF Computational Nanotechnology and Molecular Engineering Pan-American Advanced Studies Institutes (PASI) Workshop January 5-16, 2004 California Institute of Technology, Pasadena, CA Andrew S. Ichimura

2. Outline • Shortcomings of the SCF-RHF procedure • Configuration Interaction • MCSCF • Size-consistency and size-extensivity • Perturbation theory • Coupled Cluster Methods

3. What is electron correlation and why do we need it? F0 is a single determinantal wavefunction. Slater Determinant Recall that the SCF procedure accounts for electron-electron repulsion by optimizing the one-electron MOs in the presence of an average field of the other electrons. The result is that electrons in the same spatial MO are too close together; their motion is actually correlated (as one moves, the other responds). Eel.cor. = Eexact - EHF (B.O. approx; non-relativistic H)

4. RHF dissociation problem Consider H2 in a minimal basis composed of one atomic 1s orbital on each atom. Two AOs (c) leads to two MOs (f)…

5. The ground state wavefunction is: Slater determinant with two electrons in the bonding MO Expand the Slater Determinant Factor the spatial and spin parts H does not depend on spin Four terms in the AO basis Ionic terms, two electrons in one Atomic Orbital Covalent terms, two electrons shared between two AOs

6. H2 Potential Energy Surface H. + H. E At the dissociation limit, H2 must separate into two neutral atoms. 0 H H Bond stretching H H At the RHF level, the wavefunction, F, is 50% ionic and 50% covalent at all bond lengths. H2 does not dissociate correctly at the RHF level!! Should be 100% covalent at large internuclear separations.

7. RHF dissociation problem has several consequences: • Energies for stretched bonds are too large. Affects transition state structures - Ea are overestimated. • Equilibrium bond lengths are too short at the RHF level. (Potential well is too steep.) HF method ‘overbinds’ the molecule. • Curvature of the PES near equilibrium is too great, vibrational frequencies are too high. • The wavefunction contains too much ‘ionic’ character; dipole moments (and also atomic charges) at the RHF level are too large. On the bright side, SCF procedures recover ~99% of the total electronic energy. But, even for small molecules such as H2, the remaining fraction of the energy - the correlation energy - is ~110 kJ/mol, on the order of a chemical bond.

8. To overcome the RHF dissociation problem,Use a trial function that is a combination of F0 and F1 First, write a new wavefunction using the anti-bonding MO. The form is similar to F0, but describes an excited state: MO basis AO basis Ionic terms Covalent terms

9. Trial function- Linear combination of F0 and F1; two electron configurations. Ionic terms Covalent terms Three points: As the bond is displaced from equilibrium, the coefficients (a0, a1) vary until at large separations, a1 = -a0: Ionic terms disappear and the molecule dissociates correctly into two neutral atoms. Y = YCI, an example of configuration interaction. The inclusion of anti-bonding character in the wavefunction allows the electrons to be farther apart on average. Electronic motion is correlated. The electronic energy will be lower (two variational parameters).

10. Configuration Interaction - Excited Slater Determinants Since the HF method yields the best single determinant wavefunction and provides about 99% of the total electronic energy, it is commonly used as the reference on which subsequent improvements are based. As a starting point, consider as a trial function a linear combination of Slater determinants: Multi-determinant wavefunction a0 is usually close to 1 (~0.9). • M basis functions yield M molecular orbitals. • For N electrons, N/2 orbitals are occupied in the RHF wavefunction. • M-N/2 are unoccupied or virtual (anti-bonding) orbitals.

11. Generate excited Slater determinants by promoting up to N electrons from the N/2 occupied to M-N/2 virtuals: b b b 9 a,b,c… = virtual MOs a a a a 8 a,b 7 c,d c 6 k,l k 5 i i,j i i i 4 3 i,j,k… = occupied MOs j j j 2 1 … Ref. Single Double Triple Quadruple Excitation level

12. Represent the space containing all N-fold excitations byY(N). Then the COMPLETE CI wavefunction has the form Where Linear combination of Slater determinants with single excitations Doubly excitations Triples N-fold excitation The complete YCI expanded in an infinite basis yields the exact solution to the Schrödinger eqn. (Non-relativistic, Born-Oppenheimer approx.)

13. The various coefficients, , may be obtained in a variety of ways. A straightforward method is to use the Variation Principle. Expectation value of He. Energy is minimized wrt coeff In a fashion analogous to the HF eqns, the CI Schrodinger equation can be formulated as a matrix eigenvalue problem. The elements of the vector, , are the coefficients, And the eigenvalue, EK, approximates the energy of the Kth state. E1 = ECI for the lowest state of a given symmetry and spin. E2 = 1st excited state of the same symmetry and spin, and so on.

14. Some nomenclature… One-electron basis (one-particle basis) refers to the basis set. This limits the description of the one-electron functions, the Molecular Orbitals. The size of the many-electron basis (N-particle basis) refers to the number of Slater determinants. This limits the description of electron correlation. • In practice, • Complete CI (Full CI) is rarely done even for finite basis sets - too expensive. Computation scales factorialy with the number of basis functions (M!). • Full CI within a given one-particle basis is the ‘benchmark’ for that basis since 100% of the correlation energy is recovered. Used to calibrate approximate correlation methods. • CI expansion is truncated at a some excitation level, usually Singles and Doubles (CISD).

15. Configuration State Functions Consider a single excitation from the RHF reference. Both FRHF and F(1) have Sz=0, but F(1) is not an eigenfunction of S2. FRHF F(1) Linear combination of singly excited determinants is an eigenfunction of S2. Configuration State Function, CSF (Spin Adapted Configuration, SAC) Only CSFs that have the same multiplicity as the HF reference contribute to the correlation energy. Singlet CSF

16. Example H2O: (19 basis functions) CISD (~80-90%) Full CI

17. Example: Neon Atom Relative importance Ref. Singles 2 Doubles 1 Triples 4 Quadruples 3 Weight = for a given excitation level. (Frozen core approx., 5s4p3d basis - 32 functions) CISD (singles and doubles) is the only generally applicable method. For modest sized molecules and basis sets, ~80-90% of the correlation energy is recovered. CISD recovers less and less correlation energy as the size of the molecule increases.

18. Size Consistent and Size Extensive Size consistent method - the energy of two molecules (or fragments) computed at large separation (100 Å) is equal to the twice energy of the individual molecule (fragment). Only defined if the molecules are non-interacting. Ex. (ECISD of two H2 separated by 100Å) < 2(ECISD of one H2) Size extensive method - the energy scales properly with the number of particles. (Same fraction of correlation energy is recovered for CH4, C2H6, C3H8, etc.) Full CI is size consistent and extensive. All forms of truncated CI are not. (Some forms of CI, esp. MR-CI are approximately size consistent and size extensive with a large enough reference space.)

19. Multi-configuration Self-consistent Field (MCSCF) 9 Carry out Full CIandorbital optimization within a small active space. Six-electron in six-orbital MCSCF is shown. Written as [6,6]CASSCF. Complete Active Space Self-consistent Field (CASSCF) 8 7 6 H2O MOs Why? To have a better description of the ground or excited state. Some molecules are not well-described by a single Slater determinant, e.g. O3. To describe bond breaking/formation; Transition States. Open-shell system, especially low-spin. Low lying energy level(s); mixing with the ground state produces a better description of the electronic state. … 5 4 3 2 1

20. MCSCF Features: In general, the goal is to provide a better description of the main features of the electronic structure before attempting to recover most of the correlation energy. Some correlation energy (static correlation energy) is recovered. (So called dynamic correlation energy is obtained through CI and other methods through a large N-particle basis.) The choice of active space - occupied and virtual orbitals - is not always obvious. (Chemical intuition and experience help.) Convergence may be poor. CASSCF wavefunctions serve as excellent reference state(s) to recover a larger fraction of the dynamical correlation energy. A CISD calculation from a [n,m]-CASSCF reference is termed Multi-Reference CISD (MR-CISD). With a suitable active space, MRCISD approaches Full CI in accuracy for a given basis even though it is not size-extensive or consistent.

21. Examples of compounds that require MCSCF for a qualitatively correct description. zwitterionic Singlet state of twisted ethene, biradical. biradical Transition State

22. Mœller-Plesset Perturbation Theory In perturbation theory, the solution to one problem is expressed in terms of another one solved previously. The perturbation should be small in some sense relative to the known problem. Hamiltonian with pert., l Unperturbed Hamiltonian As the perturbation is turned on, W (the energy) and Y change. Use a Taylor series expansion in l.

23. Unperturbed H is the sum over Fock operators  Moller-Plesset (MP) pert th. Perturbation is a two-electron operator when H0 is the Fock operator. With the choice of H0, the first contribution to the correlation energy comes from double excitations. Explicit formula for 2nd order Moller-Plesset perturbation theory, MP2.

24. Advantages of MP’n’ Pert. Th. • MP2 computations on moderate sized systems (~150 basis functions) require the same effort as HF. Scales as M5, but in practice much less. • Size-extensive (but not variational). Size-extensivity is important; there is no error bound for energy differences. In other words, the error remains relatively constant for different systems. • Recovers ~80-90% of the correlation energy. • Can be extended to 4th order: MP4(SDQ) and MP4(SDTQ). MP4(SDTQ) recovers ~95-98% of the correlation energy, but scales as M7. • Because the computational effort is significanly less than CISD and the size-extensivity, MP2 is a good method for including electron correlation.

25. Coupled Cluster Theory Perturbation methods add all types of corrections, e.g., S,D,T,Q,..to a given order (2nd, 3rd, 4th,…). Coupled cluster (CC) methods include all corrections of a given type to infinite order. The CC wavefunction takes on a different form: Coupled Cluster Wavefunction F0 is the HF solution Exponential operator generates excited Slater determinants Cluster Operator N is the number of electrons

26. CC Theory cont. The T-operator acting on the HF reference generates all ith excited Slater Determinants, e.g. doubles Fijab. Expansion coefficients are called amplitudes; equivalent to the ai’s in the general multi-determinant wavefunction. HF ref. singles doubles triples Quadruple excitations The way that Slater determinants are generated is rather different…

27. CC Theory cont. HF reference Singly excited states Connected doubles Dis-connected doubles Connected triples, ‘true’ triples ‘Product’ Triples, disconnected triples True quadruples - four electrons interacting Product quadruples - two noninteracting pairs Product quadruples, and so on.

28. CC Theory cont. If all cluster operators up to TN are included, the method yields energies that are essentially equivalent to Full CI. In practice, only the singles and doubles excitation operators are used forming the Coupled Cluster Singles and Doubles model (CCSD). The result is that triple and quadruple excitations also enter into the energy expression (not shown) via products of single and double amplitudes. It has been shown that the connected triples term, T3, is important. It can be included perturbatively at a modest cost to yield the CCSD(T) model. With the inclusion of connected triples, the CCSD(T) model yields energies close to the Full CI in the given basis, a very accurate wavefunction.

29. Comparison of Models Accuracy with a medium sized basis set (single determinant reference): HF << MP2 < CISD < MP4(SDQ) ~CCSD < MP4(SDTQ) < CCSD(T) In cases where there is (a) strong multi-reference character and (b) for excited states, MR-CI methods may be the best option.