PHYS 705: Classical Mechanics. Small Oscillations: Example A Linear Triatomic Molecule

Similar documents
PHYS 601 HW3 Solution

PHY 5246: Theoretical Dynamics, Fall Assignment # 5, Solutions. θ = l mr 2 = l

OXFORD H i g h e r E d u c a t i o n Oxford University Press, All rights reserved.

The Spring. Consider a spring, which we apply a force F A to either stretch it or compress it

Matrix Eigenvalues and Eigenvectors September 13, 2017

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

Second degree generalized gauss-seidel iteration method for solving linear system of equations. ABSTRACT

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016

Summary: Method of Separation of Variables

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

Chapter 1. Linear Algebra

20 MATHEMATICS POLYNOMIALS

Best Approximation in the 2-norm

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Section 7.1 Integration by Substitution

dx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.

Chapter 2. Determinants

Formulae For. Standard Formulae Of Integrals: x dx k, n 1. log. a dx a k. cosec x.cot xdx cosec. e dx e k. sec. ax dx ax k. 1 1 a x.

Chapter 1 COUPLED ONE-DIMENSIONAL OSCILLATORS

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

The Wave Equation I. MA 436 Kurt Bryan

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

If deg(num) deg(denom), then we should use long-division of polynomials to rewrite: p(x) = s(x) + r(x) q(x), q(x)

1 2-D Second Order Equations: Separation of Variables

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

Laplace s equation in Cylindrical Coordinates

UCSD Phys 4A Intro Mechanics Winter 2016 Ch 4 Solutions

Polynomials and Division Theory

Math 1102: Calculus I (Math/Sci majors) MWF 3pm, Fulton Hall 230 Homework 2 solutions

PHYS 4390: GENERAL RELATIVITY LECTURE 6: TENSOR CALCULUS

CSCI 5525 Machine Learning

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

AMS 212A Applied Mathematical Methods I Lecture 06 Copyright by Hongyun Wang, UCSC. ( ), v (that is, 1 ( ) L i

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY

Chapter 3 The Schrödinger Equation and a Particle in a Box

Chapter (a) 2. At equilibrium, V = 0. x = (c) V( x) ke bkxe = 0. At equilibrium. b bx bx 2 bx. V = bke bkxe + b kxe = + = < 0.

Chapter 3 Solving Nonlinear Equations

Exponents and Powers

Multivariate problems and matrix algebra

Lesson 1: Quadratic Equations

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

Computer Graphics (CS 543) Lecture 3 (Part 1): Linear Algebra for Graphics (Points, Scalars, Vectors)

7-1: Zero and Negative Exponents

NORMALS. a y a y. Therefore, the slope of the normal is. a y1. b x1. b x. a b. x y a b. x y

In Section 5.3 we considered initial value problems for the linear second order equation. y.a/ C ˇy 0.a/ D k 1 (13.1.4)

Improper Integrals, and Differential Equations

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Integration Techniques

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Math 113 Exam 2 Practice

than 1. It means in particular that the function is decreasing and approaching the x-

Best Approximation. Chapter The General Case

Phys101 Lecture 4,5 Dynamics: Newton s Laws of Motion

Bridging the gap: GCSE AS Level

Complex, distinct eigenvalues (Sect. 7.6)

Quantum Physics II (8.05) Fall 2013 Assignment 2

Does the Order Matter?

r = cos θ + 1. dt ) dt. (1)

MAC-solutions of the nonexistent solutions of mathematical physics

Recitation 3: More Applications of the Derivative

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Lyapunov-type inequalities for Laplacian systems and applications to boundary value problems

Module 6: LINEAR TRANSFORMATIONS

Logarithmic Functions

MA 201: Partial Differential Equations Lecture - 12

Math Lecture 23

(4.1) D r v(t) ω(t, v(t))

SOLUTIONS TO CONCEPTS CHAPTER 6

7.3 Problem 7.3. ~B(~x) = ~ k ~ E(~x)=! but we also have a reected wave. ~E(~x) = ~ E 2 e i~ k 2 ~x i!t. ~B R (~x) = ~ k R ~ E R (~x)=!

Chapter Direct Method of Interpolation More Examples Chemical Engineering

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

Math Theory of Partial Differential Equations Lecture 2-9: Sturm-Liouville eigenvalue problems (continued).

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

Practice Problems Solution

Recitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

INTRODUCTION TO LINEAR ALGEBRA

MATH 573 FINAL EXAM. May 30, 2007

Unit 5. Integration techniques

Pressure Wave Analysis of a Cylindrical Drum

Math 270A: Numerical Linear Algebra

Chapter 6 Techniques of Integration

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

Math 231E, Lecture 33. Parametric Calculus

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Physics 215 Quantum Mechanics 1 Assignment 2

Fundamental Theorem of Calculus

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

Physics 137A - Quantum Mechanics - Spring 2018 Midterm 1. Mathematical Formulas

Lecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is

M A T H F A L L CORRECTION. Algebra I 2 1 / 0 9 / U N I V E R S I T Y O F T O R O N T O

Homework Problem Set 1 Solutions

u(x, y, t) = T(t)Φ(x, y) 0. (THE EQUATIONS FOR PRODUCT SOLUTIONS) Plugging u = T(t)Φ(x, y) in (PDE)-(BC) we see: There is a constant λ such that

Problem Set 2 Solutions

Transcription:

PHYS 75: Clssicl echnics Sll Oscilltions: Exple A Liner Tritoic olecule

A Liner Tritoic olecule x b b x x3 x Experientlly, one ight be interested in the rdition resulted fro the intrinsic oscilltion odes fro these tritoic olecule. The potentil energy for the two springs is, V x xb x3 x b

3 A Liner Tritoic olecule x b b x x3 x Now, we will introduce generlized coordintes reltive to their x, x, x equilibriu positions, : 3 x x j j j j,, 3 Note: x x x3 x b

4 A Liner Tritoic olecule Expnding the potentil energy bout its equilibriu position, we hve: V 3 ultiplying the squres out, we hve: x x b x x x x x x xx V 3 3 V 3 3 3 V V j j In trix for, this qudrtic for hs this for: V

5 A Liner Tritoic olecule Direct ethod to evlute using Recll, we hve V V j V j qjq V x xb x3 x b By ting the prtil derivtives directly nd evluting t, x j x j V V qq V : x xb nd q V : V q nd x x b x3 x b V qq V

6 A Liner Tritoic olecule The inetic energy for the three ss is given by: T x x x 3 Substituting our generlized coordintes or, x x j j j x j j T 3 T T j j In trix for, this qudrtic for hs this for: T

7 A Liner Tritoic olecule Cobining these two qudrtic fors into the chrcteristic eqution, VT V T Explicitly evluting this deterinnt, we hve the following eqution, 4

8 A Liner Tritoic olecule And, this hs three distinct solutions (eigenfrequencies): 3 Note: The solution ens tht the corresponding norl coordinte will hve the following trivil ODE: j or (unifor trnsltionl otion) The entire olecule will siply ove uniforly to the right or left; no oscilltions (not quite interesting otion by itself)

9 A Liner Tritoic olecule Now, we find the eigenvectors for ech solved eigenvlues using: V r T r r r r r r 3r Norl ode #: r 3 3

A Liner Tritoic olecule Norl ode #: r 3 3 3 Solving st nd 3 rd equtions, we hve nd Putting the together, we hve 3 (Note tht the solution lso stisfies the nd eqution.) 3 So, the eigenvector for is:

A Liner Tritoic olecule Norl ode #: r The eigenvector needs to be norlized with respect to T: This gives: T Finlly, we hve the norlized eignvector for :

A Liner Tritoic olecule Norl ode #: r 3 3 3

3 A Liner Tritoic olecule Norl ode #: r 3 Agin, we need to norlized with respect to T: T This gives:

4 A Liner Tritoic olecule Norl ode #3: r 3 3 3 3 3 3 3 33 First, let try to siply the trix eleents first: 3 3

5 A Liner Tritoic olecule Norl ode #3: r 3 Putting these vlues bc into the trix, we hve 3 3 33 3 3 3 3 33 333 3 3 3 3 3 33 33

6 A Liner Tritoic olecule Norl ode #3: r 3 3 3 33 3 3 3 Agin, we need to norlized with respect to T: T 3 3 3

7 A Liner Tritoic olecule Norl ode #3: r 3 3 4 3 3 This gives: 3 4 3

8 A Liner Tritoic olecule Then, our generl solution is given by: j jr r nd Ce Ce irt * irt r r r (the coplex coefficient will be deterined by IC) C r 3 3 3 3 3

9 Longitudinl Norl odes So, if single norl ode is ctive, the otion of the three generlized j r coordintes will loos lie the following,,, 3 :,, :,, 3: n,, 3 3 3

Longitudinl Norl odes nd Anitions for the two non-rigid-trnsltionl odes: 3 http://cis.poly.edu/~leung/cs4744/f3/ch4/l inole/linole.htl (fro Polytechnics Institute of NY Univ: K. ing Leung)

Sury. Pic generlized coordintes nd find nd. Expnd T nd V bout equilibriu q j Tq ( ) V( q j ) This gives two rel syetric qudrtic fors: with q q j j j j T ( ) V ( ) 3. Clculte eigenfrequencies fro chrcteristic eqution r det Vj r Tj 4. Clculte eigenvectors for ech eigenfrequencies using 5. Norlize eigenvectors with respect to T: 6. Generl solutions re in ters of the norl odes r * 7. Originl coordintes re relted bc through t C e C e r i t irt r r r Tjjrr j j j j V T j r j jr C deterined by IC r j jr r