Der Hamiltonsche kanonische Formalismus

Aus PhysikWiki
Version vom 17. August 2010, 18:02 Uhr von Schubotz (Diskussion | Beiträge) (Die Seite wurde neu angelegt: „==Der Hamiltonsche Kanonische Formalismus== ====Motivation==== Die Lagrange- Theorie benutzt als dynamische Variablen die verallgemeinerten Koordinaten qk und d…“)
(Unterschied) ← Nächstältere Version | Aktuelle Version (Unterschied) | Nächstjüngere Version → (Unterschied)
Zur Navigation springen Zur Suche springen

Der Hamiltonsche Kanonische Formalismus

Motivation

Die Lagrange- Theorie benutzt als dynamische Variablen die verallgemeinerten Koordinaten qk und deren Geschwindigkeiten:


L(q1,...,qf,q˙1,...,q˙f,t)ddtLq˙kLqk=0 k=1,..,f

Wir erhalten f DGL 2. Ordnung für qk(t) im Lagrangeformalismus

Bei gewissen Problemstellungen, wenn es beispielsweise zyklische Variablen gibt:


Lqk=0Lq˙k=const


oder auch bei bestimmten Erweiterungen der Theorie ( Quantenmechanik, statistische Mechanik)

ist es vorteilhaft, statt qk und deren Geschwindigkeiten qk und die zu qk konjugierten Impulse zu benutzen.

Die zu den verallgemeinerten Koordinaten konjugierten Impulse lauten:


pk:=Lq˙k


Die erforderliche Variablentransformation


(qk,q˙k,t)(qk,pk,t)


leistet die sogenannte Legendre- Transformation.

Im Hamiltonformalismus ergeben sich nun 2f DGL 1. Ordnung für

qk(t) und pk(t)

Legendre- Transformation und Hamiltonfunktion

Mathematisches Problem:

Ausgehend von einer Funktion y=f(x) wobei f(x)L,xq˙k


Soll statt x die Variable u=ddxf(x)


verwendet werden.

Also die Steigung von f(x) am Punkt x.

Das bedeutet, wir wollen eine Funktion, die nicht von bestimmten Koordinaten x abhängt sondern nur von der Steigung der Funktion selbst, die sie in Abhängigkeit von x an diesen Stellen hätte.


u=ddxf(x)x=ϕ(u)


Mit der Voraussetzung


d2dx2f(x)0 Ansonsten: u=ddxf(x)=const. und damit nicht umkehrbar.

Die Substitution in f führt dann auf:


y=f(x)=f(ϕ(u)):=F(u)


Bei dieser Trafo geht jedoch Information verloren.

Das heißt: Aus


y=F(u) ist y=f(x) nicht mehr eindeutig rekonstruierbar, weil alle Funktionen


y=f(x)+a mit beliebigem, konstanten a wegen


u=ddxy=ddx(f(x)+a)=dfdx


dann auf die selbe Funktion y=F(u) führen:

Alle Funktionen, die die gleiche Steigung u bei x haben, führen auf das selbe F(u). y=F(u) ist also nicht umkehrbar.

Auf das selbe F(u) führt jeweils die gesamte lineare Kurvenschar aller Funktionen f(x) zzgl. eines konstanten Parameters.

deshalb wird die Legendre- Transformation, eine andere, umkehrbare, Transformation, eingeführt:


xuyz=xuf(x)=ϕ(u)uf(ϕ(u))=g(u)


Statt der gerade genannten einfachen Variante


y=f(ϕ(u))


Die Trafo


y(x)g(u) heißt LEGENDRE- TRANSFORMATION

Graphische Veranschaulichung von g(u)=xdfdxf(x)

Das bedeutet: Die ursprüngliche Funktion y=f(x) wird nach der Trafo x-> u durch die Steigung u von f(x) und den (negativen) Achsenabschnitt charakterisiert.

Da der Achsenabschnitt mit berücksichtigt ist, wird nicht mehr die gesamte Kurvenschar auf die gleiche g(u) abgebildet. Die Abbildung ist bijektiv und damit eindeutig.

Die Werte (u,-z) bestimmen die Schar der Einhüllenden von y=f(x):

Man spricht deshalb von einer Legendre - oder auch Berührungstransformation.

Anwendung auf die Lagrangefunktion

yLxq˙up=Lq˙zq˙pL:=H(q,p,t)


Die Legendretransformierte H(q,p,t) heißt Hamiltonfunktion.

Die Variablen q und t werden nicht geändert.

Wichtig ist jedoch bei mehreren Variablen q1,...,qf:


L(q1,...,qf,q˙1,...,q˙f,t)pk:=Lq˙kH(q1,...,qf,p1,...,pf,t)=k=1fq˙kpkL


Die mathematischen Voraussetzungen für diese Prozedur sind:


L(q1,...,qf,q˙1,...,q˙f,t)C2det2Lq˙kq˙l0 damit pk=Lq˙k nach q˙k auflösbar

Die Hamiltonschen Gleichungen

Ziel: Auch hier natürlich sollen Bewegungsgleichungen für die pk,qk gefunden werden.

Die Ableitung einer Bewegungsgleichung für qk aus der Lagrangegleichung 2. Art ist bereits bekannt:

Eine Variable:

Differenziale:


L=L(q,q˙,t):dL=Lqdq+Lq˙dq˙+LtdtH=H(q,p,t):dH=Hqdq+Hpdp+HtdtH=q˙pLdH=q˙dp+pdq˙dL=q˙dp+pdq˙LqdqLq˙dq˙Ltdt


