Der Lagrange-Formalismus ist in der Physik eine 1788 von Joseph-Louis Lagrange eingeführte Formulierung der klassischen Mechanik, in der die Dynamik eines Systems durch eine einzige skalare Funktion, die Lagrange-Funktion, beschrieben wird. Der Formalismus ist (im Gegensatz zu der newtonschen Mechanik, die a priori nur in Inertialsystemen gilt) auch in beschleunigten Bezugssystemen gültig. Der Lagrange-Formalismus ist invariant gegen Koordinatentransformationen.[1] Aus der Lagrange-Funktion lassen sich die Bewegungsgleichungen mit den Euler-Lagrange-Gleichungen der Variationsrechnung aus dem Prinzip der kleinsten Wirkung bestimmen. Diese Betrachtungsweise vereinfacht viele physikalische Probleme, da sich, im Gegensatz zu der newtonschen Formulierung der Bewegungsgesetze, im Lagrange-Formalismus Zwangsbedingungen relativ einfach durch das explizite Ausrechnen der Zwangskräfte oder die geeignete Wahl generalisierter Koordinaten berücksichtigen lassen. Aus diesem Grund wird der Lagrange-Formalismus verbreitet bei Mehrkörpersystemen (MKS) eingesetzt. Er lässt sich auch auf den relativistischen Fall übertragen und ist auch in der relativistischen Quantenfeldtheorie zur Formulierung von Modellen von Elementarteilchen und ihrer Wechselwirkungen weit verbreitet.
Für Systeme mit einem generalisierten Potential und holonomen Zwangsbedingungen lautet die Lagrange-Funktion
wobei $ T $ die kinetische Energie und $ V $ die potentielle Energie des betrachteten Systems bezeichnen. Man unterscheidet sogenannte Lagrange-Gleichungen erster und zweiter Art. Im engeren Sinn versteht man unter dem Lagrange-Formalismus und den Lagrange-Gleichungen aber die zweiter Art, die häufig einfach als Lagrange-Gleichungen bezeichnet werden:
Dabei sind $ q_{i} $ generalisierte Koordinaten und $ {\dot {q}}_{i} $ deren Zeitableitungen.
Mit den Lagrange-Gleichungen erster Art lassen sich die Zwangskräfte berechnen. Sie sind äquivalent zu den Gleichungen, die sich aus dem D’Alembertschen Prinzip ergeben. Wir betrachten $ N $ Punktteilchen im $ \ \mathbb {R} ^{3} $ mit den Ortsvektoren $ \ \mathbf {r} _{i} $, $ i\in \{1,...,N\} $, deren Koordinaten durch $ s $ voneinander unabhängige (holonome) Zwangsbedingungen der Form $ F_{k}(\mathbf {r} _{1},\ldots ,\mathbf {r} _{N},t)=0 $ mit $ k\in \{1,\ldots ,s\} $ eingeschränkt sind (eine explizite Zeitabhängigkeit ist erlaubt). Dadurch werden die Lagen der Teilchen auf eine $ (3N-s) $-dimensionale Mannigfaltigkeit eingeschränkt ($ f=3N-s $ ist die Anzahl der Freiheitsgrade).
Die auf ein Teilchen $ i $ wirkenden Zwangskräfte sind proportional zum Gradienten $ \nabla F_{k} $, die Gesamt-Zwangskraft $ \mathbf {Z} _{i} $ ist daher
Wenn man annimmt, dass sich die äußeren Kräfte aus einem Potential ableiten lassen, kann man die Bewegungsgleichung schreiben (Lagrange-Gleichung 1. Art):[2]
Die $ m_{i} $ sind die Massen der $ N $ Punktteilchen, $ V $ ist die potentielle Energie. Dies, zusammen mit den Zwangsbedingungen $ F_{k}(\mathbf {r} _{1},\ldots ,\mathbf {r} _{N},t)=0 $, sind $ 3N+s $ unabhängige Gleichungen für die $ 3N $ Koordinaten der $ \mathbf {r} _{i} $ sowie für die $ s $ Lagrange-Multiplikatoren $ \lambda _{k} $. Somit ist die Lösung des Gleichungssystems eindeutig.
Bemerkung: Hier wurden nur holonome Zwangsbedingungen behandelt. Der Formalismus lässt sich aber auch auf Zwangsbedingungen der Form $ \sum _{k}a_{k}\delta q_{k}=0\, $ anwenden, die z. B. bei nicht-holonomen Zwangsbedingungen zwischen den Geschwindigkeiten der Teilchen folgen.[3] Diese Zwangsbedingungsgleichungen lassen sich im Gegensatz zu holonomen Zwangsbedingungen nicht als vollständiges Differential einer Funktion darstellen, das heißt, zwischen den Koeffizientenfunktionen gilt nicht $ {\tfrac {\partial a_{i}}{\partial q_{k}}}={\tfrac {\partial a_{k}}{\partial q_{i}}} $.
Im Fall von holonomen Zwangsbedingungen kann man neue Koordinaten $ q_{i} $ einführen, die diese implizit enthalten, sogenannte generalisierte Koordinaten. Mit der kinetischen Energie
und Potentialkräften
(die auch durch generalisierte Koordinaten ausgedrückt sind und dann als generalisierte Kräfte bezeichnet werden – sie haben nicht unbedingt die Dimension einer Kraft) lassen sich die Bewegungsgleichungen auch schreiben
oder mit der Lagrange-Funktion $ L=T-V $ (Lagrange-Gleichung 2. Art):
Treten wie in diesem Fall nur aus einem Potential ableitbare Kräfte (Potentialkräfte) auf, spricht man von konservativen Kräften.
Bemerkung: Manchmal lassen sich die generalisierten Kräfte durch ein geschwindigkeitsabhängiges generalisiertes Potential $ V(q_{1},\ldots ,q_{n},{\dot {q}}_{1},\ldots ,{\dot {q}}_{n},t) $ in folgender Form schreiben
Auch dann ergeben sich die Bewegungsgleichungen
mit der Lagrange-Funktion $ L $:
Das System ist dann aber nicht mehr im üblichen Sinn konservativ. Ein Beispiel ist das elektromagnetische Feld (siehe unten).
Manchmal hat man aber noch nicht-konservative Kräfte $ Q_{i}^{*} $, so dass sich die Gleichungen schreiben:
Ein Beispiel sind Systeme mit nicht-holonomen Zwangsbedingungen (siehe oben) oder Reibungskräften (zum Beispiel Rayleighsche Dissipationsfunktion).
Die Lagrange-Gleichungen zweiter Art ergeben sich als sogenannte Euler-Lagrange-Gleichungen[4] eines Variationsproblems und liefern die Bewegungsgleichungen, wenn die Lagrange-Funktion gegeben ist. Sie folgen aus der Variation des mit der Lagrange-Funktion gebildeten Wirkungsintegrals im Hamiltonschen Prinzip. Dazu betrachtet man alle möglichen Bahnkurven $ q(t) $ im Raum der generalisierten Koordinaten zwischen festen Anfangs- und Endpunkten. Man betrachtet die Änderung des Wirkungsintegrals bei Variation der Bahnkurven
Das hamiltonsche Prinzip besagt, dass für die klassische Bahn das Wirkungsintegral stationär unter Variation der Bahnkurven ist:
Eine Näherung in erster Ordnung lautet für eine gewöhnliche Funktion $ f(x,y) $
also
In erster Ordnung ergibt sich die Variation des Integrals also zu
Nun führt man eine partielle Integration in dem Term aus, der die Ableitung nach der Zeit enthält:
Hierbei wird benutzt, dass
ist, da Anfangs- und Endpunkt festgehalten werden. Daher gilt für die Randterme
Damit resultiert schließlich
Da nun $ \delta q $ als Faktor des gesamten Integrals auftritt und beliebig ist, kann das Integral nur dann nach dem Variationsprinzip verschwinden, wenn der Integrand selbst verschwindet. Es folgen die Lagrange-Gleichungen oder Lagrange-Gleichungen zweiter Art (die Euler-Lagrange-Gleichungen des hier betrachteten Variationsproblems):
Für jede generalisierte Koordinate $ q_{i} $ (und die zugehörige generalisierte Geschwindigkeit $ {\dot {q}}_{i} $) gibt es eine solche Gleichung. Die Lagrange-Gleichungen bilden ein System gewöhnlicher Differentialgleichungen zweiter Ordnung bezüglich der Zeitableitung. Wie viele Differentialgleichungen das im Endeffekt sind, weiß man erst, wenn die Zahl der Freiheitsgrade des Systems berechnet wurde.
Wenn die Lagrange-Funktion $ L $ nicht von einer Koordinate $ q $ abhängt, sondern nur von der zugehörigen Geschwindigkeit $ {\dot {q}} $, dann nennt man $ q $ zyklisch, zyklische Koordinate oder zyklische Variable. Der zur zyklischen Variablen $ q $ konjugierte Impuls
ist eine Erhaltungsgröße; sein Wert ändert sich nicht während der Bewegung, wie gleich gezeigt wird: Wenn die Lagrange-Funktion nicht von $ q $ abhängt, gilt
Dann folgt aber aus der Euler-Lagrange-Gleichung, dass die Zeitableitung des zugehörigen konjugierten Impulses verschwindet und er somit zeitlich konstant ist:
Allgemeiner gehört nach dem Noether-Theorem zu jeder kontinuierlichen Symmetrie der Wirkung eine Erhaltungsgröße. Bei einer zyklischen Variablen ist die Wirkung invariant unter der Verschiebung von $ q $ um eine beliebige Konstante, $ q\rightarrow q+c\,. $
In der Feldtheorie ergibt sich die Bewegungsgleichung aus dem hamiltonschen Prinzip für Felder zu
wobei $ \phi =\phi (x,y,z,t) $ das betrachtete Feld und $ {\mathcal {L}}={\mathcal {L}}\left(\phi ,{\frac {\partial \phi }{\partial x}},{\frac {\partial \phi }{\partial y}},{\frac {\partial \phi }{\partial z}},{\frac {\partial \phi }{\partial t}},x,y,z,t\right) $ die Lagrange-Dichte sind.
Man kann dies in Kurzform auch schreiben als
mit der so definierten Variationsableitung $ {\frac {\delta {\mathcal {L}}}{\delta \phi }}:={\frac {\partial {\mathcal {L}}}{\partial \phi _{i}}}-{\partial _{\mu }}\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{\mu }\phi _{i})}}\right) $.
Hinweis: Der Lagrange-Formalismus ist auch der Ausgangspunkt vieler Formulierungen der Quantenfeldtheorie.
In der relativistischen Mechanik kann die Lagrange-Funktion eines freien Teilchens aus dem hamiltonschen Prinzip abgeleitet werden, indem für die Wirkung der einfachste Fall eines relativistischen Skalars angenommen wird:
wobei $ \mathrm {d} s=c\,\mathrm {d} \tau =c\,\mathrm {d} t{\sqrt {1-{\frac {v^{2}}{c^{2}}}}} $ das zur Eigenzeit proportionale relativistische Linienelement ist und ein konstanter Faktor $ \,(-mc) $ gewählt wurde.
Die Lagrange-Funktion eines freien Teilchens ist hier nicht mehr mit der kinetischen Energie identisch (manchmal spricht man deshalb auch von kinetischer Ergänzungsenergie T in der Lagrange-Funktion). Die relativistische kinetische Energie eines Körpers mit der Masse $ m $ und Geschwindigkeit $ v={\dot {\mathbf {x} }} $ ohne Zwangsbedingungen beträgt
wohingegen für die Lagrange-Funktion die kinetische Ergänzungsenergie
maßgeblich ist. Die Lagrange-Funktion für ein Teilchen in einem Potential V ergibt sich dann zu
Für ein $ N $-Teilchensystem ist die Lagrange-Funktion mit den generalisierten Koordinaten
wobei $ n=3N-s $ die Anzahl der Freiheitsgrade und $ s $ die Anzahl der holonomen Zwangsbedingungen ist.
Für kleine Geschwindigkeiten $ |{\dot {\mathbf {x} }}|\ll c $ kann man die Wurzel bis zur ersten Ordnung entwickeln $ {\sqrt {1-x}}=1-x/2 $:
Die nullte Ordnung der Entwicklung ist eine Konstante, die negative Ruheenergie. Da die Lagrange-Gleichungen invariant sind unter Addition einer Konstanten zur Lagrange-Funktion, kann man den konstanten ersten Term weglassen und man erhält wieder die klassische kinetische Energie:
Richard Feynman hat als Erster diese Herangehensweise auch konsequent für die Herleitung der Gleichungen der Quantenmechanik verwendet. In der klassischen Physik ergeben sich die oben beschriebenen Lagrange-Gleichungen aus der Forderung, dass das Wirkungsintegral stationär wird. In Feynmans Pfadintegral-Formalismus ist die quantenmechanische Wahrscheinlichkeitsamplitude, dass ein System zwischen Anfangs- und Endbedingungen einen bestimmten Pfad einschlägt, proportional zu $ e^{\frac {iW}{\hbar }} $ mit dem Wirkungsintegral $ W $. Pfade in der Umgebung des klassischen Weges, für den die Variation von $ W $ verschwindet, liefern dabei meist die Hauptbeiträge, da sich in ihrer Umgebung die Beiträge mit fast gleichen Phasenfaktoren addieren.
Eine Masse $ m $ sei über zwei Federn mit Federkonstante $ c $ und festen Randbedingungen verbunden. Grundvoraussetzung zur Beschreibung des Problems im Lagrange-Formalismus ist das Aufstellen der Lagrange-Funktion, indem man die Terme für kinetische Energie $ T $ und potentielle Energie $ V $ aufstellt:
Die Lagrange-Funktion lautet daher:
Die Lagrange-Funktion wiederum wird zur analytischen Beschreibung des physikalischen Problems in die Euler-Lagrange-Gleichung eingesetzt, was dann auf Gleichungen führt, die den Bewegungsgleichungen in der Newtonschen Mechanik entsprechen. In diesem Beispiel lautet die generalisierte Koordinate $ x $, die Euler-Lagrange-Gleichung
Dies führt mit obigen Formeln für $ L $ auf
und damit auf die Bewegungsgleichung des Systems:
Die allgemeine Lösung dieser Differentialgleichung ist $ x(t)=A\cos(\omega t+\varphi ) $, $ t $ ist die Zeit, $ \textstyle \omega ={\sqrt {c/m}} $ die Kreisfrequenz. Die konstante Amplitude $ A $ und Phase $ \varphi $ können aus den Anfangsbedingungen bestimmt werden.
Eine Punktladung $ q $ mit Masse $ m $ bewege sich im elektromagnetischen Feld. Die generalisierten Koordinaten entsprechen den kartesischen Koordinaten in 3 Raumdimensionen.
Die Felder (Magnetfeld $ \mathbf {B} $ und elektrisches Feld $ \mathbf {E} $) werden über das Skalarpotential $ \phi $ und das Vektorpotential $ \mathbf {A} $ bestimmt:
Die kinetische Energie des Teilchens ist klassisch:
Das „Potential“ ist hier allerdings geschwindigkeitsabhängig, man spricht deshalb wie oben dargestellt von einem generalisierten Potential:
Somit ist die Lagrange-Funktion eines geladenen Teilchens im elektromagnetischen Feld:
Die Euler-Lagrange-Gleichung $ {\frac {d}{dt}}\nabla _{\dot {\mathbf {x} }}L-\nabla _{\mathbf {x} }L=0 $ führt auf die Bewegungsgleichung, auf deren rechter Seite die Lorentzkraft steht:
Die Achse einer Aufzugtrommel wird durch ein Drehmoment $ M $ angetrieben. Die Masse der Last beträgt $ m $, das Massenträgheitsmoment der Trommel ist $ J $. Der Radius der Trommel ist $ r $.
Zwischen den Koordinaten $ x $ und $ \varphi $ besteht folgende Beziehung:
$ \Rightarrow \;{\dot {x}}=r{\dot {\varphi }} $
$ \Rightarrow \;\delta x=r\delta \varphi $
Die kinetische Energie ist:
Die virtuelle Arbeit der eingeprägten Kräfte ist
$ \Rightarrow \;Q=-mgr+M $
Daraus folgt schließlich die Bewegungsgleichung
Die Auflösung dieser Gleichung nach der Winkelbeschleunigung ergibt
Bei der Atwoodschen Fallmaschine betrachtet man zwei Punktmassen im Gravitationsfeld der Erde, die über eine Rolle in der Höhe h aufgehängt und durch ein Seil der Länge l verbunden seien. Die Zwangsbedingung lautet in diesem Fall:
Wird das Seil berücksichtigt, das auf der Rolle (Rollenradius r) liegt, dann ergibt sich:
Die potentielle Energie V berechnet sich zu:
Für die Gradienten erhält man
Dies führt auf das System der Lagrange-Gleichungen 1. Art:
Dies kann man auflösen und erhält z. B. für bekannte Anfangsbedingungen:
Die 1. Masse ($ m_{1} $) ist auf einer dünnen Platte durch ein Loch in der Mitte der Platte durch ein Seil mit konstanter Länge ($ l $) mit einer 2. Masse ($ m_{2} $) verbunden, die sich nur in z-Richtung bewegt (die z-Achse zeige in Richtung Erdmittelpunkt).
Die Zwangsbedingungen $ g_{i} $ lauten:
Aus 4 Zwangsbedingungen bei 2 Massen im $ \mathbb {R} ^{3} $ ergeben sich $ 6-4=2 $ Freiheitsgrade. Für dieses Problem empfiehlt es sich aufgrund der Azimutalsymmetrie Zylinderkoordinaten zu verwenden. So können die generalisierten Koordinaten einfach bestimmt werden.
In Zylinderkoordinaten können die beiden generalisierten Koordinaten nun als
gewählt werden, wobei mittels der 4. Zwangsbedingung auch die Bewegung der $ m_{2} $ durch $ r_{1} $ beschrieben wird;
Die kinetische Energie des Systems lautet nun
Da $ z_{1}=0 $ und sich damit nur die potentielle Energie bei der 2. Masse verändert, lautet sie
Daraus folgt dann die Lagrangefunktion
Da bei dieser Problemstellung zwei generalisierte Koordinaten vorliegen, folgt jeweils eine Bewegungsgleichung für $ r_{1} $ und $ \phi _{1} $:
Aus der Gleichung für $ \phi $ zeigt sich in der Form
die Existenz einer Erhaltungsgröße, des Drehimpulses in $ z $-Richtung $ L_{1}^{z}=m_{1}r_{1}^{2}{\dot {\phi }}_{1} $, der nach dem Noether-Theorem aus der Unabhängigkeit der Lagrangefunktion von der Variablen $ q_{2}=\phi _{1} $ folgt.
In der allgemeinen Relativitätstheorie durchlaufen frei fallende Teilchen Weltlinien längster Zeit: Zwischen zwei (genügend nah beieinander liegenden) Ereignissen $ A $ und $ B $ vergeht auf einer mitgeführten Uhr auf der Weltlinie frei fallender Teilchen mehr Zeit als auf allen anderen Weltlinien durch diese Ereignisse. Sei $ s $ ein entlang des Pfades monoton wachsender Laufparameter, so ergibt sich die verstrichene Zeit zu
mit der Lagrange-Funktion
Dabei sind $ g_{mn}(x) $ die Komponentenfunktionen der Metrik (sowohl Raum- als auch Zeitkomponenten). Wir rechnen einfachheitshalber in Maßsystemen, in denen die Lichtgeschwindigkeit dimensionslos ist und den Wert $ c=1 $ hat, und verwenden die Einsteinsche Summenkonvention.
Der zu $ x^{k} $ konjugierte Impuls ist
und die Euler-Lagrange-Gleichungen lauten
Verwenden wir hier als Abkürzung das Christoffel-Symbol
so erweist sich die Weltlinie längster Dauer als Gerade: Die Richtung der Tangente an die Weltlinie
ändert sich nicht bei Parallelverschiebung längs der Weltlinie
Die Parametrisierung wird nicht festgelegt. Verfügen wir so über sie, dass der Tangentialvektor überall gleich lang ist, dann ist $ {\sqrt {g_{mn}\,{\dot {x}}^{m}\,{\dot {x}}^{n}}} $ konstant und der Tangentialvektor geht beim Durchlaufen der Weltlinie in sich über. Sie erfüllt die Geodätengleichung
Dies ist die allgemein-relativistische Form der Bewegungsgleichung eines frei fallenden Teilchens. Die Gravitation ist in den $ \Gamma _{rs}^{l} $ voll berücksichtigt.
Der Lagrange-Formalismus wird in vielen ein- und weiterführenden Lehrbüchern der klassischen Mechanik behandelt.
Literatur zu Pfadintegralen.