HartreeFock equations BornOppenheimer Approximation TN0 VNNct Core integral

Hartree-Fock equations Born-Oppenheimer Approximation TN=0, VNN=ct

- Core integral (one-electron integral - Coulombian integral (two-electron integral) - Exchange integral (two-electron integral) Looking for the spin-orbitals which give the best EHF under the orthonormalization conditions and using the method of the Lagrange multipliers one obtains (see Jensen pg. 62 -63)

Hartree-Fock equations Lagrangean multipliers are eigenvalues of the Fock operator: - one-electron operator -Hartree-Fock potential -The average repulsive potential experienced by the i-th electron due to the remaining N-1 electrons -Replaces the complicated 1/r 12 two-electron operator (electron-electron repulsion is taken into account only in an average way)

With: - represents the average local potential at position x 1 arising from an electron located in j Operating with Jj(x 1) we obtain: depends only on the value of i(x 1) Jj is a local operator The second term of v. HF(x 1) is the exchange contribution to the Hartree. Fock potential: Kj(x 1) exchanges the variables of the two spinorbitals depends on the value of i on all points in space because it depends on the integrating variable Kj(x 1) is a non-local operator Jij Kij 0 Jii=Kii no self interaction in Hartree-Fock approximation

Since fi depends on i HF equations must be solved iteratively (SCF procedure) No Initial guess v. HF, f { i} Conv. ? Yes Molecular orbital energies: Total energy: SD, EHF

For a closed-shell system in RHF formalism, the total energy and molecular orbital energies are given by (see Szabo and Ostlund, pag. 83): Ø Each occupied spin-orbital i contributes a term Hi to the energy Ø Each unique pair of electrons (irrespective of their spins) in spatial orbitals i and j contributes the term Jij to the energy Ø Each unique pair of electrons with parallel spins in spatial orbitals ψ i and ψj contributes the term –Kij to the energy Examples: 3 2 1 a) a) b) c) d) b) c) E=2 H 1+J 12 E=H 1+H 2+J 12 -K 12 E=H 1+H 2+J 12 E=H 1+2 H 2+H 3+2 J 12+J 22+J 13+2 J 23 -K 12 -K 13 -K 23 d)

Koopmans’ Theorem What is the physical meaning of orbital energies? Total electronic energy: Summing all εi: why? because εi includes the Coulombian and exchange interactions with all the other electrons (also with εj). Similarly, εj will include also the interactions with εi, so that the interaction between the i-th and j-th electrons will be counted twice in the sum over orbital energies So, what εi is? Consider the ionization of the molecule (one electron removed from orbital number k) and suppose that no change of MO’s occurs during the ionization process

Koopmans’ theorem vertical ionization potential (IP) photoelectron spectroscopy (obtained without optimizing the cation geometry) adiabatic ionization potential is obtained when the geometry of the cation is optimized physical meaning of εi: the i-th orbital ionization energy!

Hartree-Fock-Roothaan Equations LCAO-MO i=1, 2, . . . , K { μ} – a set of known functions The more complete set { μ}, the more accurate representation of the exact MO, the more exact the eigenfunctions of the Fock operator The problem of calculating HF MO the problem of calculating the set cμi LCAO coefficients matrix equation for the cμi coefficients Multiplying by μ*(r 1) on the left and integrating we get: - Fock matrix (Kx. K Hermitian matrix) - overlap matrix (Kx. K Hermitian matrix) - Roothaan equations

More compactly: FC=SC where -the matrix of the expansion coefficients (its columns describes the molecular orbitals) The requirement that the molecular orbitals be orthonormal in the LCAO approximation demands that: The problem of finding the molecular orbitals { i} and orbital energies i involves solving the matrix equation FC=SC. For this, we need an explicit expression for the Fock matrix

Charge density For a closed shell molecule, described by a single determinant wave-function The integral of this charge density is just the total number of electrons: Inserting the molecular orbital expansion into the expression for the charge density we get:

