Der Hamiltonsche kanonische Formalismus

Der Hamiltonsche Kanonische Formalismus

Motivation

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

$\begin{array}{rl}& L(q_1,...,q_f,{\dot{q}}_1,...,{\dot{q}}_f,t)\\ & \Rightarrow \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_k}-\frac{\partial L}{\partial q_k}=0\\ \end{array}$ k=1,..,f

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

Bei gewissen Problemstellungen, wenn es beispielsweise zyklische Variablen gibt:

$\frac{\partial L}{\partial q_k}=0\Rightarrow \frac{\partial L}{\partial {\dot{q}}_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:

$p_k:=\frac{\partial L}{\partial {\dot{q}}_k}$

Die erforderliche Variablentransformation

$(q_k,{\dot{q}}_k,t)\to (q_k,p_k,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)\cong L,x\cong {\dot{q}}_k$

Soll statt x die Variable $u=\frac{d}{dx}f(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=\frac{d}{dx}f(x)\Rightarrow x=\varphi (u)$

Mit der Voraussetzung

$\frac{d^2}{d{x}^2}f(x)\ne 0$ Ansonsten: $u=\frac{d}{dx}f(x)=const.$ und damit nicht umkehrbar.

Die Substitution in f führt dann auf:

$y=f(x)=f(\varphi (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=\frac{d}{dx}y=\frac{d}{dx}(f(x)+a)=\frac{df}{dx}$

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:

$\begin{array}{rl}& x\to u\\ & y\to z=xu-f(x)=\varphi (u)u-f(\varphi (u))=g(u)\\ \end{array}$

Statt der gerade genannten einfachen Variante

$y=f(\varphi (u))$

Die Trafo

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

Graphische Veranschaulichung von $g(u)=x\frac{df}{dx}-f(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

$\begin{array}{rl}& y\cong L\\ & x\cong \dot{q}\\ & u\cong p=\frac{\partial L}{\partial \dot{q}}\\ & z\cong \dot{q}p-L:=H(q,p,t)\\ \end{array}$

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:

$\begin{array}{rl}& L(q_1,...,q_f,{\dot{q}}_1,...,{\dot{q}}_f,t)\\ & p_k:=\frac{\partial L}{\partial {\dot{q}}_k}\\ & H(q_1,...,q_f,p_1,...,p_f,t)=\sum _{k=1}^f{\dot{q}}_k{p}_k-L\\ \end{array}$

Die mathematischen Voraussetzungen für diese Prozedur sind:

$\begin{array}{rl}& L(q_1,...,q_f,{\dot{q}}_1,...,{\dot{q}}_f,t)\in C^2\\ & \det \frac{{\partial }^{2}L}{\partial {\dot{q}}_k\partial {\dot{q}}_l}\ne 0\\ \end{array}$ damit $p_{k=}\frac{\partial L}{\partial {\dot{q}}_k}$ nach ${\dot{q}}_k$ auflösbar

Die Hamiltonschen Gleichungen

Ziel: Auch hier natürlich sollen Bewegungsgleichungen für die $p_k,q_k$ gefunden werden.

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

Eine Variable:

Differenziale:

$\begin{array}{rl}& L=L(q,\dot{q},t):dL=\frac{\partial L}{\partial q}dq+\frac{\partial L}{\partial \dot{q}}d\dot{q}+\frac{\partial L}{\partial t}dt\\ & H=H(q,p,t):dH=\frac{\partial H}{\partial q}dq+\frac{\partial H}{\partial p}dp+\frac{\partial H}{\partial t}dt\\ & H=\dot{q}p-L\Rightarrow dH=\dot{q}dp+pd\dot{q}-dL=\dot{q}dp+pd\dot{q}-\frac{\partial L}{\partial q}dq-\frac{\partial L}{\partial \dot{q}}d\dot{q}-\frac{\partial L}{\partial t}dt\\ \end{array}$

wegen

$\begin{array}{rl}& \frac{\partial L}{\partial \dot{q}}=p\\ & \Rightarrow dH=\frac{\partial H}{\partial q}dq+\frac{\partial H}{\partial p}dp+\frac{\partial H}{\partial t}dt=\dot{q}dp-\frac{\partial L}{\partial q}dq-\frac{\partial L}{\partial t}dt\\ \end{array}$

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

$\begin{array}{rl}& \frac{\partial H}{\partial q}=-\frac{\partial L}{\partial q}\\ & \frac{\partial H}{\partial p}=\dot{q}\\ & \frac{\partial H}{\partial t}=-\frac{\partial L}{\partial t}\\ \end{array}$

Mit Hilfe der Lagrange Bewegungsgleichung

$\begin{array}{rl}& \frac{d}{dt}\frac{\partial L}{\partial \dot{q}}=\frac{\partial L}{\partial q};\frac{\partial L}{\partial \dot{q}}=p;\frac{\partial L}{\partial q}=-\frac{\partial H}{\partial q}\\ & \Rightarrow \dot{p}=-\frac{\partial H}{\partial q}\\ & \frac{\partial H}{\partial p}=\dot{q}\\ \end{array}$

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

$\begin{array}{rl}& L(q_1,...,q_f,{\dot{q}}_1,...,{\dot{q}}_f,t)\\ & p_k:=\frac{\partial L}{\partial {\dot{q}}_k}\\ & H(q_1,...,q_f,p_1,...,p_f,t)=\sum _{k=1}^f{\dot{q}}_k{p}_k-L\\ & \\ \end{array}$

$\begin{array}{rl}& dH=\sum _k(\frac{\partial H}{\partial q_k}d{q}_k+\frac{\partial H}{\partial p_k}d{p}_k)+\frac{\partial H}{\partial t}dt=\sum _k{\dot{q}}_{k}d{p}_k+\sum _k{p}_{k}d{\dot{q}}_k-\sum _k\frac{\partial L}{\partial {\dot{q}}_k}d{\dot{q}}_k-\sum _k\frac{\partial L}{\partial q_k}d{q}_k-\frac{\partial L}{\partial t}dt\\ & =\sum _k{\dot{q}}_{k}d{p}_k-\sum _k\frac{\partial L}{\partial q_k}d{q}_k-\frac{\partial L}{\partial t}dt\\ & \Rightarrow \frac{\partial H}{\partial q_k}=-\frac{\partial L}{\partial q_k}\\ & \frac{\partial L}{\partial q_k}=\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_k}=\frac{d}{dt}{p}_k\\ & \frac{\partial H}{\partial p_k}={\dot{q}}_k;\frac{\partial H}{\partial t}=-\frac{\partial L}{\partial t}\\ \end{array}$

Somit folgen hier die Hamiltonschen Gleichungen (Kanonische Gleichungen)

$\begin{array}{rl}& {\dot{p}}_k=-\frac{\partial H}{\partial q_k}\\ & \frac{\partial H}{\partial p_k}={\dot{q}}_{k}k=1,...,f\\ \end{array}$

Der 2f- dimensionale Raum

$\mathrm{\Gamma }:=\{q_1,...,q_f,{\dot{q}}_1,...,{\dot{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 $\frac{\partial }{\partial t}{\overline{r}}_i(q_1,...,q_f)=0und\frac{\partial }{\partial t}L=0$

Dann nämlich ist

$\sum _{k=1}^f\frac{\partial T}{\partial {\dot{q}}_k}{\dot{q}}_k=2T$

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

${\dot{q}}_k$ .

Somit:

$\sum _{k=1}^f\frac{\partial L}{\partial {\dot{q}}_k}{\dot{q}}_k-L=\sum _{k=1}^f\frac{\partial T}{\partial {\dot{q}}_k}{\dot{q}}_k-L=2T-T+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:

$\begin{array}{rl}& \frac{dH}{dt}=\frac{d}{dt}(\sum _{k=1}^f\frac{\partial L}{\partial {\dot{q}}_k}{\dot{q}}_k-L)=\sum _{k=1}^f(\frac{\partial H}{\partial q_k}{\dot{q}}_k+\frac{\partial H}{\partial p_k}{\dot{p}}_k)+\frac{\partial H}{\partial t}=\sum _{k=1}^f(\frac{\partial H}{\partial q_k}\frac{\partial H}{\partial p_k}-\frac{\partial H}{\partial p_k}\frac{\partial H}{\partial q_k})-\frac{\partial L}{\partial t}=0\\ & wegen\frac{\partial L}{\partial t}=0\\ \end{array}$ 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=-m{q}^2{\omega }^2$ .

Aus $\frac{\partial L}{\partial t}=0$ folgt dann ohnehin wieder ein Erhaltungssatz: H=const.

Typisches Anwendungsschema des Hamilton- Formalismus:

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

$\overline{q}=(q_1,...,q_f)$

1. Transformation des Radiusvektors

$\begin{array}{rl}& {\overline{r}}_i={\overline{r}}_i(q_1,...,q_f,t)\\ & {\dot{\overline{r}}}_i={\dot{\overline{r}}}_i(\overline{q},\dot{\overline{q}},t)\\ \end{array}$

1. Aufstellung der Lagrangegleichung:

$L(\overline{q},\dot{\overline{q}},t)=T-V=\frac{1}{2}m\sum _i{{\dot{\overline{r}}}_i}^2-V$

1. Bestimmung der generalisierten Impulse:

$\begin{array}{rl}& p_k:=\frac{\partial L}{\partial {\dot{q}}_k}\Rightarrow p_k=p_k(\overline{q},\dot{\overline{q}},t)\\ & Umkehrung:{\dot{q}}_k={\dot{q}}_k(\overline{q},\overline{p},t)\\ \end{array}$

1. Anschließend Legendre Trafo:

$H(q_1,...,q_f,p_1,...,p_f,t)=\sum _{k=1}^f{\dot{q}}_k{p}_k-L$

1. Aufstellung und Integration der kanonischen Gleichungen:

$\begin{array}{rl}& {\dot{p}}_k=-\frac{\partial H}{\partial q_k}\\ & \frac{\partial H}{\partial p_k}={\dot{q}}_{k}k=1,...,f\\ \end{array}$

Beispiele:

Teilchen in Zylinderkoordinaten ganz ohne Zwnagsbedingungen

1. q1=3, q2=Phi, q3 = z

$\begin{array}{rl}& x=r\cos \varphi ,\dot{x}=\dot{r}\cos \varphi -r\dot{\varphi }\sin \varphi \\ & y=r\sin \varphi ,\dot{y}=\dot{r}\sin \varphi +r\dot{\varphi }\cos \varphi \\ & z=z\\ \end{array}$

$\begin{array}{rl}& T=\frac{1}{2}m({\dot{x}}^2+{\dot{y}}^2+{\dot{z}}^2)=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2+{\dot{z}}^2)\\ & V=V(r,\varphi ,z)\\ & L=L(r,\varphi ,z,\dot{r},\dot{\varphi },\dot{z})=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2+{\dot{z}}^2)-V\\ \end{array}$

1. Generalisierte Impulse:

$\begin{array}{rl}& p_k=\frac{\partial L}{\partial {\dot{q}}_k}\\ & p_r=m\dot{r}\\ & p_{\varphi }=m{r}^2\dot{\varphi }\\ & p_z=m\dot{z}\\ \end{array}$

                                                                                            Radialimpuls, z-Komponente des Drehimpulses und z-Komponente des Impulses

1. Aufstellung der Legendretrafo:

$\begin{array}{rl}& H=m{\dot{r}}^2+m{r}^2{\dot{\varphi }}^2+m{\dot{z}}^2=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2+{\dot{z}}^2)+V\\ & H=\frac{1}{2m}({p_r}^2+\frac{{p_{\varphi }}^2}{r^2}+{p_z}^2)+V(r,\varphi ,z)\\ \end{array}$

1. Kanonische Gleichungen:

$\begin{array}{rl}& {\dot{p}}_k=-\frac{\partial H}{\partial q_k}\\ & \frac{\partial H}{\partial p_k}={\dot{q}}_{k}k=1,...,f\\ & \dot{r}=\frac{\partial H}{\partial p_r}=\frac{p_r}{m},\dot{\varphi }=\frac{\partial H}{\partial p_{\varphi }}=\frac{p_{\varphi }}{m{r}^2},\dot{z}=\frac{\partial H}{\partial p_z}=\frac{p_z}{m}\\ & {\dot{p}}_r=-\frac{\partial H}{\partial r}=\frac{{p_{\varphi }}^2}{m{r}^3}-\frac{\partial V}{\partial r},{\dot{p}}_{\varphi }=-\frac{\partial H}{\partial \varphi }=-\frac{\partial V}{\partial \varphi },{\dot{p}}_z=-\frac{\partial H}{\partial z}=-\frac{\partial V}{\partial z}\\ \end{array}$

Interessant ist das Ergebnis der Zentrifugalkraft ( Scheinkraft):

F(Zentrifugal)= $\frac{{p_{\varphi }}^2}{m{r}^3}$ , die den radialen Impuls ändert.

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

${\dot{p}}_r=-\frac{\partial H}{\partial r}=\frac{{p_{\varphi }}^2}{m{r}^3}-\frac{\partial V}{\partial r},{\dot{p}}_{\varphi }=0,{\dot{p}}_z=0$

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

$z,\varphi$ sind zyklische Variablen

$\begin{array}{rl}& p_z=const.=o.B.d.A.=0\\ & p_{\varphi }=const.\\ \end{array}$ oBdA: ebene Bewegung, Drehimpulserhaltung in der Ebene

Beispiel: eindimensionaler harmonischer Oszi:

Das System ist skleronom wegen $\frac{\partial L}{\partial t}=0$ , also folgt Energieerhaltung: E=H=T+V

$\frac{1}{2}m({\dot{q}}^2+{{\omega }_o}^2{q}^2)=E=\frac{1}{2}m(\frac{p^2}{m^2}+{{\omega }_o}^2{q}^2)\Rightarrow \frac{p^2}{2mE}+\frac{q^2}{\left(\frac{2E}{m{{\omega }_o}^2}\right)}=1$

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

Die Halbachsen sind:

$\sqrt{2mE}=a,b=\sqrt{\frac{2E}{m{{\omega }_o}^2}}$ ( bestimmt durch 1. Integral).

Als kanonische Gleichungen ergibt sich:

$\begin{array}{rl}& {\dot{p}}_{}=-\frac{\partial H}{\partial q}=-m{{\omega }_o}^{2}q\\ & \dot{q}=\frac{\partial H}{\partial p}=\frac{p}{m}\\ \end{array}$

Daraus folgt dann gerade die Bewegungsgleichung $\begin{array}{rl}& \ddot{q}=\frac{d}{dt}\frac{\partial H}{\partial q}=\frac{\dot{p}}{\stackrel{´}{}m}=-{{\omega }_o}^{2}q\\ & \ddot{q}+{{\omega }_o}^{2}q=0\\ \end{array}$

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(\overline{q},\dot{\overline{q}},t)=T-V=\frac{1}{2}m\sum _i{{\dot{\overline{r}}}_i}^2-V=\frac{m}{2}{\dot{\overline{q}}}^2+e(\dot{\overline{q}}\cdot \overline{A}(\overline{q},t)-\mathrm{\Phi }(\overline{q},t))$

die kanonischen konjugierten Impulse lauten:

$\begin{array}{rl}& p_k=\frac{\partial L(\overline{q},\dot{\overline{q}},t)}{\partial {\dot{q}}_k}=m{\dot{q}}_{\stackrel{´}{}k}+e{A}_k(\overline{q},t)\\ & \Rightarrow {\dot{q}}_k=\frac{1}{m}(p_k-e{A}_k)\\ & H=\sum _{k=1}^3{p}_k{\dot{q}}_k-L=\sum _{k=1}^3{p}_k\frac{1}{m}(p_k-e{A}_k)-\frac{1}{2m}\sum _{k=1}^3{(p_k-e{A}_k)}^2-\sum _{k=1}^3\frac{e}{m}(p_k-e{A}_k)A_k+e\mathrm{\Phi }\\ & H(\overline{q},\overline{p},t)=\frac{1}{2m}{({\overline{p}}_{}-e\overline{A}{(\overline{q},t)}_{})}^2+e\mathrm{\Phi }(\overline{q},t)\\ \end{array}$

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

$m\dot{\overline{q}}=\overline{p}-e\overline{A}$ als kinetischer Impuls ( der auch tatsächlich mit der Geschwindigkeit verknüpft ist).

$p_k=\frac{\partial L}{\partial {\dot{q}}_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:

$\overline{q}=(q_1,...,q_f)\to \overline{Q}=(Q_1,...,Q_f)$

Dabei gilt dann:

$\overline{L}(\overline{Q},\dot{\overline{Q}},t)=L(\overline{q}(\overline{Q},t),\dot{\overline{q}}(\overline{Q},\dot{\overline{Q}},t),t)$

Nun kann man sich fragen, unter welchen Transformationen $(\overline{q},\overline{p})\to (\overline{Q},\overline{P})$

die Hamiltonfunktionen forminvariant sind, also:

mit $\begin{array}{rl}& {\dot{p}}_k=-\frac{\partial H}{\partial q_k}\\ & {\dot{q}}_k=\frac{\partial H}{\partial p_k}\\ \end{array}$ soll auch $\begin{array}{rl}& {\dot{P}}_k=-\frac{\partial \overline{H}}{\partial Q_k}\\ & {\dot{Q}}_k=\frac{\partial \overline{H}}{\partial P_k}\\ \end{array}$

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:

$\frac{\partial L}{\partial q_j}=0\Rightarrow \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_j}=\frac{d}{dt}{p}_j=0\Rightarrow p_j=\frac{\partial L}{\partial {\dot{q}}_j}=const$

Allerdings ist damit keine Aussage über ${\dot{q}}_j$ gemacht. Diese muss natürlich weiter als Variable behandelt werden.

Hamilton-Gleichungen:

In $H(q_1,...,q_f,p_1,...,p_f,t)$ heißt $q_j$ zyklisch, wenn

$\frac{\partial H}{\partial q_j}=0\Rightarrow -\frac{\partial L}{\partial q_j}=\frac{d}{dt}{p}_j=0\Rightarrow p_j:={\alpha }_j=const$

Das bedeutet nun, dass $q_j$ in H gar nicht auftritt. $p_j$ kann dagegen durch die Bewegungskonstante

${\alpha }_j$ ersetzt werden:

$H(q_1,...,q_{j-1},q_{j+1},...,q_f,p_1,...,p_{j-1},{\alpha }_j,p_{j+1},...,p_f,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 $(\overline{q},\overline{p})\to (\overline{Q},\overline{P})$

einführt, bis alle 

$\overline{Q}$ zyklisch sind:

$H=H(P_1,...,P_f,t)$ mit $P_k={\alpha }_k=const.$

$\begin{array}{rl}& {\dot{Q}}_k=\frac{\partial H}{\partial P_k}=:v_k(t)\\ & \Rightarrow Q_k=\underset{t_0}{\overset{t}{\int }}{v}_k(t\stackrel{´}{})dt\stackrel{´}{}+{\beta }_k\\ \end{array}$

Insgesamt finden sich 2f Konstanten der Bewegung:

${\alpha }_k,{\beta }_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:

$\begin{array}{rl}& T=\frac{1}{2}m({\dot{x}}^2+{\dot{y}}^2+{\dot{z}}^2)=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2)\\ & V=V(r)\\ & L=L(r,\dot{r},\varphi ,\dot{\varphi },)=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2)-V\\ \end{array}$

$\varphi$ ist zyklisch: $\frac{\partial L}{\partial \varphi }=0\Rightarrow \frac{\partial L}{\partial \dot{\varphi }}=m{r}^2\dot{\varphi }=l=cons$

Die Hamiltonschen Gleichungen lauten:

$\begin{array}{rl}& {\dot{p}}_k=-\frac{\partial H}{\partial q_k}\\ & \frac{\partial H}{\partial p_k}={\dot{q}}_{k}k=1,...,f\\ & p_r=\frac{\partial L}{\partial \dot{r}}=m\dot{r}\Rightarrow \dot{r}=\frac{\partial H}{\partial p_r}=\frac{p_r}{m}\\ & p_{\varphi }=\frac{\partial L}{\partial \dot{\varphi }}=m{r}^2\dot{\varphi }\Rightarrow \dot{\varphi }=\frac{\partial H}{\partial p_{\varphi }}=\frac{p_{\varphi }}{m{r}^2}\\ \end{array}$

$\begin{array}{rl}& H=p_r\dot{r}+p_{\varphi }\dot{\varphi }-L=m{\dot{r}}^2+m{r}^2{\dot{\varphi }}^2-L=\frac{m}{2}({\dot{r}}^2+r^2{\dot{\varphi }}^2)+V(r)\\ & H=\frac{{p_r}^2}{2m}+\frac{{p_{\varphi }}^2}{2m{r}^2}+V(r)\\ \end{array}$

$\frac{\partial H}{\partial \varphi }=0\Rightarrow p_{\varphi }={\alpha }_{\varphi }=l=cons$

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

$H=\frac{{p_r}^2}{2m}+\frac{l^2}{2m{r}^2}+V(r)$

Definition der kanonischen Transformationen

Kanonische Transformationen sind diffeomorphe Transformationen (umkehrbar eindeutig und zweimal stetig diffbar): $(\overline{q},\overline{p})\to (\overline{Q},\overline{P})$

$H(\overline{q},\overline{p},t)\to \overline{H}(\overline{Q},\overline{P},t)$ , die die Hamilton- Gleichungen forminvariant lassen.

Bedingung für eine kanonische Transformation:

Die Hamiltonschen Gleichungen folgen aus dem Hamiltonschen Prinzip:

$\begin{array}{rl}& \delta W=0\\ & \delta W=\delta \underset{t_1}{\overset{t_2}{\int }}dtL=\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{p}_k{\dot{q}}_k(t)-H(\overline{q},\overline{p},t)\}\\ \end{array}$ (Legendre Trafo)

Ganz entsprechend muss für das System $(\overline{Q},\overline{P},\overline{H})$ gelten:

$\begin{array}{rl}& \delta W=0\\ & \delta W=\delta \underset{t_1}{\overset{t_2}{\int }}dtL=\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)\}=0\\ \end{array}$

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

