University:
Massachusetts Institute of Technology
Course:
2.141 | Modeling and Simulation of Dynamic Systems
Academic year:

2019
Views:

396

Pages:

4
Author:

tubofib

CANONICAL TRANSFORMATION THEORY A canonical transformation may express new displacements and momenta as functions of both the original displacements and momenta, but is restricted such that it preserves the Hamiltonian form of the differential equations. Original variables ("old coordinates") H = H(q,p) dq/dt = ∂H(q,p)/∂p – f(q,p,t) dp/dt = –∂H(q,p)/∂q + e(q,p,t) The functions f(q,p,t) and e(q,p,t) contain the non-conservative and forcing terms and are known as "canonical forces" (an unfortunate historical misnomer). ⎡dq/dt⎤ ⎡ 0 ⎢ ⎥ ⎢ ⎢ ⎥ =⎢ ⎣dp/dt⎦ ⎣ –1 1 ⎤⎡∂H/∂q⎤ ⎡–f⎤ ⎥⎢ ⎥ ⎢ ⎥ ⎥⎢ ⎥ +⎢ ⎥ 0 ⎦⎣∂H/∂p⎦ ⎣ e ⎦ Symplectic notation: ⎡q⎤ r = ⎢⎣p⎥⎦ H = H(r) ⎡–f⎤ G = ⎢ e⎥ ⎣ ⎦ ⎡ 0 1 ⎤ J = ⎢ –1 0 ⎥ ⎣ ⎦ The 2n x 2n matrix J is the so-called symplectic matrix. Because J J = –1 it is analogous to the complex number j (the square root of –1). It is also analogous to a 90° rotation matrix. dr/dt = J ∂H(r)/∂r + G(r,t) Change variables ("new coordinates") q = q(q*,p*) p = p(p*,q*) Transform the Hamiltonian by substitution H(q,p) = H(q(q*,p*),p(p*,q*)) = H(q*,p*) In symplectic notation r = r(r*) H(r*) = H(r*(r)) = H(r) Differentiate the transformation from r to r* dr/dt = T dr*/dt where T = ∂r/∂r* is the Jacobian of the "reverse" transformation. Note that T is, in general, a function of r*. By design, the differential equations in the new variables are to have a Hamiltonian form. dr*/dt = J ∂H/∂r* + G* Substitute for dr*/dt dr/dt = T J ∂H/∂r* + T G* The gradient ∂H/∂r* may be expressed as follows ∂H/∂r* = Tt ∂H/∂r To justify pre-multiplication by the transpose, expand the partial derivative and note that Tij = ∂ri/∂r*j ∂H ∂r*i = ∑ j ∂H ∂rj ∂rj ∂r*i = ∑ ∂H ∂rj Tji j Therefore dr/dt = T J Tt ∂H/∂r + T G* The Hamiltonian form is to be preserved even if the system is unforced and conservative (i.e. G = 0, G* = 0) so the symplectic condition for a canonical transformation is T J Tt = J In the general case the relation between the "canonical forces" is G = T G* We need a transformation from G* to G. We could write G* = T-1 G but it is possible to avoid inverting the transformation by rearranging the symplectic condition as follows: Post-multiply by T-t -t TJ=JT Pre-multiply by J, post-multiply by –J and use J J = –1 -t JT=T J Pre-multiply by Tt Tt J T = J This is an alternative form of the symplectic condition. Apply it to the transformation of canonical forces by pre-multiplying by Tt J and using J-1 = –J Tt J G = Tt J T G* = J G* G* = – J Tt J G The symplectic notation is compact and elegant, but tends to disguise some of the structure of the equations. It can be instructive to express the results in the usual notation. The "reverse" Jacobian T is d ⎡⎢q⎤⎥ ⎡⎢ ∂q/∂q* ∂q/∂p* ⎤⎥ d ⎡⎢q*⎤⎥ dt ⎢⎣p⎥⎦ = ⎢⎣ ∂p/∂q* ∂p/∂p* ⎥⎦ dt ⎢⎣p*⎥⎦ To express the transformation of canonical forces, note that ⎡ 0 J G* = ⎢⎢ ⎣ –1 ⎡ 0 J G = ⎢⎢ ⎣ –1 1 ⎤ ⎡–f*⎤ ⎡e*⎤ ⎥⎢ ⎥ ⎢ ⎥ ⎥⎢ ⎥ = ⎢ ⎥ 0 ⎦ ⎣ e* ⎦ ⎣ f* ⎦ 1 ⎤ ⎡–f⎤ ⎡e⎤ ⎢ ⎥ ⎥⎢ ⎥ ⎥⎢ ⎥ = ⎢ ⎥ 0 ⎦ ⎣ e⎦ ⎣f⎦ Thus the transformation is ⎡ [∂q/∂q*]t [∂p/∂q*]t ⎤ ⎡e⎤ ⎡e*⎤ ⎡e⎤ ⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥ t ⎢ ⎥ =T ⎢ ⎥ =⎢ ⎥⎢ ⎥ ⎣f⎦ ⎣ f* ⎦ ⎣ [∂q/∂p*]t [∂p/∂p*]t ⎦ ⎣ f ⎦ Expand e* = [∂q/∂q*]t e + [∂p/∂q*]t f f* = [∂q/∂p*]t e + [∂p/∂p*]t f Writing the summations explicitly e*i = ∑(ej ∂qj/∂q*i + fj ∂pj/∂q*i) j f*i = ∑(ej ∂qj/∂p*i + fj ∂pj/∂p*i) j

Canonical Transformation Theory

Related Documents

Recommended Documents

Description

Related Documents

Free up your schedule!

Take 5 seconds to unlock