Chapter 5 Vibrational Motion

Size: px
Start display at page:

Download "Chapter 5 Vibrational Motion"

Transcription

1 Fall 4 Chapter 5 Vibratioal Motio Potetial Eergy Surfaces, Rotatios ad Vibratios Harmoic Oscillator Geeral Solutio for H.O.: Operator Techique Vibratioal Selectio Rules... 7 Polyatomic Molecules Beyod the harmoic oscillator approximatio Chapter 5 Vibratioal Motio Potetial Eergy Surfaces, Rotatios ad Vibratios Suppose we assume the uclei of a molecule are fixed, the we ca solve the Schrodiger equatio for the electros H ˆ ( r, r,... r ) = E ( r, r,... r ) N N This is a complicated problem, that will be discussed later (Chapter 7 ad beyod) We would get the (groud state) eergy at a particular uclear cofiguratio { R } Hece, assume we ca solve this We ca fit a curve through the poits V({ R }) This V({ R }) is called the potetial eergy surface (PES) ad is a crucial cocept i chemistry, eg. Chapter 5 Vibratioal Motio 65

2 Fall 4 Miima o the PES we associate with differet isomers, for example with the reactio A+ B C + D. Saddle poits o the PES we associate with trasitio states Obtaiig the miima (+ curvature) we ca get thermochemical iformatio. From TS we get ito rates of reactios (Chem54, Chem 35) Ideally, oce we have obtaied PES, we should solve for the Quatum eergies ivolvig the uclei. This is very complicated, ad today ca oly be doe for small molecules (up to 5 atoms). What ca be doe is solvig for molecular rotatios + vibratios i the harmoic oscillator + Rigid Rotor approximatio The crux is to make a quadratic approximatio to the PES aroud each miimum ( = molecular species) V (! R) V (R e ) + αβ k αβ ) α ) β Chapter 5 Vibratioal Motio 66

3 Fall 4 A molecule with N atoms has 3 N degrees of freedom 3 overall traslatios 3 overall rotatios or rotatios for liear molecules (3N 6) Vibratios, [3N 5 for liear molecule] k αβ : α,β =...3N All of the iformatio V (! R e ),! R e ad k αβ (! R e ) is obtaied from a electroic structure program (like Gaussia, Gamess, Turbomole ) Oce we have these, oe ca use exact calculatios to solve for Harmoic oscillator + Rigid Rotor eergies. The harmoic (quadratic) potetial is a approximatio, but ofte works well, certaily for stiff molecules/ degrees of freedom. What we will do ext is: discuss harmoic oscillator for diatomic molecules Groud state of harmoic oscillator All excited states usig operator techique Geeralize to polyatomic molecules More accurate solutio for diatomic (Beyod H.O.) Selectio rules, vibratioal spectroscopy Harmoic Oscillator I the followig chapter we will go o to discuss rotatios For a geeral polyatomic molecule we ca defie H.O. Hamiltoia as H = α For a diatomic this ca be reduced to! P M α + α a,b=...3n Ĥ =! d µ + kx k ab ) a ) b Here µ = M M is the [kg] reduced mass M + M x = ( R R e ) is deviatio [m] from equilibrium k is the force costat [Nm ] Chapter 5 Vibratioal Motio 67

4 Fall 4 I McQuarrie this is derived usig classical equatios of motio. I will post o the Website a geeral derivatio to get the H.O. form for a geeral polyatomic, but will ot discuss i class. Here I will simply use the form for H. The groud state solutio has the form o I will discuss the full solutio. e α x /. Let us determie the costat α, ad the eergy. Later! d e α x / αx / = α xe d x / x / x / e α e α x e α = α + α d µ + kx e α x / =! µ α! µ α x + kx e α x / = E e α x / E =! µ α k! α µ = α =! µk E =! µ µk! =! k µ =!ω ω = k µ ; α = µ! ω Geeral Solutio for H.O.: Operator Techique (see appedix 5 i McQuarrie) Defie q = x α e α x / Solutio Ĥ =! m d + kx e α x / α = µω /! = µk /! e q / q =! µk x x =! µk q Chapter 5 Vibratioal Motio 68