$\sum _{k=1}^f{p}_k{\dot{q}}_k(t)-H(\overline{q},\overline{p},t)=\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)+\frac{d}{dt}{M}_1$

Mit einer beliebigen Funktion

$M_1(\overline{q},\overline{Q},t)$ , die "Erzeugende der kanonischen Trafo" genannt wird.

M1 ist dabei eine Verallgemeinerung der Eichfunktion $M(\overline{q},t)$ aus dem Kapitel Eichtrafo der Lagrangefunktion (2.3)

Beweis:

$\frac{d}{dt}{M}_1=\sum _{k=1}^f(\frac{\partial M_1}{\partial q_k}{\dot{q}}_k(t)+\frac{\partial M_1}{\partial Q_k}{\dot{Q}}_k(t))+\frac{\partial M_1}{\partial t}$

Es folgt dann aus

$\sum _{k=1}^f{p}_k{\dot{q}}_k(t)-H(\overline{q},\overline{p},t)=\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)+\frac{d}{dt}{M}_1$ , dass

$\sum _{k=1}^f(p_k-\frac{\partial M_1}{\partial q_k}){\dot{q}}_k(t)=\sum _{k=1}^f(P_k+\frac{\partial M_1}{\partial Q_k}){\dot{Q}}_k(t)+H(\overline{q},\overline{p},t)-\overline{H}(\overline{Q},\overline{P},t)+\frac{\partial M_1}{\partial t}$

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

