Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Size: px
Start display at page:

Download "Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1"

Transcription

1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems that are boune. It allows us to etermine statistical properties of systems. The other is the Poisson bracket, which has allows us to etermine properties of systems without solving the equations of motion, an is connecte to the commutator in quantum mechanics. 1.1 The Virial Theorem Let s consier the variable Q p i q i, which remains boune if both p i an q i are finite. From Hamilton s equations of motion we have for the total erivative of Q Q p i q i q i ṗ i p i q i (1) so that Q(t) Q(0) 1 T Qt 1 T ( ) t p i q i (2) T T 0 T 0 If the system is perioic an we integrate over one perio, then left han sie is zero. If the system is boune, then as T the left han sie goes to zero. In these case, we have p i q i (3) where the brackets enote the time average values of the enclose quantities. The quantity V q i (4) is a generalization of the virial function introuce into the kinetic theory of gases by Clausius 1. Thus the time average virial function is V p i (5) For a system of N particles, the right-han sie becomes p i q i 2 T i 2N T (6) where the bar enotes an average over time an kinetic energy of each particle. The virial function for this system can be also expresse as follows V ṗ i q i I P r ˆnS (7) 1 This is just the total work one on the gas, incluing external an internal forces. The Virial Theorem an The Poisson Bracket-1

2 where the left ha sie is the work one by the intermolecular forces, an the secon term is the work one in containing the gas with in the container. Using the ivergence theorem, we arrive at V I V P r I 3P V (8) where V is the total volume. Combining the two expression for the average virial function we arrive at P V 2 3 N T I 3 Nkθ I 3 where θ is the absolute temperature from the expression T 3 2kθ. If I 0, this is the ieal gas law. Another application of the virial theorem is the case of a particle in a central potential V (r). In this case the average virial function is given by V p r p r (10) which is reuce to For a potential of the form V (r) Kr n, this expression becomes 2 T r V (r) (11) T 1 n V (r) (12) 2 For the gravitational or Coulomb forces n 1. This means that the average kinetic energy is 1/2 the average potential energy. 1.2 The Poisson Bracket The Poisson bracket appears often in classical mechanics an translates itself to quantum mechanics as the commutator. Let s consier the total time erivative of a function u that epens on the canonical variables q an p, an can also epen on the time t u t ( u q i u ) ṗ i u Using the canonical equations of motion, the total time erivative becomes u t ( u u ) u which can be written using the following shorthan notation where we efine the Poisson bracket as [u, H] qi,p i u t [u, H] q i,p i u ( u u ) (9) (13) (14) (15) (16) The Virial Theorem an The Poisson Bracket-2

3 Let s consier the case of u not epening explicitly on time, an the Poisson bracket being zero. In this case the quantity u(q, p) is a constant of the motion. This is the same as we have in quantum mechanics. If the commutator of an operator with the Hamiltonian is zero an oes not epen explicitly on time, then the operator correspons to a variable that is a constant of the motion. As an example of how this works, let s consier a simple harmonic oscillator in two imensions. Using the Poisson bracket, we will show that the angular momentum is conserve (a constant of the motion). The Hamiltonian for this system is H 1 2 (p2 1 p 2 2) 1 2 (q2 1 q 2 2) (17) where for simplicity we take m k 1. The canonical equations of motion are ṗ i q i The angular momentum associate with this system is The Poisson bracket is q i p i with i 1, 2 (18) J q 1 p 2 q 2 p 1 (19) J t J J (p 2 p 1 p 1 p 2 ) (q 2 q 1 q 1 q 2 ) 0 (20) therefore the angular momentum is a constant of the motion. Base on the efinition of the Poisson bracket, the canonical equations of motion can be expresse in a compact form. Let s start by calculating q i using the Poisson bracket q i [q i, H] q i (21) Similarly the canonical equation of motion for ṗ i can be erive ṗ i [p i, H] ṗ i (22) At this point we efine the Poisson bracket of any two functions that epen on the canonical variables ( u v [u, v] v ) u (23) If the functions we use are the canonical variables themselves, we get the following funamental brackets [q i, q j ] 0 [q i, p j ] δ ij [p i, p j ] 0 (24) Notice how closely these resemble the commutators in quantum mechanics. Just replace δ ij with i hδ ij. Before going further, let s consier a few properties of the Poisson bracket. Assume that we have the ynamical variables u(q i, p i ), v(q i, p i ), w(q i, p i ). The Poisson bracket has the following properties The Virial Theorem an The Poisson Bracket-3

