Hamiltonsche Mechanik

Die Hamiltonsche Mechanik ist eine Formulierung der klassischen Mechanik, die von dem irischen Mathematiker und Physiker William Rowan Hamilton entwickelt wurde. Sie beschreibt die Bewegung physikalischer Systeme – etwa von Planeten, Pendeln oder Teilchen – mithilfe von Größen wie Ort und Impuls.

Im Unterschied zur Newtonschen Mechanik, die unmittelbar mit Kräften und Beschleunigungen arbeitet, betrachtet die Hamiltonsche Mechanik den Zustand eines Systems als Punkt in einem abstrakten Raum aller möglichen Orte und Impulse, dem sogenannten Phasenraum. Die Zeitentwicklung dieses Punktes wird durch die Hamilton-Funktion mittels der Hamiltonschen Bewegungsgleichungen bei Vernachlässigung von Reibung bestimmt. Hamilton-Funktion hängt meist eng mit der Gesamtenergie des Systems zusammen.

Die Bewegungsgleichungen wurden 1834 von William Rowan Hamilton angegeben. Die Hamiltonsche Mechanik setzt voraus, dass das betrachtete System durch eine Lagrange-Funktion beschrieben werden kann. Sie ermöglicht somit keine Behandlung zusätzlicher Problemklassen, verfügt jedoch über eine besonders elegante mathematische Struktur.<ref>Klaus Lichtenegger: Schlüsselkonzepte zur Physik : Von den Newton-Axiomen bis zur Hawking-Strahlung. 1. Auflage. Springer Spektrum, Berlin 2015, ISBN 978-3-8274-2384-9, S. 34. </ref> Diese hat zwei entscheidende Vorteile:

Zum einen besagt der Satz von Liouville, dass die Bewegung im Phasenraum volumentreu ist. Daraus folgt, dass es bei der Bewegung im Phasenraum keine Wirbel und Staupunkte gibt, vergleichbar dem Fluss einer inkompressiblen Flüssigkeit.
Zum anderen besitzen die Hamiltonschen Bewegungsgleichungen eine große Gruppe von Transformationen, die kanonischen Transformationen, die es gestatten, sie in andere, manchmal lösbare Hamiltonsche Gleichungen zu transformieren.

Mit den Hamiltonschen Bewegungsgleichungen untersucht man insbesondere integrable und chaotische Bewegungen. Sie werden in der statistischen Physik verwendet. Sie spielen auch eine zentrale Rolle in vielen Bereichen der modernen Physik, darunter die Himmelsmechanik, die Quantenmechanik und die Relativitätstheorie.

Grundidee: Hamilton-Funktion und Bewegungsgleichungen[Bearbeiten]

Die Hamiltonsche Mechanik beschreibt den Zustand eines mechanischen Systems durch die verallgemeinerten Koordinaten <math>q_i</math> und die zugehörigen kanonischen Impulse <math>p_i</math>. Sie geht von einer Lagrange-Funktion

aus, aus der die kanonischen Impulse durch

<math>p_i=\frac{\partial L}{\partial \dot q_i}</math>

definiert werden.<ref>Leonard Susskind, George E. Hrabovsky: Klassische Mechanik - Das Theoretische Minimum : Alles, was Sie brauchen, um Physik zu treiben. 1. Auflage. Springer, Berlin 2019, ISBN 978-3-662-60333-8, S. 121. </ref> Lassen sich diese Beziehungen nach den Geschwindigkeiten <math>\dot q_i</math> auflösen, so kann man die Geschwindigkeiten durch die Impulse ersetzen und die Dynamik in den Variablen <math>q_i</math> und <math>p_i</math> formulieren.

Die zentrale Größe der Hamiltonschen Mechanik ist die Hamilton-Funktion <math>H</math>. Sie wird durch eine Legendre-Transformation der Lagrange-Funktion definiert:<ref>L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Band 1, Mechanik -. 9. Auflage. Akademie Verlag, Berlin 1979, S. 161. </ref>

wobei die Geschwindigkeiten als Funktionen von <math>q_i</math>, <math>p_i</math> und <math>t</math> aufzufassen sind. Die Hamilton-Funktion entspricht bei vielen mechanischen Systemen der Gesamtenergie und besitzt daher häufig eine unmittelbare physikalische Bedeutung.

Die zeitliche Entwicklung des Systems wird durch die Hamiltonschen Bewegungsgleichungen beschrieben:<ref>L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Band 1, Mechanik -. 9. Auflage. Akademie Verlag, Berlin 1979, S. 162. </ref>

<math>\dot q_i=\frac{\partial H}{\partial p_i},

\qquad \dot p_i=-\frac{\partial H}{\partial q_i}.</math>

Im Gegensatz zu den Lagrange-Gleichungen, die gewöhnlich als Differentialgleichungen zweiter Ordnung formuliert sind, erhält man damit ein System gekoppelter Differentialgleichungen erster Ordnung. Die Kenntnis der Anfangswerte aller Koordinaten und Impulse genügt, um die weitere Bewegung des Systems eindeutig festzulegen.

Teilchen im Potential[Bearbeiten]

Bei einem Teilchen der Masse <math>m</math>, das sich nichtrelativistisch in einem Potential <math>V</math> bewegt, setzt sich die Lagrange-Funktion aus kinetischer und potentieller Energie zusammen:<ref>Leonard Susskind, George E. Hrabovsky: Klassische Mechanik - Das Theoretische Minimum : Alles, was Sie brauchen, um Physik zu treiben. 1. Auflage. Springer, Berlin 2019, ISBN 978-3-662-60333-8, S. 90. </ref>

