Detailed derivation of Gaussian orbital based matrix elements in electron structure calculations.
|
|
- Jodie Potter
- 6 years ago
- Views:
Transcription
1 Detailed derivation of Gaussian orbital based matri elements in electron structure calculations. T Petersson and B Hellsing Department of Physics, Väjö University, SE Väjö, Sweden Department of Physics, Göteborg University, SE-4 96 Göteborg, Sweden tomas.petersson@vu.se Abstract. A detailed derivation of the analytic solutions is presented for overlap, kinetic, nuclear attraction and electron repulsion integrals involving Cartesian Gaussian type orbitals. It is demonstrated how s-type orbitals can be used to evaluate integrals with higher angular momentum via the properties of Hermite polynomials and differentiation with respect to non-integration variables. PACS numbers:.3.rz, 3.5-p
2 Detailed derivation of Gaussian orbital based matri elements. Introduction In quantum chemistry, the Hartree-Fock (HF) approach is quite common in calculations of the electronic structure of molecules. Usually, one constructs the molecular orbitals by a linear combination of atomic orbitals, i.e., the molecular orbital is a sum of atomic orbitals, with different centers. The most common basis functions are the Slater type and the Gaussian-type functions. The use of Gaussian type functions was proposed by Boys [] since these are simpler to integrate, but the drawback is that one needs a larger basis set to achieve the same accuracy as with Slater type functions. In a HF calculation, the solution is found iteratively in a self-consistent procedure and various integrals have to be calculated in order to solve the equations. Here, we show how analytic formulas for these integrals can be obtained. Books in quantum chemistry [, 3] or computational physics [4] describing HF calculations seldom give any eamples beyond s-type functions, which are the easiest to calculate, and therefore contain no formulas for higher angular orbitals. Clementi and Davis [5] published the analytical result of the integrals involved, but for a novice reader it is hard to see through printing errors and the article does not present any derivations. Our derivations differ from those of Clementi and Davis. We believe the path we follow to be a slightly more pedagogical one, and that students in a more advanced quantum mechanics or computational physics course, in which the many electron problem is attacked computationally, may find this article useful. A Cartesian Gaussian type orbital (GTO), centered on atom A with orbital eponent α and non-negative integers l, m and n, is defined as φ A (r) = N( A ) l (y A y ) m (z A z ) n e α(r A) () with the normalization constant ( ) 3/4 α (8α) l+m+n / l! m! n! N =. () π (l)! (m)! (n)! When l = m = n =, we have an s-type GTO, when l + m + n = we have a p-type GTO, and so on. The word Gaussian function is used for a s-type GTO. A general idea to solve the integrals is given in section. Section 3 describes how the analytic results are derived, and the complete formulas for the integrals are written out in section 4. A suggestion for how to derive other solutions is found in section 5.. General outline of the problem It is simple to integrate a Gaussian function, but it is much more difficult to do it for a p-type or higher GTO. The idea is to generate matri elements involving orbitals with higher order GTO s from Gaussian functions without involving integration variables. This means that we can first integrate a product of Gaussian functions, and then apply an operator to get the desired result. The appropriate choice is to use Hermite polynomials since they are intimately connected to Gaussian functions. The Hermite
3 Detailed derivation of Gaussian orbital based matri elements 3 polynomials can be epanded as H n () = n! [ n] m= ( ) m () n m m! (n m)! the bracket-notation [k] is used to denote the largest integer k. The epansion of a monomial in Hermite polynomial reads n = n! [ n] H n i () n i! (n i)!. (4) i= With the use of Rodrigues formula H n () = ( )n e d n d n e the -component of an arbitrary GTO may be written as This relation gives us ( A ) l e α( A) = l! l d ( A ) l e α( A) = l! l [ l] i= [ l] i= i! (l i)! α l i l i l i i! (l i)! α l i A l i A l i (3) (5) e α( A). (6) d e α( A). (7) By taking the derivative over atomic position, A, instead of the usual position, r, the goal is achieved to epress an arbitrary GTO into a form the integration is only performed over a Gaussian function, i.e. s-type GTO. 3. The derivations We show now how one derives the analytic result for the overlap, kinetic, nuclear attraction and electron repulsion integrals for arbitrary GTO s from integrations over s-type GTO s. 3.. The overlap integral The overlap integral is defined as A B = d 3 r φ A (r)φ B (r) (8) φ A (r) and φ B (r) are defined as in (). As stated in (7), only Gaussian functions need to be integrated. reduces to solve d 3 r e α (r A) e α (r B) = e ηp(a B) d 3 r e γp(r P ) The problem = π 3/ γ 3/ p e ηp(a B) (9)
4 Detailed derivation of Gaussian orbital based matri elements 4 with η p = α α α + α γ p = α + α P = γ p (α A + α B) () and a special case of the Gaussian product theorem has been used, that is e α (r A) e α (r B) = e ηp(a B) e γp(r P ). () With the aid of (6), the -component of the unnormalized analytic result becomes A B = = l! l! l +l d ( A ) l e α ( A ) ( B ) l e α ( B ) [ l ] [ l ] i = i = l i l i A l i B l i [ i! i! (l i )! (l i )! α l i α l i ] e ηp(a B) () and similarly for the y- and z-components. It remains to establish the derivatives with respect to the atom positions, which is easily done with the help of Rodrigues formula (5), and the epansion of the Hermite polynomials (3). For the -component one obtains j A j j B j e ηp(a B) = ( ) j η j [ p H j ηp (A B ) ] e ηp(a B) [ j] = ( ) j j! r= ( ) r j r η j r p (A B ) j r e ηp(a B) (3) r! (j r)! j = j + j. (4) 3.. The kinetic integral The kinetic integral is defined as A B = d 3 r φ A (r) [ ] φ B (r). (5) In order to determine this integral, we are interested in the result of applying the Laplace operator to φ B. We restrict ourselves to the -component: [ ( B ) l e α ( B ) = α (l + )( B ) l α( B ) l + ] l (l )( B ) l e α ( B ). (6) The problem reduces to a sum of three overlap integrals with different angular momenta.
5 Detailed derivation of Gaussian orbital based matri elements The nuclear attraction integral The nuclear attraction integral is defined as A Z c r c B = Z c d 3 r φ A(r)φ B (r). (7) r r c The denominator can be taken care of by using, for eample, a Fourier transform, but a shorter way to the same result is the relation d e a = π (8) a which gives r r c = du e u (r r c). (9) π In order to solve the nuclear attraction integral with Gaussian functions, the Gaussian product theorem () is applied Z c d 3 r r r c e α (r A) e α (r B) = Z c d 3 r π = Z c e ηp(a B) d 3 r π du e α (r A) α (r B) u (r r c) = Z c e ηp(a B) du e γpp u rc π du e γp(r P ) u (r r c) d 3 r e (γp+u )r +(γ pp +u r c) r () the variables η p, γ p and P are as in (). In order to solve the last integral in (), the relation ( ) π dt e (at +bt+c) b = erf e (b ac)/a () 4a a is used to arrive at dt e (at bt) = π a eb /a. () We identify a = γ p + u and b = γ p P + u r c. The result for the nuclear attraction integral with Gaussian functions, after integration over r, will be With a change of variables πz c e ηp(a B) du (γ p + u ) 3/ e γpu (P r c) /(γ p+u ). (3) t = u (γ p + u ) du = (γ p + u ) 3/ γ p dt, (4) the end result for the nuclear attraction integral with Gaussian functions becomes πz c e ηp(a B) γ p dt e γp(p rc) t. (5)
6 Detailed derivation of Gaussian orbital based matri elements 6 The result is similar to (), but with two functions to differentiate. The solution is to apply Leibniz s theorem for differentiation of a product ( ) n n (uv) = n r u n r v. (6) n r r n r r= The net step is to establish the derivatives with respect to the atom positions. The result from differentiation of e ηp(a B) is found in (3). For the -component, the differentiation of the integral becomes, j A j j B j [ j] dt e γp(p rc) t = α j α j ( ) j+u j u (P r c ) j u j! u! (j u)! γp u u= dt t (j u) e γp(p rc) t (7) j = j + j. (8) The same integral in the right-hand side of (7) will also appear for the y- and z- components. The integral can be evaluated by or by the recursion F ν (u) = dt t ν e ut [ π = (ν)! ν! 4 ν u erf u e u ν+/ ν k= ] (ν k)! 4 k (ν k)! u k+ F ν (u) = uf ν+(u) + e u (3) ν + which is only stable numerically if it is performed downward due to round-off errors. Thus, (7) can be written j A j j B j [ j] dt e γp(p rc) t = α j α j ( ) j+u j u (P r c ) j u j! u! (j u)! γp u 3.4. The electron repulsion integral u= The electron repulsion integral is defined as A, C r B, D = (9) F j u ( γp (P r c ) ). (3) d 3 r d 3 φ A (r )φ B (r )φ C (r )φ D (r ) r. (3) r r We take care of the denominator in a similar way as in (9). With the use of the Gaussian product theorem () and some rewriting in order to take advantage of (),
7 Detailed derivation of Gaussian orbital based matri elements 7 the electron repulsion integral with Gaussian functions will be d 3 r d 3 r r r e α (r A) e α (r B) e α 3(r C) e α 4(r D) = e ηp(a B) e ηq(c D) d 3 r d 3 r du π e u r γpp γ q(r Q) e [(γp+u )r (γpp +u r ) r ] (33) η p = α α η q = α 3α 4 α + α α 3 + α 4 γ p = α + α γ q = α 3 + α 4 P = γ p (α A + α B) Q = γ q (α 3 C + α 4 D). (34) We identify a = γ p + u and b = γ p P + u r, use () and get for the electron repulsion integral with Gaussian functions, after integration over r πe ηp(a B) e ηq(c D) d 3 r du (γ p + u ) 3/ e (γpp +u r ) /(γ p+u ) e u r γpp γ q(r Q). (35) We rearrange the terms in order to integrate over r, identifies the terms a = [γ p γ q + (γ p + γ q )u ]/(γ p + u ), b = [γ p u P + (γ p + u )γ q Q]/(γ p + u ) to get the form of (). The result for the electron repulsion integral with Gaussian functions then becomes π 5/ e ηp(a B)e ηq(c D) du (η + u ) 3/ e ηu (P Q) /(η+u ) (36) (γ p + γ q ) 3/ η = With a change of variables t = γ pγ q γ p + γ q. (37) u du = (η + u ) 3/ η + u η dt (38) the end-result for the electron repulsion integral with Gaussian functions becomes π 5/ e ηp(a B)e ηq(c D) γ p γ q (γ p + γ q ) / dt e η(p Q) t. (39) The result from differentiation of e ηp(a B) is found in (3), and the differentiation of e ηq(c D) is similar. The differentiation of the integral becomes, for the -component, j A j j B j j 3 C j 3 = αj α j α j 3 3 α j 4 4 γ j +j p γ j 3+j 4 q j 4 D j 4 [ j] j! u= dt e η(p Q) t ( ) u+j +j (P Q ) j u η j u ( F ) j u η(p Q) (4) u! (j u)!
8 Detailed derivation of Gaussian orbital based matri elements 8 j = j + j + j 3 + j 4. (4) The function F ν (η(p Q) ) of (4) will also show up for the y- and z-components. 4. The complete analytic solutions for the integrals We now write the complete analytical result for each integral considered, by using the results from the last section. The Cartesian GTO used here are defined as φ A (r) = N ( A ) l (y A y ) m (z A z ) n e α (r A) (4) φ B (r) = N ( B ) l (y B y ) m (z B z ) n e α (r B) (43) φ C (r) = N 3 ( C ) l 3 (y C y ) m 3 (z C z ) n 3 e α 3(r C) (44) φ D (r) = N 4 ( D ) l 4 (y D y ) m 4 (z D z ) n 4 e α 4(r D) (45) N i = ( αi π (46) ) 3/4 [ (8αi ) l i+m i +n i ] / l i! m i! n i! (47) (l i )! (m i )! (n i )! r R = ( R )î + (y R y )ĵ + (z R z )ˆk. (48) The bracket-notation [k] for summations, is used to denote the largest integer k. 4.. The overlap integral The complete analytic solution to the overlap integral now reads ( ) 3/ π A B = N N e ηp(a B) S γ p S z (49) i,i,o j,j,p S y k,k,q γ p = α + α η p = α α γ p (5) i,i,o S = ( )l l! l! γ l +l p i γ(i +i )+o p (A B ) Ω o (l i )! (l i )! (Ω o)! i o ( ) o Ω! α l i i o α l i i o 4 i +i +o i! i! o! (5) Ω = l + l (i + i ). (5) S y and S z are similarly defined in terms of the y- and z-components. The summations are over the following ranges: i = l i = l o = Ω (53) and analogously for the y- and z-components.
9 Detailed derivation of Gaussian orbital based matri elements The kinetic integral The complete analytic solution to the kinetic integral is A B = N N [α (4(l + m + n ) + 6) A B 4α { A B, l + + A B, m + + A B, n + } l (l ) A B, l m (m ) A B, m ] n (n ) A B, n (54) B, l + = ( B ) l + (y B y ) m (z B z ) n e α (r B) (55) and similarly for the other terms. unnormalized functions The nuclear attraction integral The overlap integrals should be performed with The complete analytic solution to the nuclear attraction integral is A Z c B = Z cn N π e ( ηp(a B) A A y A z F ν γp (P r c ) ) (56) r c γ p i,i o,o r,u j,j p,p s,v k,k q,q t,w γ p = α + α η p = α + α γ p P = γ p (α A + α B) (57) A = ( ) l +l l! l! i i,i o,o r,u ( ) o+r (o + o )! 4 i +i +r i! i! o! o! r! i α o i r α o i r (A B ) o +o r (l i o )! (l i o )! (o + o r)! ( ) u µ! (P r c ) µ u 4 u u! (µ u u)! γ o +o r+u p µ = l + l (i + i ) (o + o ) o o r (58) ν = µ + µ y + µ z (u + v + w). (59) A y and A z are similarly defined in terms of the y- and z-components. The summations are over the following ranges: i = l i = l o = l i o = l i r = (o + o ) u = µ (6) and similarly for the y- and z-components. The function F ν (γ p (P r c ) ) is evaluated as in (9).
10 Detailed derivation of Gaussian orbital based matri elements 4.4. The electron repulsion integral The complete analytic solution to the electron repulsion integral is A, C B, D = N N N 3 N 4 π 5/ e ηp(a B) e ηq(c D) r γ p γ q γp + γ q ( J z F ) ν η(p Q) (6) J J y i,i,i 3,i 4 o,o,o 3,o 4 r,r,u j,j,j 3,j 4 p,p,p 3,p 4 s,s,v k,k,k 3,k 4 q,q,q 3,q 4 t,t,w γ p = α + α γ q = α 3 + α 4 P = γ p (α A + α B) i,i,i 3,i 4 o,o,o 3,o 4 r,r,u η p = α + α η q = α 3 + α 4 Q = (α 3 C + α 4 D) γ p γ q γ q η = γ pγ q (6) γ p + γ q J = ( )l +l l! l! γ l +l p i i o o ( )o+r (o + o )! 4 i +i +r i! i r! o! o! r! αo i r α o i r γ (i +i )+r p (A B ) o +o r (l i o )! (l i o )! (o + o r )! l 3! l 4! γ l 3+l 4 q i 3 i 4 o 3 o 4 ( )o3+r (o 3 + o 4 )! 4 i 3+i 4 +r i3! i r 4! o 3! o 4! r! αo 4 i 3 r 3 α o 3 i 4 r 4 γ (i 3+i 4 )+r q (C D ) o 3+o 4 r (l 3 i 3 o 3 )! (l 4 i 4 o 4 )! (o 3 + o 4 r )! ( ) u µ! η µ u (P Q ) µ u 4 u u! (µ u u)! µ = l + l + l 3 + l 4 (i + i + i 3 + i 4 ) (o + o + o 3 + o 4 ) (63) ν = µ + µ y + µ z (u + v + w). (64) J y and J z are similarly defined in terms of the y- and z-components. The summations are over the following ranges: i = l i = l i 3 = l 3 i 4 = l 4 o = l i o = l i o 3 = l 3 i 3 o 4 = l 4 i [ 4 ] r = (o + o ) r = (o 3 + o 4 ) u = µ (65)
11 Detailed derivation of Gaussian orbital based matri elements and similarly for the y- and z-components. The function F ν (η(p Q) ) is evaluated as in (9). 5. Summary and conclusions In a general electron structure calculation, matri elements between quantum state orbitals have to be calculated. We have presented a detailed derivation of the fundamental matri elements in this type of calculations, namely the overlap, kinetic, nuclear attraction and electron repulsion integrals for Gaussian type orbitals. One approach is to use the Gaussian product theorem before epanding the polynomials in Hermite polynomials, and then take the derivatives with respect to atomic positions in order to get the result presented in the article of Clementi and Davis. We have presented a different approach in this article. By epanding a polynomial in Hermite polynomials and using Rodrigues formula, an arbitrary Gaussian type orbital can be epressed as a sum over derivatives of Gaussian functions. The derivatives with respect to atomic positions are taken in order to integrate the Gaussian functions before the derivatives are taken. The derived recursion formulas are not restricted to electron structure calculations; they can be used in, e.g., electrostatistics when the charge density is epanded in terms of Cartesian Gaussian-type functions. References [] Boys S F 95 Proc. R. Soc. A [] Levine I N 999 Quantum Chemistry Fifth Edition (Upper Saddle River: Prentice Hall) [3] Szabo A and Ostlund N S 98 Modern Quantum Chemistry (New York: Macmillan Publishing Company) [4] Thijssen J M 999 Computational Physics (Cambridge: Cambridge University Press) [5] Clementi E and Davis D R 966 J. Comp. Phys. 3 44
Molecular Integral Evaluation
Molecular Integral Evaluation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 13th Sostrup Summer School Quantum Chemistry and
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 informationWorksheet Week 7 Section
Worksheet Week 7 Section 8.. 8.4. This worksheet is for improvement of your mathematical writing skill. Writing using correct mathematical epression and steps is really important part of doing math. Please
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES 3.2 The Product and Quotient Rules In this section, we will learn about: Formulas that enable us to differentiate new functions formed from old functions by
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 informationPre-Algebra 8 Notes Exponents and Scientific Notation
Pre-Algebra 8 Notes Eponents and Scientific Notation Rules of Eponents CCSS 8.EE.A.: Know and apply the properties of integer eponents to generate equivalent numerical epressions. Review with students
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 informationConservation of Linear Momentum
Conservation of Linear Momentum Once we have determined the continuit equation in di erential form we proceed to derive the momentum equation in di erential form. We start b writing the integral form of
More informationPolynomials and Factoring
7.6 Polynomials and Factoring Basic Terminology A term, or monomial, is defined to be a number, a variable, or a product of numbers and variables. A polynomial is a term or a finite sum or difference of
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationINTEGRATION OF RATIONAL FUNCTIONS BY PARTIAL FRACTIONS
INTEGRATION OF RATIONAL FUNCTIONS BY PARTIAL FRACTIONS. Introduction It is possible to integrate any rational function, constructed as the ratio of polynomials by epressing it as a sum of simpler fractions
More informationIntroduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 27, December 5, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 27, December 5, 2014 (Some material in this lecture has been adapted from Cramer, C. J.
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More 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 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 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 informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationA Simply Regularized Derivation of the Casimir Force
EJTP 3, No. 13 (2006) 121 126 Electronic Journal of Theoretical Physics A Simply Regularized Derivation of the Casimir Force H. Razmi Habibollah Razmi, 37185-359, Qom, I. R. Iran Department of Physics,
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 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 informationCharge renormalization at the large-d limit for N-electron atoms and weakly bound systems
Charge renormalization at the large-d limit for N-electron atoms and weakly bound systems S. Kais and R. Bleil Department of Chemistry, Purdue University, West Lafayette, Indiana 47907 Received 25 January
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 informationLecture 3: methods and terminology, part I
Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,
More informationMore On Carbon Monoxide
More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations of CO 1 / 26 More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationLesson 50 Integration by Parts
5/3/07 Lesson 50 Integration by Parts Lesson Objectives Use the method of integration by parts to integrate simple power, eponential, and trigonometric functions both in a mathematical contet and in a
More informationPHYS 3313 Section 001 Lecture # 22
PHYS 3313 Section 001 Lecture # 22 Dr. Barry Spurlock Simple Harmonic Oscillator Barriers and Tunneling Alpha Particle Decay Schrodinger Equation on Hydrogen Atom Solutions for Schrodinger Equation for
More informationLecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel
More informationTHE DISTRIBUTIVE LAW. Note: To avoid mistakes, include arrows above or below the terms that are being multiplied.
THE DISTRIBUTIVE LAW ( ) When an equation of the form a b c is epanded, every term inside the bracket is multiplied by the number or pronumeral (letter), and the sign that is located outside the brackets.
More informationMathematics for Chemistry: Exam Set 1
Mathematics for Chemistry: Exam Set 1 June 18, 017 1 mark Questions 1. The minimum value of the rank of any 5 3 matrix is 0 1 3. The trace of an identity n n matrix is equal to 1-1 0 n 3. A square matrix
More informationNext topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization"
Next topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization" Outline 1 Bosons and Fermions 2 N-particle wave functions ("first quantization" 3 The method of quantized
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 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 informationLecture 9: Molecular integral. Integrals of the Hamiltonian matrix over Gaussian-type orbitals
Lecture 9: Molecular integral evaluation Integrals of the Hamiltonian matrix over Gaussian-type orbitals Gaussian-type orbitals The de-facto standard for electronic-structure calculations is to use Gaussian-type
More informationEuler-Maclaurin summation formula
Physics 4 Spring 6 Euler-Maclaurin summation formula Lecture notes by M. G. Rozman Last modified: March 9, 6 Euler-Maclaurin summation formula gives an estimation of the sum N in f i) in terms of the integral
More informationAlgebra Final Exam Review Packet
Algebra 1 00 Final Eam Review Packet UNIT 1 EXPONENTS / RADICALS Eponents Degree of a monomial: Add the degrees of all the in the monomial together. o Eample - Find the degree of 5 7 yz Degree of a polynomial:
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules CHEM3006P or similar background knowledge is required for this course. This course has two parts: Part 1: Quantum Chemistry techniques for simulations of molecular
More informationTitleAnalytical Expression of the Hartre Author(s) Mukoyama, Takeshi; Yasui, Jun Citation Bulletin of the Institute for Chemi University (1992), 70(4): 385-391 Issue Date 1992-11-30 URL http://hdl.handle.net/2433/77474
More informationE nuc rep = Z AZ B r AB m =0.499 hartree =13.60 ev. E nuc rep (H 2 )= m/a m =0.714 hartree =19.43 ev.
Chemistr 31 Phsical Chemistr Homework Assignment # 7 1. Sketch qualitative contour maps of the following molecular orbitals for a diatomic molecule AB. Identif each as a bonding or an anti-bonding molecular
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
More informationLesson 53 Integration by Parts
5/0/05 Lesson 53 Integration by Parts Lesson Objectives Use the method of integration by parts to integrate simple power, eponential, and trigonometric functions both in a mathematical contet and in a
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
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 informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
More informationThe Quantum Harmonic Oscillator
The Classical Analysis Recall the mass-spring system where we first introduced unforced harmonic motion. The DE that describes the system is: where: Note that throughout this discussion the variables =
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationLecture Outlines Chapter 9. Physics, 3 rd Edition James S. Walker
Lecture Outlines Chapter 9 Physics, 3 rd Edition James S. Walker 2007 Pearson Prentice Hall This work is protected by United States copyright laws and is provided solely for the use of instructors in teaching
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 informationALGEBRAIC EXPRESSIONS AND POLYNOMIALS
MODULE - ic Epressions and Polynomials ALGEBRAIC EXPRESSIONS AND POLYNOMIALS So far, you had been using arithmetical numbers, which included natural numbers, whole numbers, fractional numbers, etc. and
More informationA summary of factoring methods
Roberto s Notes on Prerequisites for Calculus Chapter 1: Algebra Section 1 A summary of factoring methods What you need to know already: Basic algebra notation and facts. What you can learn here: What
More informationMini Lecture 9.1 Finding Roots
Mini Lecture 9. Finding Roots. Find square roots.. Evaluate models containing square roots.. Use a calculator to find decimal approimations for irrational square roots. 4. Find higher roots. Evaluat. a.
More information7.2 Multiplying Polynomials
Locker LESSON 7. Multiplying Polynomials Teas Math Standards The student is epected to: A.7.B Add, subtract, and multiply polynomials. Mathematical Processes A.1.E Create and use representations to organize,
More informationPart 2 Electronic Structure Theory
HS 2013: Vorlesung Physikalische Chemie III Part 2 Electronic Structure Theory Lecture 6 Dr. Mar'n O. Steinhauser Fraunhofer Ins'tute for High- Speed Dynamics, Ernst- Mach- Ins'tut, EMI, Freiburg Email:
More informationChapter 6: Vector Analysis
Chapter 6: Vector Analysis We use derivatives and various products of vectors in all areas of physics. For example, Newton s 2nd law is F = m d2 r. In electricity dt 2 and magnetism, we need surface and
More informationSOLUTIONS. Math 130 Midterm Spring True-False: Circle T if the statement is always true. Otherwise circle F.
SOLUTIONS Math 13 Midterm Spring 16 Directions: True-False: Circle T if the statement is always true. Otherwise circle F. Partial-Credit: Please show all work and fully justify each step. No credit will
More informationExact and heuristic minimization of Boolean functions
Eact and heuristic minimization of Boolean functions Quine a McCluskey Method a) Generation of all prime implicants b) Selection of a minimum subset of prime implicants, which will represent the original
More informationThe Sommerfeld Polynomial Method: Harmonic Oscillator Example
Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic
More informationOutline. Escaping the Double Cone: Describing the Seam Space with Gateway Modes. Columbus Workshop Argonne National Laboratory 2005
Outline Escaping the Double Cone: Describing the Seam Space with Gateway Modes Columbus Worshop Argonne National Laboratory 005 Outline I. Motivation II. Perturbation Theory and Intersection Adapted Coordinates
More informationarxiv: v1 [hep-ph] 30 Dec 2015
June 3, 8 Derivation of functional equations for Feynman integrals from algebraic relations arxiv:5.94v [hep-ph] 3 Dec 5 O.V. Tarasov II. Institut für Theoretische Physik, Universität Hamburg, Luruper
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 7 I. A SIMPLE EXAMPLE OF ANGULAR MOMENTUM ADDITION Given two spin-/ angular momenta, S and S, we define S S S The problem is to find the eigenstates of the
More informationA One-Slide Summary of Quantum Mechanics
A One-Slide Summary of Quantum Mechanics Fundamental Postulate: O! = a! What is!?! is an oracle! operator wave function (scalar) observable Where does! come from?! is refined Variational Process H! = E!
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationThe Product and Quotient Rules
The Product and Quotient Rules In this section, you will learn how to find the derivative of a product of functions and the derivative of a quotient of functions. A function that is the product of functions
More information5 3w. Unit 2 Function Operations and Equivalence Standard 4.1 Add, Subtract, & Multiply Polynomials
Unit Function Operations and Equivalence This document is meant to be used as an eample guide for each of the skills we will be holding students accountable for with Standard 4.1. This document should
More informationHigher-order C n dispersion coefficients for hydrogen
Higher-order C n dispersion coefficients for hydrogen J Mitroy* and M W J Bromley Faculty of Technology, Charles Darwin University, Darwin NT 0909, Australia Received 2 November 2004; published 11 March
More informationMultiplying a Polynomial by a Monomial
Lesson -3 Multiplying a Polynomial by a Monomial Lesson -3 BIG IDEA To multiply a polynomial by a monomial, multiply each term of the polynomial by the monomial and add the products. In earlier chapters,
More informationInfinitesimal Rotations
Universit of Connecticut DigitalCommons@UConn Chemistr Education Materials Department of Chemistr Januar 007 Infinitesimal Rotations Carl W. David Universit of Connecticut, Carl.David@uconn.edu Follow
More informationQuantization of the Spins
Chapter 5 Quantization of the Spins As pointed out already in chapter 3, the external degrees of freedom, position and momentum, of an ensemble of identical atoms is described by the Scödinger field operator.
More informationtakes the values:, +
rof. Dr. I. Nasser hys- (T-) February, 4 Angular Momentum Coupling Schemes We have so far considered only the coupling of the spin and orbital momentum of a single electron by means of the spin-orbit interaction.
More informationMolecular Integral Evaluation. The 9th Sostrup Summer School. Quantum Chemistry and Molecular Properties June 25 July 7, 2006.
1 Molecular Integral Evaluation The 9th Sostrup Summer School Quantum Chemistry and Molecular Properties June 25 July 7, 2006 Trygve Helgaker Department of Chemistry, University of Oslo, Norway Overview
More informationarxiv:quant-ph/ v1 20 Apr 1995
Combinatorial Computation of Clebsch-Gordan Coefficients Klaus Schertler and Markus H. Thoma Institut für Theoretische Physik, Universität Giessen, 3539 Giessen, Germany (February, 008 The addition of
More informationPOLYNOMIALS CHAPTER 2. (A) Main Concepts and Results
CHAPTER POLYNOMIALS (A) Main Concepts and Results Meaning of a Polynomial Degree of a polynomial Coefficients Monomials, Binomials etc. Constant, Linear, Quadratic Polynomials etc. Value of a polynomial
More information0.1. Linear transformations
Suggestions for midterm review #3 The repetitoria are usually not complete; I am merely bringing up the points that many people didn t now on the recitations Linear transformations The following mostly
More informationFINAL EXAMINATION PAPER # 64 TITLE PAGE COURSE: MATH 1210
TITLE PAGE FAMILY NAME: (Print in ink) GIVEN NAME(S): (Print in ink) STUDENT NUMBER: SEAT NUMBER: SIGNATURE: (in ink) (I understand that cheating is a serious offense) Please indicate your instructor by
More informationExpanding brackets and factorising
CHAPTER 8 Epanding brackets and factorising 8 CHAPTER Epanding brackets and factorising 8.1 Epanding brackets There are three rows. Each row has n students. The number of students is 3 n 3n. Two students
More informationarxiv:math-ph/ v1 10 Jan 2005
Asymptotic and eact series representations for the incomplete Gamma function arxiv:math-ph/5119v1 1 Jan 5 Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo 34, Colima, Colima,
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationUnit 2: Rational Expressions
Rational Epressions Pure Math 0 Notes Unit : Rational Epressions -: Simplifing Rational Epressions Rational Epressions: - fractions with polnomials as numerator and / or denominator. To Simplif (Reduce)
More informationQuantum Mechanics: Vibration and Rotation of Molecules
Quantum Mechanics: Vibration and Rotation of Molecules 8th April 2008 I. 1-Dimensional Classical Harmonic Oscillator The classical picture for motion under a harmonic potential (mass attached to spring
More informationCh. 7.1 Polynomial Degree & Finite Differences
Ch. 7.1 Polynomial Degree & Finite Differences Learning Intentions: Define terminology associated with polynomials: term, monomial, binomial & trinomial. Use the finite differences method to determine
More informationNotes on Ewald summation techniques
February 3, 011 Notes on Ewald summation techniques Adapted from similar presentation in PHY 71 he total electrostatic potential energy of interaction between point charges {q i } at the positions {r i
More informationGravitational radiation
Lecture 28: Gravitational radiation Gravitational radiation Reading: Ohanian and Ruffini, Gravitation and Spacetime, 2nd ed., Ch. 5. Gravitational equations in empty space The linearized field equations
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
More informationNotes on Density Functional Theory
Notes on Density Functional Theory Rocco Martinazzo E-mail: rocco.martinazzo@unimi.it Contents 1 Introduction 1 Density Functional Theory 7 3 The Kohn-Sham approach 11 1 Introduction We consider here a
More information4.8 Partial Fraction Decomposition
8 CHAPTER 4. INTEGRALS 4.8 Partial Fraction Decomposition 4.8. Need to Know The following material is assumed to be known for this section. If this is not the case, you will need to review it.. When are
More informationWhy use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationProblem 1: Step Potential (10 points)
Problem 1: Step Potential (10 points) 1 Consider the potential V (x). V (x) = { 0, x 0 V, x > 0 A particle of mass m and kinetic energy E approaches the step from x < 0. a) Write the solution to Schrodinger
More informationQuantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 21 Square-Integrable Functions
Quantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 21 Square-Integrable Functions (Refer Slide Time: 00:06) (Refer Slide Time: 00:14) We
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 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 informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationAlgebra. Robert Taggart
Algebra Robert Taggart Table of Contents To the Student.............................................. v Unit 1: Algebra Basics Lesson 1: Negative and Positive Numbers....................... Lesson 2: Operations
More informationHarmonic Oscillator Eigenvalues and Eigenfunctions
Chemistry 46 Fall 217 Dr. Jean M. Standard October 4, 217 Harmonic Oscillator Eigenvalues and Eigenfunctions The Quantum Mechanical Harmonic Oscillator The quantum mechanical harmonic oscillator in one
More informationEXPONENT REVIEW!!! Concept Byte (Review): Properties of Exponents. Property of Exponents: Product of Powers. x m x n = x m + n
Algebra B: Chapter 6 Notes 1 EXPONENT REVIEW!!! Concept Byte (Review): Properties of Eponents Recall from Algebra 1, the Properties (Rules) of Eponents. Property of Eponents: Product of Powers m n = m
More informationMathematics for Chemistry: Exam Set 1
Mathematics for Chemistry: Exam Set 1 March 19, 017 1 mark Questions 1. The maximum value of the rank of any 5 3 matrix is (a (b3 4 5. The determinant of an identity n n matrix is equal to (a 1 (b -1 0
More information1 Summary of Chapter 2
General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical
More information