4 1. Linearity: 2. Antisymmetric: 3. The prouct rule: 4. The Jacobi ientity: [αu βw, v] α [u, v] β [w, v] (25) [u, w] [w, u] (26) [u, wv] [u, w] v w [u, v] (27) [w, [u, v]] [u, [v, w]] [v, [w, u]] 0 (28) 5. The chain rule: if the variables are expresse in terms of a new set P i an Q i that are functions of the ol variables p i an qi, then the Poisson bracket of the ynamical variables is [w, u] w η i [η i, η j ] u η j (29) where η i Q i an η in P i for i 1 n an we have a total of 2n variables. 6. The funamental brackets: 1.3 Poisson s Theorem [q i, q j ] [p i, p j ] 0 [q i, p j ] δ ij (30) Assume that we have a set of ynamical variables an we woul like to know whether their time evolution is generate by a Hamiltonian or not. Poisson s theorem states that for the time evolution of a set of variables to be generate by a Hamiltonian, the variables must satisfy ] [ ] [Ṙ, t [R, S] S R, Ṡ Assume that the time evolution is generate by H(p i, q i, t), then accoring to Eq. 15 we have The first term can be reuce using the Jacobi ientity (31) t [R, S] [[R, S], H] [R, S] (32) [[R, S], H] [[R, H], S] [R, [S, H]] (33) The secon term can be expane as follows [R, S] ( R S R ) S ( 2 R S 2 R S ) ( R 2 S R 2 R ) [ R, S ] [ R, S ] (34) The Virial Theorem an The Poisson Bracket-4

5 Next we combine the two sets [R, S] t [ [R, H] R, S ] [ R, [S, H] S ] ] [ ] [Ṙ, S R, Ṡ which is the statement of the theorem Notice that this theorem also implies that if R an S are constants of the motion, then so is [R, S]. Let s consier two examples: The first we take the one-imensional simple harmonic oscillator. The two ynamical variables are q an p. These have the following equations of motion where we take m k 1. Applying Poisson s theorem we get (35) ṗ q q p (36) [q, p] 0 t (37) [ṗ, q] [p, q] [q, q] [p, p] 0 (38) which implies that the time evolution of the ynamical variables is generate by a Hamiltonian as expecte. As a secon example, consier the same two variables but with the following equations of motion Again we apply Poisson s theorem q pq ṗ pq (39) [q, p] 0 t (40) [ q, p] [q, ṗ] [pq, p] [q, pq] [q, p] p q [p, p] [q, p] q [q, q] p p q (41) which oes not satisfy Poisson s theorem, therefore the time evolution of these variables is not generate by a Hamiltonian. To show that a transformation is canonical, we take a set of ynamical variables R an S that are functions of q i an p i an transform the variables to a new set Q i an P i. Then if the Poisson bracket for both sets of phase coorinates are equal the transformation is canonical We start by taking the time erivative of both sies then expan this out [R, S] qi,p i [R, S] Qi,P i (42) t [R, S] q i,p i t [R, S] Q i,p i (43) ] t [R, S] Qi,P i [Ṙ, S Imposing the Poisson bracket conition, we get ] t [R, S] Qi,P i [Ṙ, S q i,p i Q i,p i [ ] R, Ṡ (44) q i,p i [ ] R, Ṡ (45) Q i,p i which accoring to Poisson s theorem, the time evolution of the variables are generate by a Hamiltonian. The Virial Theorem an The Poisson Bracket-5

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10 Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Markov Chains in Continuous Time

Markov Chains in Continuous Time Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in

More information

221A Lecture Notes Notes on Classica Mechanics I

221A Lecture Notes Notes on Classica Mechanics I 1A Lecture Notes Notes on Classica Mechanics I 1 Precursor: Fermat s Principle in Geometric Optics In geometric optics, you talk about how light rays go. In homogeneous meiums, the light rays go straight.

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems D. D. Holm M3-4-5A16 Assesse Problems # 1 Due 1 Nov 2012 1 1 M3-4-5A16 Assesse Problems # 1: Do 4 out of 5 problems Exercise 1.1 (Poisson brackets for the Hopf map) Figure 1: The Hopf map. In coorinates

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N Physics 15 Solution Set Winter 018 1. At time t = 0, the wave function of a free particle moving in a one-imension is given by, ψ(x,0) = N where N an k 0 are real positive constants. + e k /k 0 e ikx k,

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

Darboux s theorem and symplectic geometry

Darboux s theorem and symplectic geometry Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus

More information

15 Functional Derivatives

15 Functional Derivatives 15 Functional Derivatives 15.1 Functionals A functional G[f] is a map from a space of functions to a set of numbers. For instance, the action functional S[q] for a particle in one imension maps the coorinate

More information

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3 13.021 Marine Hyroynamics, Fall 2004 Lecture 3 Copyright c 2004 MIT - Department of Ocean Engineering, All rights reserve. 1.2 - Stress Tensor 13.021 Marine Hyroynamics Lecture 3 Stress Tensor τ ij:. The

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Math 115 Section 018 Course Note

Math 115 Section 018 Course Note Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

Interpolated Rigid-Body Motions and Robotics

Interpolated Rigid-Body Motions and Robotics Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Physics 5153 Classical Mechanics. Canonical Transformations-1

Physics 5153 Classical Mechanics. Canonical Transformations-1 1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Kinematics of Rotations: A Summary

Kinematics of Rotations: A Summary A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are

More information

Kramers Relation. Douglas H. Laurence. Department of Physical Sciences, Broward College, Davie, FL 33314

Kramers Relation. Douglas H. Laurence. Department of Physical Sciences, Broward College, Davie, FL 33314 Kramers Relation Douglas H. Laurence Department of Physical Sciences, Browar College, Davie, FL 333 Introuction Kramers relation, name after the Dutch physicist Hans Kramers, is a relationship between

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x =

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x = 6.4 Integration using tan/ We will revisit the ouble angle ientities: sin = sin/ cos/ = tan/ sec / = tan/ + tan / cos = cos / sin / tan = = tan / sec / tan/ tan /. = tan / + tan / So writing t = tan/ we

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

Control Volume Derivations for Thermodynamics

Control Volume Derivations for Thermodynamics Control olume Derivations for Thermoynamics J. M. Powers University of Notre Dame AME 327 Fall 2003 This ocument will give a summary of the necessary mathematical operations necessary to cast the conservation

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like)

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like) Stockton Univeristy Chemistry Program, School of Natural Sciences an Mathematics 101 Vera King Farris Dr, Galloway, NJ CHEM 340: Physical Chemistry II Solving the Schröinger Equation for the 1 Electron

