MOLECULAR VIBRATIONS

Size: px
Start display at page:

Download "MOLECULAR VIBRATIONS"

Transcription

1 MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal harmoc oscllator s gve by the sum of the ketc ad potetal eergy. Suppose a partcle of mass m s movg a potetal gve by V (x) = 1 kx oe-dmeso, the, the Hamltoa s gve by H = 1 mẋ + 1 kx (1) where the frst term s the ketc eergy ad the secod the potetal eergy. Classcally, the partcle obeys Newto s Law, whch states that: mẍ = dv () dx gvg rse to the followg equato of moto upo substtutg the form of the potetal V (x) = 1 kx : mẍ = kx (3) whch solves to gve oscllatory solutos havg the form of s, cos or expoetal. We could choose the cose form: x(t) = A cos(ωt + φ) (4) where A ad φ are costats determed by tal codtos (sce Eq. (3) s a secod-order dfferetal equato, we expect two costats here) ad ω s the agular frequecy gve by k ω = (5) m Usg Eq. (5), we could wrte the Hamltoa as H = 1 mẋ + 1 mω x (6) 1

2 . Datomc molecule Now, cosder a datomc molecule havg atoms wth masses m 1 ad m teractg through a harmoc potetal descrbed wth a force costat k. Assumg that the moto of the two ucle s cofed to oe dmeso, the Hamltoa of the problem becomes H = 1 m x + 1 k(x 1 x ) (7) where x 1 ad x are the dsplacemet coordates relatve to the equlbrum geometry of the molecule (F. 1). Fgure 1. Dsplacemet coordates of a datomc molecule. We frst troduce a mass-weghted coordate: q = m x (8) to remove ay explct mass depedece from Eq. (7): H = 1 q + 1 ( k q1 q ) m1 m = + 1 ( k q k 1 q 1 q + k ) q m 1 m1 m m = 1 q + 1 q K j q j (9) where K s the dyamcal or Hessa matrx gve, ths case, by k k m K = 1 m1 m k k m1 m m j We ca fd the egevalues ad egevectors of the dyamcal matrx K by solvg The solutos ca be wrtte as j (10) K λi = 0 (11) K j j = λ (1) c Xglog Zhag 017

3 where the superscrpt refers to the th ormalsed egevector, ad λ the correspodg egevalue. Sce the matrx K s Hermta (deed K s symmetrc sce all ts compoets are real), ther egevectors are orthogoal, we ca also ormalsed them to gve orthoormal egevectors: c (m) = δ m (13) Now we wsh to fd ormal mode coordates, Q, terms of (a lear combato of) scaled atomc dsplacemets q, such that these ormal modes are decoupled from each other ad ca be aalysed separately. Expressed more formally, we wsh to be able to wrte the Hamltoa Eq. (9) a dagoal form: H = ( 1 Q + 1 ω Q ). (14) Sce the Hamltoa Eq. (9) cotas the atomc dsplacemets q, we cosder the followg coordate trasform: q = Q (15).e., wrtg the scaled coordates as a lear combato of the ormal mode coordates, where the coeffcets s the th compoet of the th ormalsed egevector from Eq. (1). Substtutg ths to the Hamltoa of the system Eq. (9), we have H = 1 q + 1 q K j q j j = 1 [ ] [ ] Q c (m) Q m + 1 [ ] [ ] Q K j c (m) Q m m j m = 1 [ ] Q c (m) Q m + 1 K j c (m) Q m m m j = 1 Q δ m Q m + 1 λ m δ m Q m m m = 1 Q + 1 λ Q (16) where achevg the secod last le, we used the equatos from Eq. (1) ad (13). Ths ow deed has the dagoal form of Eq. (14) whch we detfy the egevalues from solvg the dyamcal matrx as λ = ω (17) c Xglog Zhag 017 3