<math>L(\mathbf q,\dot{\mathbf q})=\tfrac12 m\dot{\mathbf q}^2-V(\mathbf q)</math>

Mit dem kanonischen Impuls

<math>\mathbf p=m\dot{\mathbf q}=\quad\Rightarrow\quad \quad\dot{\mathbf q}(\mathbf p)=\frac{\mathbf p}{m}</math>

lautet dann die Hamilton-Funktion:

<math>\mathcal H(\mathbf q,\mathbf p)=\mathbf p\cdot\mathbf q(\mathbf p)-L(\mathbf q,\dot{\mathbf q}(\mathbf p))=\frac{\mathbf p^2}{2\,m}+V(\mathbf q)</math>

Die zugehörigen Hamiltonschen Bewegungsgleichungen

<math>\dot q_k=\frac{\partial \mathcal H}{\partial p_k} =\frac{p_k}{m},\quad \dot{p}_k=-\frac{\partial \mathcal H}{\partial q_k} = - \frac{\partial V}{\partial q_k}</math>

sind Newtons Gleichungen für die Bewegung in einem konservativen Kraftfeld,

<math>m\,\ddot q_k = F_k = - \frac{\partial V}{\partial q_k}\,.</math>

Beispielsweise ist die potentielle Energie eines eindimensionalen <math>(n=1)</math> harmonischen Oszillators nach dem Hookeschen Gesetz <math>V(q)=\tfrac 1 2 \,m\,\omega^2 \,q^2</math>. Die Lagrange-Funktion ist

<math>L(q,\dot{q})=\tfrac12 m\dot{q}^2-\tfrac 1 2 \,m\,\omega^2 \,q^2</math>

und die Hamilton-Funktion mit <math>\dot{q}=p/m</math> lautet:<ref name="Nolting 108">Wolfgang Nolting: Grundkurs Theoretische Physik 2 - Analytische Mechanik. 8. Auflage. Springer Verlag, Heidelberg 2011, ISBN 978-3-642-12949-0, S. 108. </ref>

<math>\mathcal H(q,p)=p\, q(p)-L( q,\dot{q}( p))=\frac{p^2}{2\,m}+\tfrac 1 2 \,m\,\omega^2 \,q^2</math>

Eine zeitunabhängige Hamilton-Funktion ist eine Erhaltungsgröße <math>\mathcal H(q,p)=E</math> und die Bewegung des Oszillators verläuft entlang von Ellipsen im Phasenraum

<math>E=\frac{p^2}{2\,m}+\tfrac 1 2 \,m\,\omega^2 \,q^2\quad\text{bzw.}\quad 1=\frac{p^2}{2\,m\,E}+\frac{q^2}{2\,E/\,m\,\omega^2}</math>

mit den Halbachsen <math>\sqrt{2\,m\,E}</math> und <math>\sqrt{2\,E/\,m}\omega</math>.<ref name="Nolting 108"/>

Das Zeitverhalten liefern die Hamiltonschen Bewegungsgleichungen

<math>\dot q=\frac{\partial \mathcal H}{\partial p} =\frac{p}{m},\quad \dot{p}=-\frac{\partial \mathcal H}{\partial q} = - \frac{\partial V}{\partial q}=- m\,\omega^2\, q,</math>

denn sie entsprechen der Newtonschen Gleichung für die Bewegung in einem konservativen Kraftfeld:

<math>m\,\ddot q = - m\,\omega^2\, q</math>

und man erhält, dass die Bahn um die Ruhelage <math>q=0</math> schwingt:

<math>q(t) = A\, \cos \bigl(\omega\,(t-t_0)\bigr)\,.</math>

Dabei ist <math>A</math> die Amplitude und <math>t_0</math> eine Zeit, zu der diese maximale Auslenkung durchlaufen wird.

Freies relativistisches Teilchen[Bearbeiten]

Für ein relativistisches, freies Teilchen mit der Energie-Impuls-Beziehung <math>E^2-\mathbf p^2\,c^2=m^2\,c^4</math> ist die Hamilton-Funktion<ref name="LL" details="S. 32.">L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Band 2, Klassische Feldtheorie -. 8. Auflage. Akademie Verlag, Berlin 1981. </ref>

<math>\mathcal H(\mathbf q,\mathbf p)=\sqrt{m^2\,c^4+ \mathbf p^2\,c^2}.</math>

Die Hamiltonschen Bewegungsgleichungen besagen, wie die Geschwindigkeit mit dem Impuls zusammenhängt und dass sich der Impuls nicht mit der Zeit ändert:

<math>\dot q_k = \frac{\partial \mathcal H }{\partial p_k} = \frac{p_k\,c^2}{\sqrt{m^2\,c^4 + \mathbf p^2\,c^2}} \,,\quad k=1, 2, 3 </math>

<math>\dot p_k =-\frac{\partial \mathcal H}{\partial q_k} =0 </math>

Wenn die Hamilton-Funktion wie in diesen Beispielen nicht von der Zeit abhängt, behält das System von Teilchen seine anfängliche Energie, sie ist dann eine Erhaltungsgröße.

Relativistische Hamilton-Funktion eines geladenen Teilchens im elektromagnetischen Feld[Bearbeiten]

