Ab initio

2. Ab initio Reactant – Transition State Structure – Product • Selection of the theoretical model • Geometry optimization • Frequency calculation • Energy calculation • Refining the theoretical model

3. Description of the theoretical model QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) Energy calculation Geometry optimisation MP2/6-311G(d,p) method basis set

4. Basis Set

5. Introduction The most important result is the ENERGY !!! Schrödinger equation : HΨ=EΨ Aim to adequate molecular energy Problem:- energy is notavailable • Approximations: • -computer capacity • -CPU time • -size of molecule

6. Model chemistry: theoretical method and basis set Goal: select the most accurate calculation that is computationally feasible for a given molecular system

8. Model Chemistries - three areas of consideration • Basis sets • Theoretical methods QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) Energy calculation Geometry optimisation MP2/6-311G(d,p) method basis set

9. Five statements for demythologization: Nº 1.: The term “orbital” is a synonym for the term “One-Electron” Function (OEF) Nº 2.: A single centered OEF is synonymous with “Atomic Orbital”. A multi centered OEF is synonymous with “Molecular Orbital”. Orbital == OEF

10. Five statements …. • Nº 3.: 3 ways to express a mathematical function: • Explicitly in analytical form • (hydrogen-like AOs) • As a table of numbers • (Hartree-Fock type AOs for numerous atoms) • In the form of an expansion • (expression of an MO in terms of a set of AO)

11. Five statements …. Nº 4.: The generation of MOs (f-s) from AOs (h-s) is equivalent to the transformation of an N-dimensional vector space where {h}is the original set of non-orthogonal functions. After orthogonalization of the non-orthogonal AO basis set {h} the orthogonal set {c} is rotated to the another orthogonal set{f}. O

12. Five statements …. Nº 5.: There are certain differences between the shape of numerical Hartree-Fock atomic orbitals (HF-AO), the analytic Slater type orbitals (STO) and the analytic Gaussian type functions (GTF). However , these differences are irrelevant to the final results as the MO can be expanded in terms of any of these complete sets of functions to any desired degree of accuracy. Atomic orbital basis sets

13. Basis set • Basis functions approximate orbitals of atoms in molecule • Linear combination of basis functions approximates total electronic wavefunction • Basis functions are linear combinations of gaussian functions • Contracted gaussians • Primitive gaussians

14. STOs v. GTOs • Slater-type orbitals (J.C. Slater) • Represent electron density well in valence region and beyond (not so well near nucleus) • Evaluating these integrals is difficult • Gaussian-type orbitals (F. Boys) • Easier to evaluate integrals, but don’t represent electron density well • Overcome this by using linear combination of GTOs (Sμ)=

17. Minimal basis set • One basis function for every atomic orbital required to describe the free atom • Most-common: STO-3G • Linear combination of 3 Gaussian-type orbitals fitted to one Slater-type orbital • CH4: H(1s); C(1s,2s,2px,2py,2pz)

24. More basis functions per atom • Split valence basis sets • Double-zeta: 2 “sizes” of basis functions for each valence atomic orbital • 3-21G CH4: H(1s,1s'), C(1s,2s,2s',2px,2py,2pz,2px',2py',2pz') • Triple-zeta: 3 “sizes” of basis functions for each valence atomic orbital • 6-311G CH4: H(1s,1s',1s''), C(1s,2s,2s',2s'',2px,2py,2pz, 2px',2py',2pz',2px'',2py'',2pz'')

25. More basis functions per atom • Split valence basis sets • Double-zeta: • Triple-zeta:

26. 22 Total 36

27. Total 42 22

29. Ways to increase a basis set • Add more basis functions per atom • allow orbitals to “change size” • Add polarization functions • allow orbitals to “change shape” • Add diffuse functions for electrons with large radial extent • Add high angular momentum functions

30. Add polarization functions • Allow orbitals to change shape • Add p orbitals to H • Add d orbitals to 2nd row atoms • Add f orbitals to transition metals • 6-31G(d) - d functions per heavy atoms • Also denoted: 6-31G* • 6-31G(d,p) - d functions per heavy atoms and p functions to H atoms • Also deonoted: 6-31G**

36. Add diffuse functions • “Large” s and p orbitals for “diffuse electrons” • Lone pairs, anions, excited states, etc. • 6-31+G - diffuse functions per heavy atom • 6-31++G - diffuse functions both per heavy atom and per H atom

38. High angular momentum functions • “Custom-made” basis sets • 6-31G(2d) - 2d functions per heavy atom • 6-311++G(3df,3pd) • Triple-zeta valence • Diffuse functions on heavy atoms, H atoms • 3d, 1f functions per heavy atom; 3p, 1d functions per H atom

49. Minimal basis sets A common naming convention for minimal basis sets is STO-XG, where X is an integer. This X value represents the number of Gaussian primitive functions comprising a single basis function. In these basis sets, the same number of Gaussian primitives comprise core and valence orbitals. Minimal basis sets typically give rough results that are insufficient for research-quality publication, but are much cheaper than their larger counterparts. Here is a list of commonly used minimal basis sets: STO-2G STO-3G STO-6G STO-3G* - Polarized version of STO-3G

50. Split-valence basis sets During most molecular bonding, it is the valence electrons which principally take part in the bonding. In recognition of this fact, it is common to represent valence orbitals by more than one basis function, (each of which can in turn be composed of a fixed linear combination of primitive Gaussian functions). The notation for these split-valence basis sets is typically X-YZg. In this case, X represents the number primitive Gaussians comprising each core atomic orbital basis function. The Y and Z indicate that the valence orbitals are composed of two basis functions each Here is a list of commonly used split-valence basis sets: 3-21g 3-21g* - Polarized 3-21+g - Diffuse functions 3-21+g* - With polarization and diffuse functions 6-31g 6-31g* 6-31+g* 6-31g(3df, 3pd) 6-311g 6-311g* 6-311+g* SV(P) SVP