Where: - elements of the density matrix The integral of (r) is By means of the last equation, the electronic charge distribution may be decomposed into contributions associated with the various basis functions of the LCAO expansion. Off-diagonal elements -the electronic population of the atomic overlap distribution -give an indication of contributions to chemical binding when and centered on different atoms Diagonal elements - the net electronic charges residing in orbital

Population analysis ≡ allocate the electrons in the molecule in a fractional manner, among the various parts of the molecule (atoms, bonds, basis functions) - Mulliken population analysis Substituting the basis set expansion we get: Basis set functions are normalized Sμμ=1 Pμμ - number of electrons associated with a particular BF - net population of φμ Qμ = 2 Pμ Sμ (μ≠ ) overlap population associated with two basis functions which may be on the same atom or on two different atoms Total electronic charge in a molecule consists of two parts: first term is associated with individual BF second term is associated with pairs of BF

- gross population for φμ the net charge associated with the atom A total overlap population between atoms A and B

Formaldehyde Mulliken population analysis #P RHF/STO-3 G scf(conventional) Iop(3/33=6) Extralinks=L 316 Noraff Symm=Noint Iop(3/33=1) pop(full) Basis functions:

The summation is over occupied molecular orbitals Example

= sum over the line (or column) corresponding to the C(1 s) basis function = sum over the line (or column) corresponding to the O(2 px) basis function

Atomic populations (AP) Total atomic charges (Q=Z-AP) 1 O 8. 186789 1 O -0. 186789 2 C 5. 926642 2 C 0. 073358 3 H 4 H 0. 943285 3 H 4 H 0. 056715

Molecular orbitals of formaldehyde (RHF/STO-3 G)

Basis Sets with: { μ} – a set of known functions for UHF wave-functions two sets of coefficients are needed: if μ AO LCAO-MO if μ AO LCBF-MO

Basis functions • mathematical functions designed to give the maximum flexibility to the molecular orbitals • must have physical significance • their coefficients are obtained variationally

Slater Type Orbitals (STO) - similar to atomic orbitals of the hydrogen atom - more convenient (from the numerical calculation point of view) than AO, especially when n-l≥ 2 (radial part is simply r 2, r 3, . . . and not a polinom) STO – are labeled like hydrogen atomic orbitals and their normalized form is: STO provide reasonable representations of atomic orbitals however they are not well suited to numerical (fast) calculations of especially two-electron integrals their use in practical molecular orbital calculations has been limited

STO Advantages: • Physically, the exponential dependence on distance from the nucleus is very close to the exact hydrogenic orbitals. • Ensures fairly rapid convergence with increasing number of functions. Disadvantages: • Three and four center integrals cannot be performed analytically. • No radial nodes. These can be introduced by making linear combinations of STOs. Practical Use: • Calculations of very high accuracy, atomic and diatomic systems. • Semi-empirical methods where 3 - and 4 -center integrals are neglected.

Gaussian Type Orbitals (GTO) -introduced by Boys (1950) -powers of x, y, z multiplied by -α is a constant (called exponent) that determines the size (radial extent) of the function or: N - normalization constant f - scaling factor scale all exponents in the related gaussians in molecular calculations l, m, n are not quantum numbers L=l+m+n - used analogously to the angular momentum quantum number for atoms to mark functions as s-type (L=0), p-type (L=1), d-type (L=2), etc (shells)

The absence of rn-1 pre-exponential factor restricts single gaussian primitives to approximate only 1 s, 2 p, 3 d, 4 f, . . . orbitals. However, combinations of gaussians are able to approximate correct nodal properties of atomic orbitals GTO – uncontracted gaussian function (gaussian primitive) - contracted gaussian function (gaussian contraction) STO= GTOs are inferior to STOs in three ways: GTO’s behavior near the nucleus is poorly represented. At the nucleus, the GTO has zero slope; the STO has a cusp. GTOs diminish too rapidly with distance. The ‘tail’ behavior is poorly represented. Extra d-, f-, g-, etc. functions (from Cart. rep. )may lead to linear dependence of the basis set. They are usually dropped when large basis sets are used. Advantage: GTOs have analytical solutions. Use a linear combination of GTOs to overcome these deficiencies.

