An Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine

Size: px
Start display at page:

Download "An Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine"

Transcription

1 An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine

2 This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical problems?

3 What do quantum chemists do? Approximate the electronic structures of molecules Extract useful chemical information from them

4 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R.

5 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance

6 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

7 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

8 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

9 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

10 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

11 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy Distance

12 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy F = de dr Distance

13 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy F = de dr Distance

14 Potential energy surface (PES) defined The PES, E(R), is the energy of a molecule as a function of the locations of its nuclei, R. Distance Energy F = de dr Distance

15 Using the PES Many important chemical question can be rephrased in terms of the PES We rarely calculate the whole PES, but instead explore and describe it in a variety of ways The PES can be calculated in many ways which vary in Accuracy Computational cost The portions of the PES which they are capable of describing

16 The PES is approximate The PES idea comes from the approximate solution of the full molecular time-independent Schrodinger equation Hˆ Ψ ( R, r) = EΨ( R, r)

17 The PES is approximate The PES idea comes from the approximate solution of the full molecular time-independent Schrodinger equation Hˆ Ψ ( R, r) = EΨ( R, r) Nuclear coordinates Electronic coordinates

18 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1

19 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1 I and J index nuclei i and j index electrons

20 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1 m is represents mass r represents distance q represents charge

21 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1 Nuclear Kinetic Energy Electronic Kinetic Energy

22 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1 Nuclear Kinetic Energy Electronic Kinetic Energy Nuclear Repulsion Electron- Nuclear Attraction Electron- Electron Repulsion

23 The PES is approximate The molecular Hamiltonian is ˆ q q q H = ( R) ( r) I J I 2mI I 2 i I i I J rij I i rij i j rij 1 Atomic Units ħ = 1 e = 1 m e = 1 1/(4 πε ) = 1 0

24 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Ψ ( R, r) = χ ( R) ψ ( r; R) Nuclear Wavefunction Electronic Wavefunction (which depends parametrically on the nuclear coordinates)

25 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Ψ ( R, r) = χ ( R) ψ ( r; R) Fast Electrons Slow Nucleus

26 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Ψ ( R, r) = χ ( R) ψ ( r; R) Hˆ ψ ( r; R) = V ( R) ψ ( r; R) elec q q ˆ 1 2 I J I elec = 2 i ( r) + + i I J rij I j rij i j rij H q 1

27 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Ψ ( R, r) = χ ( R) ψ ( r; R) Hˆ elecψ ( r; R) = V ( R) ψ ( r; R) q q ˆ 1 2 I J I elec = 2 i ( r) + + i I J rij I j rij i j rij H V(R) is the (adiabatic) potential energy surface q 1

28 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Ψ ( R, r) = χ ( R) ψ ( r; R) Hˆ elecψ ( r; R) = V ( R) ψ ( r; R) q q ˆ 1 2 I J I elec = 2 i ( r) + + i I J rij I j rij i j rij H V(R) is the (adiabatic) potential energy surface This is the electronic structure problem q 1

29 The PES is approximate Separate the nuclear and electronic problems (Born-Oppenheimer Approx.) Hˆ ψ ( r; R) = V ( R) ψ ( r; R) Hˆ χ ( R) = Eχ ( R) nuc ˆ 1 2 ( ) ( ) nuc = 2m I I R + V R I H Ψ ( R, r) = χ ( R) ψ ( r; R) elec q q ˆ 1 2 I J I elec = 2 i ( r) + + i I J rij I j rij i j rij H q 1

30 The PES is approximate The Born-Oppenheimer Approximation (BOA) neglects couplings between electronic and nuclear motions These missing terms are called nonadiabatic or vibronic coupling terms The BOA is a good approximation when electronic states are not near each other in energy (ground states of stable molecules with simple electronic structures) Electronic states are near each other in energy sometimes: Some transition states (nonadiabatic reaction dynamics) Excited states/photochemistry Jahn-Teller systems The BOA is valid for the majority of problems you see in the literature and in this course

31 Adiabatic vs Diabatic PES Adiabatic PES are the solution to the electronic Schrodinger equation Diabatic PES are defined according to state character (e.g. covalent vs ionic)

32 Adiabatic vs Diabatic PES Adiabatic PES are the solution to the electronic Schrodinger equation Diabatic PES are defined according to state character (e.g. covalent vs ionic) Ionic Covalent Diabatic

33 Adiabatic vs Diabatic PES Adiabatic PES are the solution to the electronic Schrodinger equation Diabatic PES are defined according to state character (e.g. covalent vs ionic) Ionic Covalent Diabatic Adiabatic