wegen


Lq˙=pdH=Hqdq+Hpdp+Htdt=q˙dpLqdqLtdt


Dies gilt fuer beliebige Differenziale in q, p und t. Somit kann die Gleichung nur erfüllt werden für


Hq=LqHp=q˙Ht=Lt


Mit Hilfe der Lagrange Bewegungsgleichung


ddtLq˙=Lq;Lq˙=p;Lq=Hqp˙=HqHp=q˙


Die Hamiltonschen Gleichungen sind also beide gefunden.

Es handelt sich um 2 DGLn 1. Ordnung für q und p statt 1 DGL 2. Ordnung für q(t)

Mehrere Variablen

L(q1,...,qf,q˙1,...,q˙f,t)pk:=Lq˙kH(q1,...,qf,p1,...,pf,t)=k=1fq˙kpkL


dH=k(Hqkdqk+Hpkdpk)+Htdt=kq˙kdpk+kpkdq˙kkLq˙kdq˙kkLqkdqkLtdt=kq˙kdpkkLqkdqkLtdtHqk=LqkLqk=ddtLq˙k=ddtpkHpk=q˙k;Ht=Lt


Somit folgen hier die Hamiltonschen Gleichungen (Kanonische Gleichungen)


p˙k=HqkHpk=q˙kk=1,...,f


Der 2f- dimensionale Raum


Γ:={q1,...,qf,q˙1,...,q˙f}


heißt Phasenraum.

Er findet besonders in der klassischen statistischen Mechanik Anwendung. Dabei b4trachtet man Wahrscheinlichkeitsverteilungen auf dem Phasenraum.

Physikalische Bedeutung der Ham- Funktion

  • wegen L= T-V bei holonomen Zwangsbed. und konservativen Kräften
  • und wegen p(d/dt q)= 2T folgt: H = T+V

Dies gilt bei zeitlicher Translationsinvarianz ( skleronome Zwangsbed. ):

mit tr¯i(q1,...,qf)=0undtL=0


Dann nämlich ist