The first ten normalized gaussian primitives are:

There are 6 possible d-type cartesian gaussians while there are only 5 linearly independent and orthogonal d orbitals The gs, gx, gy and gz primitives have the angular symmetries of the four corresponding AO. The 6 d-type gaussian primitives may be combined to obtain a set of 5 d-type functions: gxy dxy gxz dxz gyz dyz The 6 -th linear combination gives an s-type function: In a similar manner, the 10 f-type gaussian primitives may be combined to obtain a set of 7 f-type functions

GTOs are less satisfactory than STOs in describing the AOs close to the nucleus. The two type functions substantially differ for r=0 and also, for very large values of r. cusp condition: for STO: [d/dr e-ξr]r ≠ 0 for GTO: With GTO the two-electron integrals are more easily evaluated. The reason is that the product of two gaussians, each on different centers, is another gaussian centered between the two centers: where: KAB=(2αβ/[(α+β)π])3/4 exp(-αβ/(α+β)|RA-RB|2] The exponent of the new gaussian centered at Rp is: p=α+β and the third center P is on line joining the centers A and B (see the Figure below) RP=(αRA+βRB)/(α+β)

The product of two 1 s gaussian is a third 1 s gaussian allow a more rapidly and efficiently calculation of the two-electron integrals GTO have different functional behavior with respect to known functional behavior of AOs. contractions (CGF or CGTO) L – the length of the contraction dpμ – contraction coefficients

How the gaussian primitives are derived? by fitting the CGF to an STO using a least square method varying the exponents in quantum calculations on atoms in order to minimize the energy Example STO-3 G basis set for H 2 molecule Each BF is approximated by a STO, which in turn, is fitted to a CGF of 3 primitives hydrogen 1 s orbital in STO-3 G basis set For molecular calculations, first we need a BF to describe the H 1 s atomic orbital then: MO(H 2) = LCBF 3 gaussian primitives: exponent 0. 222766 0. 405771 0. 109818 If we use a scaling factor: coefficient 0. 154329 0. 535328 0. 444636

βi=αif 2 ! Using normalized primitives we do not need a normalization factor for the whole contraction If the primitives are not normalized, we have to obtain a normalization factor. For this we use the condition: S=F 2[I 1+I 2+I 3+2 I 4+2 I 5+2 I 6]

But: so that: Analogously: and thus:

Now, Imposing that S=1 we obtain: In the general case of a contraction of dimension n, the above expression become:

Summary The 1 s hydrogen orbital in STO-3 G basis set will be: with: - normalization factors for primitives - normalization factor for the whole contraction (when unnormalized primitives or segmented contractions are used) αi βi ci Ni 2. 227660 3. 425250 0. 154329 1. 794441 0. 276934 0. 405771 0. 623913 0. 535328 0. 500326 0. 267839 0. 109818 0. 168856 0. 444635 0. 187736 0. 083474 N=1. 0000002

Explicitly: If the exponents are not scaled:

STO-3 G basis set example http: //www. chem. utas. edu. au/staff/yatesb/honours/modules/mod 5/c_sto 3 g. html This is an example of the STO-3 G basis set for methane in the format produced by the "gfinput" command in the Gaussian computer program. The first atom is carbon. The other four are hydrogens. Standard basis: STO-3 G (5 D, 7 F) Basis set in the form of general basis input: 1 0 S 3 1. 00 . 7161683735 D+02 . 1543289673 D+00 . 1304509632 D+02 . 5353281423 D+00 . 3530512160 D+01 . 4446345422 D+00 SP 3 1. 00 . 2941249355 D+01 -. 9996722919 D-01 . 1559162750 D+00 . 6834830964 D+00 . 3995128261 D+00 . 6076837186 D+00 . 2222899159 D+00 . 7001154689 D+00 . 3919573931 D+00 **** 2 0 S 3 1. 00 . 3425250914 D+01 . 1543289673 D+00 . 6239137298 D+00 . 5353281423 D+00 . 1688554040 D+00 . 4446345422 D+00 **** 3 0 S 3 1. 00 . 3425250914 D+01 . 1543289673 D+00 . 6239137298 D+00 . 5353281423 D+00 . 1688554040 D+00 . 4446345422 D+00 **** 4 0 S 3 1. 00 . 3425250914 D+01 . 1543289673 D+00 . 6239137298 D+00 . 5353281423 D+00 . 1688554040 D+00 . 4446345422 D+00 **** 5 0 S 3 1. 00 . 3425250914 D+01 . 1543289673 D+00 . 6239137298 D+00 . 5353281423 D+00 . 1688554040 D+00 . 4446345422 D+00 ****