4 At ths stage, after solvg for the egevalues ad the egevectors of the dyamcal matrx, we are fally able to wrte the ormal mode coordates as a lear combato of the (scaled) atomc dsplacemets. We have the coordate trasform Eq. (15), left multply by traspose of aother egevector ad usg the orthoormalty codto of Eq. (13), we have ] c (m) q = Q = c (m) Q = [ c (m) Q = δ m Q = Q m q (free relabellg of m ) (18) Now, the above equato wrte the ormal mode coordates as a lear combato of the mass-scaled atomc dsplacemets so that each ormal mode s decoupled from the other ormal modes (Ref Eq. 16). To solve the system of molecular motos for the datomc system, oe s the requred to solve the dagoalsato of the dyamcal matrx Eq. (10). Ths s left to the readers as a smple exercse. 3. Geeralsato to polyatomcs Cosder a molecule cosstg of N atoms. Classcally, the Hamltoa for molecular vbratos has the geeral form H = 1 3N m Ṙ + V (R) (19) where R s the posto vector for atom, m s the mass of the atom wth compoets R ( = 1,, 3N). The postos of all ucle ca be wrtte as The quatum expresso s gve stead by Ĥ = 1 R {R 1, R,, R N } (0) 1 ˆp + V (R) = 1 1 m m R + V (R) (1) However, we wll ote that usg the mass-weghted coordate Eq. (8), we ca get the dyamcal matrx usg the potetal term V (R) aloe. Suppose the molecule s at equlbrum geometry, deoted by R (0) {R (0) 1, R (0),, R (0) N }, () we ca defe ay uclear motos relatve to the equlbrum geometry by R = R (0) + x (3) c Xglog Zhag 017 4

5 where x s a collecto of 3 coordate dsplacemets for ucleus : x = (x 1, x, x 3 ). (4) The codto for equlbrum geometry for a molecule s gve by ( ) V = 0 (5) x where the subscrpt 0 deotes that the partal dervatves are evaluated at equlbrum geometry. For small-ampltude dsplacemets aroud the equlbrum, we ca Taylor-expad the potetal eergy as ) V (R) = V ( V j x x j 0 x x j + cubc terms x x j x k + (6) 0 where V 0 = V (R (0) ). Note that the lear term vashes due to equlbrum codto. I the harmoc approxmato, the Taylor expaso above s trucated after the quadratc terms. Comparg the quadratc term above wth the potetal eergy term (the secod term) Eq. (9), ad usg the mass-weghted coordates Eq. (8), we ca see that 1 the dyamcal matrx s gve by ( 1 ) V K j = (7) m m j x x j 0 ( ) V = (8) q q j We ow have a geeral method for solvg for ormal modes of a system wth a potetal eergy V. We ca costruct the dyamcal matrx from t usg Eq. (8). We the solve for ts egevalues ad egevectors to gve the (square) of the vbratoal frequeces (Eq. 17) ad the ormal mode coordates (Eq. 18). Ths s deed the same approach we use solvg the Hückel secular equatos. Of course, symmetry could be used to smplfy the solvg of these matrces. 4. Parallels betwee Hückel theory ad molecular vbratos aalyss Although the two areas appear to be very dfferet, actual fact, the ways to solve these problems are very smlar. We ote the followg resemblace betwee Hückel theory ad molecular vbratos aalyss: 0 1 you could do a chage of varable dfferetato by usg cha rule c Xglog Zhag 017 5

6 Hückel Theory atomc orbtals, φ molecular orbtals, ψ Hückel matrx, H j orbtal eerges, ɛ symmetry orbtals, θ Γ mass-weghted coordates, q ormal mode coordates, Q dyamcal matrx, K j squared frequeces, ω symmetry-adapted lear combato of dsplacemets, Q Γ Table 1. Close resemblace betwee Hückel theory ad molecular vbratos. c Xglog Zhag 017 6

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices. 4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal

More information

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh

More information

Chapter 9 Jordan Block Matrices

Chapter 9 Jordan Block Matrices Chapter 9 Jorda Block atrces I ths chapter we wll solve the followg problem. Gve a lear operator T fd a bass R of F such that the matrx R (T) s as smple as possble. f course smple s a matter of taste.

More information

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix Assgmet 7/MATH 47/Wter, 00 Due: Frday, March 9 Powers o a square matrx Gve a square matrx A, ts powers A or large, or eve arbtrary, teger expoets ca be calculated by dagoalzg A -- that s possble (!) Namely,

More information

Lecture 07: Poles and Zeros

Lecture 07: Poles and Zeros Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto

More information

Ahmed Elgamal. MDOF Systems & Modal Analysis

Ahmed Elgamal. MDOF Systems & Modal Analysis DOF Systems & odal Aalyss odal Aalyss (hese otes cover sectos from Ch. 0, Dyamcs of Structures, Al Chopra, Pretce Hall, 995). Refereces Dyamcs of Structures, Al K. Chopra, Pretce Hall, New Jersey, ISBN

More information

TESTS BASED ON MAXIMUM LIKELIHOOD

TESTS BASED ON MAXIMUM LIKELIHOOD ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal

More information

Numerical Analysis Formulae Booklet

Numerical Analysis Formulae Booklet Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx

More information

Lecture Note to Rice Chapter 8

Lecture Note to Rice Chapter 8 ECON 430 HG revsed Nov 06 Lecture Note to Rce Chapter 8 Radom matrces Let Y, =,,, m, =,,, be radom varables (r.v. s). The matrx Y Y Y Y Y Y Y Y Y Y = m m m s called a radom matrx ( wth a ot m-dmesoal dstrbuto,

More information

Chapter 3. Linear Equations and Matrices

Chapter 3. Linear Equations and Matrices Vector Spaces Physcs 8/6/05 hapter Lear Equatos ad Matrces wde varety of physcal problems volve solvg systems of smultaeous lear equatos These systems of lear equatos ca be ecoomcally descrbed ad effcetly

More information

4 Inner Product Spaces

4 Inner Product Spaces 11.MH1 LINEAR ALGEBRA Summary Notes 4 Ier Product Spaces Ier product s the abstracto to geeral vector spaces of the famlar dea of the scalar product of two vectors or 3. I what follows, keep these key

More information

DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS

DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS Course Project: Classcal Mechacs (PHY 40) Suja Dabholkar (Y430) Sul Yeshwath (Y444). Itroducto Hamltoa mechacs s geometry phase space. It deals

More information

Engineering Vibration 1. Introduction

Engineering Vibration 1. Introduction Egeerg Vbrato. Itroducto he study of the moto of physcal systems resultg from the appled forces s referred to as dyamcs. Oe type of dyamcs of physcal systems s vbrato, whch the system oscllates about certa

More information

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution: Chapter 4 Exercses Samplg Theory Exercse (Smple radom samplg: Let there be two correlated radom varables X ad A sample of sze s draw from a populato by smple radom samplg wthout replacemet The observed

More information

Non-degenerate Perturbation Theory

Non-degenerate Perturbation Theory No-degeerate Perturbato Theory Proble : H E ca't solve exactly. But wth H H H' H" L H E Uperturbed egevalue proble. Ca solve exactly. E Therefore, kow ad. H ' H" called perturbatos Copyrght Mchael D. Fayer,

More information

ENGI 4421 Propagation of Error Page 8-01

ENGI 4421 Propagation of Error Page 8-01 ENGI 441 Propagato of Error Page 8-01 Propagato of Error [Navd Chapter 3; ot Devore] Ay realstc measuremet procedure cotas error. Ay calculatos based o that measuremet wll therefore also cota a error.

More information

III-16 G. Brief Review of Grand Orthogonality Theorem and impact on Representations (Γ i ) l i = h n = number of irreducible representations.

III-16 G. Brief Review of Grand Orthogonality Theorem and impact on Representations (Γ i ) l i = h n = number of irreducible representations. III- G. Bref evew of Grad Orthogoalty Theorem ad mpact o epresetatos ( ) GOT: h [ () m ] [ () m ] δδ δmm ll GOT puts great restrcto o form of rreducble represetato also o umber: l h umber of rreducble

More information

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace

More information

ESS Line Fitting

ESS Line Fitting ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here

More information

x y exp λ'. x exp λ 2. x exp 1.

x y exp λ'. x exp λ 2. x exp 1. egecosmcd Egevalue-egevector of the secod dervatve operator d /d hs leads to Fourer seres (se, cose, Legedre, Bessel, Chebyshev, etc hs s a eample of a systematc way of geeratg a set of mutually orthogoal

More information

Investigation of Partially Conditional RP Model with Response Error. Ed Stanek

Investigation of Partially Conditional RP Model with Response Error. Ed Stanek Partally Codtoal Radom Permutato Model 7- vestgato of Partally Codtoal RP Model wth Respose Error TRODUCTO Ed Staek We explore the predctor that wll result a smple radom sample wth respose error whe a

More information

CHAPTER 4 RADICAL EXPRESSIONS

CHAPTER 4 RADICAL EXPRESSIONS 6 CHAPTER RADICAL EXPRESSIONS. The th Root of a Real Number A real umber a s called the th root of a real umber b f Thus, for example: s a square root of sce. s also a square root of sce ( ). s a cube

More information

1 Mixed Quantum State. 2 Density Matrix. CS Density Matrices, von Neumann Entropy 3/7/07 Spring 2007 Lecture 13. ψ = α x x. ρ = p i ψ i ψ i.

1 Mixed Quantum State. 2 Density Matrix. CS Density Matrices, von Neumann Entropy 3/7/07 Spring 2007 Lecture 13. ψ = α x x. ρ = p i ψ i ψ i. CS 94- Desty Matrces, vo Neuma Etropy 3/7/07 Sprg 007 Lecture 3 I ths lecture, we wll dscuss the bascs of quatum formato theory I partcular, we wll dscuss mxed quatum states, desty matrces, vo Neuma etropy

More information

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat

More information

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs

More information

1 Lyapunov Stability Theory

1 Lyapunov Stability Theory Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may

More information

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions. Ordary Least Squares egresso. Smple egresso. Algebra ad Assumptos. I ths part of the course we are gog to study a techque for aalysg the lear relatoshp betwee two varables Y ad X. We have pars of observatos

More information

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

We have already referred to a certain reaction, which takes place at high temperature after rich combustion. ME 41 Day 13 Topcs Chemcal Equlbrum - Theory Chemcal Equlbrum Example #1 Equlbrum Costats Chemcal Equlbrum Example #2 Chemcal Equlbrum of Hot Bured Gas 1. Chemcal Equlbrum We have already referred to a

More information

QR Factorization and Singular Value Decomposition COS 323

QR Factorization and Singular Value Decomposition COS 323 QR Factorzato ad Sgular Value Decomposto COS 33 Why Yet Aother Method? How do we solve least-squares wthout currg codto-squarg effect of ormal equatos (A T A A T b) whe A s sgular, fat, or otherwse poorly-specfed?

More information

Functions of Random Variables

Functions of Random Variables Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,

More information

7.0 Equality Contraints: Lagrange Multipliers

7.0 Equality Contraints: Lagrange Multipliers Systes Optzato 7.0 Equalty Cotrats: Lagrage Multplers Cosder the zato of a o-lear fucto subject to equalty costrats: g f() R ( ) 0 ( ) (7.) where the g ( ) are possbly also olear fuctos, ad < otherwse

More information

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Numercal Computg -I UNIT SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Structure Page Nos..0 Itroducto 6. Objectves 7. Ital Approxmato to a Root 7. Bsecto Method 8.. Error Aalyss 9.4 Regula Fals Method

More information

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class) Assgmet 5/MATH 7/Wter 00 Due: Frday, February 9 class (!) (aswers wll be posted rght after class) As usual, there are peces of text, before the questos [], [], themselves. Recall: For the quadratc form

More information

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory ROAD MAP... AE301 Aerodyamcs I UNIT C: 2-D Arfols C-1: Aerodyamcs of Arfols 1 C-2: Aerodyamcs of Arfols 2 C-3: Pael Methods C-4: Th Arfol Theory AE301 Aerodyamcs I Ut C-3: Lst of Subects Problem Solutos?

More information

Manipulator Dynamics. Amirkabir University of Technology Computer Engineering & Information Technology Department

Manipulator Dynamics. Amirkabir University of Technology Computer Engineering & Information Technology Department Mapulator Dyamcs mrkabr Uversty of echology omputer Egeerg formato echology Departmet troducto obot arm dyamcs deals wth the mathematcal formulatos of the equatos of robot arm moto. hey are useful as:

More information

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg

More information

Lecture Notes 2. The ability to manipulate matrices is critical in economics.

Lecture Notes 2. The ability to manipulate matrices is critical in economics. Lecture Notes. Revew of Matrces he ablt to mapulate matrces s crtcal ecoomcs.. Matr a rectagular arra of umbers, parameters, or varables placed rows ad colums. Matrces are assocated wth lear equatos. lemets

More information

ρ < 1 be five real numbers. The

ρ < 1 be five real numbers. The Lecture o BST 63: Statstcal Theory I Ku Zhag, /0/006 Revew for the prevous lecture Deftos: covarace, correlato Examples: How to calculate covarace ad correlato Theorems: propertes of correlato ad covarace

More information

α1 α2 Simplex and Rectangle Elements Multi-index Notation of polynomials of degree Definition: The set P k will be the set of all functions:

α1 α2 Simplex and Rectangle Elements Multi-index Notation of polynomials of degree Definition: The set P k will be the set of all functions: Smplex ad Rectagle Elemets Mult-dex Notato = (,..., ), o-egatve tegers = = β = ( β,..., β ) the + β = ( + β,..., + β ) + x = x x x x = x x β β + D = D = D D x x x β β Defto: The set P of polyomals of degree

More information

Solutions to problem set ); (, ) (

Solutions to problem set ); (, ) ( Solutos to proble set.. L = ( yp p ); L = ( p p ); y y L, L = yp p, p p = yp p, + p [, p ] y y y = yp + p = L y Here we use for eaple that yp, p = yp p p yp = yp, p = yp : factors that coute ca be treated

More information

Lecture IV : The Hartree-Fock method

Lecture IV : The Hartree-Fock method Lecture IV : The Hartree-Fock method I. THE HARTREE METHOD We have see the prevous lecture that the may-body Hamltoa for a electroc system may be wrtte atomc uts as Ĥ = N e N e N I Z I r R I + N e N e

More information

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK Far East Joural of Appled Mathematcs Volume, Number, 2008, Pages Ths paper s avalable ole at http://www.pphm.com 2008 Pushpa Publshg House ANALYSIS ON THE NATURE OF THE ASI EQUATIONS IN SYNERGETI INTER-REPRESENTATION

More information

Point Estimation: definition of estimators

Point Estimation: definition of estimators Pot Estmato: defto of estmators Pot estmator: ay fucto W (X,..., X ) of a data sample. The exercse of pot estmato s to use partcular fuctos of the data order to estmate certa ukow populato parameters.

More information

Stationary states of atoms and molecules

Stationary states of atoms and molecules Statoary states of atos ad olecules I followg wees the geeral aspects of the eergy level structure of atos ad olecules that are essetal for the terpretato ad the aalyss of spectral postos the rotatoal

More information

CIS 800/002 The Algorithmic Foundations of Data Privacy October 13, Lecture 9. Database Update Algorithms: Multiplicative Weights

CIS 800/002 The Algorithmic Foundations of Data Privacy October 13, Lecture 9. Database Update Algorithms: Multiplicative Weights CIS 800/002 The Algorthmc Foudatos of Data Prvacy October 13, 2011 Lecturer: Aaro Roth Lecture 9 Scrbe: Aaro Roth Database Update Algorthms: Multplcatve Weghts We ll recall aga) some deftos from last tme:

More information

Decomposition of Hadamard Matrices

Decomposition of Hadamard Matrices Chapter 7 Decomposto of Hadamard Matrces We hae see Chapter that Hadamard s orgal costructo of Hadamard matrces states that the Kroecer product of Hadamard matrces of orders m ad s a Hadamard matrx of

More information

Ideal multigrades with trigonometric coefficients

Ideal multigrades with trigonometric coefficients Ideal multgrades wth trgoometrc coeffcets Zarathustra Brady December 13, 010 1 The problem A (, k) multgrade s defed as a par of dstct sets of tegers such that (a 1,..., a ; b 1,..., b ) a j = =1 for all

More information

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law Module : The equato of cotuty Lecture 5: Coservato of Mass for each speces & Fck s Law NPTEL, IIT Kharagpur, Prof. Sakat Chakraborty, Departmet of Chemcal Egeerg 2 Basc Deftos I Mass Trasfer, we usually

More information

GG313 GEOLOGICAL DATA ANALYSIS

GG313 GEOLOGICAL DATA ANALYSIS GG33 GEOLOGICAL DATA ANALYSIS 3 GG33 GEOLOGICAL DATA ANALYSIS LECTURE NOTES PAUL WESSEL SECTION 3 LINEAR (MATRIX ALGEBRA OVERVIEW OF MATRIX ALGEBRA (or All you ever wated to kow about Lear Algebra but

More information

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy Bouds o the expected etropy ad KL-dvergece of sampled multomal dstrbutos Brado C. Roy bcroy@meda.mt.edu Orgal: May 18, 2011 Revsed: Jue 6, 2011 Abstract Iformato theoretc quattes calculated from a sampled

More information

X X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then

X X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then Secto 5 Vectors of Radom Varables Whe workg wth several radom varables,,..., to arrage them vector form x, t s ofte coveet We ca the make use of matrx algebra to help us orgaze ad mapulate large umbers

More information

Solutions to Homework Problems for the Complexity Explorer Course on Random Walks

Solutions to Homework Problems for the Complexity Explorer Course on Random Walks Solutos to Homework Problems for the Complexty Explorer Course o Radom Walks. Dsplacemet of a radom walk. Cosder the Pearso radom walk ay spatal dmeso whch the legth of each step has the fxed value a,

More information

Lecture 3 Probability review (cont d)

Lecture 3 Probability review (cont d) STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto

More information

MEASURES OF DISPERSION

MEASURES OF DISPERSION MEASURES OF DISPERSION Measure of Cetral Tedecy: Measures of Cetral Tedecy ad Dsperso ) Mathematcal Average: a) Arthmetc mea (A.M.) b) Geometrc mea (G.M.) c) Harmoc mea (H.M.) ) Averages of Posto: a) Meda