Für ein relativistisches Teilchen der Ruhemasse <math>m_0</math> und der Ladung <math>e</math> im elektromagnetischen Feld mit Vektorpotential <math>\vec A</math> und Skalarpotential <math>\phi</math> lautet die Lagrange-Funktion:<ref>L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Band 2, Klassische Feldtheorie -. 8. Auflage. Akademie Verlag, Berlin 1981, S. 55. </ref>

Der zugehörige verallgemeinerte Impuls ergibt sich aus der Ableitung der Lagrange-Funktion nach dem Geschwindigkeitsvektor:

<math>\vec{\cal P}=\frac{\partial{ L}}{\partial\vec v}=\frac{m_0\vec v}{\sqrt{1-v^2/c^2}}+e\vec{A}</math>

Die Hamilton-Funktion erhält man durch die Legendre-Transformation

<math>{\cal H}=\vec v\cdot\frac{\partial{ L}}{\partial\vec v}-{ L}=\vec v\cdot\left( \frac{m_0\vec v}{\sqrt{1-v^2/c^2}}+e\vec{A} \right)-\left(-m_0c^2\sqrt{1-\frac{v^2}{c^2}}+e\vec{A}\cdot\vec v-e\phi\right)</math>

mit dem Ergebnis:

Um die Hamilton-Funktion ausschließlich durch die kanonischen Variablen auszudrücken, ist die Geschwindigkeit durch den verallgemeinerten Impuls zu eliminieren.

<math>(\vec{\cal P}-e\vec{A})^2=\frac{m_0^2v^2}{1-v^2/c^2}\quad\Rightarrow\quad \sqrt{1-v^2/c^2}=\frac{m_0c}{\sqrt{m_0^2c^2+(\vec{\cal P}-e\vec{A})^2}}</math>

Einsetzen in die Hamilton-Funktion ergibt schließlich<ref>Eckhard Rebhan: Theoretische Physik: Quantenmechanik. 1. Auflage. Spektrum Akademischer Verlag, Berlin 2012, ISBN 978-3-8274-2314-6, S. 133. </ref>

Geometrische Formulierung[Bearbeiten]

Ein Hamiltonsches mechanisches System <math>(M,\omega,\mathcal{H})</math> besteht

aus einer <math>2n</math>-dimensionalen glatten Mannigfaltigkeit <math>M</math> genannt Phasenraum (zum Beispiel <math>\mathbb{R}^{2n}</math>) und einer symplektischen Struktur <math>\omega</math>, was <math>(M,\omega)</math> zu einer symplektischen Mannigfaltigkeit macht.
einer glatten Funktion <math>\mathcal{H}\colon M\to \mathbb{R}</math> genannt Hamilton-Funktion.<ref>Wladimir Igorewitsch Arnold: Mathematical Methods of Classical Mechanics. Hrsg.: Springer. Ägypten 1989, S. 161. </ref>

Die Nichtausgeartetheit der symplektischen Form gewährleistet, dass zu jeder Hamilton-Funktion eindeutig ein Vektorfeld <math>X_{\mathcal H}</math> existiert, das durch das innere Produkt des Vektors <math>X_{\mathcal H}</math> mit der Differentialform <math>\omega</math> <ref name="AbrahamMarsden187">Ralph Abraham, Jerrold E. Marsden: Foundations of Mechanics. 2. Auflage. The Benjamin / Cummings Publishing Company, Menlo Park / Reading 1978, ISBN 0-8053-0102-X, S. 187. </ref>

<math>\iota_{X_{\mathcal H}}\omega=\mathrm{d}\mathcal H\quad\text{bzw.}\quad \omega(X_{\mathcal H},Y)=\mathrm{d}\mathcal H\cdot Y\quad \forall Y</math>

definiert ist. Dieses Vektorfeld <math>X_{\mathcal H}</math> wird als Hamiltonsches Vektorfeld bezeichnet. Die zeitliche Entwicklung des Systems wird durch die Integralkurven von <math>X_{\mathcal{H}}</math> beschrieben. Ist

<math>\gamma(t)\in M</math>

eine solche Integralkurve, so gilt

<math>\dot{\gamma}(t)=X_{\mathcal{H}}(\gamma(t)).</math>

Damit wird die Dynamik vollständig durch die Geometrie des Phasenraums und die Hamilton-Funktion <math>\mathcal H(t,q,p)</math> festgelegt. Die Hamilton-Funktion eines Systems von Teilchen hängt von den (verallgemeinerten) Ortskoordinaten <math>q=(q_1, q_2, \dots, q_n)</math> und von den (verallgemeinerten) Impulskoordinaten <math>p=(p_1, p_2, \dots, p_n)</math> der Teilchen ab und kann auch von der Zeit <math>t</math> abhängen.

Die Zahl <math>n</math> der Koordinaten und Impulse nennt man die Zahl der Freiheitsgrade. Der Phasenraum ist <math>2n</math>-dimensional. Für die symplektischen Struktur

<math>\omega=\sum_{i=1}^n \mathrm{d}q_i\wedge \mathrm{d}p_i</math>

auf den kanonischen Koordinaten <math>(q_1, q_2, \dots, q_n,p_1, p_2, \dots, p_n)</math> nimmt das Vektorfeld <math>X_{\mathcal H}</math> die Form

<math>X_{\mathcal H}=\left(\frac{\partial \mathcal H}{\partial p_k}, -\frac{\partial \mathcal H}{\partial q_k}\right)=J\cdot\text{d}\mathcal H</math>

an, mit <math>J=\big(\begin{smallmatrix}

 0 & I\\
 -I & 0

\end{smallmatrix}\big)</math> und <math>I</math> als <math>n</math>-dimensionale Einheitsmatrix.<ref name="AbrahamMarsden187"/>