$\begin{array}{rl}& p_k=\frac{\partial M_1}{\partial q_k}\\ & P_k=-\frac{\partial M_1}{\partial Q_k}\\ & \overline{H}(\overline{Q},\overline{P},t)=H(\overline{q},\overline{p},t)+\frac{\partial M_1}{\partial t}\\ \end{array}$

Das bedeutet jedoch, dass die kanonische Transformation durch $M_1(\overline{q},\overline{Q},t)$ eindeutig bestimmt ist:

$\begin{array}{rl}& p_k=\frac{\partial M_1(\overline{q},\overline{Q},t)}{\partial q_k}\Rightarrow Q_j(\overline{q},\overline{p},t)\\ & Bedingung:\det \left(\frac{{\partial }^2{M}_1}{\partial q_k\partial Q_j}\right)\ne 0\\ & P_k=-\frac{\partial M_1(\overline{q},\overline{Q},t)}{\partial Q_k}=-\frac{\partial M_1(\overline{q},\overline{Q}(\overline{q},\overline{p},t),t)}{\partial Q_k}=P_k(\overline{q},\overline{p},t)\\ \end{array}$

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: $P_j=-\frac{\partial M_1}{\partial Q_j}in{p}_k=\frac{\partial M_1}{\partial q_k}$ liefert

