Semiempirical Molecular Orbital Theory
|
|
- Bennett Walters
- 6 years ago
- Views:
Transcription
1 Semiempirical Molecular Orbital Theory Slater- determinantal wave func:ons and artree- Fock theory Video III.iv
2 Constructing a 1-Electron Wave Function The units of the wave function are such that its square is electron per volume. As electrons are quantum particles with non-point distributions, sometimes we say density or probability density instead of electron per volume (especially when there is more than one electron, since they are indistinguishable as quantum particles) For instance, a valid wave function in cartesian coordinates for one electron might be: "( x, y,z; Z) = 2Z 5 / 2 81 # % ' 6 $ Z x 2 + y 2 + z 2 & ( * ye $Z x 2 +y 2 +z 2 / 3 ) normalization factor radial phase factor cartesian directionality (if any) ensures square integrability
3 Constructing a 1-Electron Wave Function A valid wave function in cartesian coordinates for one electron might be: "( x, y,z; Z) = 2Z 5 / 2 81 # % ' 6 $ Z x 2 + y 2 + z 2 & ( * ye $Z x 2 +y 2 +z 2 / 3 ) normalization factor radial phase factor cartesian directionality (if any) ensures square integrability x 1 x x 2 P y 1 y y 2 = x 2 y 2 z 2 φ 2 dx dy dz x 1 y 1 z 1 z 1 z z 2 This permits us to compute the probability of finding the electron within a particular cartesian volume element (normalization factors are determined by requiring that P = 1 when all limits are infinite, i.e., integration over all space)
4 Constructing a 1-Electron Wave Function To permit additional flexibility, we may take our wave function to be a linear combination of some set of common basis functions, e.g., atomic orbitals (LCAO). Thus ( ) = a i ϕ( r) φ r N i=1 For example, consider the wave function for an electron in a C bond. It could be represented by s and p functions on the atomic positions, or s functions along the bond axis, or any other fashion convenient. C
5 Constructing a 1-Electron Wave Function To optimize the coefficients in our LCAO expansion, we use the variational principle, which says that E = a i ϕ i a j ϕ j dr i j a i ϕ i a j ϕ j dr i j = ij ij a i a j a i a j ϕ i ϕ j dr ϕ i ϕ j dr = ij ij a i a j a i a j ij S ij which is minimized by one root E of the secular equation (the other roots are excited states) each value of E that satisfies the secular equation determines all of the a i values and thus the shape of the molecular orbital wave function (MO) 11 ES ES 12 1N ES 1N 21 ES ES 22 2N ES 2 N = 0 N1 ES N1 N2 ES N2 NN ES NN
6 What Are These Integrals? The electronic amiltonian includes kinetic energy, nuclear attraction, and, if there is more than one electron, electron-electron repulsion ij = ϕ i 1 nuclei electrons 2 Z ϕ j ϕ k φ 2 i ϕ j + ϕ i 2 k r k n r n ϕ j The final term is problematic. Solving for all electrons at once is a many-body problem that has not been solved even for classical particles. An approximation is to ignore the correlated motion of the electrons, and treat each electron as independent, but even then, if each MO depends on all of the other MOs, how can we determine even one of them? The artree-fock approach accomplishes this for a many-electron wave function expressed as an antisymmetrized product of one-electron MOs (a so-called Slater determinant)
7 Choose a basis set Choose a molecular geometry q (0) Compute and store all overlap, one-electron, and two-electron integrals Guess initial density matrix P (0) The artree- Fock procedure Construct and solve artree- Fock secular equation Replace P (n 1) with P (n) Choose new geometry according to optimization algorithm no no Optimize molecular geometry? yes Does the current geometry satisfy the optimization criteria? yes Output data for optimized geometry Construct density matrix from occupied MOs Is new density matrix P (n) sufficiently similar to old density matrix P (n 1)? no yes Output data for unoptimized geometry F 11 ES 11 F 12 ES 12 F 1N ES 1N F 21 ES 21 F 22 ES 22 F 2N ES 2N F N1 ES N1 F N2 ES N2 F NN ES NN nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ ( µν λσ) = φ µ ( 1)φ ν 1 ( ) 1 k ( µν λσ) 1 2 ( µλ νσ) occupied P λσ = 2 a λi φ r λ ( 2)φ σ ( 2)dr ( 1)dr 2 12 i = 0 a σi ( ) One MO per root E
8 Cement Your Understanding F µν = µ ν Z k µ 1 ν r k + P λσ λσ nuclei k ( µν λσ) 1 2 ( µλ νσ) Try to write out, neatly and carefully, in full mathematical notation, every generic integral in the Fock matrix for Slater-type orbital basis functions with all functions and operators expressed in Cartesian coordinates (which is how they are most typically evaluated in actual practice)
9 Semiempirical Molecular Orbital Theory Semiempirical simplifica:ons of the Fock matrix (early) Video III.v
10 Choose a basis set Choose a molecular geometry q (0) Compute and store all overlap, one-electron, and two-electron integrals Guess initial density matrix P (0) The artree- Fock procedure Construct and solve artree- Fock secular equation Replace P (n 1) with P (n) Choose new geometry according to optimization algorithm no no Optimize molecular geometry? yes Does the current geometry satisfy the optimization criteria? yes Output data for optimized geometry Construct density matrix from occupied MOs Is new density matrix P (n) sufficiently similar to old density matrix P (n 1)? no yes Output data for unoptimized geometry F 11 ES 11 F 12 ES 12 F 1N ES 1N F 21 ES 21 F 22 ES 22 F 2N ES 2N F N1 ES N1 F N2 ES N2 F NN ES NN nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ ( µν λσ) = φ µ ( 1)φ ν 1 ( ) 1 k ( µν λσ) 1 2 ( µλ νσ) occupied P λσ = 2 a λi φ r λ ( 2)φ σ ( 2)dr ( 1)dr 2 12 i = 0 a σi ( ) One MO per root E
11 nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ k ( µν λσ) 1 2 ( µλ νσ) occupied P λσ = 2 a λi φ i = a µi ϕ µ µ i a σi F 11 ES 11 F 12 ES 12 F 1N ES 1N F 21 ES 21 F 22 ES 22 F 2N ES 2N = 0 F N1 ES N1 F N2 ES N2 F NN ES NN Consider only first- and second- row atoms; consider a minimal basis set ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2
12 nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ ( µν λσ) = ϕ µ ( 1)ϕ ν 1 k ( µν λσ) 1 2 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 Things to consider in a semiempirical model: ( µλ νσ) occupied P λσ = 2 a λi φ i = a µi ϕ µ µ i a σi 1) What shall I choose for a basis set? 2) What shall I do for overlap integrals (in the secular determinant)? 3) What shall I do for Fock matrix elements? 4) As part of (3), what shall I do about electron repulsion integrals? Essen:ally all semiempirical models adopt a minimal basis set. The AO basis func:ons cover only the valence s, p, d, and f orbitals (p, d, and f only as needed), one func:on per orbital (so, 1 s, 3 p, 5d, and 7f orbitals). The basis func:ons themselves are taken to be Slater orbitals. ( ϕ STO ( r,θ,φ;ζ,n,l,m) = 2ζ ) n +1/ 2 ( 2n)! [ ] 1/ 2 rn 1 e ζr Y l ( ) m θ,φ
13 nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ ( µν λσ) = ϕ µ ( 1)ϕ ν 1 k ( µν λσ) 1 2 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 ( µλ νσ) occupied P λσ = 2 a λi φ i = a µi ϕ µ µ i a σi General choice for overlap integrals, S µν : The overlap between any two STOs can be computed analy:cally (i.e., it is a simple arithme:c func:on of posi:on and values of ζ, n, l, and m). Overlap integrals find use in compu:ng some Fock matrix elements. owever, when forming the secular determinant, the simplifica:on S µν = δ µν is adopted. ϕ STO ( ) n +1/ 2 ( r,θ,φ;ζ,n,l,m) = 2ζ ( 2n)! [ ] 1/ 2 rn 1 e ζr Y l ( ) m θ,φ
14 nuclei F µν = µ ν Z k µ 1 ν r k + P λσ λσ ( µν λσ) = ϕ µ ( 1)ϕ ν 1 k ( µν λσ) 1 2 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 ( µλ νσ) occupied P λσ = 2 a λi φ i = a µi ϕ µ µ i a σi There are 3 types of Fock matrix integrals, F µν : A Fock matrix element can be diagonal (i.e., µ = ν), or not. When the matrix element is off- diagonal, it can be monatomic (i.e., µ ν, but µ and ν are on the same atom), or it can be diatomic (the only remaining possibility) As every Fock matrix element includes ERIs, first let s consider the ERIs (we ll get to one- electron integrals later)
15 COMPLETE NEGLECT OF DIFFERENTIAL OVERLAP Reducing the N 4 scaling problem: ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 CNDO: ( µν λσ) = δ µν δ λσ ( µµ λλ) only ERIs not equal to zero are Coulomb integrals between 2 AO basis func:ons, which may or may not be on the same atom ( µµ λλ) = γ AB γ is a number or func:on that depends only on the atoms A and B on which µ and λ are found, respec:vely. γ AA = IP A EA A interac:on between electrons on same atom taken as difference between ioniza:on poten:al and electron affinity γ AB = γ AA + γ BB ( ) 2 + r AB γ AA + γ BB func:on above has proper limits at r AB = 0 (since then A must be B, the same atom) and at r AB very large (just 1/r AB, which is Coulomb s law in atomic units) M + M IP EA M + M one extra e-e repulsion on atom M
16 COMPLETE NEGLECT OF DIFFERENTIAL OVERLAP Reducing the N 4 scaling problem: ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 CNDO: ( µν λσ) = δ µν δ λσ ( µµ λλ) only ERIs not equal to zero are Coulomb integrals between 2 AO basis func:ons, which may or may not be on the same atom ( µµ λλ) = γ AB γ is a number or func:on that depends only on the atoms A and B on which µ and λ are found, respec:vely. no dis4nc4on in 2- electron repulsions: C! CC C! CC N : :! NN : : N! NN
17 Semiempirical Molecular Orbital Theory Semiempirical simplifica:ons of the Fock matrix (subsequent) Video III.vi
18 COMPLETE NEGLECT OF DIFFERENTIAL OVERLAP Reducing the N 4 scaling problem: ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 CNDO: ( µν λσ) = δ µν δ λσ ( µµ λλ) only ERIs not equal to zero are Coulomb integrals between 2 AO basis func:ons, which may or may not be on the same atom ( µµ λλ) = γ AB γ is a number or func:on that depends only on the atoms A and B on which µ and λ are found, respec:vely. no dis4nc4on in 2- electron repulsions:! CC C! CC C
19 INTERMEDIATE NEGLECT OF DIFFERENTIAL OVERLAP Reducing the N 4 scaling problem: ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 INDO: ( µν λσ) = δ µν δ λσ ( µµ λλ) s:ll taken to be true unless µ and λ are on the same atom ( ss ss) = G ss ( ss pp) = G sp ( pp pp) = G pp ( pp pʹ pʹ ) = G pp ʹ ( sp sp) = L sp For an atom with only s and p valence func:ons, 5 unique integrals are kept (they are simply values, in energy units, that can be determined from spectroscopy or taken as free parameters). Note that p refers to a p orbital different from a reference one, e.g., ( pp p p ) might be (p x p x p y p y ), i.e., the through- space repulsive interac:on between electrons in two perpendicular p orbitals This extra flexibility in the atomic ERIs tends to improve rela:ve energies for electronic states that differ by occupa:on of orbitals on single- centers. That is, for instance, excited states that are well characterized as one- electron local excita:ons. INDO s:ll finds substan:al use for spectroscopy.
20 NEGLECT OF DIATOMIC DIFFERENTIAL OVERLAP Reducing the N 4 scaling problem: ( µν λσ) = ϕ µ ( 1)ϕ ν 1 ( ) 1 r 12 ϕ λ 2 ( )ϕ σ ( 2)dr 1 dr 2 NDDO: ( µν λσ) = δ A µ X ν δ Bλ Y σ ( µν λσ) so, the only integrals that survive are those where both µ and ν are on the same atom, ( ss ss) = G ss ( ss pp) = G sp ( pp pp) = G pp ( pp pʹ pʹ ) = G pp ʹ ( sp sp) = L sp and λ and σ are on the same atom, but neither the two atoms need to be the same, nor do the two func:ons on the individual atoms need to be the same. For monatomic ERIs, the INDO- like values are retained. For diatomic ERIs, there are 100 possible integrals for s and p basis func:ons (try and verify that ) The diatomic integrals are not easily evaluated for Slater type orbitals. Instead, NDDO methods model the integrals as classical mul:pole interac:ons between the two atoms. Whether the mul:pole is a point charge (ss), a dipole (sp), or a quadrupole (pp or pp ) depends on the nature of the orbitals, and the orienta:on of the higher mul:poles does, too. The magnitudes of the dipoles and quadrupoles depend on the exponents of the Slater basis func:ons.
21 ONE- ELECTRON INTEGRALS µ 1 nuclei 2 2 ν Z k µ 1 ν r k Consider diagonal term, µ = ν: k µ µ Z k µ µ 1 nuclei µ Z r k µ 1 µ k µ r k nuclei = U µ Z k k k µ ( µµ s k s k ) k k µ where U is a parameter that should be nearly equal to the ioniza:on poten:al of an electron in atomic orbital µ and the alrac:on to a nucleus at the posi:ons of other atoms k is Z k :mes the repulsion with an electron in an s orbital on the same atom.
22 ONE- ELECTRON INTEGRALS µ 1 nuclei 2 2 ν Z k µ 1 ν r k k Consider off- diagonal term, µ ν, both on same atom µ ν Z k µ µ 1 nuclei ν Z r k µ 1 ν k µ r k nuclei = Z k k k µ ( µν s k s k ) k k µ where the first two terms are zero for s and p orbitals by symmetry considera:ons and the alrac:on to a nucleus at the posi:ons of other atoms k is Z k :mes the repulsion with an electron in an s orbital on the same atom.
23 ONE- ELECTRON INTEGRALS µ 1 nuclei 2 2 ν Z k µ 1 ν r k k Consider off- diagonal term, µ ν, each on a different atom ( = β + β µ ν) S 2 µν where β is our old friend from ückel theory, the resonance integral (a numerical parameter) and S is the overlap matrix computed from the Slater type orbitals. Note that the β values depend on both the atomic number and the orbital (i.e., carbon has a different β value for an s orbital than for a p orbital) What would a table of parameters look like for an NDDO model?
24 Semiempirical Molecular Orbital Theory Semiempirical parameteriza:ons and modern usage Video III.vii
25 Big Personali:es John Pople Michael Dewar
26 FINAL MOLECULAR- MECANICS- LIKE TOUCUP First adopted in MNDO, later in AM1 and PM3, etc. V N ( A,B) = Z A Z ( B s A s A s B s ) B 4 i=1 + Z AZ B a r A,i e b A,i r AB c A,i AB ( ) 2 + a B,i e b B,i ( r AB c B,i ) 2 where the first term says that nuclei repel each other propor:onal to the degree that their s electrons repel one another, but the following sum effec:vely creates gaussian shaped ripples, posi:ve or nega:ve, in the poten:al energy surface surrounding the atoms. Such an approach permits a fine tuning of atomic separa:ons that must be considered rather ad hoc, but not necessarily inconsistent with semiempirical principles. (Later extended by Jorgensen and co- workers using Pairwise Distance Directed Gaussians (PDDG) for MNDO and PM3.)
27 Chemical Intui:on vs Maths Z = Observables i Occurrences j ( calc i, j expt i, j ) 2 2 w i 1/2 Michael Dewar MNDO AM1 AMPAC Jimmy Stewart MNDO AM1 PM3 MOPAC
28 Errors in Computed eats of Forma:on (kcal/mol)
29 MODERN USAGE? SPECIFIC RANGE PARAMETERS Parameter AM1 AM1-SRP C U s C O O + + C 2 2 O U p ! s ! p O U s U p ! s ! p Source "E rxn D e (C ) D e ( ) AM1 AM1-SRP Expt
A One-Slide Summary of Quantum Mechanics
A One-Slide Summary of Quantum Mechanics Fundamental Postulate: O! = a! What is!?! is an oracle! operator wave function (scalar) observable Where does! come from?! is refined Variational Process H! = E!
More informationOne-Electron Hamiltonians
One-Electron Hamiltonians Hartree-Fock and Density Func7onal Theory Christopher J. Cramer @ChemProfCramer 2017 MSSC, July 10, 2017 REVIEW A One-Slide Summary of Quantum Mechanics Fundamental Postulate:
More informationThe successful wavefunction can be written as a determinant: # 1 (2) # 2 (2) Electrons. This can be generalized to our 2N-electron wavefunction:
T2. CNDO to AM1: The Semiempirical Molecular Orbital Models The discussion in sections T2.1 T2.3 applies also to ab initio molecular orbital calculations. T2.1 Slater Determinants Consider the general
More informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
More informationIntroduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 27, December 5, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 27, December 5, 2014 (Some material in this lecture has been adapted from Cramer, C. J.
More informationSame idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
More informationSemi-Empirical Methods CHEM 430
Semi-Empirical Methods CHEM 430 Cost, Hartree%Fock, scales,as,n 4,(N=#, basis,funcfons), Due,to,two% electron, integrals, within,fock, matrix, Semi%empirical,cut, cost,by,reducing, number,of, integrals,
More informationSemi-Empirical MO Methods
Semi-Empirical MO Methods the high cost of ab initio MO calculations is largely due to the many integrals that need to be calculated (esp. two electron integrals) semi-empirical MO methods start with the
More informationAPPENDIX 2. Semiempirical PM3 (Parametrization Method #3) When analytical solutions cannot be obtained, Schrödinger equation may be solved:
WRNING NOTICE: The experiments described in these materials are potentially hazardous and require a high level of safety training, special facilities and equipment, and supervision by appropriate individuals.
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationElectronega+vity Review
Electronega+vity Review Remember from the first course that electronega+vity is an es+mate of how atoms pull electrons towards themselves in a molecule. The higher the electron affinity, the more the element
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationLecture 9 Electronic Spectroscopy
Lecture 9 Electronic Spectroscopy Molecular Orbital Theory: A Review - LCAO approximaton & AO overlap - Variation Principle & Secular Determinant - Homonuclear Diatomic MOs - Energy Levels, Bond Order
More informationExercise 1: Structure and dipole moment of a small molecule
Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small
More informationBasis Sets and Basis Set Notation
Chemistry 46 Fall 215 Dr. Jean M. Standard November 29, 217 Basis Sets and Basis Set Notation Using the LCAO-MO approximation, molecular orbitals can be represented as linear combinations of atomic orbitals,
More informationPAPER:2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY. Module No. 34. Hückel Molecular orbital Theory Application PART IV
Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any), PHYSICAL -II QUANTUM Hückel Molecular orbital Theory Module No. 34 TABLE OF CONTENTS 1. Learning outcomes. Hückel Molecular Orbital (HMO) Theory
More informationComputational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica
Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron
More informationI. CSFs Are Used to Express the Full N-Electron Wavefunction
Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationChemistry 2. Your lecturers. Revision H 2 22/04/14. Lecture 1 Quantum Mechanics in Chemistry
Chemistry Lecture 1 Quantum Mechanics in Chemistry Your lecturers 8am Asaph Widmer-Cooper Room 316 asaph.widmer-cooper@sydney.edu.au 1pm Adam Bridgeman Room 543A adam.bridgeman@sydney.edu.au Revision H
More informationVALENCE Hilary Term 2018
VALENCE Hilary Term 2018 8 Lectures Prof M. Brouard Valence is the theory of the chemical bond Outline plan 1. The Born-Oppenheimer approximation 2. Bonding in H + 2 the LCAO approximation 3. Many electron
More informationσ u * 1s g - gerade u - ungerade * - antibonding σ g 1s
One of these two states is a repulsive (dissociative) state. Other excited states can be constructed using linear combinations of other orbitals. Some will be binding and others will be repulsive. Thus
More informationElectronic structure calculations: fundamentals George C. Schatz Northwestern University
Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions
More informationSection 2 Simple Molecular Orbital Theory
Section 2 Simple Molecular Orbital Theory In this section, the conceptual framework of molecular orbital theory is developed. Applications are presented and problems are given and solved within qualitative
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
More informationBe H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2
The Delocalized Approach to Bonding: The localized models for bonding we have examined (Lewis and VBT) assume that all electrons are restricted to specific bonds between atoms or in lone pairs. In contrast,
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More informationCHEMISTRY 4021/8021 MIDTERM EXAM 1 SPRING 2014
CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationSUPPLEMENTARY 3. H, B, C, N, O, F, Na, Mg, Al, Si, P, S, Cl, ZN, Ge, Br, Sn, I, and Hg.
SUPPLEMENTARY 3 Generate the Molecular Orbitals for CH 4 (T d ), CH 4 (D 4h ) and Cyclopropane using semiempirical module with the AM1 (Austin Model 1) Hamiltonian using HyperChem7.5. Please take the help
More informationMolecular Bonding. Molecular Schrödinger equation. r - nuclei s - electrons. M j = mass of j th nucleus m 0 = mass of electron
Molecular onding Molecular Schrödinger equation r - nuclei s - electrons 1 1 W V r s j i j1 M j m i1 M j = mass of j th nucleus m = mass of electron j i Laplace operator for nuclei Laplace operator for
More informationChemistry 3502/4502. Exam III. All Hallows Eve/Samhain, ) This is a multiple choice exam. Circle the correct answer.
B Chemistry 3502/4502 Exam III All Hallows Eve/Samhain, 2003 1) This is a multiple choice exam. Circle the correct answer. 2) There is one correct answer to every problem. There is no partial credit. 3)
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationElectronic Structure Models
Electronic Structure Models Hückel Model (1933) Basic Assumptions: (a) One orbital per atom contributes to the basis set; all orbitals "equal" (b) The relevant integrals involving the Hamiltonian are α
More informationMODELING MATTER AT NANOSCALES
MODELING MATTER AT NANOSCALES 6. The theory of molecular orbitals for the description of nanosystems (part II) 6.0. Ab initio methods. Basis functions. Luis A. Monte ro Firmado digitalmente por Luis A.
More informationThe Hartree-Fock approximation
Contents The Born-Oppenheimer approximation Literature Quantum mechanics 2 - Lecture 7 November 21, 2012 Contents The Born-Oppenheimer approximation Literature 1 The Born-Oppenheimer approximation 2 3
More informationDiatomic Molecules. 14th May Chemical Bonds in Diatomic Molecules: Overlaps and Delocalization of Electrons
Diatomic Molecules 14th May 2009 1 Chemical Bonds in Diatomic Molecules: Overlaps and Delocalization of Electrons 1.1 H + 2 Molecule Consider the process where 2 atomic nuclei and associated electron (1
More informationMolecular-Orbital Theory
Prof. Dr. I. Nasser atomic and molecular physics -551 (T-11) April 18, 01 Molecular-Orbital Theory You have to explain the following statements: 1- Helium is monatomic gas. - Oxygen molecule has a permanent
More informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 17, March 1, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 17, March 1, 2006 (Some material in this lecture has been adapted from Cramer, C. J.
More informationChapter 4 Symmetry and Chemical Bonding
Chapter 4 Symmetry and Chemical Bonding 4.1 Orbital Symmetries and Overlap 4.2 Valence Bond Theory and Hybrid Orbitals 4.3 Localized and Delocalized Molecular Orbitals 4.4 MX n Molecules with Pi-Bonding
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationChemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 25, November 5, 2003
Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer Lecture 25, November 5, 2003 (Some material in this lecture has been adapted from Cramer, C.
More informationThe symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then
1 The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then filled with the available electrons according to
More informationNH 3 H 2 O N 2. Why do they make chemical bonds? Molecular Orbitals
N 2 NH 3 H 2 O Why do they make chemical bonds? 5 Molecular Orbitals Why do they make chemical bonds? Stabilization Bond energy Types of Chemical Bonds Metallic Bond Ionic Bond Covalent Bond Covalent Bond
More informationInstructor background for the discussion points of Section 2
Supplementary Information for: Orbitals Some fiction and some facts Jochen Autschbach Department of Chemistry State University of New York at Buffalo Buffalo, NY 14260 3000, USA Instructor background for
More informationChem 1140; Molecular Modeling
P. Wipf 1 Chem 1140 $E = -! (q a + q b )" ab S ab + ab!k
More informationMolecular Bond Theory
Molecular Bond Theory Short comings of the localized electron model: electrons are not really localized so the concept of resonance was added no direct information about bond energies Molecular Orbital
More informationMolecular Geometry and Bonding Theories. Chapter 9
Molecular Geometry and Bonding Theories Chapter 9 Molecular Shapes CCl 4 Lewis structures give atomic connectivity; The shape of a molecule is determined by its bond angles VSEPR Model Valence Shell Electron
More informationQUANTUM MECHANICS AND MOLECULAR STRUCTURE
6 QUANTUM MECHANICS AND MOLECULAR STRUCTURE 6.1 Quantum Picture of the Chemical Bond 6.2 Exact Molecular Orbital for the Simplest Molecule: H + 2 6.3 Molecular Orbital Theory and the Linear Combination
More informationNew Methods for Computational Drug Design
Jimmy Charnley Kromann New Methods for Computational Drug Design M A S T E R S T H E S I S Department of Chemistry University of Copenhagen Abstract This thesis describes the work that has been carried
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationChemistry 2000 Lecture 8: Valence bond theory
Chemistry 000 Lecture 8: Valence bond theory Marc R. Roussel January 9, 08 Marc R. Roussel Valence bond theory January 9, 08 / 5 MO theory: a recap A molecular orbital is a one-electron wavefunction which,
More informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More informationPHYSICAL CHEMISTRY I. Chemical Bonds
PHYSICAL CHEMISTRY I Chemical Bonds Review The QM description of bonds is quite good Capable of correctly calculating bond energies and reaction enthalpies However it is quite complicated and sometime
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 4 Molecular orbitals C.-K. Skylaris Learning outcomes Be able to manipulate expressions involving spin orbitals and molecular orbitals Be able to write down
More information1s atomic orbital 2s atomic orbital 2s atomic orbital (with node) 2px orbital 2py orbital 2pz orbital
Atomic Orbitals 1s atomic orbital 2s atomic orbital 2s atomic orbital (with node) 2px orbital 2py orbital 2pz orbital Valence Bond Theory and ybridized Atomic Orbitals Bonding in 2 1s 1s Atomic Orbital
More informationPractical Aspects of Electronic Structure. Electronelectron. Nuclear-nuclear repulsion. repulsion. 2 α β
1. Born-Oppenheimer Approximation Need to elaborate our models to handle molecules. More particles, more variables to keep track of. For H, for instance, would have to solve for energy E and wavefunction
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More information5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.61 Physical Chemistry Exam III (1) PRINT your name on the cover page. (2) It is suggested that you READ THE ENTIRE EXAM before
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 31, April 12, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 31, April 12, 2006 (Some material in this lecture has been adapted from Cramer, C.
More informationCHEMISTRY Topic #1: Bonding What Holds Atoms Together? Spring 2012 Dr. Susan Lait
CHEMISTRY 2000 Topic #1: Bonding What Holds Atoms Together? Spring 2012 Dr. Susan Lait Why Do Bonds Form? An energy diagram shows that a bond forms between two atoms if the overall energy of the system
More informationLecture 4: methods and terminology, part II
So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,
More informationDUANE E. WILLIAMS UNIVERSITY OF FLORIDA
PARAMETERIZATION OF SEMIEMPIRICAL QUANTUM MECHANICAL METHODS FOR THE PREDICTION OF NUCLEAR MAGNETIC RESONANCE CHEMICAL SHIFTS IN BIOLOGICALLY RELEVANT SYSTEMS By DUANE E. WILLIAMS A DISSERTATION PRESENTED
More informationuse the backs of pages as needed
CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy
More informationTheoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics
Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics Problem 1 Draw molecular orbital diagrams for O 2 and O 2 +. E / ev dioxygen molecule, O 2 dioxygenyl cation, O 2 + 25
More informationAlkali metals show splitting of spectral lines in absence of magnetic field. s lines not split p, d lines split
Electron Spin Electron spin hypothesis Solution to H atom problem gave three quantum numbers, n,, m. These apply to all atoms. Experiments show not complete description. Something missing. Alkali metals
More informationChemistry 543--Final Exam--Keiderling May 5, pm SES
Chemistry 543--Final Exam--Keiderling May 5,1992 -- 1-5pm -- 174 SES Please answer all questions in the answer book provided. Make sure your name is clearly indicated and that the answers are clearly numbered,
More informationChapter 9. Covalent Bonding: Orbitals
Chapter 9. Covalent onding: Orbitals Models to explain the structures and/or energies of the covalent molecules Localized Electron (LE) onding Model Lewis Structure Valence Shell Electron Pair Repulsion
More informationModule:32,Huckel Molecular Orbital theory- Application Part-II PAPER: 2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY
Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any) 2, PHYSICAL -I QUANTUM Hückel Molecular orbital Theory Application PART II Module No. 32 TABLE OF CONTENTS 1. Learning outcomes 2. Hückel Molecular
More informationMolecular Orbital Approach to Bonding
Molecular Orbital Approach to Bonding Chemistry 362; spring 2019 Marcetta Y. Darensbourg, Professor Xuemei Yang, Graduate Assistant Kyle Burns, Graduate Assistant The following slides were modified from
More informationCovalent Bonding: Orbitals
Hybridization and the Localized Electron Model Covalent Bonding: Orbitals A. Hybridization 1. The mixing of two or more atomic orbitals of similar energies on the same atom to produce new orbitals of equal
More informationWe will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include,
1 Common Notation We will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include, µ, ν, λ, σ: Primary basis functions, i.e., atomic orbitals (AOs) {φ µ ( r)}. The size
More informationE nuc rep = Z AZ B r AB m =0.499 hartree =13.60 ev. E nuc rep (H 2 )= m/a m =0.714 hartree =19.43 ev.
Chemistr 31 Phsical Chemistr Homework Assignment # 7 1. Sketch qualitative contour maps of the following molecular orbitals for a diatomic molecule AB. Identif each as a bonding or an anti-bonding molecular
More informationLecture B6 Molecular Orbital Theory. Sometimes it's good to be alone.
Lecture B6 Molecular Orbital Theory Sometimes it's good to be alone. Covalent Bond Theories 1. VSEPR (valence shell electron pair repulsion model). A set of empirical rules for predicting a molecular geometry
More informationLecture 33: Intermolecular Interactions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 5.61 Physical Chemistry I Fall, 2017 Professors Robert W. Field Lecture 33: Intermolecular Interactions Recent Lectures Non-degenerate Perturbation Theory vs. Variational
More information7: Hückel theory for polyatomic molecules
7: ückel theory for polyatomic molecules Introduction Approximate treatment of π electron systems in organic molecules: 1 Approximations 3 4 5 6 1. π and σ frameworks completely separated. Trial wavefunctions
More informationHomework Problem Set 5 Solutions. E e + H corr (a.u.) a.) Determine the bond dissociation enthalpy of ethane in kcal/mol (1 a.u. = kcal/mol).
Chemistry 380.37 Dr. Jean M. Standard Homework Problem Set 5 Solutions 1. Given below are the sum of electronic and thermal enthalpies, E e + H corr, from Hartree-Fock calculations using a 6-31G(d) basis
More informationApplication of electrostatic Solvent Models to the Electronic Spectrum of Mesityl Oxide
Journal of Basrah Researches ((Sciences Volume 38. Number 1. ((2012 vailable online at: www.basra-science journal.org 2695 1817 ISSN pplication of electrostatic Solvent Models to the Electronic Spectrum
More informationOn the Uniqueness of Molecular Orbitals and limitations of the MO-model.
On the Uniqueness of Molecular Orbitals and limitations of the MO-model. The purpose of these notes is to make clear that molecular orbitals are a particular way to represent many-electron wave functions.
More informationChapter IV: Electronic Spectroscopy of diatomic molecules
Chapter IV: Electronic Spectroscopy of diatomic molecules IV.2.1 Molecular orbitals IV.2.1.1. Homonuclear diatomic molecules The molecular orbital (MO) approach to the electronic structure of diatomic
More informationHückel Molecular orbital Theory Application PART III
Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any) 2, PHYSICAL -I QUANTUM Hückel Molecular orbital Theory Module No. 33 PAPER: 2, PHYSICAL -I TABLE OF CONTENTS 1. Learning outcomes 2. Hückel
More informationPhysical Chemistry II Recommended Problems Chapter 12( 23)
Physical Chemistry II Recommended Problems Chapter 1( 3) Chapter 1(3) Problem. Overlap of two 1s orbitals recprobchap1.odt 1 Physical Chemistry II Recommended Problems, Chapter 1(3), continued Chapter
More information1 r A. r B. 2m e. The potential energy of the electron is. r A and r B are the electron s distances from the nuclei A and B. This expression can be
Introduction to Molecular Structure The Born-Oppenheimer approximation The Born-Oppenheimer approximation supposes that the nuclei, being so much heavier than the electron, move relatively slow and may
More informationLecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.
Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem
More informationMolecular Shape and Molecular Polarity. Molecular Shape and Molecular Polarity. Molecular Shape and Molecular Polarity
Molecular Shape and Molecular Polarity When there is a difference in electronegativity between two atoms, then the bond between them is polar. It is possible for a molecule to contain polar bonds, but
More informationJack Simons Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationChemistry 4681 Module: Electronic Structure of Small Molecules Background Handout
Chemistry 4681 Module: Electronic Structure of Small Molecules Background Handout C. David Sherrill Last Revised: 6 January 2000 1 Computational Chemistry The term computational chemistry is used to mean
More informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 21, November 12, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 204 Laura Gagliardi Lecture 2, November 2, 204 (Some material in this lecture has been adapted from Cramer, C. J. Essentials
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationGeneral Physical Chemistry II
General Physical Chemistry II Lecture 13 Aleksey Kocherzhenko October 16, 2014" Last time " The Hückel method" Ø Used to study π systems of conjugated molecules" Ø π orbitals are treated separately from
More informationNov 12, 2014 LB273 Prof. Vash6 Sawtelle. Today s Topics: SHM. h<ps:// v=6sas9n3omum&feature=youtu.be&t=16s
Nov 12, 2014 LB273 Prof. Vash6 Sawtelle Today s Topics: SHM h
More information