Die Hamilton-Funktion bestimmt die zeitliche Entwicklung der Teilchenorte und Teilchenimpulse durch die Hamiltonschen Bewegungsgleichungen:

<math>

\dot q_k =\frac{\mathrm d}{\mathrm dt}q_k =\frac{\partial \mathcal H}{\partial p_k}\,, \quad \dot p_k =\frac{\mathrm d}{\mathrm dt}p_k =-\frac{\partial \mathcal H}{\partial q_k} \,,\quad k=1, 2, \dots, n </math>

Dies ist ein System gewöhnlicher Differentialgleichungen erster Ordnung für die <math>2n</math> unbekannten Funktionen der Zeit, <math>q(t),p(t)\,.</math>

Ist die Lagrange-Funktion zeitunabhängig, so ist die Hamilton-Funktion eine Erhaltungsgröße. Dann schneiden sich die Lösungskurven nicht und es geht durch jeden Punkt des Phasenraums eine Lösungskurve.

Die Hamilton-Funktion ist die Gesamtenergie als Funktion des Phasenraums, wenn<ref>Friedhelm Kuypers: Klassische Mechanik - mit über 300 Beispielen und Aufgaben mit Lösungen. 9. Auflage. WILEY-VCH, Weinheim 2010, ISBN 978-3-527-40989-1, S. 98. </ref>

die Zwangsbedingungen skleronom, holonom sind
das System konservativ ist
das Bezugssystem ruht.

Bei zeitabhängigen <math>\mathcal H(t,q,p)</math> kann man die Zeit als einen zusätzlichen Freiheitsgrad <math>t=q_0</math> mit zugehörigem Impuls <math>p_0</math> und der zeitunabhängigen Hamilton-Funktion <math>\hat{\mathcal H}(q_0,q,p_0,p) = \mathcal H(q_0,q,p)+p_0</math> auffassen. Daher beschränken wir uns im Folgenden auf zeitunabhängige Hamilton-Funktionen. Allerdings ist die Funktion <math>\mathcal H(q_0,q,p)+p_0</math> nicht nach unten beschränkt und die Hyperfläche konstanter Energie <math>\hat{\mathcal H}=E </math> ist nicht, wie bei einigen Überlegungen vorausgesetzt, kompakt.

Harmonischer Oszillator in der geometrischen Formulierung[Bearbeiten]

Der Phasenraum des eindimensionalen harmonischen Oszillators ist die zweidimensionale Mannigfaltigkeit

<math>M=\mathbb{R}^2,</math>

deren Punkte durch die Koordinaten <math>(q,p)</math> beschrieben werden. Die moderne symplektisch-geometrische Formulierung der Hamiltonschen Mechanik beruht auf der Gleichung

<math>\iota_{X_{\cal H}}\omega=\mathrm{d}{\cal H},</math>

wobei <math>\omega</math> die symplektische Form, <math>{\cal H}</math> die Hamilton-Funktion und <math>X_{\cal H}</math> das zugehörige Hamiltonsche Vektorfeld bezeichnet. Wir verwenden für den zweidimensionalen Phasenraum <math>M</math> die kanonischen Koordinaten <math>q</math> und <math>p</math>. Als symplektische Form wählt man

<math>\omega = \mathrm{d}q \wedge \mathrm{d}p.</math>

Die Hamilton-Funktion des harmonischen Oszillators lautet

<math>{\cal H}(q,p)=\frac{p^2}{2m}+\frac{m\omega_0^2}{2}q^2.</math>

Das zugehörige Hamiltonsche Vektorfeld <math>X_{\cal H}</math> ist durch die Gleichung

<math>\iota_{X_{\cal H}}\omega=\mathrm{d}{\cal H}</math>

definiert. Aus der Hamilton-Funktion ergibt sich zunächst

<math>\mathrm{d}{\cal H}

= \frac{p}{m}\,\mathrm{d}p+m\omega_0^2 q\,\mathrm{d}q.</math>

Setzt man allgemein

<math>X_{\cal H}

= a(q,p)\frac{\partial}{\partial q} + b(q,p)\frac{\partial}{\partial p},</math>

so folgt mit der Definition des inneren und äußeren Produkts

<math>\iota_{X_{\cal H}}\omega(\cdot)=\omega(X_{\cal H},\cdot)

= (\mathrm{d}q\wedge \mathrm{d}p)(X_{\cal H},\cdot) = \mathrm{d}q(X_{\cal H})\,\mathrm{d}p(\cdot)-\mathrm{d}p(X_{\cal H})\,\mathrm{d}q(\cdot) = a\,\mathrm{d}p(\cdot)-b\,\mathrm{d}q(\cdot).</math>

Durch Vergleich mit <math>\mathrm{d}{\cal H}</math> erhält man

<math>a=\frac{p}{m},

\qquad b=-m\omega_0^2 q.</math>

Damit lautet das Hamiltonsche Vektorfeld für den harmonischen Oszillator

<math>X_{\cal H}

= \frac{p}{m}\frac{\partial}{\partial q} - m\omega_0^2 q\,\frac{\partial}{\partial p}.</math>

Eine Integralkurve des Vektorfeldes <math>X_{\cal H}</math> ist eine Kurve

<math>\gamma(t)=\bigl(q(t),p(t)\bigr)</math>

im Phasenraum, deren Tangentialvektor in jedem Punkt mit dem Wert des Vektorfeldes an diesem Punkt übereinstimmt. Mathematisch wird dies durch