5 Fall 4 Ĥ =! µ x = q x µk! q =! µk q q + k! µk q =! k µ q + q =!ω d dq + q q are called dimesioless coordiates (Check that q =! µk x is ideed dimesioless!) Commutatio Relatio: d q, = dq d d q q f( q) = f( q) dq dq Defie ew operators: ˆ d b + = q dq The!ωb + b =!ω =!ω ˆ d b = q+ dq q d dq q + d dq q d dq + q, d dq = Ĥ!ω Hˆ ˆ ω b + = h b+ Commutatio Relatios, betwee b operators: bb ˆ, ˆ = b ˆ +, bˆ + = ˆˆ+ d d bb, = q+, q dq dq d d d d = [ qq, ] +, q q,, dq dq dq dq Chapter 5 Vibratioal Motio 69

6 Fall 4 Hece = ( + ( ) ) = bb ˆ, ˆ + = ˆ+ b, bˆ = Usig this form of the Hamiltoia ad the commutatio relatios we ca derive the eigevalues ad eigefuctios of H.O.!! a) If ( q ) is eigefuctio with eigevalue E the Proof: i) ˆ + b ( q) is eigefuctio with E + = E +!ω ii) bˆ ( q) is eigefuctio with E = E!ω i) Ĥ ˆb + ii) ( ) =!ω ˆb + ˆb + ˆb + =!ω ˆb + ˆb, ˆb+ + ˆb + ˆb+ ˆb + = ˆb ( +!ω + Ĥ ) = (E +!ω ) ˆb + ( ) =!ω ˆb + ˆb + Ĥ ˆb ˆb =!ω ˆb +, ˆb + ˆb ˆb + + ˆb =!ω ˆb + ˆb ˆb + ˆb + = ˆb (!ω + Ĥ ) = E!ω ( ) ˆb ˆb + ad ˆb are called the raisig ad lowerig operators, or ladder operators Chapter 5 Vibratioal Motio 7

7 Fall 4 What about the groud state? ˆb = E!ω ( ) ˆb Still true, but E!ω < E!! Oly way out: b ˆ ( q ) = d q+ ( q) dq = Differetial equatio with solutio ( ) q q = e d q q q+ e = ( q q) e = dq Puttig it all together 4 q ( q) = e ormalized π ( q ) ( ˆ ) ( )! b + = q ormalized ˆ d b + = q dq E = + h ω,,,3... = First couple of uormalized fuctios i terms of q : Chapter 5 Vibratioal Motio 7

8 Fall 4 ( ) q q = e ( ) d q q q e = qe dq q / ( ) d q / / ( q q q qe = q ) e dq ( ) d 3 ( q q q q ) qe dq = ( q 6 q) e Etc. We ca obtai all eigefuctios by differetiatio / 3 q / ( x ) substitute x = q α + ormalize The wave fuctios take the form H ( q) e q H ( ) q are Hermite polyomials, they are either odd or eve fuctios. (Each polyomial cotais oly odd or oly eve terms) Vibratioal Selectio Rules Later we will discuss more geerally the radiatio process. For ow the trasitio dipole momet is used to defie the stregth of a spectroscopic traslatio *( x ) µ ( x ) m( x ) m d µ µ ( x) = µ ( x) + x+... x µ is the dipole at equilibrium distace, d µ is the chage i dipole with chagig x (iterucleardistace) Eg. N : d µ = HF: d CO: d *( x) µ ( x) ( x) µ large µ small m = µ *( x) m( x) (if m) Chapter 5 Vibratioal Motio 7

