Beyond the Hartree-Fock Approximation: Configuration Interaction
|
|
- Job Gibbs
- 6 years ago
- Views:
Transcription
1 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 a given basis set, the method provides the (variationally) lowest energy possible (as the expectation value of the Hamiltonian) and hence the best wavefunction possible. The method provides a fairly good description of the ground state of most molecules near the equilibrium geometry, and the Hartree-Fock method forms the basis for most electronic structure calculations in the literature up to 1995 or so. Nevertheless, no matter how large and flexible the basis set is, the HF method can only provide an approximation to the exact (true) wavefunction -- the exact wavefunction cannot be described as a single determinant (i.e. single electronic configuration). A closed shell (all electrons paired) singlet configuration can always be written as a single determinant. Open shell singlets require at least a two determinant description. Non-singlets (doublets, triplets, etc...) require special methods. The unrestricted HF (UHF) method uses one determinant for the up electrons and one (different) determinant for the down electrons. That effectively doubles the computation time relative to a closed shell calculation but does in many cases provide a decent description of the unpaired spin distribution. Restricted - open shell HF (ROHF) pairs as many electrons as possible in doubly occupied orbitals and then couples the remaining unpaired electrons to this closed shell using specialized techniques. A Hartree-Fock calculation with the largest, most flexible basis set imaginable produces an energy, which is called the Hartree-Fock limit. This energy is always higher (less negative) than the exact, non-relativistic energy for the system. The difference between the energy at the Hartree-Fock Limit, and the exact, nonrelativistic energy is called the Correlation Energy: E(Correlation) = E(exact, nonrelativistic) - E(Hartree-Fock Limit). The main deficiency in the Hartree-Fock method is that it does not account properly for the "correlated" motion of the electrons. We talk about "dynamical" vs. "non-dynamical" correlation. The latter arises from near-degeneracies (or real
2 degeneracies) of states in the system more than one electronic configuration of similar energy (transition metals, for example). Dynamical correlation refers more directly to the energy associated with the motion of the electrons being correlated and not independent of each other. Correlation of the motions of electrons with parallel spin is partially enforced by the Pauli principle and accounted for via the determinantal description of the wavefunction, but there is no correlation of the motion of electrons with opposite spin in the Hartree-Fock wavefunction. The Hartree-Fock method, based on the variational principle, produces a wavefunction which is too "averaged" over the positions of the electrons --- not sufficient account is made of the "instantaneous" interactions among the electrons. In the HF method, the electrons get "too close" to each other; <r 12 > is too small, i.e. <1/r 12 > is too large. There is too much e-e repulsion in a HF wavefunction. The absence of correlation among electrons of opposite spin lead to a number of undesirable features of the Hartree-Fock approximation. Chief among them is the inability of most closed-shell Hartree-Fock wavefunctions to properly describe dissociation and hence potential energy curves at longer bond lengths. There are problems with even the simplest of all molecules, H 2. The ground state of H 2 has a singlet wavefunction (symmetric space part, antisymmetric spin part) Ψ = φ 1 (1)φ 1 (2)[α(1)β(2) - β(1)α(2)](2) -1/2 with φ 1 (1) = (2+2S AB ) -1/2 {1s A (1) + 1s B (1)} and φ 1 (2) = (2+2S AB ) -1/2 {1s A (2) + 1s B (2)} within a minimal basis set framework. Ignoring all square roots and the spin functions, we find (by multiplying out the terms) that the spatial part of the wavefunction takes the form Ψ ~ φ 1 (1)φ 1 (2) = {1s A (1) + 1s B (1)}{1s A (2) + 1s B (2)} = {1s A (1)1s A (2) + 1s A (1)1s B (2) + 1s B (1)1s A (2) +1s B (1)1s B (2)} The first term in this expression has both electrons on atom A (corresponding to the electronic configuration H A - ); the last term has both electrons on atom B 2
3 - (corresponding to H B ); the two central terms have one electron on A and one on B. The two central terms are covalent -- but both the first and the last terms are ionic. Thus, this wavefunction has 50% covalent and 50% ionic character -- doesn't seem right for a "covalently" bound molecule like H 2. Furthermore, at large internuclear separations, the wavefunction for H 2 must represent two neutral hydrogen atoms and hence have the "purely covalent" form Ψ cov = {1s A (1)1s B (2) + 1s B (1)1s A (2)}[α(1)β(2) - β(1)α(2)](2) -1/2 which cannot be written as a single determinant! So the basic HF wavefunction (single determinant) can never be correct and good at long distances, so it will always have problems describing bond breaking processes correctly. Although inclusion of some ionic character to the H 2 ground state wavefunction does lower the energy relative to a purely covalent description, a 50% admixture is way too much and contributes to the Hartree-Fock energy always being too high (and binding energies consequently too low). In most cases, the equilibrium geometries of molecules are produced well by Hartree-Fock theory, but dissociation energies are underestimated. Energies of transition state structures, which often involve partially broken bonds and unusual bonding arrangements, are typically overestimated by Hartree-Fock theory leading to overestimation of activation energy barriers. Vibrational frequencies, in particular bond stretches, are mostly overestimated by Hartree-Fock theory, often by 10-15%. Configuration Interaction Improvements to the theory may be achieved from the use of a linear combination of Slater determinants, each Slater determinant describing an electronic configuration different from the (ground state) Hartree-Fock determinant. 3
4 Let us assume we have carried out a Hartree-Fock calculation for the electronic ground state of some closed-shell (molecular) system with N basis functions and n electrons. We thus have the n/2 molecular orbitals of lowest energy doubly occupied and (N - n/2) molecular orbitals (the virtuals) unoccupied. In the following, the occupied orbitals will be denoted by indices i,j,k.. etc., and the virtuals with indices a,b,c.. etc. We construct new determinants by replacing occupied orbitals with (formerly) unoccupied orbitals, corresponding to configurations which are electronically excited relative to the ground state. A single substitution, i -> a, produces a singly excited configuration state function, Ψ a i ; two substitutions, i -> a and j -> b, produce a doubly excited configuration, Ψ ab ij, and so on. We now write an improved wavefunctions as Ψ = Σ s' a s' Ψ s' = a 0 Ψ 0 + Σ s a s Ψ s = a 0 Ψ HF + Σa i a Ψi a + Σaij ab Ψij ab +... where Ψ 0 =Ψ HF is the original ground state Hartree-Fock wavefunction and Ψ s (s = 1,2,3...) is a Slater determinant for an excited configuration. The unknown coefficients, a s, are determined by the linear variation principle from the solution of a set of secular equations Σ s (Hst - E i δ st )a s = 0 t = 0,1,2,3... and hence a secular determinant H st - E i δ st = 0 The various configurations are mutually orthogonal (<Ψ 0 Ψ i a > = <Ψ0 Ψ ij ab > = <Ψ i a Ψjk bc > = 0; <Ψi a Ψj b > = δij δ ab and <Ψ ij ab Ψkl cd > = δik δ jl δ ac δ bd and so on), 4
5 hence the presence of the delta functions in the secular equations. H st is an "interaction" matrix element between two configurations H st = <Ψ s H Ψ t > Since the Hamiltonian includes only one- and two-electron operators, all matrix elements H st involving two configurations Ψ s and Ψ t which differ by three molecular orbitals or more is automatically zero. The lowest root of the determinant provides the (improved) energy of the electronic ground state; the second lowest root is an estimate of the energy of the lowest electronically excited state; and so on. Substitution of a particular root back into the secular equations provides the appropriate set of determinantal coefficients for the state in question. The improved wavefunction is a linear combination of configuration state functions which have been allowed to interact with each other, a configuration interaction wavefunction. If we consider the form of the configuration interaction Hamiltonian matrix, we find that along the diagonal we have the energies of the pure configurations, E 0 = E HF = <Ψ 0 H Ψ 0 >, E a i = <Ψi a H Ψi a >, Eij ab = <Ψij ab H Ψij ab >, and so on. The off-diagonal elements are the interaction elements between the ground state configuration, Ψ 0, and the singly excited configurations, Ψ a i ; between the ground state configuration, Ψ 0, and the doubly excited configurations, Ψ ab ij ; between singly, Ψ a i, and singly Ψj b, excited configurations; between singly and doubly excited configurations; between doubly and doubly excited configurations, and so on. The number of possible configurations increases rapidly as the degree of electron excitation (substitution) increases. In practice, it is customary to truncate the extent of substitution and limit the calculation to explicitly include only configurations singly and doubly excited relative to the ground state determinant (CISD, QCISD, or CCSD). Sometimes, estimates are made regarding what triple excitations might do (CISD(T); QCISD(T); CCSD(T)) 5
6 [Example: Methanol in a 6-31G(d) basis set has 38 bsf and 18 electrons. Distributing 18 electrons over 38 orbitals, corresponding to full CI, can be done in approximately different ways!] Inclusion of only singly excited configurations leads to the CIS (configuration interaction, singles) method (also called SECI (singly excited CI) or SCI). However, since H 0s = <Ψ 0 H Ψ a i > = <ΨHF H Ψ a i > = 0 always (a result known as "Brillouin's Theorem"), there is no mixing of the singly excited states with the (closed shell) ground state and hence no lowering of the ground state energy from this type of configuration interaction. H st = <Ψ s H Ψ t > = <Ψ a i H Ψj b > is generally not zero, so this type of calculation provides energies and wavefunctions of the electronically excited states (singlets and/or triplets), but it does not improve the ground state wavefunction or energy beyond its Hartree-Fock description. All triplet states are of course orthogonal to the closed shell ground state due to the difference in spin multiplicities, so configuration interaction calculations can be carried out for one spin multiplicity at a time. Inclusion of only doubly excited configurations leads to the CID (configuration interaction, doubles) method (also called DCI). The doubly excited configurations have nonzero matrix elements with the ground state and with each other (i.e., H 0s = <Ψ 0 H Ψ ab ij > = <ΨHF H Ψ ab ij > and Hst = <Ψ s H Ψ t > = <Ψ ab ij H Ψkl cd > are generally not equal to zero), so there will be mixing of the configurations. In particular, there will be a mixing of the doubly excited states into the ground state. The inclusion of doubly substituted determinants in the CI leads to an improved ground state wavefunction and a considerable lowering of the ground state energy. When both singly and doubly excited configurations are included in the calculation, we have the CISD (configuration interaction, singles and doubles) method (also called SDCI). The singly excited configurations do interact with the doubly excited configurations (i.e., H st = <Ψ s H Ψ t > = <Ψ a i H Ψjk bc > are generally not equal to 0) and do hence influence, and lower, the energy of the ground state and improve the ground state wavefunction. It is an indirect effect and thus not 6
7 nearly as large as the effect of the doubles alone. A CISD calculation may easily include a million (or more) configurations and special techniques have been developed to extract the lowest eigenvalues and vectors for matrices of that size ("Davidson's method"). One of the largest configuration interaction calculations carried out to date ( date ~ year 2000) was a full configuration interaction calculation (i.e. all levels of substitutions possible were considered) on water in a double-zeta basis set, which included more than one billion (10 9 ) configurations. The number of configurations included in a calculation may be reduced via the use of a "window". For example, a "10up-20down" CISD calculation would include all single and double excitations involving the 10 lowest virtual and 20 highest (doubly) occupied orbitals. Such reductions in excitation space are popular in semi-empirical methods ("ZINDO") but should be approached with caution, however, since the calculated results may depend strongly on the size of the window, in particular when small windows are applied. It is customary, though, in many CI applications to ignore excitations from the core (frozen core approximation) and thus only include excitations of the valence electrons into the virtual orbitals. The most time-consuming step of a CISD calculation scales approximately as n 2 N 4 (n = number of electrons, N = number of basis functions). The "rule of thumb" is that configuration interaction calculations should not be performed unless the computational facilities permit the use of basis sets of at least double-zeta plus polarization function quality ---- triple-zeta plus diffuse functions on the "atoms of importance" are, of course, even better. [Radial correlation (in-out); angular correlation (left-right)] Thus, configuration interaction calculations on real molecules are always very demanding in terms of computer memory, cpu-speed, and disk space. Differential correlation energy contributions can strongly influence reaction and, in particular, activation energies. Also, relative stabilities of isomers may be strongly affected by electron correlation, when the isomers show very different bonding patterns (such as single vs. multiple bonds; lone pairs vs. single bonds). It is still relatively rare that geometries are optimized at the CISD level for a molecule 7
8 containing more than 3-4 atoms. Instead, the correlation energy calculations are carried out at the (fixed) molecular geometry obtained from a lower level calculation such as Hartree-Fock (good basis set). Such calculations are called "single-point" calculations. Pople developed a simple nomenclature (involving "slashes", /) for reporting the level of a calculation and the geometry used in the calculation. The format is Method2/Basis2//Method1/Basis1 which means " a single point calculation using Method2 and Basis2 at the geometry optimized using Method1 and Basis1". For example, CISD/ G**//HF/6-311G** implies performing one configuration interaction calculation including all singly and doubly excited determinants with the G** basis set (valence triple-zeta + polarization + diffuse function quality) on a geometry, that was obtained by geometry optimization at the Hartree-Fock level with the 6-311G** basis set (valence triple-zeta plus polarization quality). An estimate of the effect of triple substitutions is available through the QCISD(T) method ("Quadratic CI including Single and Double substitutions with the effect of Triple substitutions included non-iteratively (i.e., not included to variational convergence)). Other popular methods based on configuration interaction include the CCSD ("Coupled Cluster Single and Double substitutions") and CCSD(T) methods. CISD is comparable to CCSD and QCISD(T) is comparable to CCSD(T). CCSD(T) calculations with large basis sets (~ triple-zeta quality plus several sets of polarization functions) are generally considered to be "state of the art" at this moment --- next year may be different. A "simple" correction formula has been developed by Davidson ('Davidson correction", S. R. Langhoff and E. R. Davidson, Int. J. Quantum. Chem., 8, 61 (1974)) to account for the energetic effect of all higher (triple, quadruple,..) substitutions, E(correction) = (1-a 2 0 ) ECISD where E CISD is the correlation energy obtained at the CISD level (i.e. E CISD -E HF ) and a 0 is the coefficient of the Hartree-Fock configuration in the CISD wavefunction. 8
Performance of Hartree-Fock and Correlated Methods
Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular
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 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 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 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 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 informationMethods for Treating Electron Correlation CHEM 430
Methods for Treating Electron Correlation CHEM 430 Electron Correlation Energy in the Hartree-Fock approximation, each electron sees the average density of all of the other electrons two electrons cannot
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 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 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 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 informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
More informationConsequently, the exact eigenfunctions of the Hamiltonian are also eigenfunctions of the two spin operators
VI. SPIN-ADAPTED CONFIGURATIONS A. Preliminary Considerations We have described the spin of a single electron by the two spin functions α(ω) α and β(ω) β. In this Sect. we will discuss spin in more detail
More informationLec20 Fri 3mar17
564-17 Lec20 Fri 3mar17 [PDF]GAUSSIAN 09W TUTORIAL www.molcalx.com.cn/wp-content/uploads/2015/01/gaussian09w_tutorial.pdf by A Tomberg - Cited by 8 - Related articles GAUSSIAN 09W TUTORIAL. AN INTRODUCTION
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 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 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 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 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 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 informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
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 informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (Tensor Contraction Engine) Why CC is important?! Correlation effects are important!! CC is size-extensive theory: can be used to describe dissociation processes.! Higher-order
More informationQUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 1, 2017 QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE OUTLINE In this project, you will carry out quantum mechanical calculations of
More information4 Post-Hartree Fock Methods: MPn and Configuration Interaction
4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the
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 informationLUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G
Q1: Display all the MOs for N2 in your report and classify each one of them as bonding, antibonding or non-bonding, and say whether the symmetry of the orbital is σ or π. Sketch a molecular orbital diagram
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 information( ) dσ 1 dσ 2 + α * 2
Chemistry 36 Dr. Jean M. Standard Problem Set Solutions. The spin up and spin down eigenfunctions for each electron in a many-electron system are normalized and orthogonal as given by the relations, α
More information( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)
Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
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 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 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 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 informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State
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 informationThis is called a singlet or spin singlet, because the so called multiplicity, or number of possible orientations of the total spin, which is
9. Open shell systems The derivation of Hartree-Fock equations (Chapter 7) was done for a special case of a closed shell systems. Closed shell means that each MO is occupied by two electrons with the opposite
More informationMulticonfigurational Quantum Chemistry. Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden
Multiconfigurational Quantum Chemistry Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden April 20, 2009 1 The Slater determinant Using the spin-orbitals,
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
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 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 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 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 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 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 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 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 informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
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 informationIntroduction to Electronic Structure Theory
CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,
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 informationCOPYRIGHTED MATERIAL. Quantum Mechanics for Organic Chemistry &CHAPTER 1
&CHAPTER 1 Quantum Mechanics for Organic Chemistry Computational chemistry, as explored in this book, will be restricted to quantum mechanical descriptions of the molecules of interest. This should not
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 informationComputational Chemistry: Molecular Simulations with Chemical and Biological Applications. Prof. M. Meuwly, Dr. PA Cazade, Dr. M.-W.
Computational Chemistry: Molecular Simulations with Chemical and Biological Applications Prof. M. Meuwly, Dr. PA Cazade, Dr. M.-W. Lee Overview 1. Electronic Structure of Molecules 1.1 The Electronic Problem
More informationTHE QUANTUM CHEMISTRY OF OPEN-SHELL SPECIES
4 THE QUANTUM CHEMISTRY OF OPEN-SHELL SPECIES Anna I. Krylov Department of Chemistry, University of Southern California, Los Angeles, CA, United States INTRODUCTION AND OVERVIEW The majority of stable
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 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 informationPRACTICE PROBLEMS Give the electronic configurations and term symbols of the first excited electronic states of the atoms up to Z = 10.
PRACTICE PROBLEMS 2 1. Based on your knowledge of the first few hydrogenic eigenfunctions, deduce general formulas, in terms of n and l, for (i) the number of radial nodes in an atomic orbital (ii) the
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 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 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 informationRelativistic and correlated calculations on the ground, excited, and ionized states of iodine
Relativistic and correlated calculations on the ground, excited, and ionized states of iodine W. A. de Jong, L. Visscher, a) and W. C. Nieuwpoort Laboratory for Chemical Physics and Materials Science Centre,
More informationHighly accurate quantum-chemical calculations
1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationLec20 Wed 1mar17 update 3mar 10am
564-17 Lec20 Wed 1mar17 update 3mar 10am Figure 15.2 Shows that increasing the diversity of the basis set lowers The HF-SCF energy considerably, but comes nowhere near the exact experimental energy, regardless
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 informationMolecular Term Symbols
Molecular Term Symbols A molecular configuration is a specification of the occupied molecular orbitals in a molecule. For example, N : σ gσ uπ 4 uσ g A given configuration may have several different states
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 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 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 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 informationChemistry 433 Computational Chemistry Fall Semester 2002 Dr. Rainer Glaser
Chemistry 433 Computational Chemistry Fall Semester 2002 Dr. Rainer Glaser Second 1-Hour Examination Electron Correlation Methods Wednesday, November 6, 2002, 11:00-11:50 Name: Answer Key Question 1. Electron
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 informationIntroduction to Heisenberg model. Javier Junquera
Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics
More informationChapter 9. Molecular Geometry and Bonding Theories
Chapter 9. Molecular Geometry and Bonding Theories PART I Molecular Shapes Lewis structures give atomic connectivity: they tell us which atoms are physically connected to which atoms. The shape of a molecule
More informationFROM MOLECULAR STRUCTURE TO FUNCTION
CEM421 FROM MOLECULAR STRUCTURE TO FUNCTION PARMACOLOGY MEDICINE MATERIAL SCIENCE BIOCEMISTRY CEMICAL ENGINEERING ORGANIC CEMISTRY ORGANOMETALLIC CEMISTRY INORGANIC CEMISTRY MOLECULAR STRUCTURE STRUCTURE
More informationExchange Mechanisms. Erik Koch Institute for Advanced Simulation, Forschungszentrum Jülich. lecture notes:
Exchange Mechanisms Erik Koch Institute for Advanced Simulation, Forschungszentrum Jülich lecture notes: www.cond-mat.de/events/correl Magnetism is Quantum Mechanical QUANTUM MECHANICS THE KEY TO UNDERSTANDING
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 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 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 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 informationELECTRONIC STRUCTURE CALCULATIONS OF OPEN-SHELL MOLECULES. Lyudmila V. Slipchenko. A Dissertation Presented to the FACULTY OF THE GRADUATE SCHOOL
ELECTRONIC STRUCTURE CALCULATIONS OF OPEN-SHELL MOLECULES by Lyudmila V. Slipchenko A Dissertation Presented to the FACULTY OF THE GRADUATE SCHOOL UNIVERSITY OF SOUTHERN CALIFORNIA In Partial Fulfillment
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 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 informationPotential energy, from Coulomb's law. Potential is spherically symmetric. Therefore, solutions must have form
Lecture 6 Page 1 Atoms L6.P1 Review of hydrogen atom Heavy proton (put at the origin), charge e and much lighter electron, charge -e. Potential energy, from Coulomb's law Potential is spherically symmetric.
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 informationNature of the Chemical Bond with applications to catalysis, materials science, nanotechnology, surface science, bioinorganic chemistry, and energy
Lecture 12, October 21, 2016 Transition metals Heme-Fe Nature of the Chemical Bond with applications to catalysis, materials science, nanotechnology, surface science, bioinorganic chemistry, and energy
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 informationChemistry 483 Lecture Topics Fall 2009
Chemistry 483 Lecture Topics Fall 2009 Text PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon A. Background (M&S,Chapter 1) Blackbody Radiation Photoelectric effect DeBroglie Wavelength Atomic
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 information3/30/2015. Third energy level. Second energy level. Energy absorbed. First energy level. Atomic nucleus. Energy released (as light)
Chapter 2 An Introduction Chemistry Lecture 2: Energy Levels and Chemical Bonding Electrons are always moving Outside the nucleus in atomic orbitals Maybe usually Average distance from nucleus (size of
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 informationOther methods to consider electron correlation: Coupled-Cluster and Perturbation Theory
Other methods to consider electron correlation: Coupled-Cluster and Perturbation Theory Péter G. Szalay Eötvös Loránd University Institute of Chemistry H-1518 Budapest, P.O.Box 32, Hungary szalay@chem.elte.hu
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationValence bond theory accounts, at least qualitatively, for the stability of the covalent bond in terms of overlapping atomic orbitals.
Molecular Orbital Theory Valence bond theory accounts, at least qualitatively, for the stability of the covalent bond in terms of overlapping atomic orbitals. Using the concept of hybridization, valence
More information