<math>\dot{\gamma}(t)=X_{\cal H}\bigl(\gamma(t)\bigr)</math>

ausgedrückt. Der Tangentialvektor der Kurve <math>\gamma(t)</math> ist

<math>\dot{\gamma}(t)

= \dot q(t)\,\frac{\partial}{\partial q} + \dot p(t)\,\frac{\partial}{\partial p}.</math>

Andererseits ergibt sich für das Hamiltonsche Vektorfeld entlang der Kurve

<math>X_{\cal H}\bigl(\gamma(t)\bigr)

= \frac{p(t)}{m}\,\frac{\partial}{\partial q} - m\omega_0^2 q(t)\,\frac{\partial}{\partial p}.</math>

Da beide Ausdrücke denselben Vektor darstellen, müssen die Koeffizienten der Basisvektoren <math>\partial/\partial q</math> und <math>\partial/\partial p</math> übereinstimmen. Daraus folgt unmittelbar

sowie

<math>\dot p(t)=-m\omega_0^2 q(t).</math>

Dies sind genau die Hamiltonschen Bewegungsgleichungen des harmonischen Oszillators.

Wirkungsprinzip[Bearbeiten]

Die Hamiltonschen Bewegungsgleichungen folgen aus dem Hamiltonschen Prinzip der stationären Wirkung. Von allen denkbaren Bahnen im Phasenraum,

<math>\Gamma\colon t \mapsto \bigl(q(t),p(t)\bigr)\,,</math>

die anfänglich zur Zeit <math>\underline{t}</math> durch den Anfangspunkt

<math>\bigl(\underline q, \underline p\bigr)=\bigl(q(\underline t),p(\underline t)\bigr)</math>

und schließlich zur Zeit <math>\overline{t}</math> durch den Endpunkt

<math>\bigl(\overline q, \overline p\bigr)=\bigl(q(\overline t),p(\overline t)\bigr)</math>

laufen, ist die physikalisch durchlaufene Bahn diejenige, auf der die Wirkung

<math>W[\Gamma]=\int\limits_{\underline t}^{\overline t}\mathrm d t\,\left(\sum_{i=1}^n p_i(t)\,

\frac{\mathrm d q_i(t)}{\mathrm d t} - \mathcal H(q(t),p(t)) \right)</math> stationär ist.

Betrachtet man nämlich eine einparametrige Schar von Kurven

<math>\Gamma_\alpha\colon t \mapsto \bigl(q(t,\alpha),p(t,\alpha)\bigr)\,,</math>

die anfänglich zur Zeit <math>\underline{t}</math> durch den Anfangspunkt

<math>\bigl(\underline q, \underline p\bigr)=\bigl(q(\underline t,\alpha),p(\underline t,\alpha)\bigr)</math>

und schließlich zur Zeit <math>\overline{t}</math> durch den Endpunkt

<math>\bigl(\overline q, \overline p\bigr)=\bigl(q(\overline t,\alpha),p(\overline t,\alpha)\bigr)</math>

laufen, so ist die Wirkung <math>W[\Gamma_\alpha]</math> für <math>\alpha=0</math> extremal, falls dort die Ableitung nach <math>\alpha</math> verschwindet.

Wir bezeichnen diese Ableitung als Variation der Wirkung

<math>\delta W = \frac{\partial W[\Gamma_\alpha]}{\partial \alpha}_{|_{\alpha=0}}\,.</math>

Ebenso ist

<math>\delta q_i = \frac{\partial q_i(t,\alpha)}{\partial \alpha}_{|_{\alpha=0}}</math>

die Variation des Ortes und

<math>\delta p_i = \frac{\partial p_i(t,\alpha)}{\partial \alpha}_{|_{\alpha=0}}</math>

die Variation des Impulses.

Die Variation der Wirkung ist nach der Kettenregel

<math>

\delta W=\sum_{i=1}^n \int\limits_{\underline t}^{\overline t}\mathrm d t \left( \delta p_i(t)\,\frac{\mathrm d q_i(t)}{\mathrm d t}+ p_i(t)\,\frac{\mathrm d \delta q_i(t)}{\mathrm d t} - \delta q_i(t)\,\frac{\partial \mathcal H}{\partial q_i}_{|_{(q(t),p(t))}} - \delta p_i(t)\,\frac{\partial \mathcal H}{\partial p_i}_{|_{(q(t),p(t))}}\right)\,. </math>

Den zweiten Term schreiben wir als vollständige Zeitableitung und einen Term, bei dem <math>\delta q_i</math> ohne Zeitableitung auftritt:

<math>

p_i(t)\,\frac{\mathrm d \delta q_i(t)}{\mathrm d t}= \frac{\mathrm d \bigl(p_i(t)\, \delta q_i(t)\bigr)}{\mathrm d t} - \frac{\mathrm d p_i(t)}{\mathrm d t}\, \delta q_i(t) </math> Das Integral über die vollständige Ableitung ergibt <math>\bigl(p_i(t)\, \delta q_i(t)\bigr)</math> zur Anfangs- und Endzeit und verschwindet, weil dann <math>\delta q_i</math> verschwindet, denn es gehen alle Kurven der Schar durch dieselben Anfangs- und Endpunkte. Fassen wir schließlich die Terme mit <math>\delta q_i</math> und <math>\delta p_i</math> zusammen, so beträgt die Variation der Wirkung

<math>