More information

Physics 114 Exam 2 Fall Name:

Physics 114 Exam 2 Fall Name: Physcs 114 Exam Fall 015 Name: For gradg purposes (do ot wrte here): Questo 1. 1... 3. 3. Problem Aswer each of the followg questos. Pots for each questo are dcated red. Uless otherwse dcated, the amout

More information

16 Homework lecture 16

16 Homework lecture 16 Quees College, CUNY, Departmet of Computer Scece Numercal Methods CSCI 361 / 761 Fall 2018 Istructor: Dr. Sateesh Mae c Sateesh R. Mae 2018 16 Homework lecture 16 Please emal your soluto, as a fle attachmet,

More information

Harley Flanders Differential Forms with Applications to the Physical Sciences. Dover, 1989 (1962) Contents FOREWORD

Harley Flanders Differential Forms with Applications to the Physical Sciences. Dover, 1989 (1962) Contents FOREWORD Harley Fladers Dfferetal Forms wth Applcatos to the Physcal Sceces FORWORD Dover, 989 (962) Cotets PRFAC TO TH DOVR DITION PRFAC TO TH FIRST DITION.. xteror Dfferetal Forms.2. Comparso wth Tesors 2.. The

More information

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION Bars Erkus, 4 March 007 Itroducto Ths docuet provdes a revew of fudaetal cocepts structural dyacs ad soe applcatos hua-duced vbrato aalyss ad tgato of