k=1fTq˙kq˙k=2T

 ( nach dem Eulerschen Satz: T ist quadratische , homogene Funktion der 

q˙k .

Somit:


k=1fLq˙kq˙kL=k=1fTq˙kq˙kL=2TT+V=T+V


beschreibt die Gesamtenergie des Systems: Nur bei skleronomen Zwangsbedingungen und konservativen Kräften !

Nach dem Noether- Theorem, speziell unter dem Kapitel ZEITLICHE TRanslationsinvarianz

folgt dann Gesamtenergieerhaltung.

Dies läßt sich leicht nachweisen:


dHdt=ddt(k=1fLq˙kq˙kL)=k=1f(Hqkq˙k+Hpkp˙k)+Ht=k=1f(HqkHpkHpkHqk)Lt=0wegenLt=0 Dies gilt also nur für skleronome Zwangsbedingungen. Bei rheonomen Zwangsbed. ist im Allgemeinen H nicht T+V !!

Beispiel: Perle an starrem rotierendem Draht:

Eine Perle der Masse m sei auf einem starren Draht, der in der -y- Ebene rotiert ( Reibung durch Erdpotenzial zu vernachlässigen): Generalisierte Koordinaten q ist der Abstand der Perle vom Mittelpunkt:


Man kann sich H=T+V denken. Dabei gilt das effektive Potenzial mit V=mq2ω2 .

Aus Lt=0 folgt dann ohnehin wieder ein Erhaltungssatz: H=const.

Typisches Anwendungsschema des Hamilton- Formalismus:

  1. Zunächst sind die generalisierten Koordinaten zu wählen:

q¯=(q1,...,qf)

  1. Transformation des Radiusvektors

r¯i=r¯i(q1,...,qf,t)r¯˙i=r¯˙i(q¯,q¯˙,t)

  1. Aufstellung der Lagrangegleichung:

L(q¯,q¯˙,t)=TV=12mir¯˙i2V

  1. Bestimmung der generalisierten Impulse:

pk:=Lq˙kpk=pk(q¯,q¯˙,t)Umkehrung:q˙k=q˙k(q¯,p¯,t)

  1. Anschließend Legendre Trafo:

H(q1,...,qf,p1,...,pf,t)=k=1fq˙kpkL

  1. Aufstellung und Integration der kanonischen Gleichungen:

p˙k=HqkHpk=q˙kk=1,...,f


Beispiele:

Teilchen in Zylinderkoordinaten ganz ohne Zwnagsbedingungen
  1. q1=3, q2=Phi, q3 = z

x=rcosϕ,x˙=r˙cosϕrϕ˙sinϕy=rsinϕ,y˙=r˙sinϕ+rϕ˙cosϕz=z

T=12m(x˙2+y˙2+z˙2)=12m(r˙2+r2ϕ˙2+z˙2)V=V(r,ϕ,z)L=L(r,ϕ,z,r˙,ϕ˙,z˙)=12m(r˙2+r2ϕ˙2+z˙2)V

  1. Generalisierte Impulse:

pk=Lq˙kpr=mr˙pϕ=mr2ϕ˙pz=mz˙

                                                                                            Radialimpuls, z-Komponente des Drehimpulses und z-Komponente des Impulses
  1. Aufstellung der Legendretrafo:

H=mr˙2+mr2ϕ˙2+mz˙2=12m(r˙2+r2ϕ˙2+z˙2)+VH=12m(pr2+pϕ2r2+pz2)+V(r,ϕ,z)

  1. Kanonische Gleichungen:

p˙k=HqkHpk=q˙kk=1,...,fr˙=Hpr=prm,ϕ˙=Hpϕ=pϕmr2,z˙=Hpz=pzmp˙r=Hr=pϕ2mr3Vr,p˙ϕ=Hϕ=Vϕ,p˙z=Hz=Vz

Interessant ist das Ergebnis der Zentrifugalkraft ( Scheinkraft):

F(Zentrifugal)= pϕ2mr3 , die den radialen Impuls ändert.

Bekannt aus dem Keplerproblem ist uns bereits der Fall V®, ein Zentralpotenzial bei ebener Bewegung:


p˙r=Hr=pϕ2mr3Vr,p˙ϕ=0,p˙z=0


Somit sind Drehimpuls in der Ebene und z-Impuls des Systems erhalten.


z,ϕ sind zyklische Variablen


pz=const.=o.B.d.A.=0pϕ=const. oBdA: ebene Bewegung, Drehimpulserhaltung in der Ebene

Beispiel: eindimensionaler harmonischer Oszi:

Das System ist skleronom wegen Lt=0 , also folgt Energieerhaltung: E=H=T+V


12m(q˙2+ωo2q2)=E=12m(p2m2+ωo2q2)p22mE+q2(2Emωo2)=1


Also ist die Lösung der Phasenraumkurve eine Ellipse. Die Ellipsengröße variiert je nach Energie:

Die Halbachsen sind:


2mE=a,b=2Emωo2 ( bestimmt durch 1. Integral).

Als kanonische Gleichungen ergibt sich:


p˙=Hq=mωo2qq˙=Hp=pm


Daraus folgt dann gerade die Bewegungsgleichung q¨=ddtHq=p˙´m=ωo2qq¨+ωo2q=0


Diese definiert ein Richtungsfeld im Phasenraum

Geladenes Teilchen im elektromagnetischen Feld:

Aus dem Kapitel Eichtransformation der Lagrangefunktion ist das nötige Handwerkszeugs bereits bekannt:


L(q¯,q¯˙,t)=TV=12mir¯˙i2V=m2q¯˙2+e(q¯˙A¯(q¯,t)Φ(q¯,t))


die kanonischen konjugierten Impulse lauten:


pk=L(q¯,q¯˙,t)q˙k=mq˙´k+eAk(q¯,t)q˙k=1m(pkeAk)H=k=13pkq˙kL=k=13pk1m(pkeAk)12mk=13(pkeAk)2k=13em(pkeAk)Ak+eΦH(q¯,p¯,t)=12m(p¯eA¯(q¯,t))2+eΦ(q¯,t)


Dabei begegnen uns die feinen Unterschiede im Impuls, nämlich


mq¯˙=p¯eA¯ als kinetischer Impuls ( der auch tatsächlich mit der Geschwindigkeit verknüpft ist).


pk=Lq˙k ist kanonischer Impuls

Kanonische Transformationen

Wir wissen bereits, dass die Wahl der verallgemeinerten Koordinaten nicht eindeutig ist ( Kapitel 2.4: Forminvarianz der Lagrangegleichungen).

Dabei haben wir gesehen, dass die Lagrangegleichungen 2. Art forminvariant bleiben unter beliebigen diffeomorphen Transformationen der Koordinaten:


q¯=(q1,...,qf)Q¯=(Q1,...,Qf)


Dabei gilt dann:


L¯(Q¯,Q¯˙,t)=L(q¯(Q¯,t),q¯˙(Q¯,Q¯˙,t),t)


Nun kann man sich fragen, unter welchen Transformationen (q¯,p¯)(Q¯,P¯)


die Hamiltonfunktionen forminvariant sind, also:

mit p˙k=Hqkq˙k=Hpk soll auch P˙k=H¯QkQ˙k=H¯Pk

gelten !

Nebenbemerkungen:

  • die Klasse der erlaubten Transformationen muss größer sein als beim Lagrangeformalismus, da jetzt die pk neben den qk als UNABHÄNGIGE Variablen betrachtet werden, die ebenfalls und vor allem völlig unabhängig transformiert werden können.
  • Die neuen Qk und Pk haben unter Umständen gar nicht mehr den Charakter von Orts- und Impulsvariablen.

In den Lagrangegleichungen der 2. Art heißt qj zyklisch, wenn:


Lqj=0ddtLq˙j=ddtpj=0pj=Lq˙j=const


Allerdings ist damit keine Aussage über q˙j gemacht. Diese muss natürlich weiter als Variable behandelt werden.

Hamilton-Gleichungen:

In H(q1,...,qf,p1,...,pf,t) heißt qj zyklisch, wenn


Hqj=0Lqj=ddtpj=0pj:=αj=const


Das bedeutet nun, dass qj in H gar nicht auftritt. pj kann dagegen durch die Bewegungskonstante


αj ersetzt werden:


H(q1,...,qj1,qj+1,...,qf,p1,...,pj1,αj,pj+1,...,pf,t)


Damit jedoch hat das kanonische System nur noch f-1 Freiheitsgrade.

Idee ist es nun, die Hamiltongleichungen zu lösen, indem man Schritt für Schritt zyklische Variablen durch geeignete Trafos der (q¯,p¯)(Q¯,P¯)

einführt, bis alle 

Q¯ zyklisch sind:


H=H(P1,...,Pf,t) mit Pk=αk=const.


Q˙k=HPk=:vk(t)Qk=t0tvk(t´)dt´+βk


Insgesamt finden sich 2f Konstanten der Bewegung:


αk,βk k=1,...,f

Als Beispiel ( Vergl. Kapitel 3.5) betrachten wir das reduzierte 2-Körper-Problem in der Ebene senkrecht zum Drehimpuls l:


T=12m(x˙2+y˙2+z˙2)=12m(r˙2+r2ϕ˙2)V=V(r)L=L(r,r˙,ϕ,ϕ˙,)=12m(r˙2+r2ϕ˙2)V


ϕ ist zyklisch: Lϕ=0Lϕ˙=mr2ϕ˙=l=cons


Die Hamiltonschen Gleichungen lauten:


p˙k=HqkHpk=q˙kk=1,...,fpr=Lr˙=mr˙r˙=Hpr=prmpϕ=Lϕ˙=mr2ϕ˙ϕ˙=Hpϕ=pϕmr2


H=prr˙+pϕϕ˙L=mr˙2+mr2ϕ˙2L=m2(r˙2+r2ϕ˙2)+V(r)H=pr22m+pϕ22mr2+V(r)


Hϕ=0pϕ=αϕ=l=cons


Somit läßt sich die Hamiltonfunktion von f=2 auf f=1 Freiheitsgrade reduzieren:


H=pr22m+l22mr2+V(r)


Definition der kanonischen Transformationen

Kanonische Transformationen sind diffeomorphe Transformationen (umkehrbar eindeutig und zweimal stetig diffbar): (q¯,p¯)(Q¯,P¯)


H(q¯,p¯,t)H¯(Q¯,P¯,t) , die die Hamilton- Gleichungen forminvariant lassen.

Bedingung für eine kanonische Transformation:

Die Hamiltonschen Gleichungen folgen aus dem Hamiltonschen Prinzip:


δW=0δW=δt1t2dtL=δt1t2dt{k=1fpkq˙k(t)H(q¯,p¯,t)} (Legendre Trafo)

Ganz entsprechend muss für das System (Q¯,P¯,H¯) gelten:


δW=0δW=δt1t2dtL=δt1t2dt{k=1fPkQ˙k(t)H¯(Q¯,P¯,t)}=0


Man kann sich leicht überzeugen, dass diese beiden Forderungen äquivalent sind, falls:


k=1fpkq˙k(t)H(q¯,p¯,t)=k=1fPkQ˙k(t)H¯(Q¯,P¯,t)+ddtM1


Mit einer beliebigen Funktion


M1(q¯,Q¯,t) , die "Erzeugende der kanonischen Trafo" genannt wird.

M1 ist dabei eine Verallgemeinerung der Eichfunktion M(q¯,t) aus dem Kapitel Eichtrafo der Lagrangefunktion (2.3)

Beweis:

ddtM1=k=1f(M1qkq˙k(t)+M1QkQ˙k(t))+M1t


Es folgt dann aus


k=1fpkq˙k(t)H(q¯,p¯,t)=k=1fPkQ˙k(t)H¯(Q¯,P¯,t)+ddtM1 , dass


k=1f(pkM1qk)q˙k(t)=k=1f(Pk+M1Qk)Q˙k(t)+H(q¯,p¯,t)H¯(Q¯,P¯,t)+M1t


Da aber q¯ und Q¯ unabhängige Variablen sind kann obige Gleichung nur für alle denkbaren unabhängigen Variablen erfüllt werden, falls


pk=M1qkPk=M1QkH¯(Q¯,P¯,t)=H(q¯,p¯,t)+M1t


Das bedeutet jedoch, dass die kanonische Transformation durch M1(q¯,Q¯,t) eindeutig bestimmt ist:


pk=M1(q¯,Q¯,t)qkQj(q¯,p¯,t)Bedingung:det(2M1qkQj)0Pk=M1(q¯,Q¯,t)Qk=M1(q¯,Q¯(q¯,p¯,t),t)Qk=Pk(q¯,p¯,t)


Somit kann der Impuls durch die alten Koordinaten Ort,Impuls und zeit ausgedrückt werden und die Abhängigkeit von zeitabhängigkeiten verschwindet. ( Der Ausdruck von Q durch q, p und t ist als Umkehrung der Bestimmung von p zu sehen).

Für die gesamte Umkehrtrafo gilt: Pj=M1Qjinpk=M1qk liefert


qk(Q¯,P¯,t)ausPj=M1Qjundpj(Q¯,P¯,t)auspk=M1qk


Äquivalenzrelation:

δW=δt1t2dtL=δt1t2dt{k=1fpkq˙k(t)H(q¯,p¯,t)}=0 (Legendre Trafo)


δW=δt1t2dtL=δt1t2dt{k=1fPkQ˙k(t)H¯(Q¯,P¯,t)}=0


Beweis:


δt1t2dt{k=1fpkq˙k(t)H(Q¯,P¯,t)}=0δt1t2dt{k=1fPkQ˙k(t)H¯(Q¯,P¯,t)+ddtM1}=0=δt1t2dt{k=1fPkQ˙k(t)H¯(Q¯,P¯,t)}+δ{M1(q(t2),Q(t2),t2)M1(q(t1),Q(t1),t1)}


Dabei gelten die Relationen:


δt1t2dt{k=1fPkQ˙k(t)H¯(Q¯,P¯,t)}=t1t2dt{k=1fδPkQ˙k(t)+PkδQ˙k(t)H¯(Q¯,P¯,t)QkδQkH¯(Q¯,P¯,t)PkδPk}


δ{M1(q(t2),Q(t2),t2)M1(q(t1),Q(t1),t1)}=k(M1qkδqk|t1t2+M1QkδQk|t1t2)mitM1qkδqk|t1t2=0undM1QkδQk|t1t20


Außerdem:


t1t2dtPkδQ˙k(t)=PkδQk|t1t2t1t2dtP˙kδQk


Nebenbemerkung: Für die Variation gilt bekanntlich: δq¯(t1)=δq¯(t2)=0 . Jedoch sind p(t1) und p(t2) beliebig. Dadurch können sich nun insbesondere die Randbedingungen für Q(q¯,p¯,t) ändern.

Unter Beachtung der obigen relationen gilt nun:


0=δt1t2dtL=k=1f(Pk+M1Qk)δQk|t1t2+t1t2dtk=1f{(Q˙k(t)H¯(Q¯,P¯,t)Pk)δPk(P˙k(t)+H¯(Q¯,P¯,t)Qk)δQk}


Aus den obigen Relationen ist bekannt:


(Pk+M1Qk)=0


Gleichzeitig sind jedoch Pi und Qi unabhängig und können demnach unabhängig variiert werden. Das bedeutet, dass δPkundδQk unabhängig sind.

Somit muss jeweils für sich gelten:


0=(Q˙k(t)H¯(Q¯,P¯,t)Pk)0=(P˙k(t)+H¯(Q¯,P¯,t)Qk)


und es sind die Hamiltonschen Gleichungen äquivalent in den neuen Koordinaten, was zu beweisen war.

Äquivalente Formen der erzeugenden Funktion

Eine Legendre- Transformation von M1 liefert:


Aus dem vorigen Beweis ist bekannt:


k=1f(pkq˙kPkQ˙k)(HH¯)=ddtM1


Außerdem gilt:


ddtM1=ddt(M2(q¯(t),P¯(t),t)kPkQk)=k=1f(M2qkq˙k+M2PkP˙kP˙kQkPkQ˙k)+M2t


So dass folgt:


k=1f(pkM2qk)q˙k+(QkM2Pk)P˙k+(PkPk)Q˙k=(HH¯)+M2t


Da dies für beliebige q˙k,P˙k gilt, kann die Summe nur allgemein identisch sein, wenn gilt:


(pkM2qk)=0pk=M2qkQk=M2PkH¯=H+M2t


Analog kann gezeigt werden, dass für


M3(p¯,Q¯,t)=M1(q¯,Q¯,t)k=1fM1qkqk


Hier folgt ( Übungsaufgabe):


(qk+M3pk)=0qk=M3pkPk=M2Qk


oder


M4(p¯,P¯,t)=M1(q¯,Q¯,t)k=1f(M1qkqk+M1QkQk)


(qk+M4pk)=0qk=M4pkQk=M4Pk


Beispiele für kanonische Transformationen

Erzeugende sei:


M1(q¯,Q¯,t)=j=1fqjQjpj=M1qj=QjPj=M1Qj=qj(q¯,p¯)(P¯,Q¯)


Bei dieser Trafo werden also Ort und Impuls vertauscht.

Beispiel 2:


M2(q¯,P¯,t)=j=1fqjPjpj=M2qj=PjQj=M2Pj=qj(q¯,p¯)(Q¯,P¯)


Dies ist also die identische Transformation

Beispiel: Harmonischer Oszillator:

H=p22m+mω22q2M1(q,Q)=mω2q2cotQp=M1q=mωqcotQP=M1Q=mω2q2sin2Qq=(2mωP)12sinQp=(2mωP)12cosQ


H=H¯(M1t)=0H=2mωPcos2Q2m+mω22P2mωsin2Q=ωP


Die Variable Q ist also zyklisch.


P˙=HQ=0P=α=constQ˙=HP=ωQ=ωt+β


Somit kann q(t) durch Integration ( 2 Integrationskonstanten !!) gefunden werden:


q(t)=(2αmω)12sin(ωt+β)


Dabei beschreibt α die Amplitude und β die Phase.

Symplektische Struktur des Phasenraums

Da die kanonischen Transformationen generalisierte Koordinaten und Impulse ineinander transformieren können, sollten q und p nicht gegeneinander ausgezeichnet sein. Um diese Symmetrie des kanonischen Formalismus auszuzeichnen, wird eine neue Notation eingeführt.

Sei zunächst f= 1


x¯:=(qp) ist Vektor im Phasenraum


H,x:=(HqHp) ist Ableitungsvektor


J:=(0110) ist Metrik im Phasenraum ( metrischer Tensor)

In diesem Fall lassen sich die kanonischen Gleichungen vereinfacht schreiben als:


x¯˙:=JH,xJx¯˙=H,xq˙=Hp,p˙=Hq


Leicht läßt sich zeigen:


J2=1J1=JT=J


Verallgemeinerung auf mehr Freiheitsgrade

x¯:=(q1...qfp1...pf)H¯x:=(Hq1...HqfHp1...Hpf)J:=(01f1f0)


Die kanonischen Gleichungen lauten


x¯˙:=JHxJx¯˙=Hx


Beispiel ist ein lineares autonomes System in einer Dimension, also der verallgemeinerte eindimensionale harmonische Oszillator:


x¯˙:=Ax¯=JHx


Diese Gleichung ist abzuleiten aus der Hamiltonfunktion:


H=12(aq2+2bqp+cp2)z.B.a=ω02,b=0,c=1


x¯˙:=(0110)(HqHp)=bq+cpaqbp


Somit ergibt sich eine Einschränkung an die Matrix A:


A=(bcab)tr(A)=0


Dies gilt für Hamiltonsche Systeme ! ( Einschränkung an die Dynamik im Phasenraum)

Kanonische Transformationen in kompakter Notation

Aus den 4 Äquivalenten Formen der Erzeugenden für kanonische Transformationen folgt:

M1(q¯,Q¯,t):pj=M1qjPj=M1QjpjQk=2M1Qkqj=Pkqj


M2(q¯,P¯,t)=M1(q¯,Q¯,t)j=1fM1QjQjpj=M2qjQj=M2PjpjPk=2M2Pkqj=Qkqj


M3(p¯,Q¯,t)=M1(q¯,Q¯,t)j=1fM1qjqjqj=M3pjPj=M3QjqjQk=2M3Qkpj=Pkpj


M4(p¯,P¯,t)=M1(q¯,Q¯,t)j=1f(M1QjQj+M1qjqj)qj=M4pjQj=M1Pj=qjqjPk=2M1Pkpj=Qkpj


Dabei sind:


x¯:=(q1...qfp1...pf)y¯:=(Q1...QfP1...Pf)


Mαβ=xαyβ(M1)αβ:=yαxβα,β=1,...,2f


Beweis:

γ=12fMαγ(M1)γβ=γ=12fxαyγyγxβ=xαxβ=δαβ


Damit läßt sich eine einheitliche Schreibweise finden für die Relationen aller Erzeugenden:


Mαβ=μ,ν=12fJαμJβν(M1)μν


Beweis:

In Matrixform lautet diese Gleichung:


M=J(JM1)T


Die linke Seite (M) lautet:


M=(qQqPpQpP)


Die rechte Seite lautet:


J(JM1)T=(0110)[(0110)(QqQpPqPp)]T=(0110)(PqPpQqQp)T=(0110)((Pq)T(Qq)T(Pp)T(Qp)T)=((Pp)T(Qp)T(Pq)T(Qq)T)


Die Matrixform für die Erzeugenden läßt sich folgendermaßen äquivalent umformen:


M=J(JM1)TJM=(JM1)T=(M1)TJTMTJM=MT(M1)TJT=(M1M)TJT=JT=JMTJM=J


Dabei ist J der metrische Tensor und M die Matrix der 2. Ableitungen der Erzeugenden der kanonischen Transformation, also die Jacobi- Matrix für die Erzeugenden der kanonischen Trafo.

Dies bedeutet jedoch nichts anderes als: Die Metrik im Phasenraum ist invariant unter kanonischen Transformationen !

J definiert dabei eine Metrik über das verallgemeinerte schiefsymmetrische Skalarprodukt:


(x¯,y¯):=x¯TJy¯=i,k=12fxiJikyk


es handelt sich dabei um eine schiefsymmetrische, nichtentartete Bilinearform

Eigenschaften:

  1. Schiefsymmetrie:

(x¯,y¯)=(y¯,x¯) , Beweis: (x¯,y¯)=x¯TJy¯=(y¯TJTx¯)T=y¯TJx¯=(y¯,x¯)

  1. bilinear:

(x¯,λ1y¯1+λ2y¯2)=λ1(x¯,y¯1)+λ2(x¯,y¯2)

  1. nichtentartet:

(x¯,y¯)=0y¯x¯=0


Nebenbemerkung: Es gilt: (x¯,x¯)=0x¯

Also Selbstorthogonalität

Beweis: x¯TJx¯=(qp)(0110)(qp)=qppq=0


Die Symplektische Struktur auf dem R2f ist von einer euklidischen Metrik grundsätzlich zu unterscheiden:


(x¯,y¯)Eu=ixiyi=x¯Tgy¯


Mit dem metrischen Tensor g, einer 2fx2f dimensionalen Einheitsmatrix !

Im Euklidischen gelten jedoch die Relationen:


(x¯,y¯)=(y¯,x¯)(x¯,x¯)0


Definition:

Die Menge der Matrizen M ( kanonische Trafo) mit


MTJM=J bildet die reelle symplektische Gruppe S über R2f .

Dies ist die Symmetriegruppe der symplektischen Struktur.

Gruppeneigenschaften

1. M1,M2SM3=M1M2S


Beweis: M3TJM3=(M1M2)TJ(M1M2)=M2TM1TJM1M2=M2TJM2=J


2. Assoziativität ( matrixmultiplikation !)

3. Einselement Einheitsmatrix !

  1. Inverse:

M1:=J1MTJ

Beweis: M1M=(J1MTJ)M=J1(MTJM)=J1J=1


Dabei gilt : MT,JS Beweis: Übungsaufgabe

  1. Weiter gilt:

detM=1 Beweis: Übungsaufgabe oder Scheck, S. 102

Fazit:

Die Invarianz der kanonischen Gleichungen x¯˙:=Ax¯=JH¯,x kann durch di symplektische Struktur des Phasenraums beschrieben werden:


y˙i=kyixkx˙ky¯˙=M1x¯˙=(J1MTJ)JH¯,xH¯yi=kH¯xkxkyiH¯,y=MTH¯,xy¯˙=(J1MTJ)J(MT)1H¯,y=J(1)MT(MT)1H¯,y=JH¯,y


Der Satz von Liouville

Lösung der Differenzialgleichung x¯˙:=Ax¯=JH¯,x


x¯(t,t0,x¯0)=:Φ¯t,t0(x¯0) Definition: Fluß im Phasenraum

to und xo beschreibt die Anfangskonfiguration und Phi den Fluß.

Der Fluß beschreibt dabei die Zeitentwicklung der Anfangskonfiguration:


Φ¯t,t0(x¯0):ΓΓx¯0(t0)x¯(t)


Dies entspricht einer Kurvenschar, die durch die Zeit parametrisiert ist:

Beispiel: eindimensionaler harmonischer Oszi:


(q,p)Γq˙=pp˙=ω02qx¯˙=Ax¯A=(01ω020)


Die Lösung lautet:


x¯(t)=Φ¯tt0(x¯0)=exp[(tt0)A]x¯0


Dies ist also gerade das Exponenzial der Matrix.

Aufschluss liefert eine Reihenentwicklung:


x¯(t)=n[(tt0)A]nn!x¯0=[1cosω0(tt0)+Aω0sinω0(tt0)]x¯0=(cosω0(tt0)1ω0sinω0(tt0)ω0sinω0(tt0)cosω0(tt0))x¯0


Beweis:


A2=(01ω020)(01ω020)=ω021A2n=(1)2nω02n1A2n+1=(1)nω02n+11ω0A


Als Ergebnis erhalten wir, dass alle Phasenpunkte mit gleicher, konstanter Winkelgeschwindigkeit wo, rotieren: Ein Ensemble von Anfangskonfigurationen Uto läuft zum Zeitpunkt Ut insbesondere nicht auseinander.

Das bedeutet, das Gebiet Uto wandert ohne Änderung der Form und Orientierung um den Nullpunkt:

Man erhält als markantes Ergebnis, dass das Phasenvolumen bei der Zeitentwicklung erhalten ist. Im Allgemeinen ändert sich zwar die Form, stets gilt jedoch der Liouvillesche Satz:

Bei der Hamiltonschen Zeitentwicklung ist das Phasenvolumen erhalten ( auch seine Orientierung). Der Fluß im Phasenraum ist also divergenzfrei.

Beweis ( integrale Form):

Gegeben sei eine Menge von Anfangskonfigurationen (to), die das Phasenraumgebiet Uto mit dem Volumen Vto ausfüllen:


Vto=Utod2fx0 Bei t: Vt=Utd2fx0=Ut0d2fx0det(xx0)=Ut0d2fx0det(DΦt,t0(x¯0))


Mit der Jacobi- Matrix:


(DΦt,t0(x¯0))ik:=Φit,t0(x¯0)x0k=xix0k


Dies kann für Zeiten nahe t0 reihenentwickelt werden:


Φt,t0(x¯0)=x¯0+F¯(x¯0,t)(tt0)+O((tt0)2)F¯(x¯0,t)=JH¯,x=(HpHq)


Somit folgt:


Φit,t0(x¯0)x0k=δik+F¯i(x¯0,t)x0k(tt0)+O((tt0)2)


Mit Hilfe det(1+Bε)=1+εtr(B)+O(ε2) folgt:


det(DΦt,t0)=|Φit,t0(x¯0)x0k|=1+(tt0)i=12fF¯i(x¯0,t)x0i(tt0)+O((tt0)2)i=12fF¯i(x¯0,t)x0i=divF¯=qHppHq=0


Der Fluß im Phasenraum ist also divergenzfrei. Dann folgt jedoch für die Jacobideterminante:


det(DΦt,t0)=|Φit,t0(x¯0)x0k|=1+O((tt0)2)1


Vt=Ut0d2fx0det(DΦt,t0(x¯0))=Ut0d2fx0(1+O(tt0)2)Vt0


Nebenbemerkung:

Der Satz von Liouville kann auch in der LOKALEN Form formuliert werden:

Für den Fluß Φt,t0 zu x¯˙:=JH¯,x ist DΦt,t0 eine symplektische Matrix, das heißt det(DΦt,t0)=1 .

Das bedeutet, das Volumenelement dx1...dx2f im Phasenraum ist unter dem Fluß invariant:


dx1...dx2f=Det(DΦ)dx01....dx02f=dx01....dx02f


Beispiel: eindimensionaler harmonischer Oszi

xi(t)=k=12Pikx0kP=(cosω0(tt0)1ω0sinω0(tt0)ω0sinω0(tt0)cosω0(tt0))detP=cos2ω0(tt0)+sin2ω0(tt0)=1dx1dx2=(detP)dx01dx02=dx01dx02


Poisson- Klammern

Jede Observable läßt sich in der klassischen Mechanik als Funktion von Ort, Impuls und Zeit darstellen:


Observable=g(q¯,p¯,t)


Die zeitliche Änderung längs der Bahn q¯(t),p¯(t) im Phasenraum Γ


dg(q¯,p¯,t)dt=i=1f(gqiq˙i+gpip˙i)+gt=i=1f(gqiHpigpiHqi)+gt=:{g,H}+gt


Definition:

Für zwei beliebige Observablen g(q¯,p¯,t) und f(q¯,p¯,t) heißt


i=1f(gqifpigpifqi)=:{g,f}


Poisson- Klammer

Eigenschaften

  1. die Poissonklammer ist eine schiefsymmetrische nicht entartete Bilinearform. Das bedeutet jedoch, sie definiert ein symplektisches Skalarprodukt im Phasenraum:


{g,f}=(f¯x,g¯x)=f¯xTJg¯x=i,k=1f(fxiJikgxk)=(fqfp)(0110)(gqgq)=(fqfp)(gpgq) Aufgrund der schiefsymmetrischen Struktur und der Bilinearität sowie der Nichtentartung und der daraus folgenden Selbstorthogonalität gilt:

1.Schiefsymmetrie: {f,g}={g,f}


2.bilinear: {f,λ1g1+λ2g2}=λ1{f,g1}+λ2{f,g2}


3.nichtentartet: (f,g)=0gf=const. (Nullelement, wegen f¯,x=0 )

Nebenbemerkung: Es gilt: {f,f}=0f

Also Selbstorthogonalität

Weiter gilt die Produktregel ( Leibnizregel): {f,gh}=g{f,h}+{f,g}h


Die Jacobi- Identität: {f,{g,h}}={{f,g},h}+{g,{f,h}}


Weiter gilt: gqk={g,pk}gpk={g,qk}


Die Poissonklammer ist invariant unter kanonischen Transformationen:

Beweis: Trafo: x->y

Die Jacobi- Determinante Mαβ=xαyβ ist symplektische Matrix,

das heißt , es gilt: MTJM=JM1J(M1)T=J , da ja M1J(M1)T=J1MTJJ(M1)T=J1MT(M1)T=J1=J


Nun muss man umrechnen von :


fxi=kfykykxi=kMki1fykf¯x=(M1)Tf¯yf¯xT=f¯yT(M1)f¯xTJg¯x=f¯yT(M1)J(M1)Tg¯y=f¯yTJg¯y


Also:


(f¯x,g¯x)=(f¯y,g¯y)


Für nicht explizit zeitabhängige Observable g(q¯,p¯) gilt:


dgdt={g,H}


g ist genau dann Bewegungskonstante, wenn gilt:


{g,H}=0


Speziallfall: g ist Koordinate der Impuls: g=qk,g=pk


q˙k={qk,H}p˙k={pk,H}


So folgen die Hamiltonschen Gleichungen

Kompakt kann geschrieben werden:


x˙k={xk,H}=i,j=1fxkxiJijHxj=j=1fJkjHxjalso:x¯˙=JH¯,x


Fundamentale Poisson- Klammern:


{qk,qj}=0{pk,pj}=0{qk,pj}=δkj


Kompakt:


{xk,xj}=l,mxkxlJlmxjxm=Jkj(0110)


Die Poissonklammer ist invariant unter kanonischen Transformationen, da


MTJM=J


Jedoch ist auch die Umkehrung richtig: ist die Transformation kanonisch, so gelten die obigen Poissonklammer- Beziehungen.

Somit:

Satz: Die Transformation (q¯,p¯)(Q¯,P¯) ist genau dann kanonisch, wenn :


{Qk,Qj}=0{Pk,Pj}=0{Qk,Pj}=δkj


Beweis: Zur Vereinfachung: Nicht explizit zeitabhängige Trafos: Mt=0H¯=H


Bewegungsgleichung:


y˙k={yk,H}=i,j=1fykxiJijHxj=j=1fJkjHxj


Wegen (f¯x,g¯x)=(f¯y,g¯y)

kann nun die Bewegungsgleichung in den alten Koordinaten gebildet werden:


y˙k={yk,H}=i,j,l=1fykxiJijH¯ylylxj=l=1fH¯yli,j=1fykxiJijylxj=l=1fH¯yl{yk,yl}


Also folgt:


y˙k={yk,H}=i,j,l=1fykxiJijH¯ylylxj=l=1fH¯yli,j=1fykxiJijylxj=l=1fH¯yl{yk,yl}y˙k=l=1fJklH¯ylJkl={yk,yl}


Mit der Bedeutung


y˙k=l=1fJklH¯yl Hamiltonsche Bewegungsgleichung in den neuen Koordinaten -> Trafo kanonisch


Jkl={yk,yl}

 fundamentale Poissonklammern in den neuen Koordinaten

Somit ergibt sich ein einfach nachprüfbares Kriterium für kanonische Transformationen !

Folgende Aussagen sind äquivalent:


x¯=(q¯p¯)y¯=(Q¯P¯) ist kanonisch


die kanonischen Gleichungen x¯˙=JH¯,x sind invariant


die Poissonklammern {f,g} sind invariant für alle f und g


die fundamentalen Poissonklammern Jkl={xk,xl} sind ivariant


die Jacobi- Matrix Mαβ=xαxβ ist symplektisch, das heißt MTJM=J


es existiert eine Erzeugende !

Bezug zur Quantenmechanik

Ein Übergang zur Quantenmechanik ist möglich:

Von der klassischen Variablen g(q¯,p¯,t) zum qm. Operator: g:H>H mit dem Hilbertraum H

Von der Poissonklammer: {f,g}1i[f,g] zum Kommutator

Aus den fundamentalen Poisson- Klammern folgen die kanonischen Vertauschiungsrelationen:


{qk,qj}=0[qk,qj]=0{pk,pj}=0[pk,pj]=0{qk,pj}=δkj[qk,pj]=iδkj


Die Hamiltonfunktion H(q,p,t) geht über zum Hamilton- Operator

Die Bewegungsgleichungen:


dg(q¯,p¯,t)dt=i=1f(gqiq˙i+gpip˙i)+gt=i=1f(gqiHpigpiHqi)+gt=:{g,H}+gtdgdt=1i[g,H]+gt


Wobei auch nur der Zusammenhang zwischen Poisson- Klammer und Kommutator recycled wurde.

Da in diesem Bild die Operatoren zeitabhängig sind haben wir es mit der Heisenbergschen bewegungsgleichung zu tun. Im Schrödingerbild ist der Operator zeitunabhängig und die Schrödingergleichung gibt eine Bewegungsgleichung für die Zustände an.