9 Fall 4 d µ + *( xx ) m( x ) d µ α *( qq ) m( qdq ) dµ α *( )( ˆ+ ˆ q b b) m( q) dq + bˆ( m )~ ˆ + m b ( m )~ m+ For Harmoic Oscillator we ca oly get allowed trasitios if Δ = + or Δ = ΔE = ±!ω Polyatomic Molecules! Ĥ = + α µ α α,β k αβ ) α ) β By some maipulatio (see otes o webpage) this ca be writte as a sum of H.O. Hamiltoias.tedious derivatios Ĥ =! ω i d i dq + i q i + =! ω ˆbi ˆbi i + i The coordiates q : are liear combiatios of atomic displacemet vectors 3 N coordiate 3 traslatio, 3 rotatio, (3 N 6) vibratio For liear molecule, rotatio aroud axis is ot a displacemet rotatios, ad (3 N 5) vibratios Eg. Water Normal modes : display symmetry of molecule Chapter 5 Vibratioal Motio 73

10 Fall 4 The eigefuctios are simply products of dimesioal H.O. fuctios ( q ) ( q ) ( q ) k l m 3 E k,l,m = k +!ω + l +!ω + m +!ω 3 Very simple solutio. Oe eeds to diagoalize the mass weighted hessia k αβ M α M β 3N 3N Matrix Reasoable approximatios to all vibratioal levels Statistical Mechaics (recall Chem 35) Beyod the harmoic oscillator approximatio The true potetial is ot quadratic. For large molecules this is ot so easy to correct (people would use low order perturbatio theory, based o a quartic force field) For smalls molecules, i particular diatomics, oe ca solve for the full vibratioal problem The exact eergies are ot equidistat. A better approximatio is E() =!ω e +!ω x + e e The eergy levels are usually called g(v), v =,vibratioal quatum umber There are a fiite umber of boud levels Aother effect is that trasitio momets ca be ozero eve whe Δ ±. This leads to the observatio of overtoes. A ofte used form for the potetial for a diatomic is the Morse potetial Chapter 5 Vibratioal Motio 74

11 Fall 4 ( x ) V( x) = D e β V x x = V D eβ x e = = k x = equilibrium geometry β = k D D ca be used from experimet ( e D measurable) e This is still a approximate potetial Today: Potetial eergy curves ca be calculated usig electroic structure programs. Chapter 5 Vibratioal Motio 75

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples: 5.6 4 Lecture #3-4 page HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT Do t restrict the wavefuctio to a sigle term! Could be a liear combiatio of several wavefuctios e.g. two terms:

More information

1. Szabo & Ostlund: 2.1, 2.2, 2.4, 2.5, 2.7. These problems are fairly straightforward and I will not discuss them here.

1. Szabo & Ostlund: 2.1, 2.2, 2.4, 2.5, 2.7. These problems are fairly straightforward and I will not discuss them here. Solutio set III.. Szabo & Ostlud:.,.,.,.5,.7. These problems are fairly straightforward ad I will ot discuss them here.. N! N! i= k= N! N! N! N! p p i j pi+ pj i j i j i= j= i= j= AA ˆˆ= ( ) Pˆ ( ) Pˆ

More information

1 Adiabatic and diabatic representations

1 Adiabatic and diabatic representations 1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular

More information

5.76 Lecture #33 5/08/91 Page 1 of 10 pages. Lecture #33: Vibronic Coupling

5.76 Lecture #33 5/08/91 Page 1 of 10 pages. Lecture #33: Vibronic Coupling 5.76 Lecture #33 5/8/9 Page of pages Lecture #33: Vibroic Couplig Last time: H CO A A X A Electroically forbidde if A -state is plaar vibroically allowed to alterate v if A -state is plaar iertial defect

More information

Vibrational Spectroscopy 1

Vibrational Spectroscopy 1 Applied Spectroscopy Vibratioal Spectroscopy Recommeded Readig: Bawell ad McCash Chapter 3 Atkis Physical Chemistry Chapter 6 Itroductio What is it? Vibratioal spectroscopy detects trasitios betwee the

More information

1. Hydrogen Atom: 3p State

1. Hydrogen Atom: 3p State 7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).