More information

Newton s Power Flow algorithm

Newton s Power Flow algorithm Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other

More information

Third handout: On the Gini Index

Third handout: On the Gini Index Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The

More information

On the construction of symmetric nonnegative matrix with prescribed Ritz values

On the construction of symmetric nonnegative matrix with prescribed Ritz values Joural of Lear ad Topologcal Algebra Vol. 3, No., 14, 61-66 O the costructo of symmetrc oegatve matrx wth prescrbed Rtz values A. M. Nazar a, E. Afshar b a Departmet of Mathematcs, Arak Uversty, P.O. Box

More information

General Method for Calculating Chemical Equilibrium Composition

General Method for Calculating Chemical Equilibrium Composition AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt

More information

Chapter 4 Multiple Random Variables

Chapter 4 Multiple Random Variables Revew for the prevous lecture: Theorems ad Examples: How to obta the pmf (pdf) of U = g (, Y) ad V = g (, Y) Chapter 4 Multple Radom Varables Chapter 44 Herarchcal Models ad Mxture Dstrbutos Examples:

More information

L5 Polynomial / Spline Curves

L5 Polynomial / Spline Curves L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a

More information

1 Convergence of the Arnoldi method for eigenvalue problems

1 Convergence of the Arnoldi method for eigenvalue problems Lecture otes umercal lear algebra Arold method covergece Covergece of the Arold method for egevalue problems Recall that, uless t breaks dow, k steps of the Arold method geerates a orthogoal bass of a