Split valence basis sets http: //www. chem. utas. edu. au/staff/yatesb/honours/modules/mod 5/split_bas. html Valence orbitals are represented by more than one basis function, (each of which can in turn be composed of a fixed linear combination of primitive Gaussian functions). Depending on the number of basis functions used for the reprezentation of valence orbitals, the basis sets are called valence double, triple, or quadruple-zeta basis sets. Since the different orbitals of the split have different spatial extents, the combination allows the electron density to adjust its spatial extent appropriate to the particular molecular environment. Split is often made for valence orbitals only, which are chemically important. 3 -21 G basis set The valence functions are split into one basis function with two GTOs, and one with only one GTO. (This is the "two one" part of the nomenclature. ) The core consists of three primitive GTOs contracted into one basis function, as in the STO-3 G basis set. 1 0 //C atom S 3 1. 00 . 1722560000 D+03 . 6176690000 D-01 . 2591090000 D+02 . 3587940000 D+00 . 5533350000 D+01 . 7007130000 D+00 SP 2 1. 00 . 3664980000 D+01 -. 3958970000 D+00 . 2364600000 D+00 . 7705450000 D+00 . 1215840000 D+01 . 8606190000 D+00 SP 1 1. 00 . 1958570000 D+00 . 1000000000 D+01 **** 2 0 //H atom S 2 1. 00 . 5447178000 D+01 . 1562850000 D+00 . 8245472400 D+00 . 9046910000 D+00 S 1 1. 00 . 1831915800 D+00 . 100000 D+01 **** 6 -311 G basis set

Extended basis sets The most important additions to basis sets are polarization functions and diffuse basis functions. Polarization basis functions The influence of the neighboring nuclei will distort or polarize the electron density near a given nucleus. In order to take this effect into account, orbitals that have more flexible shapes in a molecule than the s, p, d, etc shapes in the free atoms are used. An s orbital is polarized by using a p-type orbital A p orbital is polarized by mixing in a d-type orbital 6 -31 G(d) a set of d orbitals is used as polarization functions on heavy atoms 6 -31 G(d, p) a set of d orbitals are used as polarization functions on heavy atoms and a set of porbitals are used as polarization functions on hydrogen atoms

Diffuse basis functions For excited states and anions where the electronic density is more spread out over the molecule, some basis functions which themselves are more spread out are needed (i. e. GTOs with small exponents). These additional basis functions are called diffuse functions. They are normally added as single GTOs. 6 -31+G - adds a set of diffuse sp orbitals to the atoms in the first and second rows (Li - Cl). 6 -31++G - adds a set of diffuse sp orbitals to the atoms in the first and second rows (Li- Cl) and a set of diffuse s functions to hydrogen. Diffuse functions can also be added along with polarization functions. This leads, for example, to the 6 -31+G(d), 6 -31+G(d, p) and 6 -31++G(d, p) basis sets. Standard basis: 6 -31+G (6 D, 7 F) Basis set in the form of general basis input: 1 0 S 6 1. 00 . 3047524880 D+04 . 1834737130 D-02 . 4573695180 D+03 . 1403732280 D-01 . 1039486850 D+03 . 6884262220 D-01 . 2921015530 D+02 . 2321844430 D+00 . 9286662960 D+01 . 4679413480 D+00 . 3163926960 D+01 . 3623119850 D+00 SP 3 1. 00 . 7868272350 D+01 -. 1193324200 D+00 . 6899906660 D-01 . 1881288540 D+01 -. 1608541520 D+00 . 3164239610 D+00 . 5442492580 D+00 . 1143456440 D+01 . 7443082910 D+00 SP 1 1. 00 . 1687144782 D+00 . 1000000000 D+01 SP 1 1. 00 . 4380000000 D-01 . 1000000000 D+01 **** 2 0 S 3 1. 00 . 1873113696 D+02 . 3349460434 D-01 . 2825394365 D+01 . 2347269535 D+00 . 6401216923 D+00 . 8137573262 D+00 S 1 1. 00 . 1612777588 D+00 . 100000 D+01 ****