\delta W=\sum_{i=1}^n \int\limits_{\underline t}^{\overline t}\mathrm d t \left( -\delta q_i(t)\left(\frac{\mathrm d p_i(t)}{\mathrm d t}+ \frac{\partial \mathcal H}{\partial q_i}_{|_{(q(t),p(t))}}\right) + \delta p_i(t)\left(\frac{\mathrm d q_i(t)}{\mathrm d t}- \frac{\partial \mathcal H}{\partial p_i}_{|_{(q(t),p(t))}}\right)\right)\,. </math> Damit die Wirkung stationär ist, muss dieses Integral für alle <math>\delta q_i</math> und alle <math>\delta p_i</math> verschwinden, die anfänglich und schließlich verschwinden. Das ist genau dann der Fall, wenn die Faktoren verschwinden, mit denen sie im Integral auftreten:

<math>

0 = \frac{\mathrm d p_i(t)}{\mathrm d t}+ \frac{\partial \mathcal H}{\partial q_i}_{|_{(q(t),p(t))}} </math>

<math>

0= \frac{\mathrm d q_i(t)}{\mathrm d t}- \frac{\partial \mathcal H}{\partial p_i}_{|_{(q(t),p(t))}} </math> Die Wirkung ist also stationär, wenn die Hamiltonschen Bewegungsgleichungen gelten.

Zusammenhang zur Lagrange-Funktion[Bearbeiten]

Die Hamilton-Funktion ist die bezüglich der Geschwindigkeiten <math>\dot q=(\dot q_1,\dot q_2\dots \dot q_n)</math> Legendre-Transformierte der Lagrange-Funktion <math>\mathcal L(q,\dot q):</math>

<math>\mathcal H(q,p)=

\sum_{k=1}^n p_k\, \dot q_k(q,p) - \mathcal L(q,\dot q(q,p))</math> Dabei sind auf der rechten Seite mit den Geschwindigkeiten <math>\dot q</math> diejenigen Funktionen <math>\dot q(q,p)</math> gemeint, die man erhält, wenn man die Definition der Impulse

<math>p_k = \frac{\partial \mathcal L}{\partial \dot q_k}</math>

nach den Geschwindigkeiten auflöst.

Wenn man die Definition der Impulse invertieren und nach den Geschwindigkeiten auflösen kann, dann gelten die Hamiltonschen Bewegungsgleichungen genau dann, wenn die Euler-Lagrange-Gleichungen der Wirkung

<math>W[\Gamma]=\int_{\underline t}^{\overline t}\mathrm d t\,

\mathcal L(q(t),\dot q(t))</math> erfüllt sind. Denn die partielle Ableitung von <math>\mathcal H</math> nach den Impulsen ergibt nach der Kettenregel und der Definition der Impulse

<math>

\frac{\partial \mathcal H}{\partial p_i} = \dot q_i + \sum_j p_j \frac{\partial \dot q_j}{\partial p_i} - \sum_j \frac{\partial \dot q_j}{\partial p_i} \underbrace{\frac{\partial \mathcal L}{\partial \dot q_j}}_{p_j} = \dot q_i\,. </math> Ebenso ergibt die Ableitung nach den Ortskoordinaten

<math>

\frac{\partial \mathcal H}{\partial q_i} = \sum_j p_j \frac{\partial \dot q_j}{\partial q_i} - \frac{\partial \mathcal L}{\partial q_i} - \sum_j \frac{\partial \dot q_j}{\partial q_i} \underbrace{\frac{\partial \mathcal L}{\partial \dot q_j}}_{p_j} = -\frac{\partial \mathcal L}{\partial q_i}\,. </math> Die Euler-Lagrange-Gleichung besagt

<math>

\frac{\partial \mathcal L}{\partial q_i} = \frac{\mathrm d}{\mathrm d t}\frac{\partial \mathcal L}{\partial \dot q_i} = \dot{p}_i\,. </math> Also gelten die Hamiltonschen Bewegungsgleichungen, wenn die Euler-Lagrange-Gleichung gilt. Umgekehrt gilt die Euler-Lagrange-Gleichung, wenn die Hamiltonschen Bewegungsgleichungen gelten.

Beispielsweise hängt beim freien, relativistischen Teilchen mit der Lagrange-Funktion<ref name="LL" details="S. 30." />

<math>\mathcal L= - m\,c^2 \sqrt{1-\dot{\mathbf q}^2/c^2}</math>

der Impuls gemäß

<math>\mathbf p = \frac{\partial \mathcal L}{\partial \dot{\mathbf q}} =\frac{m \dot{\mathbf q}}{\sqrt{1-\dot{\mathbf q}^2/c^2}}</math>

von der Geschwindigkeit ab. Umgekehrt ist die Geschwindigkeit daher die Funktion

<math>\dot{\mathbf q}=\frac{\mathbf p\,c^ 2}{\sqrt{m^2\,c^4+\mathbf p^2\,c^2}}</math>

des Impulses. In die obige Gleichung für <math>\mathcal H</math> eingesetzt

<math>\begin{align}

\mathcal H(\mathbf q,\mathbf p) & = \mathbf p\cdot \dot{\mathbf q}(\mathbf q,\mathbf p) - {\mathcal L}(\mathbf q,\dot{\mathbf q}(\mathbf q,\mathbf p)) = \frac{{\mathbf p}^2 \,c^2}{\sqrt{m^2\,c^4+\mathbf p^2\,c^2}} + m\,c^2 \sqrt{1-\frac{{\mathbf p}^2\,c^2}{m^2\,c^4+\mathbf p^2\,c^2}}= \\
& = \frac{{\mathbf p}^2 \,c^2}{\sqrt{m^2\,c^4+\mathbf p^2\,c^2}} + \frac{m^2 \,c^4}{\sqrt{m^2\,c^4+\mathbf p^2\,c^2}}=\sqrt{m^2\,c^4+ \mathbf p^2\,c^2}