More information

Beam Warming Second-Order Upwind Method

Beam Warming Second-Order Upwind Method Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet

More information

Summary of the lecture in Biostatistics

Summary of the lecture in Biostatistics Summary of the lecture Bostatstcs Probablty Desty Fucto For a cotuos radom varable, a probablty desty fucto s a fucto such that: 0 dx a b) b a dx A probablty desty fucto provdes a smple descrpto of the

More information

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM Jose Javer Garca Moreta Ph. D research studet at the UPV/EHU (Uversty of Basque coutry) Departmet of Theoretcal

More information

PHYS Look over. examples 2, 3, 4, 6, 7, 8,9, 10 and 11. How To Make Physics Pay PHYS Look over. Examples: 1, 4, 5, 6, 7, 8, 9, 10,

PHYS Look over. examples 2, 3, 4, 6, 7, 8,9, 10 and 11. How To Make Physics Pay PHYS Look over. Examples: 1, 4, 5, 6, 7, 8, 9, 10, PHYS Look over Chapter 9 Sectos - Eamples:, 4, 5, 6, 7, 8, 9, 0, PHYS Look over Chapter 7 Sectos -8 8, 0 eamples, 3, 4, 6, 7, 8,9, 0 ad How To ake Phscs Pa We wll ow look at a wa of calculatg where the

