SCF calculation on HeH +
|
|
- Violet Townsend
- 5 years ago
- Views:
Transcription
1 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 two-electron system HeH +. The text follows the procedures outlined in Szabo and Ostlund in Modern Quantum Chemistry and is directly related to the computer code hehp.f provided in the lectures. 1
2 Procedures The calculations are carried out with a minimal basis set: STO-3G. For an SCF calculation one first chooses a geometry of the nuclei which are located at R 1 and R 2 with an internuclear distance R = R 1 R 2. In the program R = a 0 is used. Next, the basis functions for the H and He atoms need to be specified. Here, a contraction (sum) of 3 Gaussian orbitals approximating the Slater orbital of H and He are used. The two Slater functions are φ 1 (ξ 3 1 /π)1/2 exp ξ 1 r R 1 φ 2 (ξ 3 2 /π)1/2 exp ξ 2 r R 2 (1) (2) where φ 1 describes the electron of He and φ 2 that of H. The standard Slater exponent ξ 2 for H in a molecule is ξ 2 = 1.24 a 1 0 and for He a variational calculation together with a correction for the positive charge of the molecule yields ξ 1 = a 1 0. Both Slater functions in Eqs. 2 are now represented as a sum (i.e. contraction ) over 3 Gauss functions φ CGF = 3 p=1 d pµφp GF (α pµ, r R a ) (therefore STO-3G) with φp GF (α pµ, r R a )=(2α pµ /π) 3/4 exp α pµ r R a 2. They are for H atom (including the normalization of (2α/π) 3/4 ): φ CGF = φ GF and for the He atom: (α = ) φ GF GF (α = ) φ (α = ) (3) φ CGF = φ GF (α = ) φ GF GF (α = ) φ (α = ) In a next step the necessary integrals over the basis functions are evaluated. Here the formulae involving functions are summarized. What is required are a) overlap integrals (2- center overlap), b) electron kinetic energy integrals, c) nuclear attraction and d) two-electron (4) 2
3 repulsion. The reason why Gaussian functions are used is the fact that the product of two Gaussian functions centered at R A and R B is again a Gaussian centered at R P. With g ( r R A ) = exp α r R A 2 one obtains g ( r R A )g ( r R B ) = exp α r R A 2 exp β r R B 2 = Kg ( r R P ) (5) with K = exp [ αβ/(α+β) R A R B 2 ] (6) and R P = (α R A + β R B )/(α + β) (7) Integrals a) and b) above can be evaluated using standard integrals over Gauss functions, c) and d) use Fourier transform techniques. The standard integral required is 0 drr 2 exp pr2 = (π/p) 3/2 /(4π) (8) With this, the overlap integral a) can be evaluated as (A B) = d r 1 g ( r 1 R A )g ( r 1 R B ) = [π/(α + β)] 3/2 exp [αβ/(α+β) R A R B 2 ] (9) and the kinetic integral b) yields (A B) = αβ 2αβ [3 (α + β) (α + β) R A R B 2 π ][ (α + β) ]3/2 exp [αβ/(α+β) R A R B 2 ] (10) For integrals c) and d) the following expressions are obtained: (A Z C B) = 2π r 1C α + β Z C exp [ αβ α+β R A R B 2 ] Fo [(α + β) R P R C 2 ] (11) 3
4 where F o is the error function F o (t) = 1 2 (π t )1/2 er f(t 1/2 ) (12) and (AB CD) = 2π 5/2 /[(α + β)(γ + δ)(α + β + γ + δ) 1/2 ] (13) αβ [ exp α+β R A R B 2 ] [ γδ γ+δ R C R D 2 ] (14) F 0 [(α + β)(γ + δ)/(α + β + γ + δ) R P R Q 2 ] (15) with R P = (α R A + β R B )/(α + β) (16) and R Q = (γ R C + δ R D )/(γ + δ). (17) These are all integrals required for the STO-3G calculation and the corresponding procedures in the program are functions S(A, B, RAB2), T(A, B, RAB2), V(A, B, RAB2, RCP2, ZC) and TWOE(A, B,C, D, RAB2, RCD2, RPQ2). From these integrals the overlap matrix S, the kinetic energy matrix T, the matrix of nuclear attraction energy to nucleus 1 (the He nucleus), V (1), and the matrix of nuclear attraction to nucleus 2 (the H nucleus), V (2), are calculated. Adding H core = T +V (1) +V (2) (18) the core Hamiltonian is computed. This Hamiltonian is the correct one for a single electron in the field of the nuclear point charges. The core Hamiltonian only needs to be evaluated once during the iterative calculations. In the next step, the overlap matrix S is diagonalized and the basis functions are orthonormalized. Here, canonical orthogonalization is used. The basis functions φ µ are transformed to 4
5 φ µ by means of an orthogonal transformation matrix X φ µ = ν X µν φ ν (19) For a canonical transformation X = U/ s (20) where U is the matrix that diagonalizes S, and s is a diagonal matrix with the eigenvalues of S. A matrix element of X is therefore X i j = U i j / s j. For the present problem we obtain φ 1 = φ φ 2 (21) φ 2 = φ φ 2 (22) Now the SCF procedure can be started. First, an initial guess for the density matrix P is required. Here P = 0 and F H core. This may be a poor choice and in more complicated situations a better guess may be needed (for example from a semiempirical calculation - this is what Gaussian will do for you). Next, the Fock matrix is transformed to the orthonormalized basis set F = X FX (23) Diagonalizing this matrix (solving F C = C ε) gives a unitary matrix of coefficients C and eigenvalues ε. From C the coefficients of the original basis functions φ µ are calculated via C = XC (24) Remember: eigenvalues remain unchanged under orthogonal transformations. In the present case the coefficient matrix is C = (25) 5
6 which indicates that the lowest molecular orbital ψ 1 is mainly (0.9291) composed of φ 1 with only a little contribution (0.1398) of φ 2. Without electron repulsion (remember: G = 0 was used up to here) the electrons cluster around the nucleus with the higher charge (He). From equation (25) the first real guess of the density matrix P can be formed N/2 P µν = 2 a C µa C νa (26) which has elements P = (27) here. From P we now form G which contains the two-electron part of the Fock matrix. The explicit expression for a matrix element of G is G µν = P λσ [(µν σλ) 1 (µλ σν)] (28) 2 λσ For the present situation with two basis functions φ 1 and φ 2 there are 2 4 = 16 possible integrals (µν λσ). There are, however, only 6 unique ones here. (φ 1 φ 1 φ 1 φ 1 ) = E h (29) (φ 2 φ 1 φ 1 φ 1 ) = E h (30) (φ 2 φ 1 φ 2 φ 1 ) = E h (31) (φ 2 φ 2 φ 1 φ 1 ) = E h (32) (φ 2 φ 2 φ 2 φ 1 ) = E h (33) (φ 2 φ 2 φ 2 φ 2 ) = E h (34) (35) The one-center integrals (φ 1 φ 1 φ 1 φ 1 ) and (φ 2 φ 2 φ 2 φ 2 ) are the repulsions between an electron in 6
7 φ 1 and another electron in the same orbital φ 1 (ditto for φ 2 ). The two-center integral (φ 2 φ 2 φ 1 φ 1 ) is the repulsion between an electron in φ 1 and one in φ 2. It has an asymptotic value of 1/R for large R. The other three integrals do not have simple classical interpretations. With these integrals the matrix elements G µν in eq. (28) can be evaluated and a new Fock matrix can be formed F = H core + G (36) Now the eigenvalue problem is solved with this new Fock matrix F, and a new guess for C (eq. (24)) and P (eq. (26)) is obtained. I.e. we return to equation 23 and iterate until convergence. Convergence is determined, e.g., with respect to changes in the density matrix P. I.e. if elements P µν in P do not change anymore by more than a pre-specified threshold, the calculation is called converged. The final wavefunctions and orbital energies for the present case are C = (37) and ε = (38) The lowest orbital, ψ 1, is a bonding orbital (same sign of both coefficients). It is mainly composed of the He function φ 1. The virtual orbital, ψ 2, is antibonding (opposite sign of the coefficients). Koopman s theorem allows us to quantify ionization potential and electron affinity: the IP = E h is large as it should be for a cationic species and EA = E h is positive which means that HeH + will bind an additional electron. This does, however, not mean that HeH is stable because the dissociation products (He + H + ) will bind an additional electron much more strongly than HeH +. A Mulliken population analysis yields the partial charges on the two nuclei. For this, one considers the diagonal elements of the matrix PS. The total number of electrons N can be related 7
8 to the molecular orbitals and the matrices P and S in the following way (the second equation follows from ψ a ( r) = K µ=1 C µaφ µ ( r)): N = 2 N/2 a d r ψ a ( r) 2 = µ ν P µν S νµ = µ (PS) µµ = tr(ps) (39) It should be noted that there is no unique definition of the number of electrons associated with a given atom or nucleus. The above analysis is the Mulliken population analysis and associates a charge q A = Z A (PS) µµ (40) µ A where Z A is the charge of nucleus A. Because the trace (tr, sum of diagonal matrix elements) obeys tr(ab) = tr(ba) one also has N = µ (S 1/2 PS 1/2 ) µµ = µ P µµ (41) where P is the density matrix in a symmetrically orthogonalized basis set φ µ( r) = ν (S 1/2 ) νµ φ ν ( r). (42) The diagonal elements associated with P yield the Löwdin population analysis q A = Z A (S 1/2 PS 1/2 ) µµ. (43) µ A 8
9 Short summary 1. Specify a molecule (nuclear coordinates R A ), atomic numbers Z A, the number of electrons N and an electronic basis set φ µ. 2. Calculate all required molecular integrals S µν, H core µν and (µν λσ). 3. Diagonalize overlap matrix S and obtain transformation matrix X (here canonical orthogonalization). 4. Obtain a guess for the density matrix P (here P = H core ; better choice may be required). 5. Calculate the matrix G from P and the two-electron integrals (µν λσ) via G µν = P λσ [(µν σλ) 1 (µλ σν)] (44) 2 λσ 6. Add G to the core Hamiltonian to obtain the Fock matrix F = H core + G 7. Calculate the transformed Fock matrix F = X FX 8. Diagonalize F to obtain C and ε. 9. Calculate C = XC. 10. Form a new density matrix P from C by using N/2 P µν = 2 a C µa C νa (45) 11. Determine whether procedure has converged (e.g. compare old P and new P). If not, return to step (5) with the new P. 12. If procedure has converged, use matrices C, P, F, etc. to determine expectation values and other properties of interest. 9
Chem 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 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 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 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 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 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 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 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 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 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 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 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 informationSemiempirical Molecular Orbital Theory
Semiempirical Molecular Orbital Theory Slater- determinantal wave func:ons and artree- Fock theory Video III.iv Constructing a 1-Electron Wave Function The units of the wave function are such that its
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 informationCHEM6085: Density Functional Theory Lecture 10
CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired
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 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 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 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 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 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 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 informationHartree-Fock. It can also be viewed as the lowest order term in perturbative expansion (linear in U) with respect to interaction between electrons.
Hartree-Fock It is probably the simplest method to treat the many-particle system. The dynamic many particle problem is replaced by an effective one-electron problem: electron is moving in an effective
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 informationChemistry 6 (9 am section) Spring Covalent Bonding
Chemistry 6 (9 am section) Spring 000 Covalent Bonding The stability of the bond in molecules such as H, O, N and F is associated with a sharing (equal) of the VALENCE ELECTRONS between the BONDED ATOMS.
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 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 informationMolecular Structure Both atoms and molecules are quantum systems
Molecular Structure Both atoms and molecules are quantum systems We need a method of describing molecules in a quantum mechanical way so that we can predict structure and properties The method we use is
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 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 informationCHEMISTRY 333 Fall 2006 Exam 3, in class
CHEMISTRY 333 Fall 006 Exam 3, in class Name SHOW YOUR WORK. Mathematica & calculators. You may use Mathematica or a calculator to do arithmetic. Do not use Mathematica for any higher-level manipulations
More informationPauli Deformation APPENDIX Y
APPENDIX Y Two molecules, when isolated say at infinite distance, are independent and the wave function of the total system might be taken as a product of the wave functions for the individual molecules.
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 informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
More information( ) + ( ) + ( ) = 0.00
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 32, April 14, 2006 (Some material in this lecture has been adapted from Cramer, C.
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 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 informationIntroduction to numerical projects
Introduction to numerical projects Here follows a brief recipe and recommendation on how to write a report for each project. Give a short description of the nature of the problem and the eventual numerical
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 informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
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 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 information( R) MOLECULES: VALENCE BOND DESCRIPTION VALENCE BOND APPROXIMATION METHOD. Example: H 2. at some fixed R, repeat for many different values
5.6 004 Lecture #9 page MOLECULES: VLENCE OND DESCRIPTION VLENCE OND PPROXIMTION METHOD Example: H H = + + r r r r r R Determine E ( R) at some fixed R, repeat for many different values Electronic Schrodinger
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 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 informationLecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 10 The Dirac equation WS2010/11: Introduction to Nuclear and Particle Physics The Dirac equation The Dirac equation is a relativistic quantum mechanical wave equation formulated by British physicist
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 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 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 informationMO theory is better for spectroscopy (Exited State Properties; Ionization)
CHEM 2060 Lecture 25: MO Theory L25-1 Molecular Orbital Theory (MO theory) VB theory treats bonds as electron pairs. o There is a real emphasis on this point (over-emphasis actually). VB theory is very
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 informationHartree-Fock-Roothan Self-Consistent Field Method
Hartree-Fock-Roothan Self-Consistent Field Method 1. Helium Here is a summary of the derivation of the Hartree-Fock equations presented in class. First consider the ground state of He and start with with
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 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 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 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 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 informationGustavus Adolphus College. Lab #5: Computational Chemistry
CHE 372 Gustavus Adolphus College Lab #5: Computational Chemistry Introduction In this investigation we will apply the techniques of computational chemistry to several of the molecular systems that we
More informationX-Ray transitions to low lying empty states
X-Ray Spectra: - continuous part of the spectrum is due to decelerated electrons - the maximum frequency (minimum wavelength) of the photons generated is determined by the maximum kinetic energy of the
More information362 Lecture 6 and 7. Spring 2017 Monday, Jan 30
362 Lecture 6 and 7 Spring 2017 Monday, Jan 30 Quantum Numbers n is the principal quantum number, indicates the size of the orbital, has all positive integer values of 1 to (infinity) l is the angular
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 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 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 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 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 informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationChemistry 2. Lecture 1 Quantum Mechanics in Chemistry
Chemistry 2 Lecture 1 Quantum Mechanics in Chemistry Your lecturers 8am Assoc. Prof Timothy Schmidt Room 315 timothy.schmidt@sydney.edu.au 93512781 12pm Assoc. Prof. Adam J Bridgeman Room 222 adam.bridgeman@sydney.edu.au
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 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 informationwbt Λ = 0, 1, 2, 3, Eq. (7.63)
7.2.2 Classification of Electronic States For all diatomic molecules the coupling approximation which best describes electronic states is analogous to the Russell- Saunders approximation in atoms The orbital
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 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 informationGeneral Physical Chemistry II
General Physical Chemistry II Lecture 10 Aleksey Kocherzhenko October 7, 2014" Last time " promotion" Promotion and hybridization" [He] 2s 2 2p x 1 2p y 1 2p z0 " 2 unpaired electrons" [He] 2s 1 2p x 1
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 information1 Quantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More informationAn Iterative Technique for Solving the N-electron Hamiltonian: The Hartree-Fock method
An Iterative Technique for Solving the N-electron Hamiltonian: The Hartree-Fock method Tony Hyun Kim Abstract The problem of electron motion in an arbitrary field of nuclei is an important quantum mechanical
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 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 informationCHAPTER 11 MOLECULAR ORBITAL THEORY
CHAPTER 11 MOLECULAR ORBITAL THEORY Molecular orbital theory is a conceptual extension of the orbital model, which was so successfully applied to atomic structure. As was once playfuly remarked, a molecue
More informationCHAPTER 10 Tight-Binding Model
CHAPTER 0 Tight-Binding Model Linear Combination of Atomic Orbitals (LCAO) Application to Bands from s-levels General Features of Tight-Binding Levels Wannier Functions 6 a : S S P 3S Core FE Semicore
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 informationLecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel
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 informationA Computer Study of Molecular Electronic Structure
A Computer Study of Molecular Electronic Structure The following exercises are designed to give you a brief introduction to some of the types of information that are now readily accessible from electronic
More informationwe have to deal simultaneously with the motion of the two heavy particles, the nuclei
157 Lecture 6 We now turn to the structure of molecules. Our first cases will be the e- quantum mechanics of the two simplest molecules, the hydrogen molecular ion, H +, a r A r B one electron molecule,
More informationECE440 Nanoelectronics. Lecture 07 Atomic Orbitals
ECE44 Nanoelectronics Lecture 7 Atomic Orbitals Atoms and atomic orbitals It is instructive to compare the simple model of a spherically symmetrical potential for r R V ( r) for r R and the simplest hydrogen
More informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
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 informationCHEMISTRY Topic #1: Atomic Structure and Nuclear Chemistry Fall 2017 Dr. Susan Findlay See Exercise 5.3
CHEMISTRY 1000 Topic #1: Atomic Structure and Nuclear Chemistry Fall 2017 Dr. Susan Findlay See Exercise 5.3 Periodic Trends and Effective Nuclear Charge Imagine four atoms/ions: One has a nucleus with
More informationd 2 dx 2. H = 2 2m p 0 dξ 2 H = 2m e ma2 H = 1 2m dξ 2
CHAPTER 14 The Origins of Chemical Bonding SECTION 14.1 14.1 A proton (of mass m p ) in a 1-D box of length L = 10a 0, the system in part (a), has a Hamiltonian we write in terms of the ordinary length
More informationMean-field concept. (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1
Mean-field concept (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1 Static Hartree-Fock (HF) theory Fundamental puzzle: The
More informationCHAPTER 9 COVALENT BONDING: ORBITALS 323
APTER 9 OVALET BODIG: ORBITALS 323 2 3 2 2 2 3 3 2 2 3 2 3 O * * 2 o; most of the carbons are not in the same plane since a majority of carbon atoms exhibit a tetrahedral structure (19.5 bond angles).
More informationNucleons in Nuclei: Interactions, Geometry, Symmetries
Nucleons in Nuclei: Interactions, Geometry, Symmetries Jerzy DUDEK Department of Subatomic Research, CNRS/IN 2 P 3 and University of Strasbourg, F-67037 Strasbourg, FRANCE September 28, 2010 Mathematial
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 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 information26 Group Theory Basics
26 Group Theory Basics 1. Reference: Group Theory and Quantum Mechanics by Michael Tinkham. 2. We said earlier that we will go looking for the set of operators that commute with the molecular Hamiltonian.
More informationExercises Symmetries in Particle Physics
Exercises Symmetries in Particle Physics 1. A particle is moving in an external field. Which components of the momentum p and the angular momentum L are conserved? a) Field of an infinite homogeneous plane.
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 information