34 Adiabatic vs Diabatic PES Adiabatic PES are the solution to the electronic Schrodinger equation Diabatic PES are defined according to state character (e.g. covalent vs ionic) Ionic S 1 Covalent S 0 Diabatic Adiabatic

35 Diabatic PES Diabats are conceptually convenient There is no unique way to define diabats Diabats can be very strongly coupled over a broad range of configuration space Ionic Covalent Diabatic

36 The PES is many-dimensional Each atom is a point in a 3-D space the PES is 3N dimensional (where N is the number of atoms) The molecular geometry (configuration, structure) can be represented many ways Cartesian coordinates Internal Coordinates (e.g. Z-matrix)

37 Cartesian Coordinates Also known as XYZ coordinates y O H H x Used in various common storage formats:.xyz,.pdb

38 Z-Matrix Coordinates Define bond lengths, angles, dihedrals

39 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O

40 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1,

41 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1,

42 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2,

43 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1,

44 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1, 104.5, 3,

45 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1, 104.5, 3,

46 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1, 104.5, 3,

47 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1, 104.5, 3, 0.0 0

48 Z-Matrix Coordinates Define bond lengths, angles, dihedrals O O, 1, H, 1, 0.905, 2, H, 2, 0.905, 1, 104.5, 3,

49 Z-Matrix Coordinates Advantages 3 Chemically intuitive Non-redundant (eliminates translation, rotation) Disadvantages For large or cyclic molecules the geometry can be very sensitive to small changes in the parameters Equations of motion/geometry optimization are more nature in Cartesians

50 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES

51 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES Distance Energy Distance

52 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES Distance Energy Equilibrium Distance Distance

53 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES Distance Energy What does this energy mean? Equilibrium Distance Distance

54 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES Distance Energy What does this energy mean? NOTHING! Equilibrium Distance Distance

55 Local Minima Correspond to Stable Structures At reasonable temperatures molecules tend to spend most of their time at low energy points on the PES Distance Energy What does this energy mean? NOTHING! Equilibrium Distance Distance Dissociation Energy

56 A Good Minimum is Hard to Find In many dimensional systems we use geometry optimization algorithms Iterative process Calculate gradient of the PES (AKA forces) Update geometry Check for convergence Repeat if necessary

57 A Good Minimum is Hard to Find Energy

58 A Good Minimum is Hard to Find Energy

59 A Good Minimum is Hard to Find Energy

60 A Good Minimum is Hard to Find Energy

61 A Good Minimum is Hard to Find Energy

62 A Good Minimum is Hard to Find Energy

63 A Good Minimum is Hard to Find Energy

64 A Good Minimum is Hard to Find Energy

65 A Good Minimum is Hard to Find Convergence criteria The gradient (force) is below a threshold The change in geometry from the previous step is below a threshold The change in energy from the previous step is below a threshold

66 A Good Minimum is Hard to Find A word on forces Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

67 A Good Minimum is Hard to Find A word on forces D Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

68 A Good Minimum is Hard to Find A word on forces D+ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