More information

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES FDM: Appromato of Frst Order Dervatves Lecture APPROXIMATION OF FIRST ORDER DERIVATIVES. INTRODUCTION Covectve term coservato equatos volve frst order dervatves. The smplest possble approach for dscretzato

More information

6.867 Machine Learning

6.867 Machine Learning 6.867 Mache Learg Problem set Due Frday, September 9, rectato Please address all questos ad commets about ths problem set to 6.867-staff@a.mt.edu. You do ot eed to use MATLAB for ths problem set though

More information

Rademacher Complexity. Examples

Rademacher Complexity. Examples Algorthmc Foudatos of Learg Lecture 3 Rademacher Complexty. Examples Lecturer: Patrck Rebesch Verso: October 16th 018 3.1 Itroducto I the last lecture we troduced the oto of Rademacher complexty ad showed

More information

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i CHEMICAL EQUILIBRIA The Thermodyamc Equlbrum Costat Cosder a reversble reacto of the type 1 A 1 + 2 A 2 + W m A m + m+1 A m+1 + Assgg postve values to the stochometrc coeffcets o the rght had sde ad egatve

More information

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each 01 Log1 Cotest Roud Theta Complex Numbers 1 Wrte a b Wrte a b form: 1 5 form: 1 5 4 pots each Wrte a b form: 65 4 4 Evaluate: 65 5 Determe f the followg statemet s always, sometmes, or ever true (you may

More information

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Joural of Appled Matematcs ad Computatoal Mecacs 04, 3(4), 7-34 FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Ata Cekot, Stasław Kukla Isttute of Matematcs, Czestocowa Uversty of Tecology Częstocowa,

More information

1 Solution to Problem 6.40

1 Solution to Problem 6.40 1 Soluto to Problem 6.40 (a We wll wrte T τ (X 1,...,X where the X s are..d. wth PDF f(x µ, σ 1 ( x µ σ g, σ where the locato parameter µ s ay real umber ad the scale parameter σ s > 0. Lettg Z X µ σ we

More information

n -dimensional vectors follow naturally from the one

n -dimensional vectors follow naturally from the one B. Vectors ad sets B. Vectors Ecoomsts study ecoomc pheomea by buldg hghly stylzed models. Uderstadg ad makg use of almost all such models requres a hgh comfort level wth some key mathematcal sklls. I

More information

Unit 9. The Tangent Bundle

Unit 9. The Tangent Bundle Ut 9. The Taget Budle ========================================================================================== ---------- The taget sace of a submafold of R, detfcato of taget vectors wth dervatos at

More information

MATH 247/Winter Notes on the adjoint and on normal operators.

MATH 247/Winter Notes on the adjoint and on normal operators. MATH 47/Wter 00 Notes o the adjot ad o ormal operators I these otes, V s a fte dmesoal er product space over, wth gve er * product uv, T, S, T, are lear operators o V U, W are subspaces of V Whe we say

More information

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006 Mult-Layer Refracto Problem Rafael Espercueta, Bakersfeld College, November, 006 Lght travels at dfferet speeds through dfferet meda, but refracts at layer boudares order to traverse the least-tme path.

More information

Mu Sequences/Series Solutions National Convention 2014

Mu Sequences/Series Solutions National Convention 2014 Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed

More information

AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET

AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET Abstract. The Permaet versus Determat problem s the followg: Gve a matrx X of determates over a feld of characterstc dfferet from

More information

Transforms that are commonly used are separable

Transforms that are commonly used are separable Trasforms s Trasforms that are commoly used are separable Eamples: Two-dmesoal DFT DCT DST adamard We ca the use -D trasforms computg the D separable trasforms: Take -D trasform of the rows > rows ( )

More information

CHAPTER 3 POSTERIOR DISTRIBUTIONS

CHAPTER 3 POSTERIOR DISTRIBUTIONS CHAPTER 3 POSTERIOR DISTRIBUTIONS If scece caot measure the degree of probablt volved, so much the worse for scece. The practcal ma wll stck to hs apprecatve methods utl t does, or wll accept the results

More information

Arithmetic Mean and Geometric Mean

Arithmetic Mean and Geometric Mean Acta Mathematca Ntresa Vol, No, p 43 48 ISSN 453-6083 Arthmetc Mea ad Geometrc Mea Mare Varga a * Peter Mchalča b a Departmet of Mathematcs, Faculty of Natural Sceces, Costate the Phlosopher Uversty Ntra,

More information

Dynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load

Dynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load Dyamc Aalyss of Axally Beam o Vsco - Elastc Foudato wth Elastc Supports uder Movg oad Saeed Mohammadzadeh, Seyed Al Mosayeb * Abstract: For dyamc aalyses of ralway track structures, the algorthm of soluto

More information

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model Lecture 7. Cofdece Itervals ad Hypothess Tests the Smple CLR Model I lecture 6 we troduced the Classcal Lear Regresso (CLR) model that s the radom expermet of whch the data Y,,, K, are the outcomes. The

More information

Econometric Methods. Review of Estimation

Econometric Methods. Review of Estimation Ecoometrc Methods Revew of Estmato Estmatg the populato mea Radom samplg Pot ad terval estmators Lear estmators Ubased estmators Lear Ubased Estmators (LUEs) Effcecy (mmum varace) ad Best Lear Ubased Estmators

More information

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler Mathematcal ad Computatoal Applcatos, Vol. 8, No. 3, pp. 293-300, 203 BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS Aysegul Ayuz Dascoglu ad Nese Isler Departmet of Mathematcs,

More information

3D Geometry for Computer Graphics. Lesson 2: PCA & SVD

3D Geometry for Computer Graphics. Lesson 2: PCA & SVD 3D Geometry for Computer Graphcs Lesso 2: PCA & SVD Last week - egedecomposto We wat to lear how the matrx A works: A 2 Last week - egedecomposto If we look at arbtrary vectors, t does t tell us much.

More information

Applied Fitting Theory VII. Building Virtual Particles

Applied Fitting Theory VII. Building Virtual Particles Appled Fttg heory II Paul Avery CBX 98 38 Jue 8, 998 Apr. 7, 999 (rev.) Buldg rtual Partcles I Statemet of the problem I may physcs aalyses we ecouter the problem of mergg a set of partcles to a sgle partcle

More information

Lecture 9: Tolerant Testing

Lecture 9: Tolerant Testing Lecture 9: Tolerat Testg Dael Kae Scrbe: Sakeerth Rao Aprl 4, 07 Abstract I ths lecture we prove a quas lear lower boud o the umber of samples eeded to do tolerat testg for L dstace. Tolerat Testg We have

More information

Generalized Linear Regression with Regularization

Generalized Linear Regression with Regularization Geeralze Lear Regresso wth Regularzato Zoya Bylsk March 3, 05 BASIC REGRESSION PROBLEM Note: I the followg otes I wll make explct what s a vector a what s a scalar usg vec t or otato, to avo cofuso betwee

More information

Matricial Potentiation

Matricial Potentiation Matrcal Potetato By Ezo March* ad Mart Mates** Abstract I ths short ote we troduce the potetato of matrces of the same sze. We study some smple propertes ad some example. * Emertus Professor UNSL, Sa Lus

More information

h-analogue of Fibonacci Numbers

h-analogue of Fibonacci Numbers h-aalogue of Fboacc Numbers arxv:090.0038v [math-ph 30 Sep 009 H.B. Beaoum Prce Mohammad Uversty, Al-Khobar 395, Saud Araba Abstract I ths paper, we troduce the h-aalogue of Fboacc umbers for o-commutatve

More information

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I MMJ FINITE EEMENT METHOD Cotut requremets Assume that the fuctos appearg uder the tegral the elemet equatos cota up to (r) th order To esure covergece N must satsf Compatblt requremet the fuctos must have

More information

Lecture 3. Sampling, sampling distributions, and parameter estimation

Lecture 3. Sampling, sampling distributions, and parameter estimation Lecture 3 Samplg, samplg dstrbutos, ad parameter estmato Samplg Defto Populato s defed as the collecto of all the possble observatos of terest. The collecto of observatos we take from the populato s called

More information

Analysis of Lagrange Interpolation Formula

Analysis of Lagrange Interpolation Formula P IJISET - Iteratoal Joural of Iovatve Scece, Egeerg & Techology, Vol. Issue, December 4. www.jset.com ISS 348 7968 Aalyss of Lagrage Iterpolato Formula Vjay Dahya PDepartmet of MathematcsMaharaja Surajmal

More information

Class 13,14 June 17, 19, 2015

Class 13,14 June 17, 19, 2015 Class 3,4 Jue 7, 9, 05 Pla for Class3,4:. Samplg dstrbuto of sample mea. The Cetral Lmt Theorem (CLT). Cofdece terval for ukow mea.. Samplg Dstrbuto for Sample mea. Methods used are based o CLT ( Cetral

More information

Application of Legendre Bernstein basis transformations to degree elevation and degree reduction

Application of Legendre Bernstein basis transformations to degree elevation and degree reduction Computer Aded Geometrc Desg 9 79 78 www.elsever.com/locate/cagd Applcato of Legedre Berste bass trasformatos to degree elevato ad degree reducto Byug-Gook Lee a Yubeom Park b Jaechl Yoo c a Dvso of Iteret

More information