More information

Zeros of Polynomials

Zeros of Polynomials Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree

More information

Solution of Quantum Anharmonic Oscillator with Quartic Perturbation

Solution of Quantum Anharmonic Oscillator with Quartic Perturbation ISS -79X (Paper) ISS 5-0638 (Olie) Vol.7, 0 Abstract Solutio of Quatum Aharmoic Oscillator with Quartic Perturbatio Adelaku A.O. Departmet of Physics, Wesley Uiversity of Sciece ad Techology, Odo, Odo

More information

Lecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box

Lecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box 561 Fall 013 Lecture #5 page 1 Last time: Lecture #5: Begi Quatum Mechaics: Free Particle ad Particle i a 1D Box u 1 u 1-D Wave equatio = x v t * u(x,t): displacemets as fuctio of x,t * d -order: solutio

More information

5.80 Small-Molecule Spectroscopy and Dynamics

5.80 Small-Molecule Spectroscopy and Dynamics MIT OpeCourseWare http://ocw.mit.edu 5.80 Small-Molecule Spectroscopy ad Dyamics Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Lecture # 33 Supplemet

More information

Perturbation Theory I (See CTDL , ) Last time: derivation of all matrix elements for Harmonic-Oscillator: x, p, H. n ij n.

Perturbation Theory I (See CTDL , ) Last time: derivation of all matrix elements for Harmonic-Oscillator: x, p, H. n ij n. Perturbatio Theory I (See CTDL 1095-1107, 1110-1119) 14-1 Last time: derivatio of all matrix elemets for Harmoic-Oscillator: x, p, H selectio rules scalig ij x i j i steps of 2 e.g. x : = ± 3, ± 1 xii

More information

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to: 2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium

More information

Lecture 14 and 15: Algebraic approach to the SHO. 1 Algebraic Solution of the Oscillator 1. 2 Operator manipulation and the spectrum 4

Lecture 14 and 15: Algebraic approach to the SHO. 1 Algebraic Solution of the Oscillator 1. 2 Operator manipulation and the spectrum 4 Lecture 14 ad 15: Algebraic approach to the SHO B. Zwiebach April 5, 2016 Cotets 1 Algebraic Solutio of the Oscillator 1 2 Operator maipulatio ad the spectrum 4 1 Algebraic Solutio of the Oscillator We

More information

PHY4905: Nearly-Free Electron Model (NFE)

PHY4905: Nearly-Free Electron Model (NFE) PHY4905: Nearly-Free Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. No-degeerate eigestates

More information

Chapter 8 Approximation Methods, Hueckel Theory

Chapter 8 Approximation Methods, Hueckel Theory Witer 3 Chem 356: Itroductory Quatum Mechaics Chapter 8 Approximatio Methods, Huecel Theory... 8 Approximatio Methods... 8 The Liear Variatioal Priciple... Chapter 8 Approximatio Methods, Huecel Theory

More information

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields Hydroge (atoms, molecules) i exteral fields Static electric ad magetic fields Oscyllatig electromagetic fields Everythig said up to ow has to be modified more or less strogly if we cosider atoms (ad ios)

More information

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

More information

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Assignment 2 Solutions SOLUTION. ϕ 1  = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ. PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio

More information

Physics 2D Lecture Slides Lecture 22: Feb 22nd 2005

Physics 2D Lecture Slides Lecture 22: Feb 22nd 2005 Physics D Lecture Slides Lecture : Feb d 005 Vivek Sharma UCSD Physics Itroducig the Schrodiger Equatio! (, t) (, t) #! " + U ( ) "(, t) = i!!" m!! t U() = characteristic Potetial of the system Differet

More information

The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.

The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m 2 + V ( r,t) (1. Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 1-1 574 TIME-DEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Time-evolutio for time-idepedet Hamiltoias The time evolutio of the state of a quatum system

More information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka) 7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.

More information

Physics 232 Gauge invariance of the magnetic susceptibilty