\end{align}</math>

ergibt sich die schon angegebene Hamilton-Funktion des freien, relativistischen Teilchens.

Hängt die Lagrange-Funktion nicht explizit von der Zeit ab, dann besagt das Noether-Theorem, dass die Energie

<math>E(q,\dot q)=\sum_k \dot q_k \frac{\partial \mathcal L}{\partial \dot q_k} - \mathcal L</math>

auf den physikalischen Bahnen ihren anfänglichen Wert behält. Der Vergleich mit der Legendre-Transformation zeigt, dass es sich bei der Hamilton-Funktion um diese Energie handelt, bei der die Geschwindigkeiten als Funktion der Impulse aufzufassen sind:

<math>\mathcal H(q,p) = E(q,\dot q(q,p))</math>

Poisson-Klammer[Bearbeiten]

Der Wert einer Phasenraumfunktion <math>\Phi(t,q,p)</math> ändert sich auf Bahnen <math>(q(t),p(t))</math> mit der Zeit dadurch, dass er explizit von <math>t</math> abhängt und dadurch, dass sich der Bahnpunkt ändert:

<math>\frac{\mathrm d \Phi(t,q(t),p(t))}{\mathrm d t} =

\frac{\partial \Phi}{\partial t} + \sum_i \left( \frac{\partial \Phi}{\partial q_i}\frac{\mathrm d q_i}{\mathrm d t}+ \frac{\partial \Phi}{\partial p_i}\frac{\mathrm d p_i}{\mathrm d t}\right) \, .</math>

Die physikalisch durchlaufenen Bahnen genügen den Hamiltonschen Bewegungsgleichungen:<ref>L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Band 1, Mechanik -. 9. Auflage. Akademie Verlag, Berlin 1979, S. 166. </ref>

<math>\frac{\mathrm d \Phi(t,q(t),p(t))}{\mathrm d t} =

\frac{\partial \Phi}{\partial t} + \sum_i \left( \frac{\partial \Phi}{\partial q_i}\frac{\partial \mathcal H}{\partial p_i} - \frac{\partial \Phi}{\partial p_i}\frac{\partial \mathcal H}{\partial q_i}\right) \, .</math>

Mit der von Siméon Denis Poisson eingeführten Poisson-Klammer zweier Phasenraumfunktionen <math>\Phi</math> und <math>\Psi</math>

<math>\bigl\{\Phi, \Psi\bigr \}=

\sum_i \left( \frac{\partial \Phi}{\partial q_i}\frac{\partial \Psi}{\partial p_i} - \frac{\partial \Phi}{\partial p_i}\frac{\partial \Psi}{\partial q_i}\right)</math>

gilt also

<math>\frac{\mathrm d \Phi(t,q(t),p(t))}{\mathrm d t} =

\frac{\partial \Phi}{\partial t} + \bigl\{\Phi, \mathcal H \bigr \} \, .</math>

Mit Poisson-Klammern geschrieben gleicht das Formelbild der Hamiltonschen Bewegungsgleichungen den Heisenbergschen Bewegungsgleichungen der Quantenmechanik.<ref>Friedhelm Kuypers: Klassische Mechanik - mit über 300 Beispielen und Aufgaben mit Lösungen. 9. Auflage. WILEY-VCH, Weinheim 2010, ISBN 978-3-527-40989-1, S. 322. </ref>

Als Koordinatenfunktionen aufgefasst haben die Phasenraumkoordinaten die Poisson-Klammern

<math>\{q_i,q_j\}=0=\{p_i,p_j\}\,,\,\{q_i,p_j\}=\delta_{ij} \, .</math>

Ihnen entsprechen in der Quantenmechanik nach kanonischer Quantisierung die kanonischen Vertauschungsrelationen.

Die Poisson-Klammer ist antisymmetrisch, linear und genügt der Produktregel und der Jacobi-Identität. Für alle Zahlen <math>a</math> und <math>b</math> und alle Phasenraumfunktionen <math>\Psi, \Phi, \Lambda</math> gilt

<math>\{\Psi,\Phi\}= -\{\Phi,\Psi\} \, ,</math>
<math>\{\Psi,a\,\Phi+b\,\Lambda\}= a\,\{\Psi,\Phi\}+b\,\{\Psi,\Lambda\} \, ,</math>
<math>\{\Psi,\Phi\,\,\Lambda\}=\{\Psi,\Phi\}\,\Lambda + \Phi\,\{\Psi,\Lambda\} \, ,</math>
<math>\{\Psi,\{\Phi,\Lambda\}\} + \{\Phi,\{\Lambda,\Psi\}\} + \{\Lambda,\{\Psi,\Phi\}\} = 0 \, .</math>

Die differenzierbaren Phasenraumfunktionen bilden eine Lie-Algebra mit der Poisson-Klammer als Lie-Produkt.<ref>Norbert Straumann: Klassische Mechanik. 1. Auflage. Springer, Berlin Heidelberg 1987, ISBN 3-540-18527-5, S. 204. </ref>

Hamiltonscher Fluss[Bearbeiten]

Zu jeder (zeitunabhängigen) Phasenraumfunktion <math>\Phi</math> gehört das Vektorfeld <math>v_{\Phi}\,,</math>