$\begin{array}{rl}& q_k(\overline{Q},\overline{P},t)aus{P}_j=-\frac{\partial M_1}{\partial Q_j}\\ & und\\ & p_j(\overline{Q},\overline{P},t)aus{p}_k=\frac{\partial M_1}{\partial q_k}\\ \end{array}$

Äquivalenzrelation:

$\delta W=\delta \underset{t_1}{\overset{t_2}{\int }}dtL=\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{p}_k{\dot{q}}_k(t)-H(\overline{q},\overline{p},t)\}=0$ (Legendre Trafo)

$\iff \delta W=\delta \underset{t_1}{\overset{t_2}{\int }}dtL=\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)\}=0$

Beweis:

$\begin{array}{rl}& \delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{p}_k{\dot{q}}_k(t)-H(\overline{Q},\overline{P},t)\}=0\\ & \iff \delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)+\frac{d}{dt}{M}_1\}=0\\ & =\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)\}+\delta \{M_1(q(t_2),Q(t_2),t_2)-M_1(q(t_1),Q(t_1),t_1)\}\\ \end{array}$

Dabei gelten die Relationen:

$\delta \underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f{P}_k{\dot{Q}}_k(t)-\overline{H}(\overline{Q},\overline{P},t)\}=\underset{t_1}{\overset{t_2}{\int }}dt\{\sum _{k=1}^f\delta P_k{\dot{Q}}_k(t)+P_k\delta {\dot{Q}}_k(t)-\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial Q_k}\delta Q_k-\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial P_k}\delta P_k\}$

$\begin{array}{rl}& \delta \{M_1(q(t_2),Q(t_2),t_2)-M_1(q(t_1),Q(t_1),t_1)\}=\sum _k({\frac{\partial M_1}{\partial q_k}\delta q_k|}_{t_1}^{t_2}+{\frac{\partial M_1}{\partial Q_k}\delta Q_k|}_{t_1}^{t_2})\\ & mit{\frac{\partial M_1}{\partial q_k}\delta q_k|}_{t_1}^{t_2}=0und{\frac{\partial M_1}{\partial Q_k}\delta Q_k|}_{t_1}^{t_2}\ne 0\\ \end{array}$

Außerdem:

$\underset{t_1}{\overset{t_2}{\int }}dt{P}_k\delta {\dot{Q}}_k(t)={P_k\delta Q_k|}_{t_1}^{t_2}-\underset{t_1}{\overset{t_2}{\int }}dt{\dot{P}}_k\delta Q_k$

Nebenbemerkung: Für die Variation gilt bekanntlich: $\delta \overline{q}(t_1)=\delta \overline{q}(t_2)=0$ . Jedoch sind p(t1) und p(t2) beliebig. Dadurch können sich nun insbesondere die Randbedingungen für $Q(\overline{q},\overline{p},t)$ ändern.

Unter Beachtung der obigen relationen gilt nun:

$0=\delta \underset{t_1}{\overset{t_2}{\int }}dtL=\sum _{k=1}^f{(P_k+\frac{\partial M_1}{\partial Q_k})\delta Q_k|}_{t_1}^{t_2}+\underset{t_1}{\overset{t_2}{\int }}dt\sum _{k=1}^f\{({\dot{Q}}_k(t)-\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial P_k})\delta P_k-({\dot{P}}_k(t)+\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial Q_k})\delta Q_k\}$

Aus den obigen Relationen ist bekannt:

$(P_k+\frac{\partial M_1}{\partial Q_k})=0$

Gleichzeitig sind jedoch Pi und Qi unabhängig und können demnach unabhängig variiert werden. Das bedeutet, dass $\delta P_{k}und\delta Q_k$ unabhängig sind.

Somit muss jeweils für sich gelten:

$\begin{array}{rl}& 0=({\dot{Q}}_k(t)-\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial P_k})\\ & 0=({\dot{P}}_k(t)+\frac{\partial \overline{H}(\overline{Q},\overline{P},t)}{\partial Q_k})\\ \end{array}$

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:

$\sum _{k=1}^f(p_k{\dot{q}}_k-P_k{\dot{Q}}_k)-(H-\overline{H})=\frac{d}{dt}{M}_1$

Außerdem gilt:

$\frac{d}{dt}{M}_1=\frac{d}{dt}(M_2(\overline{q}(t),\overline{P}(t),t)-\sum _k{P}_k{Q}_k)=\sum _{k=1}^f(\frac{\partial M_2}{\partial q_k}{\dot{q}}_k+\frac{\partial M_2}{\partial P_k}{\dot{P}}_k-{\dot{P}}_k{Q}_k-P_k{\dot{Q}}_k)+\frac{\partial M_2}{\partial t}$

So dass folgt:

$\sum _{k=1}^f(p_k-\frac{\partial M_2}{\partial q_k}){\dot{q}}_k+(Q_k-\frac{\partial M_2}{\partial P_k}){\dot{P}}_k+(P_k-P_k){\dot{Q}}_k=(H-\overline{H})+\frac{\partial M_2}{\partial t}$

Da dies für beliebige ${\dot{q}}_k,{\dot{P}}_k$ gilt, kann die Summe nur allgemein identisch sein, wenn gilt:

$\begin{array}{rl}& (p_k-\frac{\partial M_2}{\partial q_k})=0\Rightarrow p_k=\frac{\partial M_2}{\partial q_k}\\ & Q_k=\frac{\partial M_2}{\partial P_k}\\ & \overline{H}=H+\frac{\partial M_2}{\partial t}\\ \end{array}$

Analog kann gezeigt werden, dass für

$M_3(\overline{p},\overline{Q},t)=M_1(\overline{q},\overline{Q},t)-\sum _{k=1}^f\frac{\partial M_1}{\partial q_k}{q}_k$

Hier folgt ( Übungsaufgabe):

$\begin{array}{rl}& (q_k+\frac{\partial M_3}{\partial p_k})=0\Rightarrow q_k=-\frac{\partial M_3}{\partial p_k}\\ & P_k=-\frac{\partial M_2}{\partial Q_k}\\ \end{array}$

oder

$M_4(\overline{p},\overline{P},t)=M_1(\overline{q},\overline{Q},t)-\sum _{k=1}^f(\frac{\partial M_1}{\partial q_k}{q}_k+\frac{\partial M_1}{\partial Q_k}{Q}_k)$

$\begin{array}{rl}& \Rightarrow (q_k+\frac{\partial M_4}{\partial p_k})=0\Rightarrow q_k=-\frac{\partial M_4}{\partial p_k}\\ & Q_k=-\frac{\partial M_4}{\partial P_k}\\ \end{array}$

Beispiele für kanonische Transformationen

Erzeugende sei:

$\begin{array}{rl}& M_1(\overline{q},\overline{Q},t)=\sum _{j=1}^f{q}_j{Q}_j\\ & \Rightarrow p_j=\frac{\partial M_1}{\partial q_j}=-Q_j\\ & P_j=-\frac{\partial M_1}{\partial Q_j}=q_j\\ & (\overline{q},\overline{p})\to (\overline{P},-\overline{Q})\\ \end{array}$

Bei dieser Trafo werden also Ort und Impuls vertauscht.

Beispiel 2:

$\begin{array}{rl}& M_2(\overline{q},\overline{P},t)=\sum _{j=1}^f{q}_j{P}_j\\ & \Rightarrow p_j=\frac{\partial M_2}{\partial q_j}=P_j\\ & Q_j=\frac{\partial M_2}{\partial P_j}=q_j\\ & (\overline{q},\overline{p})\to (\overline{Q},\overline{P})\\ \end{array}$

Dies ist also die identische Transformation

Beispiel: Harmonischer Oszillator:

$\begin{array}{rl}& H=\frac{p^2}{2m}+\frac{m{\omega }^2}{2}{q}^2\\ & M_1(q,Q)=\frac{m\omega }{2}{q}^2\cot Q\\ & \Rightarrow p=\frac{\partial M_1}{\partial q}=m\omega q\cot Q\\ & P=-\frac{\partial M_1}{\partial Q}=\frac{m\omega }{2}\frac{q^2}{{\sin }^{2}Q}\\ & q={\left(\frac{2}{m\omega }P\right)}^{\frac{1}{2}}\sin Q\\ & p={\left(2m\omega P\right)}^{\frac{1}{2}}\cos Q\\ & \\ \end{array}$

$\begin{array}{rl}& H=\overline{H}\left(\frac{\partial M_1}{\partial t}\right)=0\\ & H=\frac{2m\omega P{\cos }^{2}Q}{2m}+\frac{m{\omega }^{2}2P}{2m\omega }{\sin }^{2}Q=\omega P\\ \end{array}$

Die Variable Q ist also zyklisch.

$\begin{array}{rl}& \dot{P}=-\frac{\partial H}{\partial Q}=0\Rightarrow P=\alpha =const\\ & \dot{Q}=\frac{\partial H}{\partial P}=\omega \Rightarrow Q=\omega t+\beta \\ \end{array}$

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

$q(t)={\left(\frac{2\alpha }{m\omega }\right)}^{\frac{1}{2}}\sin (\omega t+\beta )$

Dabei beschreibt $\alpha$ die Amplitude und $\beta$ 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

$\overline{x}:=\left(\begin{array}{c}q\\ p\\ \end{array}\right)$ ist Vektor im Phasenraum

${H_{,}}_x:=\left(\begin{array}{c}\frac{\partial H}{\partial q}\\ \frac{\partial H}{\partial p}\\ \end{array}\right)$ ist Ableitungsvektor

$J:=\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)$ ist Metrik im Phasenraum ( metrischer Tensor)

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

$\dot{\overline{x}}:=J{H}_{,x}\iff -J\dot{\overline{x}}=H_{,x}\iff \dot{q}=\frac{\partial H}{\partial p},\dot{p}=-\frac{\partial H}{\partial q}$

Leicht läßt sich zeigen:

$\begin{array}{rl}& J^2=-1\\ & J^{-1}=J^T=-J\\ \end{array}$

Verallgemeinerung auf mehr Freiheitsgrade

$\begin{array}{rl}& \overline{x}:=\left(\begin{array}{c}{q}_1\\ ...\\ q_f\\ p_1\\ ...\\ p_f\\ \end{array}\right)\\ & {\overline{H}}_x:=\left(\begin{array}{c}\frac{\partial H}{\partial q_1}\\ ...\\ \frac{\partial H}{\partial q_f}\\ \frac{\partial H}{\partial p_1}\\ ...\\ \frac{\partial H}{\partial p_f}\\ \end{array}\right)J:=\left(\begin{array}{cc}0& 1_f\\ -1_f& 0\\ \end{array}\right)\\ \end{array}$

Die kanonischen Gleichungen lauten

$\dot{\overline{x}}:=J{H}_x\iff -J\dot{\overline{x}}=H_x$

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

$\dot{\overline{x}}:=A\overline{x}=J{H}_x$

Diese Gleichung ist abzuleiten aus der Hamiltonfunktion:

$H=\frac{1}{2}(a{q}^2+2bqp+c{p}^2)z.\mathbf{B}.a={{\omega }_0}^2,b=0,c=1$

$\Rightarrow \dot{\overline{x}}:=\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)\left(\begin{array}{c}\frac{\partial H}{\partial q}\\ \frac{\partial H}{\partial p}\\ \end{array}\right)=\begin{array}{c}bq+cp\\ -aq-bp\\ \end{array}$

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

$\begin{array}{rl}& A=\left(\begin{array}{cc}b& c\\ -a& -b\\ \end{array}\right)\\ & tr\left(A\right)=0\\ \end{array}$

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:

$\begin{array}{rl}& M_1(\overline{q},\overline{Q},t):\\ & \Rightarrow p_j=\frac{\partial M_1}{\partial q_j}\\ & P_j=-\frac{\partial M_1}{\partial Q_j}\\ & \Rightarrow \frac{\partial p_j}{\partial Q_k}=\frac{{\partial }^2{M}_1}{\partial Q_k\partial q_j}=-\frac{\partial P_k}{\partial q_j}\\ \end{array}$

$\begin{array}{rl}& M_2(\overline{q},\overline{P},t)=M_1(\overline{q},\overline{Q},t)-\sum _{j=1}^f\frac{\partial M_1}{\partial Q_j}{Q}_j\\ & \Rightarrow p_j=\frac{\partial M_2}{\partial q_j}\\ & Q_j=\frac{\partial M_2}{\partial P_j}\\ & \Rightarrow \frac{\partial p_j}{\partial P_k}=\frac{{\partial }^2{M}_2}{\partial P_k\partial q_j}=\frac{\partial Q_k}{\partial q_j}\\ \end{array}$

$\begin{array}{rl}& M_3(\overline{p},\overline{Q},t)=M_1(\overline{q},\overline{Q},t)-\sum _{j=1}^f\frac{\partial M_1}{\partial q_j}{q}_j\\ & \Rightarrow q_j=\frac{\partial M_3}{\partial p_j}\\ & P_j=-\frac{\partial M_3}{\partial Q_j}\\ & \Rightarrow \frac{\partial q_j}{\partial Q_k}=-\frac{{\partial }^2{M}_3}{\partial Q_k\partial p_j}=\frac{\partial P_k}{\partial p_j}\\ \end{array}$

$\begin{array}{rl}& M_4(\overline{p},\overline{P},t)=M_1(\overline{q},\overline{Q},t)-\sum _{j=1}^f(\frac{\partial M_1}{\partial Q_j}{Q}_j+\frac{\partial M_1}{\partial q_j}{q}_j)\\ & \Rightarrow q_j=-\frac{\partial M_4}{\partial p_j}\\ & Q_j=\frac{\partial M_1}{\partial P_j}=q_j\\ & \Rightarrow \frac{\partial q_j}{\partial P_k}=\frac{{\partial }^2{M}_1}{\partial P_k\partial p_j}=-\frac{\partial Q_k}{\partial p_j}\\ \end{array}$

Dabei sind:

$\overline{x}:=\left(\begin{array}{c}{q}_1\\ ...\\ q_f\\ p_1\\ ...\\ p_f\\ \end{array}\right)\overline{y}:=\left(\begin{array}{c}{Q}_1\\ ...\\ Q_f\\ P_1\\ ...\\ P_f\\ \end{array}\right)$

$\begin{array}{rl}& M_{\alpha \beta }=\frac{\partial x_{\alpha }}{\partial y_{\beta }}\\ & {\left(M^{-1}\right)}_{\alpha \beta }:=\frac{\partial y_{\alpha }}{\partial x_{\beta }}\alpha ,\beta =1,...,2f\\ \end{array}$

Beweis:

$\sum _{\gamma =1}^{2f}{M}_{\alpha \gamma }{\left(M^{-1}\right)}_{\gamma \beta }=\sum _{\gamma =1}^{2f}\frac{\partial x_{\alpha }}{\partial y_{\gamma }}\frac{\partial y_{\gamma }}{\partial x_{\beta }}=\frac{\partial x_{\alpha }}{\partial x_{\beta }}={\delta }_{\alpha \beta }$

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

$M_{\alpha \beta }=\sum _{\mu ,\nu =1}^{2f}{J}_{\alpha \mu }{J}_{\beta \nu }{\left(M^{-1}\right)}_{\mu \nu }$

Beweis:

In Matrixform lautet diese Gleichung:

$M=J{\left(J{M}^{-1}\right)}^T$

Die linke Seite (M) lautet:

$M=\left(\begin{array}{cc}\frac{\partial q}{\partial Q}& \frac{\partial q}{\partial P}\\ \frac{\partial p}{\partial Q}& \frac{\partial p}{\partial P}\\ \end{array}\right)$

Die rechte Seite lautet:

$\begin{array}{rl}& J{\left(J{M}^{-1}\right)}^T=\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right){\left[\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)\left(\begin{array}{cc}\frac{\partial Q}{\partial q}& \frac{\partial Q}{\partial p}\\ \frac{\partial P}{\partial q}& \frac{\partial P}{\partial p}\\ \end{array}\right)\right]}^T=\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right){\left(\begin{array}{cc}\frac{\partial P}{\partial q}& \frac{\partial P}{\partial p}\\ -\frac{\partial Q}{\partial q}& -\frac{\partial Q}{\partial p}\\ \end{array}\right)}^T\\ & =\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right){\left(\begin{array}{cc}{\left(\frac{\partial P}{\partial q}\right)}^T& -{\left(\frac{\partial Q}{\partial q}\right)}^T\\ {\left(\frac{\partial P}{\partial p}\right)}^T& -{\left(\frac{\partial Q}{\partial p}\right)}^T\\ \end{array}\right)}^{}=\left(\begin{array}{cc}{\left(\frac{\partial P}{\partial p}\right)}^T& -{\left(\frac{\partial Q}{\partial p}\right)}^T\\ -{\left(\frac{\partial P}{\partial q}\right)}^T& {\left(\frac{\partial Q}{\partial q}\right)}^T\\ \end{array}\right)\\ \end{array}$

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

$\begin{array}{rl}& M=J{\left(J{M}^{-1}\right)}^T\\ & \Rightarrow JM=-{\left(J{M}^{-1}\right)}^T=-{\left(M^{-1}\right)}^T{J}^T\\ & \Rightarrow M^{T}JM=-M^T{\left(M^{-1}\right)}^T{J}^T=-{\left(M^{-1}M\right)}^T{J}^T=-J^T=J\\ & M^{T}JM=J\\ \end{array}$

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:

$(\overline{x},\overline{y}):={\overline{x}}^{T}J\overline{y}=\sum _{i,k=1}^{2f}{x}_i{J}_{ik}{y}_k$

es handelt sich dabei um eine schiefsymmetrische, nichtentartete Bilinearform

Eigenschaften:

1. Schiefsymmetrie:

$(\overline{x},\overline{y})=-(\overline{y},\overline{x})$ , Beweis: $(\overline{x},\overline{y})={\overline{x}}^{T}J\overline{y}={\left({\overline{y}}^T{J}^T\overline{x}\right)}^T=-{\overline{y}}^T{J}^{}\overline{x}=-(\overline{y},\overline{x})$

1. bilinear:

$(\overline{x},{\lambda }_1{\overline{y}}_1+{\lambda }_2{\overline{y}}_2)={\lambda }_1(\overline{x},{\overline{y}}_1)+{\lambda }_2(\overline{x},{\overline{y}}_2)$

1. nichtentartet:

$(\overline{x},\overline{y})=0\mathrm{\forall }\overline{y}\Rightarrow \overline{x}=0$

Nebenbemerkung: Es gilt: $(\overline{x},\overline{x})=0\mathrm{\forall }\overline{x}$

Also Selbstorthogonalität

Beweis: ${\overline{x}}^{T}J\overline{x}=\left(\begin{array}{cc}q& p\\ \end{array}\right)\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)\left(\begin{array}{c}q\\ p\\ \end{array}\right)=qp-pq=0$

Die Symplektische Struktur auf dem $R^{2f}$ ist von einer euklidischen Metrik grundsätzlich zu unterscheiden:

${(\overline{x},\overline{y})}_{Eu}=\sum _i{x}_i{y}_i={\overline{x}}^{T}g\overline{y}$

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

Im Euklidischen gelten jedoch die Relationen:

$\begin{array}{rl}& (\overline{x},\overline{y})=(\overline{y},\overline{x})\\ & (\overline{x},\overline{x})\ge 0\\ \end{array}$

Definition:

Die Menge der Matrizen M ( kanonische Trafo) mit

$M^{T}JM=J$ bildet die reelle symplektische Gruppe S über $R^{2f}$ .

Dies ist die Symmetriegruppe der symplektischen Struktur.

Gruppeneigenschaften

1. $M_1,M_2\in S\Rightarrow M_3=M_1{M}_2\in S$

Beweis: ${M_3}^{T}J{M}_3={\left(M_1{M}_2\right)}^{T}J\left(M_1{M}_2\right)={M_2}^T{M_1}^{T}J{M}_1{M}_2={M_2}^{T}J{M}_2=J$

2. Assoziativität ( matrixmultiplikation !)

3. Einselement Einheitsmatrix !

1. Inverse:

$M^{-1}:=J^{-1}{M}^{T}J$

Beweis: $M^{-1}M=\left(J^{-1}{M}^{T}J\right)M=J^{-1}\left(M^{T}JM\right)=J^{-1}J=1$

Dabei gilt : $M^T,J\in S$ Beweis: Übungsaufgabe

1. Weiter gilt:

$\det M=1$ Beweis: Übungsaufgabe oder Scheck, S. 102

Fazit:

Die Invarianz der kanonischen Gleichungen $\dot{\overline{x}}:=A\overline{x}=J{\overline{H}}_{,x}$ kann durch di symplektische Struktur des Phasenraums beschrieben werden:

$\begin{array}{rl}& {\dot{y}}_i=\sum _k^{}\frac{\partial y_i}{\partial x_k}{\dot{x}}_k\iff \dot{\overline{y}}=M^{-1}\dot{\overline{x}}=\left(J^{-1}{M}^{T}J\right)J{\overline{H}}_{,x}\\ & \frac{\partial \overline{H}}{\partial y_i}=\sum _k^{}\frac{\partial \overline{H}}{\partial x_k}\frac{\partial x_k}{\partial y_i}\iff {\overline{H}}_{,y}=M^T{\overline{H}}_{,x}\\ & \Rightarrow \dot{\overline{y}}=\left(J^{-1}{M}^{T}J\right)J{\left(M^T\right)}^{-1}{\overline{H}}_{,y}=-J(-1)M^T{\left(M^T\right)}^{-1}{\overline{H}}_{,y}=J{\overline{H}}_{,y}\\ \end{array}$

Der Satz von Liouville

Lösung der Differenzialgleichung $\dot{\overline{x}}:=A\overline{x}=J{\overline{H}}_{,x}$

$\overline{x}(t,t_0,{\overline{x}}_0)=:{\overline{\mathrm{\Phi }}}_{t,t_0}({\overline{x}}_0)$ Definition: Fluß im Phasenraum

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

Der Fluß beschreibt dabei die Zeitentwicklung der Anfangskonfiguration:

$\begin{array}{rl}& {\overline{\mathrm{\Phi }}}_{t,t_0}({\overline{x}}_0):\mathrm{\Gamma }\to \mathrm{\Gamma }\\ & {\overline{x}}_0(t_0)\to \overline{x}(t)\\ \end{array}$

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

Beispiel: eindimensionaler harmonischer Oszi:

$\begin{array}{rl}& (q,p)\in \mathrm{\Gamma }\\ & \dot{q}=p\\ & \dot{p}=-{{\omega }_0}^{2}q\\ & \dot{\overline{x}}=A\overline{x}\\ & A=\left(\begin{array}{cc}0& 1\\ -{{\omega }_0}^2& 0\\ \end{array}\right)\\ \end{array}$

Die Lösung lautet:

$\overline{x}(t)={\overline{\mathrm{\Phi }}}_{t-t_0}({\overline{x}}_0)=\exp \left[(t-t_0)A\right]{\overline{x}}_0$

Dies ist also gerade das Exponenzial der Matrix.

Aufschluss liefert eine Reihenentwicklung:

$\begin{array}{rl}& \overline{x}(t)=\sum _n\frac{{\left[(t-t_0)A\right]}^n}{n!}{\overline{x}}_0=[1\cos {\omega }_0(t-t_0)+\frac{A}{{\omega }_0}\sin {\omega }_0(t-t_0)]{\overline{x}}_0\\ & =\left(\begin{array}{cc}\cos {\omega }_0(t-t_0)& \frac{1}{{\omega }_0}\sin {\omega }_0(t-t_0)\\ -{\omega }_0\sin {\omega }_0(t-t_0)& \cos {\omega }_0(t-t_0)\\ \end{array}\right){\overline{x}}_0\\ \end{array}$

Beweis:

$\begin{array}{rl}& A^2=\left(\begin{array}{cc}0& 1\\ -{{\omega }_0}^2& 0\\ \end{array}\right)\left(\begin{array}{cc}0& 1\\ -{{\omega }_0}^2& 0\\ \end{array}\right)=-{{\omega }_0}^{2}1\\ & A^{2n}={(-1)}^{2n}{{\omega }_0}^{2n}1\\ & A^{2n+1}={(-1)}^n{{\omega }_0}^{2n+1}\frac{1}{{\omega }_0}A\\ \end{array}$

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:

$V_{to}=\underset{U_{to}}{\overset{}{\int }}{d}^{2f}{x}_0$ Bei t: $V_t=\underset{U_t}{\overset{}{\int }}{d}^{2f}{x}_0=\underset{U_{t_0}}{\overset{}{\int }}{d}^{2f}{x}_0\det \left(\frac{\partial x}{\partial x_0}\right)=\underset{U_{t_0}}{\overset{}{\int }}{d}^{2f}{x}_0\det (D{\mathrm{\Phi }}_{t,t_0}({\overline{x}}_0))$

Mit der Jacobi- Matrix:

${(D{\mathrm{\Phi }}_{t,t_0}({\overline{x}}_0))}_{ik}:=\frac{\partial {{\mathrm{\Phi }}^i}_{t,t_0}({\overline{x}}_0)}{\partial {x_0}^k}=\frac{\partial x^i}{\partial {x_0}^k}$

Dies kann für Zeiten nahe t0 reihenentwickelt werden:

$\begin{array}{rl}& {\mathrm{\Phi }}_{t,t_0}({\overline{x}}_0)={\overline{x}}_0+\overline{F}({\overline{x}}_0,t)(t-t_0)+O({(t-t_0)}^2)\\ & \overline{F}({\overline{x}}_0,t)=J{\overline{H}}_{,x}=\left(\begin{array}{c}\frac{\partial H}{\partial p}\\ -\frac{\partial H}{\partial q}\\ \end{array}\right)\\ \end{array}$

Somit folgt:

$\frac{\partial {{\mathrm{\Phi }}^i}_{t,t_0}({\overline{x}}_0)}{\partial {x_0}^k}={\delta }_{ik}+\frac{\partial {\overline{F}}^i({\overline{x}}_0,t)}{\partial {x_0}^k}(t-t_0)+O({(t-t_0)}^2)$

Mit Hilfe $\det (1+B\epsilon )=1+\epsilon tr(B)+O({\epsilon }^2)$ folgt:

$\begin{array}{rl}& \det \left(D{\mathrm{\Phi }}_{t,t_0}\right)=\left|\frac{\partial {{\mathrm{\Phi }}^i}_{t,t_0}({\overline{x}}_0)}{\partial {x_0}^k}\right|=1+(t-t_0)\sum _{i=1}^{2f}\frac{\partial {\overline{F}}^i({\overline{x}}_0,t)}{\partial {x_0}^i}(t-t_0)+O({(t-t_0)}^2)\\ & \sum _{i=1}^{2f}\frac{\partial {\overline{F}}^i({\overline{x}}_0,t)}{\partial {x_0}^i}=div\overline{F}=\frac{\partial }{\partial q}\frac{\partial H}{\partial p}-\frac{\partial }{\partial p}\frac{\partial H}{\partial q}=0\\ \end{array}$

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

$\det \left(D{\mathrm{\Phi }}_{t,t_0}\right)=\left|\frac{\partial {{\mathrm{\Phi }}^i}_{t,t_0}({\overline{x}}_0)}{\partial {x_0}^k}\right|=1+O({(t-t_0)}^2)\cong 1$

$\Rightarrow V_t=\underset{U_{t_0}}{\overset{}{\int }}{d}^{2f}{x}_0\det (D{\mathrm{\Phi }}_{t,t_0}({\overline{x}}_0))=\underset{U_{t_0}}{\overset{}{\int }}{d}^{2f}{x}_0(1+O{(t-t_0)}^2)\cong V_{t0}$

Nebenbemerkung:

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

Für den Fluß ${\mathrm{\Phi }}_{t,t_0}$ zu $\dot{\overline{x}}:=J{\overline{H}}_{,x}$ ist $D{\mathrm{\Phi }}_{t,t_0}$ eine symplektische Matrix, das heißt $\det \left(D{\mathrm{\Phi }}_{t,t_0}\right)=1$ .

Das bedeutet, das Volumenelement $d{x}^1...d{x}^{2f}$ im Phasenraum ist unter dem Fluß invariant:

$d{x}^1...d{x}^{2f}=Det(D\mathrm{\Phi })d{x_0}^1....d{x_0}^{2f}=d{x_0}^1....d{x_0}^{2f}$

Beispiel: eindimensionaler harmonischer Oszi

$\begin{array}{rl}& x^i(t)=\sum _{k=1}^2{P}_{ik}{x_0}^k\\ & P=\left(\begin{array}{cc}\cos {\omega }_0(t-t_0)& \frac{1}{{\omega }_0}\sin {\omega }_0(t-t_0)\\ -{\omega }_0\sin {\omega }_0(t-t_0)& \cos {\omega }_0(t-t_0)\\ \end{array}\right)\\ & \det P={\cos }^2{\omega }_0(t-t_0)+{\sin }^2{\omega }_0(t-t_0)=1\\ & d{x}^{1}d{x}^2=(\det P)d{x_0}^{1}d{x_0}^2=d{x_0}^{1}d{x_0}^2\\ \end{array}$

Poisson- Klammern

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

$Observable=g(\overline{q},\overline{p},t)$

Die zeitliche Änderung längs der Bahn $\overline{q}(t),\overline{p}(t)$ im Phasenraum $\mathrm{\Gamma }$

$\frac{dg(\overline{q},\overline{p},t)}{dt}=\sum _{i=1}^f(\frac{\partial g}{\partial q_i}{\dot{q}}_i+\frac{\partial g}{\partial p_i}{\dot{p}}_i)+\frac{\partial g}{\partial t}=\sum _{i=1}^f(\frac{\partial g}{\partial q_i}\frac{\partial H}{\partial p_i}-\frac{\partial g}{\partial p_i}\frac{\partial H}{\partial q_i})+\frac{\partial g}{\partial t}=:\{g,H\}+\frac{\partial g}{\partial t}$

Definition:

Für zwei beliebige Observablen $g(\overline{q},\overline{p},t)$ und $f(\overline{q},\overline{p},t)$ heißt

$\sum _{i=1}^f(\frac{\partial g}{\partial q_i}\frac{\partial f}{\partial p_i}-\frac{\partial g}{\partial p_i}\frac{\partial f}{\partial q_i})=:\{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\}=({\overline{f}}_x,{\overline{g}}_x)={{\overline{f}}_x}^{T}J{\overline{g}}_x=\sum _{i,k=1}^f\left(\frac{\partial f}{\partial x_i}{J}_{ik}\frac{\partial g}{\partial x_k}\right)=\left(\begin{array}{cc}\frac{\partial f}{\partial q}& \frac{\partial f}{\partial p}\\ \end{array}\right)\left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)\left(\begin{array}{c}\frac{\partial g}{\partial q}\\ \frac{\partial g}{\partial q}\\ \end{array}\right)=\left(\begin{array}{cc}\frac{\partial f}{\partial q}& \frac{\partial f}{\partial p}\\ \end{array}\right)\left(\begin{array}{c}\frac{\partial g}{\partial p}\\ -\frac{\partial g}{\partial q}\\ \end{array}\right)$ 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,{\lambda }_1{g}_1+{\lambda }_2{g}_2\}={\lambda }_1\{f,g_1\}+{\lambda }_2\{f,g_2\}$

3.nichtentartet: $(f,g)=0\mathrm{\forall }g\Rightarrow f=const.$ (Nullelement, wegen ${\overline{f}}_{,x}=0$ )

Nebenbemerkung: Es gilt: $\{f,f\}=0\mathrm{\forall }f$

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: $\frac{\partial g}{\partial q_k}=\{g,p_k\}\frac{\partial g}{\partial p_k}=-\{g,q_k\}$

Die Poissonklammer ist invariant unter kanonischen Transformationen:

Beweis: Trafo: x->y

Die Jacobi- Determinante $M_{\alpha \beta }=\frac{\partial x_{\alpha }}{\partial y_{\beta }}$ ist symplektische Matrix,