More information

Problem Set II Solutions

Problem Set II Solutions Physics 31600 R. Wal Classical Mechanics Autun, 2002 Proble Set II Solutions 1) Let L(q, q; t) be a Lagrangian [where, as in class, q stans for (q 1,..., q n )]. Suppose we introuce new coorinates (Q 1

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Implicit Differentiation Using the Chain Rule In the previous section we focuse on the erivatives of composites an saw that THEOREM 20 (Chain Rule) Suppose that u = g(x) is ifferentiable

More information

Ψ(x) = c 0 φ 0 (x) + c 1 φ 1 (x) (1) Ĥφ n (x) = E n φ n (x) (2) Ψ = c 0 φ 0 + c 1 φ 1 (7) Ĥ φ n = E n φ n (4)

Ψ(x) = c 0 φ 0 (x) + c 1 φ 1 (x) (1) Ĥφ n (x) = E n φ n (x) (2) Ψ = c 0 φ 0 + c 1 φ 1 (7) Ĥ φ n = E n φ n (4) 1 Problem 1 Ψx = c 0 φ 0 x + c 1 φ 1 x 1 Ĥφ n x = E n φ n x E Ψ is the expectation value of energy of the state 1 taken with respect to the hamiltonian of the system. Thinking in Dirac notation 1 an become

More information

CHM 532 Notes on Creation and Annihilation Operators

CHM 532 Notes on Creation and Annihilation Operators CHM 53 Notes on Creation an Annihilation Operators These notes provie the etails concerning the solution to the quantum harmonic oscillator problem using the algebraic metho iscusse in class. The operators

More information

Approximate Molecular Orbital Calculations for H 2. George M. Shalhoub

Approximate Molecular Orbital Calculations for H 2. George M. Shalhoub Approximate Molecular Orbital Calculations for H + LA SALLE UNIVESITY 9 West Olney Ave. Philaelphia, PA 94 shalhoub@lasalle.eu Copyright. All rights reserve. You are welcome to use this ocument in your

More information

Calculus of variations - Lecture 11

Calculus of variations - Lecture 11 Calculus of variations - Lecture 11 1 Introuction It is easiest to formulate the problem with a specific example. The classical problem of the brachistochrone (1696 Johann Bernoulli) is the search to fin

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lecture 3 Continuous Systems an Fiels (Chapter 13) Where Are We Now? We ve finishe all the essentials Final will cover Lectures 1 through Last two lectures: Classical Fiel Theory

More information

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) = 2 2t 1 + t 2 cos x = 1 t2 We will revisit the double angle identities:

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) = 2 2t 1 + t 2 cos x = 1 t2 We will revisit the double angle identities: 6.4 Integration using tanx/) We will revisit the ouble angle ientities: sin x = sinx/) cosx/) = tanx/) sec x/) = tanx/) + tan x/) cos x = cos x/) sin x/) tan x = = tan x/) sec x/) tanx/) tan x/). = tan