das Phasenraumfunktionen <math>\Psi</math> längs der Kurven ableitet, die die Hamiltonschen Gleichungen mit <math>\mathcal H = \Phi</math> lösen.

Die Abbildung <math>\Phi_t</math> der Anfangswerte der Lösungskurven <math>(q(0),p(0))</math> auf <math>(q(t),p(t))</math> ist der zu <math>\Phi</math> gehörige Hamiltonsche Fluss.

Symplektische Struktur[Bearbeiten]

Der Phasenraum mit seiner Poisson-Klammer ist eine symplektische Mannigfaltigkeit mit der symplektischen Form

<math>\omega=\sum_i \mathrm d q_i\wedge\,\mathrm d p_i\,.</math>

Angewendet auf die zu <math>\Phi</math> und <math>\Psi</math> gehörigen Vektorfelder ergibt diese Zweiform die Poisson-Klammer der beiden Funktionen:

<math>\omega\,(v_\Phi,v_\Psi)=\{\Phi,\Psi\}</math>

Die symplektische Form ist invariant unter jedem Hamiltonschen Fluss. Dies besagt Folgendes: Ist anfänglich eine zweidimensionale Fläche <math>F</math> im Phasenraum gegeben, dann wird sie mit der Zeit durch den Hamiltonschen Fluss einer Phasenraumfunktion <math>\Phi</math> auf die Fläche <math>\Phi_t(F)</math> abgebildet. Die mit der symplektischen Form gemessene Größe der Anfangsfläche stimmt mit der Größe zu jeder späteren Zeit überein. Hamiltonscher Fluss ist flächentreu:

<math>\int_F\omega = \int_{\Phi_t(F)}\omega</math>

Da das Flächenelement <math>\omega</math> invariant ist, ist auch das Volumenelement <math>\omega^n=n!\,\mathrm d^nq\,\mathrm d^n p</math> invariant unter Hamiltonschem Fluss. Dieser Befund ist Liouvilles Theorem. Das Volumen eines Bereichs <math>B</math> des Phasenraumes ändert sich nicht bei Hamiltonscher Zeitentwicklung:

<math>\int_B\omega^n = \int_{\Phi_t(B)}\omega^n</math>

Insbesondere bleibt der Bereich, innerhalb dessen sich das System anfänglich wegen der Messfehler befindet, gleich groß. Daraus kann man allerdings nicht schließen, dass sich anfängliche Unkenntnis nicht vergrößert. Bei chaotischer Bewegung können Anfangswerte, die sich zunächst nur durch kleine Messfehler unterschieden, auf einen großen Bereich mit vielen kleinen Löchern wie Schlagsahne verteilt werden. Auch Schlagen von Sahne vergrößert ihr mikroskopisch ermitteltes Volumen nicht.

Kanonische Transformation[Bearbeiten]

Die Hamilton-Gleichungen vereinfachen sich, falls die Hamilton-Funktion von einer Variablen, beispielsweise <math>q_1\,,</math> nicht abhängt. Dann liegt eine Symmetrie vor: die Hamilton-Funktion ist invariant unter der Verschiebung von <math>q_1\,.</math> Umgekehrt können bei Vorliegen einer Symmetrie (in einer Umgebung eines Punktes, der kein Fixpunkt ist) die Orts- und Impulsvariablen so gewählt werden, dass die Hamilton-Funktion von einer Variablen <math>q_1</math> nicht abhängt. Dann ist einfach <math>p_1(t)=p_1(0)\,.</math>

Integrable Bewegung[Bearbeiten]

Die Bewegungsgleichungen sind integrabel, wenn die Hamilton-Funktion nur von den Impulsen abhängt. Dann sind die Impulse konstant und die Ableitungen der Hamilton-Funktion nach den Impulsen sind die zeitlich konstanten Geschwindigkeiten, mit denen die Koordinaten <math>q</math> linear zunehmen,

<math>\dot q_k =\frac{\partial \mathcal H}{\partial p_k}=\omega_k(p)\,,\ \dot p_k =-\frac{\partial \mathcal H}{\partial q_k}=0\,,\

p_k(t)=p_k(0)\,,\ q_k(t)=\omega_k(p)\,t+q_k(0)\,.</math> Ist zudem die Phasenraumfläche konstanter Energie <math>\mathcal{H}(q,p)=E</math> kompakt, dann handelt es sich bei den Koordinaten <math>q_k</math> um die Winkel auf einem Torus, die um <math>2\pi</math> vergrößert wieder denselben Punkt benennen,

Der Phasenraum solch eines integrablen Systems besteht aus <math>n</math>-dimensionalen Tori, um die sich die Lösungskurven der Hamiltonschen Gleichungen winden.

Zusammenhang mit der Quantenmechanik[Bearbeiten]

So wie in der Mechanik die Hamilton-Funktion die Zeitentwicklung bestimmt, so bestimmt der Hamilton-Operator die Zeitentwicklung in der Quantenmechanik. Man erhält ihn für viele quantenmechanische Systeme aus der Hamilton-Funktion des entsprechenden klassischen Systems durch kanonische Quantisierung, indem man den algebraischen Ausdruck für <math>\mathcal H(q,p)</math> als Funktion von Operatoren <math>q</math> und <math>p</math> liest, die den kanonischen Vertauschungsrelationen genügen.

Quellen[Bearbeiten]

V. I. Arnold: Mathematical Methods of Classical Mechanics. Springer-Verlag, 1989, ISBN 0-387-96890-3.

Einzelnachweise[Bearbeiten]