Lagrangian Theory. Several-body Systems
|
|
- Sherman Chapman
- 5 years ago
- Views:
Transcription
1 Lagrangan Theory of Several-body Systems Ncholas Wheeler, Reed College Physcs Department November 995 Introducton. Let the N-tuple of 3-vectors {x (t) : =, 2,..., N} descrbe, relatve to an nertal frame, the confguraton of an N-partcle system at tme t. To descrbe the dynamcs of the system we would fnd t natural to ntroduce the Lagrangan L = 2 m ẋ ẋ U(x, x 2,..., x N ) () Suppose, however, we had elected to descrbe the partcles n terms of ther Cartesan relatonshp to a reference pont X(t) n arbtrarly prescrbed moton: x = X + r. We would then have L = 2 m (Ẋ + ṙ ) (Ẋ + ṙ ) U(X + r) = 2 MẊ Ẋ + Ẋ m ṙ + 2 m ṙ ṙ U(X + r) (2) n whch the dynamcal varables formerly {x, x 2,..., x N }, presently {r, r 2,..., r N } are stll (as before, and not at all surprsngly) N n number. The equatons of moton now read m r = m Ẍ U(X + r) ( =, 2,..., N) (3) n whch we nterpret the frst term on the rght to be a fcttous force term, an artfact of the crcumstance that the orgn of the X-centered r-frame s (except when Ẍ = 0) non-nertal. None of whch s n any respect problematc.
2 2 Lagrangan theory of several body systems Now wth an eye to the algebrac smplfcaton of (2) mpose upon the r-varables ths constrant: m r = 0 (4.) Equvalently, assocate X wth the center of mass of the N-partcle system: X = m x wth M = M m (4.2) Equaton (2) can now be notated L = 2 MẊ Ẋ + 2 m ṙ ṙ U(X + r) (5) (note the dsappearance of the cross-term), whch on ts face appears to refer to a system wth dynamcal varables {X, r, r 2,..., r N } more numerous that the varables of the system wth whch we started. And (5) gves rse to equatons of moton whch are not only more numerous than but also nconsstent wth equatons (3). What s gone wrong? In (4) we see that when we subjected the varables x to a cross-term-kllng constrant we effectvely promoted X to the status of a dynamcal varable; t s, accordng to (4.2), a varable a collectve varable whose t-dependence has now to be extracted from equatons of moton, and can no longer be sad to be arbtrarly prescrbed. But the Lagrangan (5) s heedless of ths crcumstance, and gves rse to equatons of moton the solutons of whch wll, n general, stand n volaton of (4.). To obtan correct results we mght, for example, ntroduce N N r N = m N m r and ṙ N = m N m ṙ nto (5) to obtan a Lagrangan of type L(Ẋ, ṙ,..., ṙ N, X, r,..., r N ), but such a procedure bears the formal blemsh of a dscrmnatory asymmetry not natural to the physcs of the stuaton. How to proceed more symmetrcally? In place of (5) wrte L = 2 MẊ Ẋ + 2 m ṙ ṙ U(X + r) g m r (6) where g s a Lagrange multpler whch wll be accorded the formal status of a supernumerary dynamcal varable. The resultng equatons of moton read MẌ = U (7.) m r = U gm ( =, 2,..., N) (7.2) 0 = m r (7.3)
3 A smple example: the one-dmensonal A 2 molecule 3 where =. These are N + 2 equatons of moton n as many varables. The last of these the g th Lagrange equaton of moton (7.3) s smply the constrant relaton (4.), and entals m r = 0. Addng equatons (7.2) together, and subtractng the result from (7.), we obtan The equatons of moton (7) reduce therefore to g = Ẍ (8) MẌ = U (9.) m { r + Ẍ} = U ( =, 2,..., N) (9.2) whch are attractvely symmetrc (no r has been dscrmnated aganst), but redundant: addng equatons (9.2) together gves back (9.). Equatons (9.2) are consstent wth (3), of whch they are a partcularzed nstance. It s nstructve to note also that the ntroducton of (8) nto (6) yelds a Lagrangan whch s dstnct from but (by Ẍ m r = d dt [Ẋ m r ] Ẋ m ṙ ) gauge equvalent to the Lagrangan of (2). A smple example: the one-dmensonal A 2 molecule. We look now to the Lagrangan theory of what mght be called a one-dmensonal A 2 molecule. The consttuent atoms resde at x and x 2 > x, and both have mass m. The molecule tself therefore has mass M = 2m. We assume the molecule to be bound together by a sprng of natural length a and strength k, and to move n an ambent potental U(x). In natural varables the Lagrangan reads L = 2 m(ẋ2 + ẋ 2 2) U(x ) U(x 2 ) 2 k[(x 2 x ) a] 2 (0) The varables ntutvely most natural to ths smple system are the external coordnate X = M (mx + mx 2 ) = 2 (x + x 2 ) and the nternal coordnate 2s = x 2 x, whch descrbes the nstantaneous length of the molecule. Immedately } x = X s x 2 = X + s () gvng L = m(ẋ2 + ṡ 2 ) U(X s) U(X + s) 2 k(2s a)2 (3) It becomes analytcally advantageous at ths pont to ntroduce the varable q = 2s a whch descrbes molecular length relatve to the rest length of the molecule;
4 4 Lagrangan theory of several body systems then s = 2 (a + q) and we have where K = 2k and We are led thus to wrte L = mẋ2 + 2 { 2 m q2 2 Kq2 } U(X, q) U(X, q) = U(X + s) + U(X s) = {U(X) + U (X)s + 2 U (X)s } + {U(X) U (X)s + 2 U (X)s } = 2U(X) + U (X)s L = {mẋ2 2U(X)} + 2 { 2 m q2 2 Kq2 } U nteracton (X, q) (4) wth U nteracton (X, q) = 4 U (X)(a 2 + 2aq + q 2 ) +... = 2 au (X) q +... The operatve assumpton here s that the ambent potental changes lttle over the dmenson of the molecule. The strkng absence of a U (X)-term s an artfact of our assumpton that m = m 2 = m, and means that the nteracton s, n leadng approxmaton, tdal. In ths respect the physcs of A 2 molecules s dstnct from the physcs of AB molecules. Lookng now to the equatons of moton mẍ = U (X) 2 qau (X) (5.) m q + Kq = au (X) (5.2) we fnd t natural on physcal grounds to abandon the 2 nd term on the rght sde of (5.). Returnng wth X(t) a soluton of the equaton thus obtaned to (5.2), we have m q + Kq = F(t) wth F(t) = au (X(t)) (6) Evdently the tdal term on the rght sde of (5.2) serves n effect to force the nternal oscllaton of the molecule. The reader who was awatng the entry of relatve varables r and r 2 nto the preceedng dscusson wll have been struck by ther absence. Ther non-appearance can be attrbuted to the crcumstance that the asymmetry problem whch motvated our ntal dscusson does not arse n the case N = 2; t s non-dscrmnatory to speak of s = 2 (r 2 r ). In ths respect the physcs of A n molecules (n > 2) s margnally more nterestng. The physcs even of A 2 molecules becomes markedly more nterestng when we gve up the one-dmensonalty or our problem, for then the molecule can be expected to experence torques, and to tumble n ways responsve to the dervatve
5 A less smple example: the one-dmensonal A 3 molecule 5 structure of U(X); I shall, however, resst the temptaton to enter nto an mmedate dscusson of the detals. A less smple example: the one-dmensonal A 3 molecule. We now assume partcles of dentcal mass m to resde at x, x 2 > x and x 3 > x 2 and to be bound by sprngs dentcal to those encountered n the prevous example. In natural varables the Lagrangan (compare (0)) reads L = 2 m(ẋ2 + ẋ ẋ 2 3) U(x ) U(x 2 ) U(x 3 ) 2 k{[(x 2 x ) a] 2 + [(x 3 x 2 ) a] 2 } (7) Drect appropraton of (6) gves rse for such a system to L = 2 MẊ2 + 2 m(ṙ2 + ṙ ṙ 2 3) U(X + r ) U(X + r 2 ) U(X + r 3 ) 2 k[(r 2 r ) a] 2 2 k[(r 3 r 2 ) a] 2 g m(r + r 2 + r 3 ) (8) where M = 3m s the mass of the A 3 molecule, and where X = 3 (x + x 2 + x 3 ) serves to locate ts center of mass. Equatons (9) acqure therefore ths partcularzed meanng: MẌ = U (X + r ) U (X + r 2 ) U (X + r 3 ) (9.) m( r + Ẍ) = U (X + r ) + k[(r 2 r ) a] m( r 2 + Ẍ) = U (X + r 2 ) k[(r 2 r ) a] + k[(r 3 r 2 ) a] (9.2) m( r 3 + Ẍ) = U (X + r 3 ) k[(r 3 r 2 ) a] Lookng frst to (9.), we observe that RHS of (9.) = 3U (X) m U (X) 3 mr 2m U (X) 3 mr (20) The second term on the rght vanshes by defnton of the center of mass. It s nterestng n ths lght to notce that the sum encountered n the thrd term on the rght serves to defne the nstantaneous moment of nerta relatve to the center of mass (.e., the centered second moment of the molecular mass dstrbuton), and that the sums encountered n hgher-order terms defne nameless hgher moments of the mass dstrbuton. Upon the abandonment of all such neglgble terms, (9.) reduces to mẍ = U (X): the center of mass moves as a sngle atom would move n the ambent potental. Lookng now to (9.2), we are motvated by the structure of the sprng terms to ntroduce } q = (r 2 r ) a (2) q 2 = (r 3 r 2 ) a Inverson of q + a = r + r 2 q 2 + a = r 2 + r 3 0 = +r + r 2 + r 3
6 6 Lagrangan theory of several body systems gves r = 3 (2q + q 2 ) a r 2 = + 3 ( q q 2 ) r 3 = + 3 ( q + 2q 2 ) + a (22) and when we return wth (2) and (22) to (8) we by computaton obtan (compare (4)) L = 3{ 2 mẋ2 U(X)} + { 3 m( q2 + q q 2 + q 2 2) 2 k(q2 + q 2 2)} U nteracton (X, q, q 2 ) (23) where U nteracton (X, q, q 2 ) = 2m U (X) m [ 2 9 a a(q 2 q ) = 9 au (X) (q 2 q ) (q2 + q q 2 + q 2 2) ] +... To better emphasze the essentals of the stuaton as t now stands, I wrte L = L center of mass (Ẋ, X) + L nternal( q, q) + L nteracton (X, q) and assert (on the physcal grounds to whch I have already alluded, and to a more careful dscusson of whch I promse to return) that L nteracton (X, q) contrbutes essentally to the moton of the nternal varables q, but nessentally to the moton of X. We are led thus to the equatons of moton mẍ = U (X) (24.) { d } L = 0 ( =, 2) (24.2) dt q q where L = L nternal ( q, q) + L nteracton (X(t), q). I have refraned from wrtng out the detaled mplcatons of (24.2) because those equatons are not welladapted to analytcal treatment. To obtan more workable equatons we very The tedum of the computaton s much reduced f one wrtes r r 2 = G q q 2 and ṙ ṙ 2 = G q 2 q 2 wth G = r 3 a ṙ and notces that G T G = for then (r 2 + r r 2 3) and (ṙ 2 + ṙ ṙ 2 3) become qute easy to evaluate.
7 A less smple example: the one-dmensonal A 3 molecule 7 much n the sprt of standard small oscllaton theory subject L nternal ( q, q) to some preparatory massagng. Wrtng T L nternal ( q, q) = 6 m q M q 2 wth M = ( q q 2 ) T 2 k q K q ( q K = q 2 ) 0 0 we ntroduce new varables (wrte q = RQ wth R a rotaton matrx: R T R = I) to obtan T T Q L nternal ( Q, Q) = 6 m Q R Q T M R 2 Q 2 k Q R T K R 2 Q 2 ( Q Q 2 ) (25) and look to the smultaneous dagonalzaton of R T M R and R T K R. Ths, of course, s standard theory of small oscllatons methodology; the only unusual crcumstance s that here t s not be sprng matrx K but the mass matrx M whch comes to us n ntally non-dagonal form. Some prelmnary observatons: from det(m λi) = λ 2 4λ + 3 we conclude that the egenvalues of M can be descrbed λ = 2 ±, and that our assgnment, therefore, s to dscover the matrx cos ϕ sn ϕ R = sn ϕ cos ϕ such that R T M R = Our assgnment (to say the same thng another way) s to dscover the ϕ such that ( 2 cos = 2 ϕ + 3 sn 2 ) ϕ 2 sn ϕ cos ϕ 2 2 sn ϕ cos ϕ sn 2 ϕ + 3 cos 2 ϕ Immedately ϕ = 45, gvng R = 2 (26) In any event qute apart from the detals of the argument whch led us to (26) we have only to nsert (26) nto (25) to obtan L nternal ( Q, Q) = 6 m( Q Q 2 2) 2 k(q2 + Q 2 ) The modal moton of the free A 3 molecule can therefore be descrbed Q + ω 2 Q = 0 (27)
8 8 Lagrangan theory of several body systems wth ω = 3k/m ω 2 = k/m (28) For such a molecule (.e., for a one-dmensonal A 3 molecule n the total absence of an ambent potental) we have X(t) = X 0 + V 0 t q (t) = A cos(ω t + δ ) + A 2 cos(ω 2 t + δ 2 ) q 2 (t) = A cos(ω t + δ ) + A 2 cos(ω 2 t + δ 2 ) so (recall (22)) n natural varables we have the followng explct descrpton of the free moton latent n the Lagrangan (7): x (t) = X(t) 3 A cos(ω t + δ ) A 2 cos(ω 2 t + δ 2 ) x 2 (t) = X(t) A cos(ω t + δ ) x 3 (t) = X(t) 3 A cos(ω t + δ ) + A 2 cos(ω 2 t + δ 2 ) (29) These equatons make clear the sense n whch the fast mode (the mode under the control of A ) s a hp-swngng mode, and the slow mode (controlled by A 2 ) s a breather mode a dance n whch m 2 does not partcpate. So much for the free moton of the system. Now renstate the ambent potental U(x). The moton X(t) of the center of mass s now no longer unform, but accelerated as descrbed (n leadng approxmaton) by (24.). Of more partcular nterest s the fact that L nteracton = 9 au (X(t)) (q q 2 ) = F(t) Q where F(t) = 2 9 au (X(t)). The mplcaton s that tdal forces couple (n leadng approxmaton) only to the fast mode. Ths I fnd somewhat counterntutve, snce t s the slow breather mode of the A 3 molecule whch most resembles the soltary mode of the A 2 molecule a breather mode whch, as we know from prevous work, does respond to tdal forces. From L molecular ( Q, Q, t) = 6 m( Q Q 2 2) 2 k(q2 + Q 2 ) + F(t) Q we obtan Q + ω 2 Q = 3 m F(t) Q 2 + ω 2 2Q 2 = 0 (30) whch gve back (27) n the absence of tdal forces. Moton of a struck A 3 molecule. In place of (7) we now have L = 2 m(ẋ2 + ẋ ẋ 2 3) 2 k{[(x 2 x ) a] 2 + [(x 3 x 2 ) a] 2 } x F (t) the nessental assumpton here beng that (snce objects are most commonly struck on ther exposed surfaces) t s the end-partcle m whch has been struck.
Physics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationPhysics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationMolecular structure: Diatomic molecules in the rigid rotor and harmonic oscillator approximations Notes on Quantum Mechanics
Molecular structure: Datomc molecules n the rgd rotor and harmonc oscllator approxmatons Notes on Quantum Mechancs http://quantum.bu.edu/notes/quantummechancs/molecularstructuredatomc.pdf Last updated
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationThe classical spin-rotation coupling
LOUAI H. ELZEIN 2018 All Rghts Reserved The classcal spn-rotaton couplng Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 louaelzen@gmal.com Abstract Ths paper s prepared to show that a rgd
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationRigid body simulation
Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More information10. Canonical Transformations Michael Fowler
10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationLAGRANGIAN MECHANICS
LAGRANGIAN MECHANICS Generalzed Coordnates State of system of N partcles (Newtonan vew): PE, KE, Momentum, L calculated from m, r, ṙ Subscrpt covers: 1) partcles N 2) dmensons 2, 3, etc. PE U r = U x 1,
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
More informationSpring 2002 Lecture #13
44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationChapter 3. r r. Position, Velocity, and Acceleration Revisited
Chapter 3 Poston, Velocty, and Acceleraton Revsted The poston vector of a partcle s a vector drawn from the orgn to the locaton of the partcle. In two dmensons: r = x ˆ+ yj ˆ (1) The dsplacement vector
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationMEASUREMENT OF MOMENT OF INERTIA
1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationThe non-negativity of probabilities and the collapse of state
The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationOne Dimension Again. Chapter Fourteen
hapter Fourteen One Dmenson Agan 4 Scalar Lne Integrals Now we agan consder the dea of the ntegral n one dmenson When we were ntroduced to the ntegral back n elementary school, we consdered only functons
More informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationQuantum Mechanics I Problem set No.1
Quantum Mechancs I Problem set No.1 Septembe0, 2017 1 The Least Acton Prncple The acton reads S = d t L(q, q) (1) accordng to the least (extremal) acton prncple, the varaton of acton s zero 0 = δs = t
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationU.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016
U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More information10-701/ Machine Learning, Fall 2005 Homework 3
10-701/15-781 Machne Learnng, Fall 2005 Homework 3 Out: 10/20/05 Due: begnnng of the class 11/01/05 Instructons Contact questons-10701@autonlaborg for queston Problem 1 Regresson and Cross-valdaton [40
More informationHow Differential Equations Arise. Newton s Second Law of Motion
page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationIntegrals and Invariants of
Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
More informationChapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D
Chapter Twelve Integraton 12.1 Introducton We now turn our attenton to the dea of an ntegral n dmensons hgher than one. Consder a real-valued functon f : R, where the doman s a nce closed subset of Eucldean
More informationThe exponential map of GL(N)
The exponental map of GLN arxv:hep-th/9604049v 9 Apr 996 Alexander Laufer Department of physcs Unversty of Konstanz P.O. 5560 M 678 78434 KONSTANZ Aprl 9, 996 Abstract A fnte expanson of the exponental
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationr i r j 3. (2) Gm j m i r i (r i r j ) r i r j 3. (3)
N-body problem There s a lot of nterestng physcs related to the nteractons of a few bodes, but there are also many systems that can be well-approxmated by a large number of bodes that nteract exclusvely
More informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and
More informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
More informationNEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system).
EWTO S LAWS Consder two partcles. 1 1. If 1 0 then 0 wth p 1 m1v. 1 1 2. 1.. 3. 11 These laws only apply when vewed from an nertal coordnate system (unaccelerated system). consder a collecton of partcles
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationApplied Nuclear Physics (Fall 2004) Lecture 23 (12/3/04) Nuclear Reactions: Energetics and Compound Nucleus
.101 Appled Nuclear Physcs (Fall 004) Lecture 3 (1/3/04) Nuclear Reactons: Energetcs and Compound Nucleus References: W. E. Meyerhof, Elements of Nuclear Physcs (McGraw-Hll, New York, 1967), Chap 5. Among
More information2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification
E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton
More informationFirst Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.
Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act
More informationClassical Mechanics Virtual Work & d Alembert s Principle
Classcal Mechancs Vrtual Work & d Alembert s Prncple Dpan Kumar Ghosh UM-DAE Centre for Excellence n Basc Scences Kalna, Mumba 400098 August 15, 2016 1 Constrants Moton of a system of partcles s often
More informationMEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut
MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationClassical Field Theory
Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to
More information4. Laws of Dynamics: Hamilton s Principle and Noether's Theorem
4. Laws of Dynamcs: Hamlton s Prncple and Noether's Theorem Mchael Fowler Introducton: Galleo and Newton In the dscusson of calculus of varatons, we antcpated some basc dynamcs, usng the potental energy
More information6. Hamilton s Equations
6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationThree views of mechanics
Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental
More informationPY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg
PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationAPPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14
APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce
More information