MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
|
|
- Madeleine Rich
- 6 years ago
- Views:
Transcription
1 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 (TISE) for the molecule. For a diatomic molecule AB having N electrons, the Hamiltonian for the molecule may be written as H = T n + T e + V ne + V nn + V ee where T n = - [ ( 2 2 /2m A ) A + ( 2 /2m B ) 2 B ] T e = - ( 2 N 2 /2m e ) i=1 i V ne = - [(Z A e 2 N /4 0 ) i=1 (1/r Ai ) + (Z B e 2 N /4 0 ) i=1 (1/r Bi )] V nn = (Z A Z B e 2 /4 0 ) (1/r AB ) V ee = (e 2 N N /4 0 ) i=1 j>i (1/r ij ) nuclear kinetic energy electron kinetic energy nuclear-electron potential energy nuclear-nuclear potential energy electron-electron potential energy The above can be extended in a straightforward manner for use in polyatomic molecules. While writing the Hamiltonian operator is simple, finding solutions to the corresponding TISE is not. Some simplification can be obtained by making the Born-Oppenheimer approximation. This approximation is based on the observation that electrons are much less massive than nuclei. As a result, one can, to a very good first approximation, assume that the nuclei are stationary relative to the electrons, which are, due to their light mass, able to adjust to changes in nuclear geometry "instantaneously". As a consequence of this, the nuclear kinetic energy term (T n ) can be removed, while the nuclear-nuclear potential energy term becomes a constant whose value depends on the geometry chosen for the molecule (which, for a diatomic molecule, means the bond length). The more simple TISE obtained when making the Born-Oppenheimer approximation is still too difficult to solve exactly (except for a single electron system, such as in the H 2 + ion). However, there are a number of approximation methods that have been developed over the past decades that make it possible to find solutions where, within the approximations made, the error is small. Molecular orbital theory The first approximation usually made is the molecular orbital approximation. Each electron in the system is assumed to occupy a separate molecular orbital. The wavefunction for the molecule is then the product of the individual one electron orbital functions. Satisfaction of the Pauli principle (that the overall wavefunction be antisymmetric with respect to exchange of any two electrons) is obtained by writing the solution in terms of a matrix (the Slater determinant). The starting point in finding the best set of orbitals for a molecule is the Hartree-Fock (HF) approach. Molecular orbitals are found by numerically solving the Hartree-Fock equation, a relationship that is derived from our initial TISE. Solutions are often obtained in terms of a set of functions called Slater type orbitals (STOs). Since the numerical calculations required when using STOs are very time consuming, the STOs are themselves usually approximated by a set of Gaussian functions. This substitution leads to great improvements in the speed with which the required numerical calculations can be carried out. The quality of a particular molecular orbital calculation generally depends on two factors: (1) Basis set - The basis set is the set of initial functions from which the final molecular orbitals are constructed. In principle the best basis set would be a complete set, which would make it possible to (within the
2 Born-Oppenheimer and orbital approximations) find an exact solution to the TISE for a molecule. However, the time required to carry out the calculations involved scales approximately as the square of the number of basis functions being used, so as a practical matter a complete basis set cannot be used. There are two common basis sets used in molecular orbital calculations: Minimal basis set. Uses all of the core and valence atomic orbitals for each atom in the molecule. Augmented basis set. Uses all of the orbitals in the minimal basis set, along with additional unoccupied atomic orbitals, called virtual orbitals. Using more virtual orbitals improves the quality of the calculation but increases the time required to carry out the calculation. As an example, a minimal basis set for a molecular orbital calculation for the HF molecule would use a 1s atomic orbital for hydrogen, and the 1s, 2s, 2p x, 2p y, and 2p z atomic orbitals for fluorine (usually represented as STOs). This minimal basis set could be augmented by including the 2p z atomic orbital for the hydrogen atom (the p orbital that points towards the fluorine atom in the HF molecule), or by including virtual atomic orbitals from both the hydrogen and the fluorine atom. The basis set used in a particular molecular orbital calculation is usually indicated by the label used to identify the calculation. For example, minimal basis sets are often written as STO-nG, where STO indicates that Slater type orbitals are being used in the calculation, and n is the number of Gaussian functions being used to approximate the Slater type orbitals. So, for example, the STO-3G basis set uses Slater type orbitals, each of which has been approximated as a combination of three Gaussian functions. More complicated basis sets will be indicated with different types of symbols (which are often difficult to decipher for those who do not carry out molecular orbital calculations on a regular basis). (2) Corrections to the Hartree-Fock method. The Hartree-Fock method makes some approximations to find the solutions to the TISE for a molecule. There are a variety of methods that have been developed to go beyond the basic Hartree-Fock approach to improve the results that are obtained. Most corrections concentrate on accounting for the correlation of motion of the electrons in a molecule. Correlation refers to the fact that the presence of an electron at a particular point in space affects the locations of the other electrons in the molecule, due both to electron-electron repulsion and spin effects related to satisfying the Pauli principle. A common method for correcting for the effects of electron correlation is called configuration interaction (CI). Configuraton interaction calculations make use of virtual orbitals. In full-ci all possible assignments of electrons to virtual orbitals are made. However, because of the time involved for such calculations, CI is often restricted to configurations where only a few electrons are allowed to occupy virtual orbitals. For example, CI-SD restricts the CI calculations to states with one (single) or two (double) electron excitations. A second method for correcting for electron correlation is called the coupled cluster method (CC). This method is based on the theory of many body interactions in multiparticle systems. As is the case for the CI method, different levels of approximation are indicated based on the number of excitations that have been included. So, for example, CCSD includes all terms with one or two electrons in excited states. An additional correction procedure, called the Moller-Plesset method, makes use of perturbation theory to improve the quality of the molecular orbital calculation. One starts with the solution to the Hartree-Fock equation, and then treats corrections to this solution as small perturbations, which are addressed using standard techniques. The label used for this type of correction indicates the order to which corrections have been made. For example, in MP2, corrections up to second order in the energy have been carried out, while in MP3 corrections up to third order have been made. Additional methods for improving the quality of a molecular orbital calculation have also been developed. As is the case with basis sets, abbreviations to represent the types of corrections being included are usually indicated when a particular calculation is carried out. 2
3 Density functional theory As an alternative to the molecular orbital approach, a second method, called density functional theory (DFT), has been developed for theoretical calculations on molecules. Instead of attempting to solve the TISE for a molecule, the density functional approach attempts to find the probability density for electrons in the molecule, (x,y,z). In 1964, Hohenberg and Kohn proved that all of the properties of a molecule can be calculated from this probability density. Unfortunately, there is no method for finding an exact expression for (x,y,z) for a molecule. However, as is the case for molecular orbital theory, there are a variety of approximation methods that have been developed to implement the density functional approach. The quality of density functional calculations for molecules is now in many respects comparable to that found using molecular orbital theory, although there are still difficulties in modeling molecular dissociation and transition states. Density functional theory is generally the best method to use in solid state calculations. The experiment In the present experiment you will carry out several different molecular orbital calculations on the diatomic molecule you were assigned. The three calculations you will do are as follows: hf/6-31g* - A Hartree-Fock calculation using an augmented minimal basis set. The STOs are approximated using six Gaussian functions. The symbol 6-31 indicates use of a "split" basis set, which gives more flexibility to the calculation, while the symbol * indicates the use of additional polarized functions, one way of augmenting a minimal basis set. mp2/6-31g* - The same as the first calculation, with the addition of Moller-Plesset correction terms out to second order. b3 yp/6-31g* - A density functional calculation using the same basis set of functions used in the Hartree- Fock calculation. The symbol b3 yp stands for the Becke three parameter approximation to the Lee-Yang-Parr exchange correlation functional, which is essentially a method for accounting for electron correlation in the calculation. The calculations will be done using Gaussian-98, an older version of a commercial pckage of programs used to carry out molecular orbital calculations (the most recent version of this package is Gaussian-09), which is available on the computer in CP 375. To get to the screen where you will input the data for the calculation click FILE---New, and then enter the following information: % Section Leave blank Route Section Enter one of the following #hf/6-31g* opt freq #mp2/6-31g* opt freq #b3 yp/6-31g* opt freq for the Hartree-Fock calculation for the Moller-Plesset calculation for the density functional calculation (Note - Opt freq stands for "optimal frequency". This instruction tells the program to find the minimum in the potential energy for the molecule, and to then find the vibrational constant e.) Title section You need to put something here or the program will not run. So, for example, you could put your name and the kind of calculation being done (Joens-HF-O2, for example). 3
4 Charge and multiplicity For a calculation on a molecule the charge will be zero, while for an ion it will be the ion charge. The multiplicity is equal to 2S+1, where S is the total spin quantum number for the ground electronic state of the molecule. For example, for O 2, S = 1, and so the multiplicity is 3. Enter this information as the charge, a space, and the multiplicity (0 3 would be the entry for a calculation for the O 2 molecule). Molecule Specification This enters the geometric information used as a starting point in the calculation, and is very simple for a diatomic molecule. Enter the data in the same way as shown below for a calculation on the O 2 molecule. O O 1 R R=1.30 This means the following O - identifies the first atom in the diatomic molecule, using the symbol for the atom (in this case, an oxygen atom) O 1 R - identifies the second atom in the diatomic (in this case, O), and tells the program to look for a starting value for R, the bond distance. blank line R= The starting value for R to be used by the program (in Å, angstroms; 1 nm = 10 Å). The results of the calculation are not sensitive to the starting value used for R, as long as a reasonable value is used. For your trial value for R you can use the experimental value for r e from the literature search, rounded to two significant figures and converted to Å. The program will vary the geometry until it find r e, the equilibrium bond distance. After entering the data, click RUN (upper right side of the screen). The program will ask for a file name. I would suggest you use your name and an abbreviation for the type of calculation (joenshf-o2, for example). There are two pieces of information that you need to obtain from the calculation: re (equilibrium bond length) The value for r e will appear in a table labeled "optimized parameters" found approximately in the middle of the output file, and again near the bottom of the output file. It will have the following appearance. Optimization completed. -- Stationary point found ! Optimized Parameters!! (Angstroms and Degrees)! ! Name Definition Value Derivative Info.! ! R1 R(1,2) DE/DX = 0.! In the above example the value of r e is Å. 4
5 e (vibrational constant) The value for e will be given towards the bottom of the output file, and will have the following appearance. Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering activities (A**4/AMU), Raman depolarization ratios, reduced masses (AMU), force constants (mdyne/a) and normal coordinates: 1 SGG Frequencies Red. masses Frc consts IR Inten Raman Activ Depolar Atom AN X Y Z In the above example the value of e is cm -1. Note that there is a large amount of additional information given in the output file, including the distribution of electrons, the rotational constant, the dipole moment, dissociation constant, an estimate of the infrared and Raman intensity, and thermodynamic information calculated from the molecular parameters. You should make a table of results that includes the values for r e and e for each calculation. Your values for r e and e should be compared to values from the literature, which can be obtained from the NIST website. Briefly discuss your results. Focus on which method seems to best agree with experiment. 5
Electron 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 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 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 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 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 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 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
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 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 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 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 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 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 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 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 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
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 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 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 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 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 informationIFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations
Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Introduction: The objective of the second computer laboratory exercise is to get acquainted with a program for performing quantum chemical
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 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 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 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 informationPerformance 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 informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
More informationGaussian: Basic Tutorial
Input file: # hf sto-g pop=full Water - Single Point Energy 0 H.0 H.0 H 04.5 Route Section Start with # Contains the keywords Gaussian: Basic Tutorial Route Section Title Section Charge-Multiplicity Molecule
More information221B Lecture Notes Many-Body Problems II Molecular Physics
1 Molecules 221B Lecture Notes Many-Body Problems II Molecular Physics In this lecture note, we discuss molecules. I cannot go into much details given I myself am not familiar enough with chemistry. But
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 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 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 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 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 informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More 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 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 informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
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 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 informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More information2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.
.~, ~ I, sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.)1e'" 21t-ol Je C'...-------- lj-vi, J? Jr Jr \Ji 2~:J~ -ryej- r- 2 Jr A - f3 c _,~,= ~,.,w._..._.
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
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 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 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 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 informationChemistry 4681 Module: Electronic Structure of Small Molecules Background Handout
Chemistry 4681 Module: Electronic Structure of Small Molecules Background Handout C. David Sherrill Last Revised: 6 January 2000 1 Computational Chemistry The term computational chemistry is used to mean
More 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 informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
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 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 informationQuantum Chemistry. NC State University. Lecture 5. The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy
Quantum Chemistry Lecture 5 The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy NC State University 3.5 Selective absorption and emission by atmospheric gases (source:
More informationQUANTUM CHEMISTRY PROJECT 3: PARTS B AND C
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 6, 2017 QUANTUM CHEMISTRY PROJECT 3: PARTS B AND C PART B: POTENTIAL CURVE, SPECTROSCOPIC CONSTANTS, AND DISSOCIATION ENERGY OF DIATOMIC HYDROGEN (20
More informationNMR and IR spectra & vibrational analysis
Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent
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 informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationChemistry 416 Spectroscopy Fall Semester 1993 Dr. Rainer Glaser
Chemistry 416 Spectroscopy Fall Semester 1993 Dr. Rainer Glaser Second 1-Hour Examination IR/Raman Spectroscopy Monday, October 4, 1993, 8:40-9:30 Name: Question 1 35 Question 2 40 Question 3 25 Total
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 informationNPTEL/IITM. Molecular Spectroscopy Lectures 1 & 2. Prof.K. Mangala Sunder Page 1 of 15. Topics. Part I : Introductory concepts Topics
Molecular Spectroscopy Lectures 1 & 2 Part I : Introductory concepts Topics Why spectroscopy? Introduction to electromagnetic radiation Interaction of radiation with matter What are spectra? Beer-Lambert
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 SPECTROSCOPY
MOLECULAR SPECTROSCOPY First Edition Jeanne L. McHale University of Idaho PRENTICE HALL, Upper Saddle River, New Jersey 07458 CONTENTS PREFACE xiii 1 INTRODUCTION AND REVIEW 1 1.1 Historical Perspective
More informationThe Potential Energy Surface
The Potential Energy Surface In this section we will explore the information that can be obtained by solving the Schrödinger equation for a molecule, or series of molecules. Of course, the accuracy of
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 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 informationFYS-6306 QUANTUM THEORY OF MOLECULES AND NANOSTRUCTURES
i FYS-6306 QUANTUM THEORY OF MOLECULES AND NANOSTRUCTURES Credit units: 6 ECTS Lectures: 48 h Tapio Rantala, prof. Tue 10 12 SC203 SG219 8 10 SG312 FirstName.LastName@tut.fi http://www.tut.fi/~trantala/opetus/
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 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 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 informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationP. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION
P. W. Atkins and R. S. Friedman Molecular Quantum Mechanics THIRD EDITION Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1997 Introduction and orientation 1 Black-body radiation 1 Heat capacities 2 The
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 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 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 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 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 informationHartree, Hartree-Fock and post-hf methods
Hartree, Hartree-Fock and post-hf methods MSE697 fall 2015 Nicolas Onofrio School of Materials Engineering DLR 428 Purdue University nonofrio@purdue.edu 1 The curse of dimensionality Let s consider a multi
More informationChem 442 Review of Spectroscopy
Chem 44 Review of Spectroscopy General spectroscopy Wavelength (nm), frequency (s -1 ), wavenumber (cm -1 ) Frequency (s -1 ): n= c l Wavenumbers (cm -1 ): n =1 l Chart of photon energies and spectroscopies
More informationr R A 1 r R B + 1 ψ(r) = αψ A (r)+βψ B (r) (5) where we assume that ψ A and ψ B are ground states: ψ A (r) = π 1/2 e r R A ψ B (r) = π 1/2 e r R B.
Molecules Initial questions: What are the new aspects of molecules compared to atoms? What part of the electromagnetic spectrum can we probe? What can we learn from molecular spectra? How large a molecule
More informationPhysical Chemistry I Fall 2016 Second Hour Exam (100 points) Name:
Physical Chemistry I Fall 2016 Second Hour Exam (100 points) Name: (20 points) 1. Quantum calculations suggest that the molecule U 2 H 2 is planar and has symmetry D 2h. D 2h E C 2 (z) C 2 (y) C 2 (x)
More informationChemistry, Physics and the Born- Oppenheimer Approximation
Chemistry, Physics and the Born- Oppenheimer Approximation Hendrik J. Monkhorst Quantum Theory Project University of Florida Gainesville, FL 32611-8435 Outline 1. Structure of Matter 2. Shell Models 3.
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 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 information( ) electron gives S = 1/2 and L = l 1
Practice Modern Physics II, W018, Set 1 Question 1 Energy Level Diagram of Boron ion B + For neutral B, Z = 5 (A) Draw the fine-structure diagram of B + that includes all n = 3 states Label the states
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More informationeigenvalues eigenfunctions
Born-Oppenheimer Approximation Atoms and molecules consist of heavy nuclei and light electrons. Consider (for simplicity) a diatomic molecule (e.g. HCl). Clamp/freeze the nuclei in space, a distance r
More information( ) x10 8 m. The energy in a mole of 400 nm photons is calculated by: ' & sec( ) ( & % ) 6.022x10 23 photons' E = h! = hc & 6.
Introduction to Spectroscopy Spectroscopic techniques are widely used to detect molecules, to measure the concentration of a species in solution, and to determine molecular structure. For proteins, most
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 informationIntroduction to Computational Chemistry Computational (chemistry education) and/or. (Computational chemistry) education
Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools to help increase student understanding of material
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 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 information6.2 Polyatomic Molecules
6.2 Polyatomic Molecules 6.2.1 Group Vibrations An N-atom molecule has 3N - 5 normal modes of vibrations if it is linear and 3N 6 if it is non-linear. Lissajous motion A polyatomic molecule undergoes a
More informationIntroduction to Hartree-Fock calculations in Spartan
EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationChemistry 416 Spectroscopy Fall Semester 1997 Dr. Rainer Glaser
Chemistry 416 Spectroscopy Fall Semester 1997 Dr. Rainer Glaser Third 1-Hour Examination Vibrational Spectroscopy Monday, November 24, 1997, 8:40-9:30 Name: Answer Key Question 1 (Combination) 25 Question
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 information