Physics 232 Gauge invariance of the magnetic susceptibilty Physics 232 Gauge ivariace of the magetic susceptibilty Peter Youg (Dated: Jauary 16, 2006) I. INTRODUCTION We have see i class that the followig additioal terms appear i the Hamiltoia o addig a magetic

More information

Office: JILA A709; Phone ;

Office: JILA A709; Phone ; Office: JILA A709; Phoe 303-49-7841; email: weberjm@jila.colorado.edu Problem Set 5 To be retured before the ed of class o Wedesday, September 3, 015 (give to me i perso or slide uder office door). 1.

More information

Physics 2D Lecture Slides Lecture 25: Mar 2 nd

Physics 2D Lecture Slides Lecture 25: Mar 2 nd Cofirmed: D Fial Eam: Thursday 8 th March :3-:3 PM WH 5 Course Review 4 th March am WH 5 (TBC) Physics D ecture Slides ecture 5: Mar d Vivek Sharma UCSD Physics Simple Harmoic Oscillator: Quatum ad Classical

More information

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j. Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α

More information

Appendix: The Laplace Transform

Appendix: The Laplace Transform Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio

More information

Lecture 25 (Dec. 6, 2017)

Lecture 25 (Dec. 6, 2017) Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed

More information

The Born-Oppenheimer approximation

The Born-Oppenheimer approximation The Bor-Oppeheimer approximatio 1 Re-writig the Schrödiger equatio We will begi from the full time-idepedet Schrödiger equatio for the eigestates of a molecular system: [ P 2 + ( Pm 2 + e2 1 1 2m 2m m

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

EXPERIMENT OF SIMPLE VIBRATION

EXPERIMENT OF SIMPLE VIBRATION EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic

More information

PAPER : IIT-JAM 2010

PAPER : IIT-JAM 2010 MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

On the Quantum Theory of Molecules

On the Quantum Theory of Molecules O the Quatum Theory of Molecules M. Bor a, J.R. Oppeheimer b a Istitute of Theoretical Physics, Göttige b Istitute of Theoretical Physics, Göttige Abstract It will be show that the familiar compoets of

More information

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered

More information

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc. Quatum Mechaics I 1 April, 14 Assigmet 5: Solutio 1 For a particle icidet o a potetial step with E < V, show that the magitudes of the amplitudes of the icidet ad reflected waves fuctios are the same Fid

More information

TIME-CORRELATION FUNCTIONS

TIME-CORRELATION FUNCTIONS p. 8 TIME-CORRELATION FUNCTIONS Time-correlatio fuctios are a effective way of represetig the dyamics of a system. They provide a statistical descriptio of the time-evolutio of a variable for a esemble

More information

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

d dx where k is a spring constant

d dx where k is a spring constant Vorlesug IX Harmoic Oscillator 1 Basic efiitios a properties a classical mechaics Oscillator is efie as a particle subject to a liear force fiel The force F ca be epresse i terms of potetial fuctio V F

More information

On the Quantum Theory of Molecules M Born and J R Oppenheimer Ann. Physik 84, 457 (1927) January 23, 2002 Translated by S M Blinder with emendations b

On the Quantum Theory of Molecules M Born and J R Oppenheimer Ann. Physik 84, 457 (1927) January 23, 2002 Translated by S M Blinder with emendations b O the Quatum Theory of Molecules M Bor ad J R Oppeheimer A. Physik 84, 457 (197) Jauary 3, Traslated by S M Blider with emedatios by Bria Sutclie ad Wolf Geppert Abstract It will be show that the familiar

More information

The Method of Least Squares. To understand least squares fitting of data.

The Method of Least Squares. To understand least squares fitting of data. The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve

More information

Forces: Calculating Them, and Using Them Shobhana Narasimhan JNCASR, Bangalore, India

Forces: Calculating Them, and Using Them Shobhana Narasimhan JNCASR, Bangalore, India Forces: Calculatig Them, ad Usig Them Shobhaa Narasimha JNCASR, Bagalore, Idia shobhaa@jcasr.ac.i Shobhaa Narasimha, JNCASR 1 Outlie Forces ad the Hellma-Feyma Theorem Stress Techiques for miimizig a fuctio

More information

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

September 2012 C1 Note. C1 Notes (Edexcel) Copyright   - For AS, A2 notes and IGCSE / GCSE worksheets 1 September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright

More information

Physics 201 Final Exam December

Physics 201 Final Exam December Physics 01 Fial Exam December 14 017 Name (please prit): This test is admiistered uder the rules ad regulatios of the hoor system of the College of William & Mary. Sigature: Fial score: Problem 1 (5 poits)

More information

Math 128A: Homework 1 Solutions

Math 128A: Homework 1 Solutions Math 8A: Homework Solutios Due: Jue. Determie the limits of the followig sequeces as. a) a = +. lim a + = lim =. b) a = + ). c) a = si4 +6) +. lim a = lim = lim + ) [ + ) ] = [ e ] = e 6. Observe that