das heißt , es gilt: $M^{T}JM=J\iff M^{-1}J{\left(M^{-1}\right)}^T=J$ , da ja $M^{-1}J{\left(M^{-1}\right)}^T=J^{-1}{M}^{T}JJ{\left(M^{-1}\right)}^T=-J^{-1}{M}^T{\left(M^{-1}\right)}^T=-J^{-1}=J$

Nun muss man umrechnen von :

$\begin{array}{rl}& \frac{\partial f}{\partial x_i}=\sum _k\frac{\partial f}{\partial y_k}\frac{\partial y_k}{\partial x_i}=\sum _k{M_{ki}}^{-1}\frac{\partial f}{\partial y_k}\iff {\overline{f}}_x={\left(M^{-1}\right)}^T{\overline{f}}_y\iff {{\overline{f}}_x}^T={{\overline{f}}_y}^T\left(M^{-1}\right)\\ & {{\overline{f}}_x}^{T}J{\overline{g}}_x={{\overline{f}}_y}^T\left(M^{-1}\right)J{\left(M^{-1}\right)}^T{\overline{g}}_y={{\overline{f}}_y}^{T}J{\overline{g}}_y\\ \end{array}$

Also:

$({\overline{f}}_x,{\overline{g}}_x)=({\overline{f}}_y,{\overline{g}}_y)$

Für nicht explizit zeitabhängige Observable $g(\overline{q},\overline{p})$ gilt:

$\frac{dg}{dt}=\{g,H\}$

g ist genau dann Bewegungskonstante, wenn gilt:

$\{g,H\}=0$

Speziallfall: g ist Koordinate der Impuls: $g=q_k,g=p_k$

$\begin{array}{rl}& {\dot{q}}_k=\{q_k,H\}\\ & {\dot{p}}_k=\{p_k,H\}\\ \end{array}$

So folgen die Hamiltonschen Gleichungen

Kompakt kann geschrieben werden:

$\begin{array}{rl}& {\dot{x}}_k=\{x_k,H\}=\sum _{i,j=1}^f\frac{\partial x_k}{\partial x_i}{J}_{ij}\frac{\partial H}{\partial x_j}=\sum _{j=1}^f{J}_{kj}\frac{\partial H}{\partial x_j}\\ & also:\dot{\overline{x}}=J{\overline{H}}_{,x}\\ \end{array}$

Fundamentale Poisson- Klammern:

$\begin{array}{rl}& \{q_k,q_j\}=0\\ & \{p_k,p_j\}=0\\ & \{q_k,p_j\}={\delta }_{kj}\\ \end{array}$

Kompakt:

$\{x_k,x_j\}=\sum _{l,m}\frac{\partial x_k}{\partial x_l}{J}_{lm}\frac{\partial x_j}{\partial x_m}=J_{kj}\cong \left(\begin{array}{cc}0& 1\\ -1& 0\\ \end{array}\right)$

Die Poissonklammer ist invariant unter kanonischen Transformationen, da

$M^{T}JM=J$

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

Somit:

Satz: Die Transformation $(\overline{q},\overline{p})\to (\overline{Q},\overline{P})$ ist genau dann kanonisch, wenn :

$\begin{array}{rl}& \{Q_k,Q_j\}=0\\ & \{P_k,P_j\}=0\\ & \{Q_k,P_j\}={\delta }_{kj}\\ \end{array}$

Beweis: Zur Vereinfachung: Nicht explizit zeitabhängige Trafos: $\frac{\partial M}{\partial t}=0\iff \overline{H}=H$

Bewegungsgleichung:

${\dot{y}}_k=\{y_k,H\}=\sum _{i,j=1}^f\frac{\partial y_k}{\partial x_i}{J}_{ij}\frac{\partial H}{\partial x_j}=\sum _{j=1}^f{J}_{kj}\frac{\partial H}{\partial x_j}$

Wegen $({\overline{f}}_x,{\overline{g}}_x)=({\overline{f}}_y,{\overline{g}}_y)$

kann nun die Bewegungsgleichung in den alten Koordinaten gebildet werden:

${\dot{y}}_k=\{y_k,H\}=\sum _{i,j,l=1}^f\frac{\partial y_k}{\partial x_i}{J}_{ij}\frac{\partial \overline{H}}{\partial y_l}\frac{\partial y_l}{\partial x_j}=\sum _{l=1}^f\frac{\partial \overline{H}}{\partial y_l}\sum _{i,j=1}^f\frac{\partial y_k}{\partial x_i}{J}_{ij}\frac{\partial y_l}{\partial x_j}=\sum _{l=1}^f\frac{\partial \overline{H}}{\partial y_l}\{y_k,y_l\}$

Also folgt:

$\begin{array}{rl}& {\dot{y}}_k=\{y_k,H\}=\sum _{i,j,l=1}^f\frac{\partial y_k}{\partial x_i}{J}_{ij}\frac{\partial \overline{H}}{\partial y_l}\frac{\partial y_l}{\partial x_j}=\sum _{l=1}^f\frac{\partial \overline{H}}{\partial y_l}\sum _{i,j=1}^f\frac{\partial y_k}{\partial x_i}{J}_{ij}\frac{\partial y_l}{\partial x_j}=\sum _{l=1}^f\frac{\partial \overline{H}}{\partial y_l}\{y_k,y_l\}\\ & {\dot{y}}_k=\sum _{l=1}^f{J}_{kl}\frac{\partial \overline{H}}{\partial y_l}\iff J_{kl}=\{y_k,y_l\}\\ \end{array}$

Mit der Bedeutung

${\dot{y}}_k=\sum _{l=1}^f{J}_{kl}\frac{\partial \overline{H}}{\partial y_l}$ Hamiltonsche Bewegungsgleichung in den neuen Koordinaten -> Trafo kanonisch

$J_{kl}=\{y_k,y_l\}$

 fundamentale Poissonklammern in den neuen Koordinaten

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

Folgende Aussagen sind äquivalent:

$\overline{x}=\left(\begin{array}{c}\overline{q}\\ \overline{p}\\ \end{array}\right)\to \overline{y}=\left(\begin{array}{c}\overline{Q}\\ \overline{P}\\ \end{array}\right)$ ist kanonisch

$\iff$ die kanonischen Gleichungen $\dot{\overline{x}}=J{\overline{H}}_{,x}$ sind invariant

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

$\iff$ die fundamentalen Poissonklammern $J_{kl}=\{x_k,x_l\}$ sind ivariant

$\iff$ die Jacobi- Matrix $M_{\alpha \beta }=\frac{\partial x_{\alpha }}{\partial x_{\beta }}$ ist symplektisch, das heißt $M^{T}JM=J$

$\iff$ es existiert eine Erzeugende !

Bezug zur Quantenmechanik

Ein Übergang zur Quantenmechanik ist möglich:

Von der klassischen Variablen $g(\overline{q},\overline{p},t)$ zum qm. Operator: $g:H->H$ mit dem Hilbertraum H

Von der Poissonklammer: $\{f,g\}\to \frac{1}{i\hslash }[f,g]$ zum Kommutator

Aus den fundamentalen Poisson- Klammern folgen die kanonischen Vertauschiungsrelationen:

$\begin{array}{rl}& \{q_k,q_j\}=0\to [q_k,q_j]=0\\ & \{p_k,p_j\}=0\to [p_k,p_j]=0\\ & \{q_k,p_j\}={\delta }_{kj}\to [q_k,p_j]=i\hslash {\delta }_{kj}\\ \end{array}$

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

Die Bewegungsgleichungen:

$\begin{array}{rl}& \frac{dg(\overline{q},\overline{p},t)}{dt}=\sum _{i=1}^f(\frac{\partial g}{\partial q_i}{\dot{q}}_i+\frac{\partial g}{\partial p_i}{\dot{p}}_i)+\frac{\partial g}{\partial t}=\sum _{i=1}^f(\frac{\partial g}{\partial q_i}\frac{\partial H}{\partial p_i}-\frac{\partial g}{\partial p_i}\frac{\partial H}{\partial q_i})+\frac{\partial g}{\partial t}=:\{g,H\}+\frac{\partial g}{\partial t}\\ & \to \frac{dg}{dt}=\frac{1}{i\hslash }[g,H]+\frac{\partial g}{\partial t}\\ \end{array}$

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.