Number of primitives and basis functions for 1, 2 -Benzosemiquinone free radical with the STO-3 G basis set Primitives: atom C: atom H: atom O: nr. primitives = 15 x nr. atoms = 6 → 90 nr. primitives = 3 x nr. atoms = 4 → 12 nr. primitives = 15 x nr. atoms = 2 → 30 TOTAL: 132 GTO primitives Basis functions: atom C: nr. BF = 5 x nr. atoms = 6 → 30 atom H: nr. BF = 1 x nr. atoms = 4 → 4 atom O: nr. BF = 5 x nr. atoms = 2 → 10 TOTAL: 44 BF

Number of primitives and basis functions for 1, 2 -Benzosemiquinone free radical with the 6 -31+G* (6 -31+G(d)) basis set Primitives: atom C: atom H: atom O: nr. primitives = 32 x nr. atoms = 6 → 192 nr. primitives = 4 x nr. atoms = 4 → 16 nr. primitives = 32 x nr. atoms = 2 → 64 TOTAL: 272 GTO primitives Basis functions: atom C: nr. BF = 19 x nr. atoms = 6 → 114 atom H: nr. BF = 2 x nr. atoms = 4 → 8 atom O: nr. BF = 19 x nr. atoms = 2 → 38 TOTAL: 160 BF

Effective core potentials (ECPs) • Core electrons, which are not chemically very important, require a large number of basis functions for an accurate description of their orbitals. • An effective core potential is a linear combination of specially designed Gaussian functions that model the core electrons, i. e. , the core electrons are represented by a effective potential and one treats only the valence electrons explicitly. • ECP reduces the size of the basis set needed to represent the atom (but introduces additional approximations) • for heavy elements, pseudopotentials can also include of relativistic effects that otherwise would be costly to treat ECP potentials are specified as parameters of the following equation: where p is the dimension of the expansion di are the coefficients for the expansion terms, r 0 is the distance from nucleus and ξi represents the exponents for each term. • Saving computational effort • Taking care of relativistic effects partly • Important for heavy atoms, e. g. , transition metal atoms

ECP example complex. Pd 1. chk #P Opt B 3 LYP/gen pseudo=read etc. H H complex Pd v 1 02 C C H C H C etc. 8. 89318310 9. 52931379 9. 29586123 10. 52592748 10. 95942133 10. 85850598 11. 51852449 10. 20972534 9. 90388210 8. 77525770 7. 93893890 8. 89096200 8. 13380930 10. 13123090 10. 22866610 11. 23549650 6. 72569337 6. 27102032 6. 60431879 5. 30965653 4. 98695425 4. 84438728 4. 19609286 5. 34144511 4. 15752044 17. 83312399 10. 48668123 5. 63848578 17. 14049639 11. 10318367 NCOH 0 6 -31 G(d) **** Pd 0 CEP-121 G Recomendations for basis set selection • Always a compromise between accuracy and computational cost! • With the increase of basis set size, calculated energy will converge (complete basis set (CBS) limit). • Special cases (anion, transition metal, transtion state) • Use smaller basis sets for preliminary calculations and for heavy duties (e. g. , geometry optimizations), and use larger basis sets to refine calculations. • Use larger basis sets for critical atoms (e. g. , atoms directly involved in bond-breaking/forming), and use smaller basis sets for unimportant atoms (e. g. , atoms distant away from active site). • Use popular and recommended basis sets. They have been tested a lot and shown to be good for certain types of calculations.