More information

Numerical Methods in Fourier Series Applications

Numerical Methods in Fourier Series Applications Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets

More information

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial. Taylor Polyomials ad Taylor Series It is ofte useful to approximate complicated fuctios usig simpler oes We cosider the task of approximatig a fuctio by a polyomial If f is at least -times differetiable

More information

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015 ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],

More information

Recurrence Relations

Recurrence Relations Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The

More information

New Version of the Rayleigh Schrödinger Perturbation Theory: Examples

New Version of the Rayleigh Schrödinger Perturbation Theory: Examples New Versio of the Rayleigh Schrödiger Perturbatio Theory: Examples MILOŠ KALHOUS, 1 L. SKÁLA, 1 J. ZAMASTIL, 1 J. ČÍŽEK 2 1 Charles Uiversity, Faculty of Mathematics Physics, Ke Karlovu 3, 12116 Prague

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

Definition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.

Definition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4. 4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad

More information

Physics Supplement to my class. Kinetic Theory

Physics Supplement to my class. Kinetic Theory Physics Supplemet to my class Leaers should ote that I have used symbols for geometrical figures ad abbreviatios through out the documet. Kietic Theory 1 Most Probable, Mea ad RMS Speed of Gas Molecules

More information

Lecture #3. Math tools covered today

Lecture #3. Math tools covered today Toay s Program:. Review of previous lecture. QM free particle a particle i a bo. 3. Priciple of spectral ecompositio. 4. Fourth Postulate Math tools covere toay Lecture #3. Lear how to solve separable

More information

Lanczos-Haydock Recursion

Lanczos-Haydock Recursion Laczos-Haydock Recursio Bor: Feb 893 i Székesehérvár, Hugary- Died: 5 Jue 974 i Budapest, Hugary Corelius Laczos From a abstract mathematical viewpoit, the method for puttig a symmetric matrix i three-diagoal

More information

Exercises and Problems

Exercises and Problems HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle

More information

CHAPTER 5. Theory and Solution Using Matrix Techniques

CHAPTER 5. Theory and Solution Using Matrix Techniques A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL

More information

True Nature of Potential Energy of a Hydrogen Atom

True Nature of Potential Energy of a Hydrogen Atom True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial

More information

Paired Data and Linear Correlation

Paired Data and Linear Correlation Paired Data ad Liear Correlatio Example. A group of calculus studets has take two quizzes. These are their scores: Studet st Quiz Score ( data) d Quiz Score ( data) 7 5 5 0 3 0 3 4 0 5 5 5 5 6 0 8 7 0

More information

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations . Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a

More information

Algebra of Least Squares

Algebra of Least Squares October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal

More information

Orthogonal polynomials derived from the tridiagonal representation approach

Orthogonal polynomials derived from the tridiagonal representation approach Orthogoal polyomials derived from the tridiagoal represetatio approach A. D. Alhaidari Saudi Ceter for Theoretical Physics, P.O. Box 374, Jeddah 438, Saudi Arabia Abstract: The tridiagoal represetatio