69 A Good Minimum is Hard to Find A word on forces D+ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate E(D+ΔD) (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

70 A Good Minimum is Hard to Find A word on forces D-ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate E(D+ΔD) (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

71 A Good Minimum is Hard to Find A word on forces D-ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate E(D+ΔD) E(D-ΔD) (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

72 A Good Minimum is Hard to Find A word on forces D-ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate E(D+ΔD) - E(D-ΔD) F = 2ΔD (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

73 A Good Minimum is Hard to Find A word on forces D-ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Central finite difference Actually differentiate E(D+ΔD) - E(D-ΔD) F = 2ΔD (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

74 A Good Minimum is Hard to Find A word on forces D-ΔD Forces are a 3N-dimensional vector (x, y, and z force on each atom) Analytical forces Actually differentiate E(D+ΔD) - E(D) F = 2ΔD (Computational) cost is similar to calculating energy Available for most simpler methods of calculating the PES Numerical forces Forward finite difference Compute derivative by finite difference Cost is at least 3N times the cost of calculating the energy Available for all methods of calculating the PES

75 A Good Minimum is Hard to Find Energy

76 A Good Minimum is Hard to Find Energy

77 Energy A Good Minimum is Hard to Find Global Minimum Local Minimum

78 Energy A Good Minimum is Hard to Find How do you know you ve found the global minimum?

79 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy

80 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy R O Ketone R' R OH Enol R'

81 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy R O Ketone R' R OH Enol R'

82 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy R R R' trans R' gauche

83 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy ΔE R R R' trans R' gauche

84 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy H C N C N H

85 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy H C N H C N Transition State (TS) C N H

86 A Good Minimum is Hard to Find How do you know you ve found the global minimum? Which minimum do you want? Energy E activation H C N H C N Transition State (TS) C N H

87 A Good Transition State is Really A TS is defined as the highest point along the lowest energy path between reactant and product is a stable point (F=0) has a single imaginary frequency (one normal mode in which the second derivative is negative) Hard to Find Energy H C N Transition State (TS)

88 A Good Transition State is Really A TS can be found using one of the TS optimization algorithms available in most quantum chemistry codes by a series of constrained optimizations with only a single reaction coordinate is constrained Hard to Find Energy H C N Transition State (TS)

89 A Good Transition State is Really Hard to Find C H N Energy (ev) H C N C N H H-C-N Angle (deg)

90 Calculating Vibrational Frequencies Energy harmonic approx. Distance true PES Calculate the second derivative matrix (Hessian): 2 d E 2 dr at the minimum Diagonalize mass weighted hessian to get normal modes Calculate frequencies according to the well known solution of the harmonic oscillator problem

91 Energy Differences and Errors Two types of errors Constant error an error that is the same for all points on the PES Non-parallelity error an error that varies from point to point of the PES

92 Energy Differences and Errors Constant Error Exact

93 Energy Differences and Errors Constant Error Exact Approx.

94 Energy Differences and Errors Constant Error Exact Energy Calculated Exact Approx. Distance

95 Energy Differences and Errors Constant Error Exact Energy Calculated Exact Approx. Distance

96 Energy Differences and Errors Constant Error Exact Energy Calculated Exact Approx. Distance

97 Energy Differences and Errors Constant Error Exact Energy Calculated Exact Approx. Distance Constant errors in energy usually don t matter

98 Energy Differences and Errors Non-parallelity Error Calculated Energy Exact Distance

99 Energy Differences and Errors Non-parallelity Error Calculated Energy Exact Wrong dissociation energy Distance

100 Energy Differences and Errors Non-parallelity Error Calculated Energy Exact Wrong dissociation energy Distance Wrong frequencies

101 Energy Differences and Errors Non-parallelity Error Calculated Energy Exact Wrong dissociation energy Distance Wrong frequencies Non-parallelity errors matter!

102 Energy v. Free Energy Free energies, not energies, determine the probability of being in a specific state at finite temperature The energy of a molecule at a specific geometry is an energy, NOT a free energy! Some additional work is required to get a free energy estimate from frequencies RRKM theory estimate reaction rates from frequencies at the reactant and transition state molecular dynamics

103 Energy v. Free Energy

104 Energy v. Free Energy Low Energy Low Energy

105 Energy v. Free Energy Wide well

106 Energy v. Free Energy Wide well Narrow well

107 Energy v. Free Energy Wide well Higher Entropy Narrow well Lower Entropy

108 Energy v. Free Energy Wide well Higher Entropy Lower Free Energy Narrow well Lower Entropy High Free Energy

109 Energy v. Free Energy Lower Frequency Wide well Higher Entropy Lower Free Energy Higher Frequency Narrow well Lower Entropy High Free Energy

110 The PES isn t everything Electronic properties can be calculated from the wavefunction Dipole, quadrupole moments Atomic charges Polarizabilities IR intensities NMR chemical shifts

111 Summary The PES is the energy as a function of the location of the nuclei Stationary points (F=0) on the PES are representative of stable structures and transition states Non-parallelity error in the PES is bad; error cancellation is your friend We can calculate electronic properties from the wavefunction

112 The Hierarchy of Methods Cost Accuracy

113 The Hierarchy of Methods Hartree-Fock (HF) Cost HF Accuracy

114 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Cost HF Accuracy

115 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Cost MP2 HF Accuracy

116 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Cost MP2 CI HF Accuracy

117 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MP2 CI CC HF Accuracy

118 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MC MP2 CI CC Multiconfigurational Methods (MC) HF Accuracy

119 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MC MP2 CI CC Multiconfigurational Methods (MC) HF Break bonds Excited States Open Shells Accuracy

120 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MC MP2 CI CC Multiconfigurational Methods (MC) HF DFT Density Functional Theory (DFT) Accuracy

121 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MC MP2 CI CC Multiconfigurational Methods (MC) HF DFT Density Functional Theory (DFT) SE Semiempirical (SE) Accuracy

122 The Hierarchy of Methods Hartree-Fock (HF) Post-HF Moller-Plesset Perturb. Theory (MP2) Configuration Interaction (CI) Coupled Cluster (CC) Cost MC MP2 CI CC Multiconfigurational Methods (MC) HF DFT Density Functional Theory (DFT) SE Semiempirical (SE) Accuracy

123 Coming Up Next Linear Algebra Review Superposition principle Variational Principle Basis Sets and the relationship between matrices and operators Hartree-Fock Approximation What is it? How is it implemented in the computer? When does it fail and why?

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University

AN 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 information

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )

MO 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 information

Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland

Example 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 information

Chemistry 4560/5560 Molecular Modeling Fall 2014

Chemistry 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 information

CHEM3023: Spins, Atoms and Molecules

CHEM3023: 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 information

Chemistry 334 Part 2: Computational Quantum Chemistry

Chemistry 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 information

VALENCE Hilary Term 2018

VALENCE Hilary Term 2018 VALENCE Hilary Term 2018 8 Lectures Prof M. Brouard Valence is the theory of the chemical bond Outline plan 1. The Born-Oppenheimer approximation 2. Bonding in H + 2 the LCAO approximation 3. Many electron

More information

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education

Session 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 information

Figure 1: Transition State, Saddle Point, Reaction Pathway

Figure 1: Transition State, Saddle Point, Reaction Pathway Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces

More information

Quantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,

Quantum 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 information

Conical Intersections. Spiridoula Matsika

Conical Intersections. Spiridoula Matsika Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)

More information

DFT calculations of NMR indirect spin spin coupling constants

DFT 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 information

CE 530 Molecular Simulation

CE 530 Molecular Simulation 1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy

More information

Hartree, Hartree-Fock and post-hf methods

Hartree, 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 information

CHEM 545 Theory and Practice of Molecular Electronic Structure. Anna I. Krylov. DON T PANIC.

CHEM 545 Theory and Practice of Molecular Electronic Structure. Anna I. Krylov.   DON T PANIC. CHEM 545 Theory and Practice of Molecular Electronic Structure Anna I. Krylov http://iopenshell.usc.edu/chem545/ DON T PANIC USC Fall 2014 Things to do: 1. Install IQmol (by this Thursday). http://iqmol.org/.

More information

23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation M I

23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation M I 23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation 1. Now we will write down the Hamiltonian for a molecular system comprising N nuclei and n electrons.

More information

Computational Methods. Chem 561

Computational 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 information

NPTEL/IITM. Molecular Spectroscopy Lectures 1 & 2. Prof.K. Mangala Sunder Page 1 of 15. Topics. Part I : Introductory concepts Topics

NPTEL/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 information

Lecture 08 Born Oppenheimer Approximation

Lecture 08 Born Oppenheimer Approximation Chemistry II: Introduction to Molecular Spectroscopy Prof. Mangala Sunder Department of Chemistry and Biochemistry Indian Institute of Technology, Madras Lecture 08 Born Oppenheimer Approximation Welcome

More information

Molecular energy levels

Molecular energy levels Molecular energy levels Hierarchy of motions and energies in molecules The different types of motion in a molecule (electronic, vibrational, rotational,: : :) take place on different time scales and are

More information

Electronic structure calculations: fundamentals George C. Schatz Northwestern University

Electronic 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 information

Born-Oppenheimer Approximation

Born-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

Electric properties of molecules

Electric 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 information

CHEM6085: Density Functional Theory Lecture 10

CHEM6085: 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

Practical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology

Practical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Practical Advice for Quantum Chemistry Computations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Choice of Basis Set STO-3G is too small 6-31G* or 6-31G** 6 probably

More information

Lecture 9. Hartree Fock Method and Koopman s Theorem

Lecture 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

Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur

Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Under the supervision of: Dr. Marcel Nooijen Associate Professor Department of Chemistry University of Waterloo

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack 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 information

Electron Correlation

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 information

A microscopic approach to nuclear dynamics. Cédric Simenel CEA/Saclay, France

A microscopic approach to nuclear dynamics. Cédric Simenel CEA/Saclay, France A microscopic approach to nuclear dynamics Cédric Simenel CEA/Saclay, France Introduction Quantum dynamics of complex systems (nuclei, molecules, BEC, atomic clusters...) Collectivity: from vibrations

More information

Diatomic Molecules. 7th May Hydrogen Molecule: Born-Oppenheimer Approximation

Diatomic Molecules. 7th May Hydrogen Molecule: Born-Oppenheimer Approximation Diatomic Molecules 7th May 2009 1 Hydrogen Molecule: Born-Oppenheimer Approximation In this discussion, we consider the formulation of the Schrodinger equation for diatomic molecules; this can be extended

More information

P. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION

P. 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 information

Yingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2

Yingwei 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 information

Session 7 Overview: Part A I. Prediction of Vibrational Frequencies (IR) Part B III. Prediction of Electronic Transitions (UV-Vis) IV.

Session 7 Overview: Part A I. Prediction of Vibrational Frequencies (IR) Part B III. Prediction of Electronic Transitions (UV-Vis) IV. Session 7 Overview: Part A I. Prediction of Vibrational Frequencies (IR) II. Thermochemistry Part B III. Prediction of Electronic Transitions (UV-Vis) IV. NMR Predictions 1 I. Prediction of Vibrational

More information

The Self Interaction Correction revisited

The Self Interaction Correction revisited The Self Interaction Correction revisited Explicit dynamics of clusters and molecules under irradiation Spectroscopic accuracy at low energy SIC problem : one electron interacts with its own mean-field!

More information

Quantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno

Quantum 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 information

Multi-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 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 information

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy.

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. What determines the UV-Vis (i.e., electronic transitions) band appearance? Usually described by HOMO LUMO electron jump LUMO

More information

Introduction to Computational Chemistry

Introduction to Computational Chemistry Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September

More information

Vibrational Spectra (IR and Raman) update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6

Vibrational Spectra (IR and Raman) update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6 Vibrational Spectra (IR and Raman)- 2010 update Tinoco has very little, p.576, Engel Ch. 18, House Ch. 6 Born-Oppenheimer approx. separate electron-nuclear Assume elect-nuclear motion separate, full wave

More information

Lecture 4: methods and terminology, part II

Lecture 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 information

σ u * 1s g - gerade u - ungerade * - antibonding σ g 1s

σ u * 1s g - gerade u - ungerade * - antibonding σ g 1s One of these two states is a repulsive (dissociative) state. Other excited states can be constructed using linear combinations of other orbitals. Some will be binding and others will be repulsive. Thus

More information

This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.

This 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 information

Condensed matter physics FKA091

Condensed matter physics FKA091 Condensed matter physics FKA091 Ermin Malic Department of Physics Chalmers University of Technology Henrik Johannesson Department of Physics University of Gothenburg Teaching assistants: Roland Jago &

More information

Teoría del Funcional de la Densidad (Density Functional Theory)

Teorí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 information

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES

More information

Performance of Hartree-Fock and Correlated Methods

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 information

The Potential Energy Surface

The 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 information

QUANTUM CHEMISTRY FOR TRANSITION METALS

QUANTUM 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 information

Chemistry 883 Computational Quantum Chemistry

Chemistry 883 Computational Quantum Chemistry Chemistry 883 Computational Quantum Chemistry Instructor Contact Information Professor Benjamin G. Levine levine@chemistry.msu.edu 215 Chemistry Building 517-353-1113 Office Hours Tuesday 9:00-11:00 am

More information

Density Functional Theory

Density 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 information

Methods for Treating Electron Correlation CHEM 430

Methods 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 information

Chemistry 2. Lecture 1 Quantum Mechanics in Chemistry

Chemistry 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 information

Introduction to Computational Chemistry

Introduction 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 information

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy.

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. What determines the UV-Vis (i.e., electronic transitions) band appearance? Usually described by HOMO LUMO electron jump LUMO

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack 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 information

Molecular Simulation I

Molecular 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 information

Principles of Molecular Spectroscopy

Principles of Molecular Spectroscopy Principles of Molecular Spectroscopy What variables do we need to characterize a molecule? Nuclear and electronic configurations: What is the structure of the molecule? What are the bond lengths? How strong

More information

Ab initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier

Ab 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 information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

OVERVIEW 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 information

Chemistry 483 Lecture Topics Fall 2009

Chemistry 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 information

Same idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)

Same 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 information

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.

Chem 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 information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 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

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 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 information

AN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS. Laren K. Mortensen

AN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS. Laren K. Mortensen AN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS by Laren K. Mortensen A senior thesis submitted to the faculty of Brigham Young University in partial fulfillment

More information

Chemistry, Physics and the Born- Oppenheimer Approximation

Chemistry, 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 information

/2Mα 2 α + V n (R)] χ (R) = E υ χ υ (R)

/2Mα 2 α + V n (R)] χ (R) = E υ χ υ (R) Spectroscopy: Engel Chapter 18 XIV 67 Vibrational Spectroscopy (Typically IR and Raman) Born-Oppenheimer approx. separate electron-nuclear Assume elect-nuclear motion separate, full wave fct. ψ (r,r) =

