Yingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
|
|
- Denis Mitchell
- 5 years ago
- Views:
Transcription
1 Purdue University CHM 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 approximation 2 4 Variational principle 4 5 Orbitals Kinds of orbitals Notations of spin orbitals Hartree product Hartree-Fock energy Two electrons Many electrons Properties Self-consistent procedure 8 8 Basis sets Basic concepts Examples Koopmans theorem 11 1
2 1 Hartree energy The hartree (symbol: Eh or Ha), also known as the Hartree energy, is the atomic unit of energy, named after the British physicist Douglas Hartree. 2 E h = m e a 2, (1.1) 0 = J, (1.2) = ev, (1.3) = kcal/mol. (1.4) Table 1: Energy levels Phenomenon Energy level nuclear physics E h ionization, excitation, dissociation, chemistry reaction E h vibration E h rotation E h spin < 10 4 E h 2 Many-body system The many body Hamiltonian is Ĥ = ˆT e + ˆT N + ˆV ee + ˆV NN + ˆV en. (2.1) Table 2: Items in many body Hamiltonian Name Notation Sign i 2 i Kinetic energy of electron ˆTe = 1 2 Kinetic energy of nuclei ˆTN = 1 Electron repulsion ˆVee = i>j Nuclear repulsion ˆVNN = A>B 1 Electron-nuclear attraction ˆV en = A,i 2M A A 2 A Negative sign comes from momentum operator p = i 1 r ij Positive sign comes from charge ( e)( e) = e 2 r AB Z A r Ai Negative sign comes from charge Z A ( e)e 3 Born-Oppenheimer approximation Main idea The core of Born-Oppenheimer approximation is that electron and nuclear degrees of freedom are separable. Electronic Schrödinger Equation After Born-Oppenheimer approximation, the electronic Hamiltonian Ĥel is defined as follows: 2
3 Table 3: Changes in Born-Oppenheimer approximation Name Notation Changes Kinetic energy of electron ˆTe = 1 2 i 2 i unchanged Kinetic energy of nuclei ˆTN = 1 2M A A 2 A 0 since M A m e Electron repulsion ˆVee = 1 i>j r ij unchanged Nuclear repulsion ˆVNN = A>B 1 r AB becomes a constant Electron-nuclear attraction ˆVeN = Z A A becomes parametric dependence, A,i r Ai instead of variable dependence. Table 4: Operators about electronic Hamiltonian Operators Notations 0-electron operator V NN = Z AZ B A>B r AB 1-electron operator ĥ(i) = i Z A A r ia 2-electron operator ˆv(i,j) = 1 r ij Electronic Hamiltonian Ĥ el = iĥ(i)+ i>j 1 r ij +V NN Exception Born-Oppenheimer approximation fails when the electron and nuclei motion is coupled. An example is malondialdehyde CH 2 (CHO) 2, where H is not only a particle but also a wave. See Fig.1. Figure 1: Malondialdehyde CH 2 (CHO) 2 3
4 4 Variational principle The variational principle has two parts: 1. The complete set of exact eigenfunctions of a Hamiltonian H define an orthogonal complete basis set for the total space of wave functions. If so, it is possible to express any function as a linear combination of the exact eigenfunctions. 2. If φ is any normalised, well-behaved function that satisfies the boundary conditions of the Hamiltonian, then the energy associated with φ is always above that of the eigenfunction of lowest energy, i.e. φ H φ E 0, (4.1) where E 0 is the true value of the lowest energy eigenvalue of H. This principle allows us to calculate an upper bound for the ground state energy by finding the trial wavefunction φ for which the integral is minimised (hence the name trial wavefunctions are varied until the optimum solution is found). 5 Orbitals 5.1 Kinds of orbitals Table 5: Orbitals Name Notation Functions Orbital a wave function for a single particle, a electron Spatial orbital φ i (r) describes the spatial distribution of an electron such that φ i (r) 2 dr is the probability of finding the particle near r Spin orbital χ i (x) = φ(r)α(ω) or φ(r)β(ω) a wave function for an electron that describes both its distribution and spin Molecular orbitals the wave functions for electrons in molecular Hartree product ψ HP = χ i (x 1 )χ j (x 2 ) χ k (x p ) Note: spin orbitals, instead of spatial obitals 5.2 Notations of spin orbitals Sometimes it is convenient to use a notation that indicates a spin orbital by its spatial part, using a bar or lack of bar to denote whether it has the β or α spin. In this notation, the Hartree Fock ground state of H 2 is ψ 0 = χ 1 χ 2 = φ 1 φ1 = 1 1, which indicates that both electrons occupy the same spatial orbital φ 1, but one has α spin and the other has β spin. 5.3 Hartree product Proposition 5.1. Define the Hartree product Ψ HP (x 1,x 2,,x N ) = χ i (x 1 )χ j (x 2 ) χ k (x N ). (5.1) 4
5 Consider the Hamiltonian in which H HP = h(i), (5.2) i=1 h(i)χ j (x i ) = ε j χ j (x i ). (5.3) Shown that the Hartree product (5.1) is an eigenfunction of the Hamiltonian (5.2) with an eigenvalue given by E = ε i +ε j + +ε k. (5.4) Proof. According to (5.3) and (5.1), we have where E is given by (5.4). h(i)ψ HP = ε i Ψ HP, h(i)ψ HP = (ε i +ε j + +ε k )Ψ HP, i=1 H HP Ψ HP = EΨ HP, 6 Hartree-Fock energy 6.1 Two electrons The Hamiltonian of system with 2-electrons is ( Ĥel N=2 = ) ( Z A + 1 r 1A A A The trail function of slater type is The corresponding energy is Z A r 2A ) + 1 r 12 (6.1) = ĥ(1)+ĥ(2)+ 1 r 12. (6.2) ψ 0 = 1 2 (χ 1 (1)χ 2 (2) χ 1 (2)χ 2 (1)). (6.3) E N=2 HF = ψ 0 ĤN=2 el ψ 0, (6.4) = d1χ 1(1)ĥ(1)χ 1(1)+ d2χ 2(2)ĥ(2)χ 2(2) + d1d2χ 1 (1)χ 2 (2) 1 χ 1 (1)χ 2 (2) d1d2χ 1 r (1)χ 2 (2) 1 χ 2 (1)χ 1 (2)(6.5) 12 r 12 = h 1 +h 2 +J 12 K 12, (6.6) where J 12 is called coulomb term and K 12 is called exchange term. 5
6 6.2 Many electrons The energy of N electron system is E HF = ψ 0 Ĥel ψ 0 (6.7) = h i + (J ij K ij ), (6.8) = i=1 h i i=1 where J ij and K ij are defined in (6.10)-(6.13). 6.3 Properties (ij kl) = ij kl = Define some operators i=1 j>i i=1 j=1 (J ij K ij ), (6.9) dr 1 dr 2 ψ i(r 1 )ψ j (r 1 )r 1 12 ψ k(r 2 )ψ l (r 2 ), (6.10) dr 1 dr 2 ψ i (r 1)ψ j (r 2)r 1 12 ψ k(r 1 )ψ l (r 2 ), (6.11) J ij = (ii jj) = ij ij, called coulomb integral (6.12) K ij = (ij ji) = ij ji, called exchange integral. (6.13) Ĵ i x j (2) = x i (1) 1 x i (1) x j (2), r 12 (6.14) ˆK i x j (2) = x i (1) 1 x j (1) x i (2), r 12 (6.15) ˆf i = ĥi + (Ĵj ˆK j ), Fock operator, (6.16) x i ˆf i x i = ĥi + j (J ij K ij ). (6.17) Proposition 6.1. Prove the following properties of coulomb and exchange integrals: j J ii = K ii, (6.18) J ij = J ij, (6.19) K ij = K ij, (6.20) J ij = J ji, (6.21) K ij = K ji. (6.22) Proof. According to the notations (6.10)-(6.13), it is easy to know the following. J ii = (ii ii) = K ii ; J ij = (ii jj) = J ij ; 6
7 K ij = (ji ij), = dr 1 dr 2 ψj (r 1)ψ i (r 1 )r12 1 ψ i (r 2)ψ j (r 2 ), = dr 1 dr 2 ψj(r 2 )ψ i (r 2 )r12 1 ψ i(r 1 )ψ j (r 1 ), = dr 1 dr 2 ψi(r 1 )ψ j (r 1 )r12 1 ψ j(r 2 )ψ i (r 2 ), = (ij ji) = K ij ; J ij = (ii jj), = dr 1 dr 2 ψi (r 1)ψ i (r 1 )r12 1 ψ j (r 2)ψ j (r 2 ), = dr 1 dr 2 ψi (r 2)ψ i (r 2 )r12 1 ψ j (r 1)ψ j (r 1 ), = dr 1 dr 2 ψj(r 1 )ψ j (r 1 )r12 1 ψ i(r 2 )ψ i (r 2 ), = (jj ii) = J ji ; K ij = (ij ji), = dr 1 dr 2 ψi (r 1)ψ j (r 1 )r12 1 ψ j (r 2)ψ i (r 2 ), = dr 1 dr 2 ψi (r 2)ψ j (r 2 )r12 1 ψ j (r 1)ψ i (r 1 ), = dr 1 dr 2 ψj(r 1 )ψ i (r 1 )r12 1 ψ i(r 2 )ψ j (r 2 ), = (ji ij) = K ji. Proposition 6.2. Verify the energies of the determinants shown in Fig.2 by inspection. a. h 11 +h 22 +J 12 K 12, b. h 11 +h 22 +J 12, c. 2h 11 +J 11, d. 2h 2 +J 22, e. 2h 11 +h 22 +J 11 +2J 12 K 12, f. 2h 22 +h 11 +J 22 +2J 12 K 12, g. 2h 11 +2h 22 +J 11 +J 22 +4J 12 2K 12. Proof. Please see Table 6. Remark 6.1. Different spin means no exchange K ij. 7
8 Figure 2: Energy determinants Table 6: Interpretation of determinantal energies Number one-electron Coulomb Exchange a h 11 +h 22 J 12 K 12 b h 11 +h 22 J 12 0 c 2h 11 J 11 0 d 2h 2 J 22 0 e 2h 11 +h 22 J 11 +2J 12 K 12 f 2h 22 +h 11 J 22 +2J 12 K 12 g 2h 11 +2h 22 J 11 +J 22 +4J 12 2K 12 7 Self-consistent procedure Nonlinear Hartree-Fock equations can be solved by a self-consistent procedure, as follows 1. An initial guess orbitals (matrix C 0 ) are generated; F(C)C = SCE (7.1) 2. Fock operator is calculated for this set of orbitals (F 0 = F(C 0 )); 3. Fock operator is diagonalized, F 0 C 1 = SC 1 E 1 (called Roothaan equations), and new orbitals C 1 replace the orbitals C 0 from previous step; 4. Steps 2-3 are repeated until the difference between orbitals from step n and step n+1 is below certain (user-specified) threshold. Remark 7.1. A typical sailing of HF method is: diagonalization K 3 and form Fock matrix K 4. 8 Basis sets 8.1 Basic concepts Slater type and Gassuian type Double zeta and triple zeta 1. Polarization functions: adding d-type basis functions to the heavy atoms (, Li-F) and p-type functions to hydrogen ( ). Note that He does not need polarization. 2. Diffusion functions: (1) Charged particles, such as HSO 4 ; (2) He 2. Remark 8.1. Polarization is more expensive than diffusion. 8
9 8.2 Examples Question: Determine the total number of basis functions and primitive basis functions for STO-3G, 6-31G, 6-31G**, and 6-311(+,+)G** calculations of methane (CH 4 ). Write the contraction scheme in general notations for each basis (Szabo & Ostlund, chapter 3.6). Assume that pure angular momentum (5d,7f etc functions) polarization functions are used. Answer: 1. The STO-3G basis set for methane consists of one 1s orbital on each H atom and a 1s, 2s, and set of three 2p orbitals on C. See Tables 7-9 for details. The contraction scheme in general notations is (6s3p/3s)[2s1p/1s]. Table 7: STO-3G: Carbon atom 1s 2s 2p(3) Total No. Basis Functions No. Gaussians for each orbital No. Primitives Table 8: STO-3G: Hydrogen atom 1s Total No. Basis Functions 1 1 No. Gaussians for each orbital 3 No. Primitives 3 3 Table 9: STO-3G: Methane molecule No. Basis Functions = 9 No. Primitives = The 6-31G basis set is a split valence double zeta basis set. For hydrogen atom, a split valence double zeta basis consists of two 1s orbitals, denoted 1s and 1s. Note that for the H atom, because the 1s electron is considered the valence shell a double zeta basis set is used. For the carbon atom, a split valence double zeta basis set consists of a single 1s orbital, along with two 2s and two each of 2px, 2py, and 2pz orbitals, denoted as 1s, 2s, 2s, 2p(3) and 2p (3). See Tables for details. The contraction scheme in general notations is (10s4p/4s)[3s2p/2s]. Table 10: 6-31G: Carbon atom 1s 2s 2s 2p(3) 2p (3) Total No. Basis Functions No. Gaussians for each orbital No. Primitives
10 Table 11: 6-31G: Hydrogen atom 1s 1s Total No. Basis Functions No. Gaussians for each orbital 3 1 No. Primitives Table 12: 6-31G: Methane molecule No. Basis Functions = 17 No. Primitives = The 6-31G** basis set is a split valence double zeta basis set with added polarization functions on both carbon and hydrogen atoms. For hydrogen atom, the split valence double zeta basis consists of two 1s orbitals, and extra polarization functions 2px, 2py, and 2pz orbitals, denoted 1s, 1s and 2p(3). For the carbon atom, a split valence double zeta basis set consists of a single 1s orbital, along with two 2s and a double set of 2p orbitals, 2p(3)and 2p. (3), andextrapolarizationfunctions d(5). See Tables13-15fordetails. The contraction scheme in general notations is (10s4p1d/4s1p)[3s2p1d/2s1p]. Table 13: 6-31G**: Carbon atom 1s 2s 2s 2p(3) 2p (3) d(5) Total No. Basis Functions No. Gaussians for each orbital No. Primitives Table 14: 6-31G**: Hydrogen atom 1s 1s 2p(3) Total No. Basis Functions No. Gaussians for each orbital No. Primitives Table 15: 6-31G**: Methane molecule No. Basis Functions = 34 No. Primitives = 55 10
11 4. The 6-311(+,+)G** basis set is a split valence triple zeta basis set with added diffuse and polarization functions on both carbon and hydrogen atoms. For hydrogen atom, the split valence triple zeta basis consists of three 1s orbitals 1s,1s,1s, and diffuse function 2s, and extra polarization functions 2px, 2py, and 2pz orbitals, 2p(3). For the carbon atom, a split valence triple zeta basis set consists of a single 1s orbital, along with three 2s orbitals 2s,2s,2s, and a triple set of 2p orbitals, 2p(3), 2p (3) and 2p (3), and diffuse functions 3s, 3p(3), and extra polarization functions 3d(5). See Tables for details. The contraction scheme in general notations is (12s6p1d/6s1p)[5s4p1d/4s1p]. Table 16: 6-311(+,+)G**: Carbon atom 1s 2s 2s 2s 2p(3) 2p (3) 2p (3) 3s 3p d(5) Total No. Basis Functions No. Gaussians for each orbital No. Primitives Table 17: 6-311(+,+)G**: Hydrogen atom 1s 1s 1s 2s 2p(3) Total No. Basis Functions No. Gaussians for each orbital No. Primitives Table 18: 6-311(+,+)G**: Methane molecule No. Basis Functions = 50 No. Primitives = 71 9 Koopmans theorem Theorem 9.1 (Koopmans). This theorem has two parts: 1. The ionization potential (IP) for removing an electron from χ c is just the negative of the orbital energy ε c, IP = N 1 E c N E 0 = ε c > 0. (9.1) 2. The electron affinity (EA) for adding an electron to the virtual spin orbital χ r is just the negative of the orbital energy of that virtual spin orbital, i.e. EA = N E 0 N+1 E r = ε r. (9.2) If ε r is negative (i.e. if N+1 Ψ r ) is more stable than N Ψ 0, the electron affinity is positive. 11
The 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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
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 informationA 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 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 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 informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
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 informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More information4πε. me 1,2,3,... 1 n. H atom 4. in a.u. atomic units. energy: 1 a.u. = ev distance 1 a.u. = Å
H atom 4 E a me =, n=,,3,... 8ε 0 0 π me e e 0 hn ε h = = 0.59Å E = me (4 πε ) 4 e 0 n n in a.u. atomic units E = r = Z n nao Z = e = me = 4πε = 0 energy: a.u. = 7. ev distance a.u. = 0.59 Å General results
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 informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationExam 4 Review. Exam Review: A exam review sheet for exam 4 will be posted on the course webpage. Additionally, a practice exam will also be posted.
Chem 4502 Quantum Mechanics & Spectroscopy (Jason Goodpaster) Exam 4 Review Exam Review: A exam review sheet for exam 4 will be posted on the course webpage. Additionally, a practice exam will also be
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 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 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 informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
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 informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
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 informationIntroduction to Computational Quantum Chemistry: Theory
Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction
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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
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 informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
More informationHartree-Fock Theory. ˆf(r)χ i (x) = ɛ i χ i (x) (1) Z A K=1 I. DERIVATION OF THE HARTREE-FOCK EQUATIONS. A. The energy of a Slater determinant
Hartree-Fock Theory The HF approximation plays a crucial role in chemistry and constitutes the starting point for more elaborate treatments of electron correlation. Furthermore, many semi-empirical methods
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationonly two orbitals, and therefore only two combinations to worry about, but things get
131 Lecture 1 It is fairly easy to write down an antisymmetric wavefunction for helium since there are only two orbitals, and therefore only two combinations to worry about, but things get complicated
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationDiatomic Molecules. 7th May Hydrogen Molecule: Born-Oppenheimer Approximation
Diatomic Molecules 7th May 2009 1 Hydrogen Molecule: Born-Oppenheimer Approximation In this discussion, we consider the formulation of the Schrodinger equation for diatomic molecules; this can be extended
More information2.4. Quantum Mechanical description of hydrogen atom
2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022
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 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 informationCh 125a Problem Set 1
Ch 5a Problem Set Due Monday, Oct 5, 05, am Problem : Bra-ket notation (Dirac notation) Bra-ket notation is a standard and convenient way to describe quantum state vectors For example, φ is an abstract
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 informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationWe also deduced that transformations between Slater determinants are always of the form
.3 Hartree-Fock The Hartree-Fock method assumes that the true N-body ground state wave function can be approximated by a single Slater determinant and minimizes the energy among all possible choices of
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 informationIntroduction to Computational Chemistry: Theory
Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationQuantum Theory of Many-Particle Systems, Phys. 540
Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
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 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 informationHandbook of Computational Quantum Chemistry
Handbook of Computational Quantum Chemistry David B. Cook Dept. of Chemistry University of Sheffield DOVER PUBLICATIONS, INC. Mineola, New York F Contents 1 Mechanics and molecules 1 1.1 1.2 1.3 1.4 1.5
More information0 belonging to the unperturbed Hamiltonian H 0 are known
Time Independent Perturbation Theory D Perturbation theory is used in two qualitatively different contexts in quantum chemistry. It allows one to estimate (because perturbation theory is usually employed
More informationLecture 10. Born-Oppenheimer approximation LCAO-MO application to H + The potential energy surface MOs for diatomic molecules. NC State University
Chemistry 431 Lecture 10 Diatomic molecules Born-Oppenheimer approximation LCAO-MO application to H + 2 The potential energy surface MOs for diatomic molecules NC State University Born-Oppenheimer approximation
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 informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationH 2 in the minimal basis
H 2 in the minimal basis Alston J. Misquitta Centre for Condensed Matter and Materials Physics Queen Mary, University of London January 27, 2016 Overview H 2 : The 1-electron basis. The two-electron basis
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 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 informationONE AND MANY ELECTRON ATOMS Chapter 15
See Week 8 lecture notes. This is exactly the same as the Hamiltonian for nonrigid rotation. In Week 8 lecture notes it was shown that this is the operator for Lˆ 2, the square of the angular momentum.
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 informationChemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron):
April 6th, 24 Chemistry 2A 2nd Midterm. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (-electron): E n = m e Z 2 e 4 /2 2 n 2 = E Z 2 /n 2, n =, 2, 3,... where Ze is
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationChapter 9: Multi- Electron Atoms Ground States and X- ray Excitation
Chapter 9: Multi- Electron Atoms Ground States and X- ray Excitation Up to now we have considered one-electron atoms. Almost all atoms are multiple-electron atoms and their description is more complicated
More informationSecond quantization. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
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 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 informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationMOLECULAR STRUCTURE. The general molecular Schrödinger equation, apart from electron spin effects, is. nn ee en
MOLECULAR STRUCTURE The Born-Oppenheimer Approximation The general molecular Schrödinger equation, apart from electron spin effects, is ( ) + + V + V + V =E nn ee en T T ψ ψ n e where the operators in
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationLecture 9: Molecular Orbital theory for hydrogen molecule ion
Lecture 9: Molecular Orbital theory for hydrogen molecule ion Molecular Orbital Theory for Hydrogen Molecule Ion We have seen that the Schrödinger equation cannot be solved for many electron systems. The
More informationStructure of diatomic molecules
Structure of diatomic molecules January 8, 00 1 Nature of molecules; energies of molecular motions Molecules are of course atoms that are held together by shared valence electrons. That is, most of each
More informationIntroduction and theoretical background
1 Introduction and theoretical background 1.1 The Schrödinger equation and models of chemistry The Schrödinger equation and its elements As early as 1929, the noted physicist P. A. M. Dirac wrote 1 The
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 informationExp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 +
Exp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 + 1. Objectives Quantum chemical solvers are used to obtain the energy and the orbitals of the simplest molecules
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
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 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 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 informationChemistry 2000 Lecture 1: Introduction to the molecular orbital theory
Chemistry 2000 Lecture 1: Introduction to the molecular orbital theory Marc R. Roussel January 5, 2018 Marc R. Roussel Introduction to molecular orbitals January 5, 2018 1 / 24 Review: quantum mechanics
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 informationQuantum Chemistry Methods
1 Quantum Chemistry Methods T. Helgaker, Department of Chemistry, University of Oslo, Norway The electronic Schrödinger equation Hartree Fock theory self-consistent field theory basis functions and basis
More informationVol. 9 COMPUTATIONAL CHEMISTRY 319
Vol. 9 COMPUTATIONAL CHEMISTRY 319 COMPUTATIONAL QUANTUM CHEMISTRY FOR FREE-RADICAL POLYMERIZATION Introduction Chemistry is traditionally thought of as an experimental science, but recent rapid and continuing
More informationChapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found.
Chapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found. In applying quantum mechanics to 'real' chemical problems, one is usually faced with a Schrödinger
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationChemistry 881 Lecture Topics Fall 2001
Chemistry 881 Lecture Topics Fall 2001 Texts PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon MATHEMATICS for PHYSICAL CHEMISTRY, Mortimer i. Mathematics Review (M, Chapters 1,2,3 & 4; M&S,
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 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 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 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 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 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 informationATOMIC STRUCTURE. Atomic Structure. Atomic orbitals and their energies (a) Hydrogenic radial wavefunctions
ATOMIC STRUCTURE Atomic orbitals and their energies (a) Hydrogenic radial wavefunctions Bundet Boekfa Chem Div, Fac Lib Arts & Sci Kasetsart University Kamphaeng Saen Campus 1 2 Atomic orbitals and their
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationLecture 3: methods and terminology, part I
Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,
More information5.61 Physical Chemistry Final Exam 12/16/09. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry
5.6 Physical Chemistry Final Exam 2/6/09 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.6 Physical Chemistry Final Examination () PRINT your name on the cover page. (2) It
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 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 information1 Rayleigh-Schrödinger Perturbation Theory
1 Rayleigh-Schrödinger Perturbation Theory All perturbative techniques depend upon a few simple assumptions. The first of these is that we have a mathematical expression for a physical quantity for which
More information