More information

PHYS-505 Parity and other Discrete Symmetries Lecture-7!

PHYS-505 Parity and other Discrete Symmetries Lecture-7! PHYS-505 Parity ad other Discrete Symmetries Lecture-7! 1 Discrete Symmetries So far we have cosidered cotiuous symmetry operators that is, operatios that ca be obtaied by applyig successively ifiitesimal

More information

x c the remainder is Pc ().

x c the remainder is Pc (). Algebra, Polyomial ad Ratioal Fuctios Page 1 K.Paulk Notes Chapter 3, Sectio 3.1 to 3.4 Summary Sectio Theorem Notes 3.1 Zeros of a Fuctio Set the fuctio to zero ad solve for x. The fuctio is zero at these

More information

Eigenvalues and Eigenfunctions of Woods Saxon Potential in PT Symmetric Quantum Mechanics

Eigenvalues and Eigenfunctions of Woods Saxon Potential in PT Symmetric Quantum Mechanics Eigevalues ad Eigefuctios of Woods Saxo Potetial i PT Symmetric Quatum Mechaics Ayşe Berkdemir a, Cüeyt Berkdemir a ad Ramaza Sever b* a Departmet of Physics, Faculty of Arts ad Scieces, Erciyes Uiversity,

More information

Math 475, Problem Set #12: Answers

Math 475, Problem Set #12: Answers Math 475, Problem Set #12: Aswers A. Chapter 8, problem 12, parts (b) ad (d). (b) S # (, 2) = 2 2, sice, from amog the 2 ways of puttig elemets ito 2 distiguishable boxes, exactly 2 of them result i oe

More information

THE KALMAN FILTER RAUL ROJAS

THE KALMAN FILTER RAUL ROJAS THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider

More information

Hilbert Space Methods Used in a First Course in Quantum Mechanics

Hilbert Space Methods Used in a First Course in Quantum Mechanics Hilbert Space Methods Used i a First Course i Quatum Mechaics Victor Poliger Physics/Mathematics Bellevue College 03/07/3-04//3 Outlie The Ifiite Square Well: A Follow-Up Timelie of basic evets Statistical

More information

Lesson 10: Limits and Continuity

Lesson 10: Limits and Continuity www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals

More information

1 Last time: similar and diagonalizable matrices

1 Last time: similar and diagonalizable matrices Last time: similar ad diagoalizable matrices Let be a positive iteger Suppose A is a matrix, v R, ad λ R Recall that v a eigevector for A with eigevalue λ if v ad Av λv, or equivaletly if v is a ozero

More information

CHAPTER I: Vector Spaces

CHAPTER I: Vector Spaces CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig

More information

Properties and Tests of Zeros of Polynomial Functions

Properties and Tests of Zeros of Polynomial Functions Properties ad Tests of Zeros of Polyomial Fuctios The Remaider ad Factor Theorems: Sythetic divisio ca be used to fid the values of polyomials i a sometimes easier way tha substitutio. This is show by

More information

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model. 5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Uiversity of Wasigto Departmet of Cemistry Cemistry 453 Witer Quarter 15 Lecture 14. /11/15 Recommeded Text Readig: Atkis DePaula: 9.1, 9., 9.3 A. Te Equipartitio Priciple & Eergy Quatizatio Te Equipartio

More information

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A

More information

SOLID MECHANICS TUTORIAL BALANCING OF RECIPROCATING MACHINERY

SOLID MECHANICS TUTORIAL BALANCING OF RECIPROCATING MACHINERY SOLID MECHANICS TUTORIAL BALANCING OF RECIPROCATING MACHINERY This work covers elemets of the syllabus for the Egieerig Coucil Exam D5 Dyamics of Mechaical Systems. O completio of this tutorial you should

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows

More information

1 Approximating Integrals using Taylor Polynomials

1 Approximating Integrals using Taylor Polynomials Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................

More information