More information

Molecular orbitals, potential energy surfaces and symmetry

Molecular orbitals, potential energy surfaces and symmetry Molecular orbitals, potential energy surfaces and symmetry mathematical presentation of molecular symmetry group theory spectroscopy valence theory molecular orbitals Wave functions Hamiltonian: electronic,

More information

The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.

The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions

More information

(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β

(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 information

we have to deal simultaneously with the motion of the two heavy particles, the nuclei

we 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 information

Intro/Review of Quantum

Intro/Review of Quantum Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the

More information

Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal

Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5

More information

2~: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.

2~: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 information

Structure of diatomic molecules

Structure 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 information

Exp. 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 + 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 information

Computational Chemistry I

Computational 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 information

Oslo node. Highly accurate calculations benchmarking and extrapolations

Oslo 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 information

EXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle:

EXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle: EXAM INFORMATION Radial Distribution Function: P() r RDF() r Br R() r B is the normalization constant., p Operator: p ^ d dx Heisenberg Uncertainty Principle: n ax n! Integrals: xe dx n1 a x p Particle

More information

Theoretical Investigation of possible Chemiluminescent Analogs of 1,2-dioxetane

Theoretical Investigation of possible Chemiluminescent Analogs of 1,2-dioxetane Theoretical Investigation of possible Chemiluminescent Analogs of 1,2-dioxetane Bachelor Thesis in Chemistry Hannes Gustafsson Degree Project C in Chemistry, 15 credits Thesis Defence: Spring 2017 Supervisors:

More information

MOLECULAR STRUCTURE. The general molecular Schrödinger equation, apart from electron spin effects, is. nn ee en

MOLECULAR 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 information

On the Uniqueness of Molecular Orbitals and limitations of the MO-model.

On 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 information

Rotations and vibrations of polyatomic molecules

Rotations and vibrations of polyatomic molecules Rotations and vibrations of polyatomic molecules When the potential energy surface V( R 1, R 2,..., R N ) is known we can compute the energy levels of the molecule. These levels can be an effect of: Rotation

More information

Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations

Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker

More information

Introduction to Hartree-Fock Molecular Orbital Theory

Introduction 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 information

When are electrons fast and are Born and Oppenheimer correct?

When are electrons fast and are Born and Oppenheimer correct? When are electrons fast and are Born and Oppenheimer correct? Wim J. van der Zande Department of Molecular and Laser Physics University of Nijmegen Han-Sur-Lesse, Winter 2003 Contents of the Lectures 0.

More information

Chemistry 2000 Lecture 1: Introduction to the molecular orbital theory

Chemistry 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 information

Quantum Mechanical Simulations

Quantum 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 information

DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory

DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,

More information

Chemistry 881 Lecture Topics Fall 2001

Chemistry 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 information

Molecules in strong magnetic fields

Molecules in strong magnetic fields Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway

More information

221B Lecture Notes Many-Body Problems II Molecular Physics

221B 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 information

The wavefunction that describes a bonding pair of electrons:

The wavefunction that describes a bonding pair of electrons: 4.2. Molecular Properties from VB Theory a) Bonding and Bond distances The wavefunction that describes a bonding pair of electrons: Ψ b = a(h 1 ) + b(h 2 ) where h 1 and h 2 are HAOs on adjacent atoms

More information

Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 17, March 1, 2006

Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 17, March 1, 2006 Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 17, March 1, 2006 (Some material in this lecture has been adapted from Cramer, C. J.

More information

Introduction to Electronic Structure Theory

Introduction 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 information

NMR and IR spectra & vibrational analysis

NMR 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 information

Wavepacket Correlation Function Approach for Nonadiabatic Reactions: Quasi-Jahn-Teller Model

Wavepacket Correlation Function Approach for Nonadiabatic Reactions: Quasi-Jahn-Teller Model Wavepacket Correlation for Nonadiabatic Reactions Bull. Korean Chem. Soc. 04, Vol. 35, No. 4 06 http://dx.doi.org/0.50/bkcs.04.35.4.06 Wavepacket Correlation Function Approach for Nonadiabatic Reactions:

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ Course overview This is a course about density functional theory (DFT)

More information