More information

arxiv:nlin/ v1 [nlin.cd] 21 Mar 2002

arxiv:nlin/ v1 [nlin.cd] 21 Mar 2002 Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Discrete Hamilton Jacobi Theory and Discrete Optimal Control

Discrete Hamilton Jacobi Theory and Discrete Optimal Control 49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin

More information

Centrum voor Wiskunde en Informatica

Centrum voor Wiskunde en Informatica Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.

More information

Can we derive Newton s F = ma from the SE?

Can we derive Newton s F = ma from the SE? 8.04 Quantum Physics Lecture XIII p = pˆ (13-1) ( ( ) ) = xψ Ψ (13-) ( ) = xψ Ψ (13-3) [ ] = x (ΨΨ ) Ψ Ψ (13-4) ( ) = xψ Ψ (13-5) = p, (13-6) where again we have use integration by parts an the fact that

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

The Natural Logarithm

The Natural Logarithm The Natural Logarithm -28-208 In earlier courses, you may have seen logarithms efine in terms of raising bases to powers. For eample, log 2 8 = 3 because 2 3 = 8. In those terms, the natural logarithm

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

Multi-component bi-hamiltonian Dirac integrable equations

Multi-component bi-hamiltonian Dirac integrable equations Chaos, Solitons an Fractals 9 (009) 8 8 www.elsevier.com/locate/chaos Multi-component bi-hamiltonian Dirac integrable equations Wen-Xiu Ma * Department of Mathematics an Statistics, University of South

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

SYMPLECTIC GEOMETRY: LECTURE 3

SYMPLECTIC GEOMETRY: LECTURE 3 SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as

More information

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives

More information

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

Lecture 2 - First order linear PDEs and PDEs from physics

Lecture 2 - First order linear PDEs and PDEs from physics 18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate

More information

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Geometric Algebra Approach to Fluid Dynamics

Geometric Algebra Approach to Fluid Dynamics Geometric Algebra Approach to Flui Dynamics Carsten Cibura an Dietmar Hilenbran Abstract In this work we will use geometric algebra to prove a number of well known theorems central to the fiel of flui

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Basic Thermoelasticity

Basic Thermoelasticity Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................

More information

QUANTUMMECHANICAL BEHAVIOUR IN A DETERMINISTIC MODEL. G. t Hooft

QUANTUMMECHANICAL BEHAVIOUR IN A DETERMINISTIC MODEL. G. t Hooft QUANTUMMECHANICAL BEHAVIOUR IN A DETERMINISTIC MODEL G. t Hooft Institute for Theoretical Physics University of Utrecht, P.O.Box 80 006 3508 TA Utrecht, the Netherlans e-mail: g.thooft@fys.ruu.nl THU-96/39

More information

Chapter 3 Definitions and Theorems

Chapter 3 Definitions and Theorems Chapter 3 Definitions an Theorems (from 3.1) Definition of Tangent Line with slope of m If f is efine on an open interval containing c an the limit Δy lim Δx 0 Δx = lim f (c + Δx) f (c) = m Δx 0 Δx exists,

More information

The Three-dimensional Schödinger Equation

The Three-dimensional Schödinger Equation The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation

More information

MATHEMATICS BONUS FILES for faculty and students

MATHEMATICS BONUS FILES for faculty and students MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department

More information

Witten s Proof of Morse Inequalities

Witten s Proof of Morse Inequalities Witten s Proof of Morse Inequalities by Igor Prokhorenkov Let M be a smooth, compact, oriente manifol with imension n. A Morse function is a smooth function f : M R such that all of its critical points

More information