Lecture 7: Properties of Random Samples

Lecture 7: Properties of Random Samples Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ

More information

Fundamental Concepts: Surfaces and Curves

Fundamental Concepts: Surfaces and Curves UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat

More information

Lecture 8: October 20, Applications of SVD: least squares approximation

Lecture 8: October 20, Applications of SVD: least squares approximation Mathematical Toolkit Autum 2016 Lecturer: Madhur Tulsiai Lecture 8: October 20, 2016 1 Applicatios of SVD: least squares approximatio We discuss aother applicatio of sigular value decompositio (SVD) of

More information

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion LETURE 4 No-liear trasverse motio Floquet trasformatio Harmoic aalysis-oe dimesioal resoaces Two-dimesioal resoaces No-liear trasverse motio No-liear field terms i the trajectory equatio: Trajectory equatio

More information

MIT Department of Chemistry 5.74, Spring 2005: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff

MIT Department of Chemistry 5.74, Spring 2005: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff MIT Departmet of Chemistry 5.74, Sprig 5: Itroductory Quatum Mechaics II Istructor: Professor Adrei Tomaoff p. 97 ABSORPTION SPECTRA OF MOLECULAR AGGREGATES The absorptio spectra of periodic arrays of

More information

Linearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation

Linearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation Liearly Idepedet Sets Bases p p c c p Review { v v vp} A set of vectors i a vector space is said to be liearly idepedet if the vector equatio cv + c v + + c has oly the trivial solutio = = { v v vp} The

More information

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + 62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of

More information

Chapter 2: Numerical Methods

Chapter 2: Numerical Methods Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,

More information

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at

More information

Review Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =

Review Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = = Review Problems ICME ad MS&E Refresher Course September 9, 0 Warm-up problems. For the followig matrices A = 0 B = C = AB = 0 fid all powers A,A 3,(which is A times A),... ad B,B 3,... ad C,C 3,... Solutio:

More information

Nonlinear regression

Nonlinear regression oliear regressio How to aalyse data? How to aalyse data? Plot! How to aalyse data? Plot! Huma brai is oe the most powerfull computatioall tools Works differetly tha a computer What if data have o liear

More information

MATH 10550, EXAM 3 SOLUTIONS

MATH 10550, EXAM 3 SOLUTIONS MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,

More information

The multiplicative structure of finite field and a construction of LRC

The multiplicative structure of finite field and a construction of LRC IERG6120 Codig for Distributed Storage Systems Lecture 8-06/10/2016 The multiplicative structure of fiite field ad a costructio of LRC Lecturer: Keeth Shum Scribe: Zhouyi Hu Notatios: We use the otatio

More information

MAT 271 Project: Partial Fractions for certain rational functions

MAT 271 Project: Partial Fractions for certain rational functions MAT 7 Project: Partial Fractios for certai ratioal fuctios Prerequisite kowledge: partial fractios from MAT 7, a very good commad of factorig ad complex umbers from Precalculus. To complete this project,

More information

Review Sheet for Final Exam

Review Sheet for Final Exam Sheet for ial To study for the exam, we suggest you look through the past review sheets, exams ad homework assigmets, ad idetify the topics that you most eed to work o. To help with this, the table give

More information

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions .03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio

More information

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Polynomial Functions. New Section 1 Page 1. A Polynomial function of degree n is written is the form:

Polynomial Functions. New Section 1 Page 1. A Polynomial function of degree n is written is the form: New Sectio 1 Page 1 A Polyomial fuctio of degree is writte is the form: 1 P x a x a x a x a x a x a 1 1 0 where is a o-egative iteger expoet ad a 0 ca oly take o values where a, a, a,..., a, a. 0 1 1.

More information

The Transition Dipole Moment

The Transition Dipole Moment The Trasitio Dipole Momet Iteractio of Light with Matter The probability that a molecule absorbs or emits light ad udergoes a trasitio from a iitial to a fial state is give by the Eistei coefficiet, B

More information