Laplace-Runge-Lenz-Vektor

Der Laplace-Runge-Lenz-Vektor (in der Literatur auch Runge-Lenz-Vektor, Lenzscher Vektor etc., nach Pierre-Simon Laplace, Carl Runge und Wilhelm Lenz) ist eine Erhaltungsgröße der Bewegung in einem 1/r-Potential (Coulomb-Potential, Gravitationspotential), d.h., er ist auf jedem Punkt der Bahn gleich (Erhaltungsgröße). Er zeigt vom Brennpunkt der Bahn (Kraftzentrum) zum nächstgelegenen Bahnpunkt (Perihel bei der Erdbahn) und hat somit eine Richtung parallel zur großen Bahnachse. Sein Betrag ist mit der Exzentrizität der Bahn verknüpft. Der Laplace-Runge-Lenz-Vektor ermöglicht daher die elegante Herleitung der Bahnkurve $r(\varphi )$ eines Teilchens (z.B. Planet im Keplerproblem, $\alpha$ -Teilchen gestreut an Atomkern) in diesem Kraftfeld, ohne eine einzige Bewegungsgleichung lösen zu müssen.

In der klassischen Mechanik wird der Vektor hauptsächlich benutzt, um die Form und Orientierung der Umlaufbahn eines astronomischen Körpers um einen anderen zu beschreiben, etwa die Bahn eines Planeten um seinen Stern.

Auch in der Quantenmechanik des Wasserstoffatoms spielt der Vektor als Laplace-Runge-Lenz- oder Laplace-Runge-Lenz-Pauli-Operator eine Rolle.

Definition

Illustration des Laplace-Runge-Lenz-Vektors an einer Ellipsenbahn für vier unterschiedliche Punkte

Sei

$V=-{\frac {k}{r}}$

ein radialsymmetrisches anziehendes Potential, das mit einer Proportionalitätskonstante umgekehrt proportional zum relativen Abstand zweier Objekte ist. Dann ist der Laplace-Runge-Lenz-Vektor ${\vec {A}}$ definiert als

${\vec {A}}={\vec {p}}\times {\vec {L}}-mk{\vec {e}}_{r}$ ,

wobei

${\vec {p}}$ den Impuls des Körpers
${\vec {L}}$ seinen Drehimpuls,
seine Masse und
${\vec {e}}_{r}$ den radialen Einheitsvektor

bezeichnet.

Beweis der Erhaltung

Direkte Rechnung

In einem System mit 1/r-Potential gilt Isotropie. Daher gilt Drehimpulserhaltung mit der Konsequenz, dass die Bewegung in einer Ebene senkrecht zum Drehimpuls stattfindet und es eine einfache Beziehung zwischen Drehimpuls und Winkelgeschwindigkeit gibt:

${\vec {L}}={\vec {r}}\times {\vec {p}}=mr^{2}{\vec {\omega }}=\mathrm {const}$

Die Winkelgeschwindigkeit bestimmt die Zeitableitung des zweiten Terms von ${\vec {A}}$ , denn ein Einheitsvektor kann sich nur durch Drehung ändern:

${\frac {\mathrm {d} }{\mathrm {d} t}}\,mk\,{\vec {e}}_{r}=mk\,{\vec {\omega }}\times {\vec {e}}_{r}$

Die Kraft ist nach dem 2. Newtonschen Gesetz die Änderungsrate des Impulses:

${\vec {F}}=-{\frac {k}{r^{2}}}\,{\vec {e}}_{r}={\frac {\mathrm {d} {\vec {p}}}{\mathrm {d} t}}$

Für den ersten Term von ${\vec {A}}$ gilt damit

${\frac {\mathrm {d} }{\mathrm {d} t}}{\vec {p}}\times {\vec {L}}=\left(-{\frac {k}{r^{2}}}{\vec {e}}_{r}\right)\times \left(mr^{2}{\vec {\omega }}\right)=mk\,{\vec {\omega }}\times {\vec {e}}_{r}.$

Durch Subtrahieren folgt nun die Konstanz des Runge-Lenz-Vektors:

${\frac {\mathrm {d} }{\mathrm {d} t}}{\vec {A}}={\vec {0}}$

Folgerung aus dem Noether-Theorem

Obwohl in der Literatur teilweise vertreten wird, es existiere zum Laplace-Runge-Lenz-Vektor keine zugehörige Symmetrietransformation der Lagrangefunktion, kann diese leicht angegeben werden. Die Lagrangefunktion eines attraktiven 1/r -Potentials lautet:

${\mathcal {L}}={\frac {1}{2}}m{\vec {v}}^{2}+{\frac {k}{r}}$

Die der Erhaltung des Laplace-Rung-Lenz-Vektors zugrunde liegende Symmetrie zeigt sich unter der Variablentransformation

$r_{i}\to r'_{i}=r_{i}+m\varepsilon _{k}\left(v_{i}r_{k}-{\frac {1}{2}}r_{i}v_{k}-{\frac {1}{2}}{\vec {r}}\cdot {\vec {v}}\delta _{ik}\right)$

mit drei infinitesimalen Parametern $\varepsilon _{k}$ . Mithilfe der Bewegungsgleichungen kann die entsprechende Transformation der Geschwindigkeiten als

$v_{i}\to v'_{i}=v_{i}+{\frac {1}{2}}m\varepsilon _{k}\left(v_{i}v_{k}-{\vec {v}}^{2}\delta _{ik}-{\frac {k}{m}}{\frac {r_{i}r_{k}}{r^{3}}}+{\frac {k}{m}}{\frac {\delta _{ik}}{r}}\right)$

identifiziert werden. Durch Einsetzen in die Lagrangefunktion und Taylor-Entwicklung bis zur Ordnung ${\mathcal {O}}(\varepsilon _{k})$ zeigt sich, dass sich diese wie

${\mathcal {L}}\to {\mathcal {L}}'={\mathcal {L}}+mk\varepsilon _{k}\left({\frac {v_{k}}{r}}-{\frac {{\vec {r}}\cdot {\vec {v}}}{r^{3}}}r_{k}\right)={\mathcal {L}}+mk\varepsilon _{k}{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {r_{k}}{r}}$

verhält, wobei der zusätzliche Term eine totale Zeitableitung ist und daher die Wirkung des Systems invariant lässt. Aus dem Noether-Theorem folgt, dass die drei Komponenten des Vektors

$A_{k}=m\left(v_{i}r_{k}-{\frac {1}{2}}r_{i}v_{k}-{\frac {1}{2}}{\vec {r}}\cdot {\vec {v}}\delta _{ik}\right){\frac {\partial {\mathcal {L}}}{\partial v_{i}}}-mk{\frac {r_{k}}{r}}=m^{2}({\vec {v}}\times ({\vec {r}}\times {\vec {v}}))_{k}-mk{\frac {r_{k}}{r}}$

erhalten sind.

Erhaltung im Hamilton-Formalismus

Mit der Hamilton-Funktion des Systems

${\mathcal {H}}={\frac {{\vec {p}}^{2}}{2m}}-{\frac {k}{r}}$

folgt für die partiellen Ableitungen der Hamilton-Funktion und des Laplace-Runge-Lenz-Vektors nach den Koordinaten und Impulsen

${\begin{aligned}&{\frac {\partial {\mathcal {H}}}{\partial r_{i}}}=k{\frac {r_{i}}{r^{3}}}&\qquad &{\frac {\partial {\mathcal {H}}}{\partial p_{i}}}={\frac {p_{i}}{m}}&\\&{\frac {\partial {\vec {A}}}{\partial r_{i}}}=p^{2}{\vec {e}}_{i}-p_{i}{\vec {p}}-mk{\frac {{\vec {e}}_{i}}{r}}+mk{\frac {r_{i}}{r^{3}}}{\vec {r}}&\qquad &{\frac {\partial {\vec {A}}}{\partial p_{i}}}=2p_{i}{\vec {r}}-({\vec {r}}\cdot {\vec {p}}){\vec {e}}_{i}-r_{i}{\vec {p}}&\end{aligned}}$

und nach den Hamiltonschen Bewegungsgleichungen

${\frac {\mathrm {d} {\vec {A}}}{\mathrm {d} t}}={\frac {\partial {\vec {A}}}{\partial t}}+\{{\vec {A}},{\mathcal {H}}\}={\frac {\partial {\vec {A}}}{\partial t}}+\left({\frac {\partial {\vec {A}}}{\partial r_{i}}}{\frac {\partial {\mathcal {H}}}{\partial p_{i}}}-{\frac {\partial {\vec {A}}}{\partial p_{i}}}{\frac {\partial {\mathcal {H}}}{\partial r_{i}}}\right)=0$

Herleitung der Bahnkurve

Hierfür ist normalerweise, d.h., wenn man das Arbeiten mit der Energie als Erhaltungsgröße vorzieht, eine aufwändige Integration mit mehreren Substitutionen nötig. Dagegen folgt aus der Multiplikation des Runge-Lenz-Vektors mit ${\vec {r}}$ nun einfach nach der Kosinusbeziehung des Skalarprodukts (pfeillose Buchstaben kennzeichnen stets die Beträge des zugehörigen Vektors):

${\vec {A}}\cdot {\vec {r}}=Ar\cos \varphi ={\vec {r}}\cdot \left({\vec {p}}\times {\vec {L}}\right)-mkr={\vec {L}}\cdot \left({\vec {r}}\times {\vec {p}}\right)-mkr=L^{2}-mkr$

Hierbei wurden die Zyklizität des Spatproduktes und die Drehimpulsdefinition genutzt. $\varphi$ bezeichnet den Winkel zwischen Runge-Lenz- und Ortsvektor.

Durch Umschreiben entsteht eine Kegelschnittgleichung in Polarkoordinaten:

$r={\frac {L^{2}/(mk)}{1+{\frac {A}{mk}}\cos \varphi }}={\frac {L^{2}/(mk)}{1+\varepsilon _{k}\cos \varphi }}$ ,

wobei der Term $A/mk$ als die numerische Exzentrizität des Kegelschnitts $\varepsilon _{k}$ , die die Bahnform Kreis ( $\varepsilon _{k}=0$ ), Ellipse ( $0<\varepsilon _{k}<1$ ), Parabel ( $\varepsilon _{k}=1$ ) oder Hyperbel ( $\varepsilon _{k}>1$ ) bestimmt, identifiziert werden kann.

Hodograph der Kepler-Bahn; die Punkte 1–4 entsprechen denen in obiger Abbildung

Weiterhin ist ebenfalls die Herleitung des Hodographe der Keplerbahn mithilfe des Laplace-Runge-Lenz-Vektors möglich. Da der Drehimpulsvektor senkrecht auf der Bewegungsebene steht, ${\vec {A}}\cdot {\vec {p}}=0$ , folgt nach

${\frac {mk}{r}}{\vec {r}}={\vec {p}}\times {\vec {L}}-{\vec {A}}$

mit der Lagrange-Identität und einer zyklischen Permutation des Spatprodukts

$m^{2}k^{2}=p^{2}L^{2}-2({\vec {A}}\times {\vec {p}})\cdot {\vec {L}}+A^{2}$ .

Bei einer Wahl des Koordinatensystems, sodass der der Drehimpuls in -Richtung zeigt, ${\vec {L}}=L{\vec {e}}_{z}$ , und der dazu orthogonale Laplace-Runge-Lenz-Vektor in -Richtung, ${\vec {A}}=A{\vec {e}}_{x}$ , folgt:

$p_{x}^{2}+\left(p_{y}-{\frac {A}{L}}\right)^{2}=\left({\frac {mk}{L}}\right)^{2}$

Der Hodograph ist somit ein Kreis mit Radius $mk/L$ , der um $A/L$ vom Zentrum der Kraft verschoben ist. Für die Schnittpunkte des Hodographen mit der -Achse $p_{0}$ gilt:

$p_{0}=\pm {\sqrt {{\frac {m^{2}k^{2}}{L^{2}}}-{\frac {A^{2}}{L^{2}}}}}=\pm {\sqrt {2m|E|}}$

Sie sind somit unabhängig vom Drehimpuls und vom Laplace-Runge-Lenz-Vektor.

Eigenschaften

Der Runge-Lenz-Vektor liegt in der Bahnebene, denn er steht senkrecht zum Drehimpulsvektor:

${\vec {L}}\cdot {\vec {A}}={\vec {L}}\cdot \left({\vec {p}}\times {\vec {L}}\right)-mk{\frac {{\vec {L}}\cdot {\vec {r}}}{r}}={\vec {p}}\cdot \left({\vec {L}}\times {\vec {L}}\right)-mk{\frac {\left({\vec {r}}\times {\vec {p}}\right)\cdot {\vec {r}}}{r}}=0$

Der Runge-Lenz-Vektor zeigt vom Kraftzentrum der Bahn (einem der beiden Brennpunkte) zum Perizentrum, d.h. zentrumnächsten Punkt der Bahn. Dies folgt sofort aus obiger Bahngleichung, da $\varphi$ den Winkel zwischen Orts- und Runge-Lenz-Vektor darstellt und minimal ist für maximalen Nenner, d.h. $\cos \varphi =1\Rightarrow \varphi =0$ .
Der Runge-Lenz-Vektor hat als Betrag das $mk$ -Fache der numerischen Exzentrizität der Bahnkurve. Dies wurde bereits bei der Herleitung derselben gezeigt.
Alle drei Komponenten des Laplace-Runge-Lenz-Vektors sind Erhaltungsgrößen. Da sein Betrag bereits durch die anderen Erhaltungsgrößen Drehimpuls und Energie und seine Lage durch die orthogonalität zum Drehimpulsvektor vorgegeben sind, liefert der Laplace-Runge-Lenz-Vektors nur eine unabhängige Erhaltungsgröße. Das Kepler-Problem hat daher fünf unabhängige Erhaltungsgrößen (Energie, 3 Komponenten des Drehimpulsvektors, Orientierung des Laplace-Runge-Lenz-Vektors) für sechs Anfangsbedingungen; es ist daher ein maximal superintegrables System.

Periheldrehung bei Abweichungen vom Kepler-Potential

Die Erhaltung des Runge-Lenz-Vektors impliziert, dass die Ellipsen der Planetenbewegung im Kepler-Potential eine feststehende Orientierung im Raum haben.

Bei kleinen Abweichungen vom 1/r-Potential, z.B. durch Anwesenheit anderer Planeten im Sonnensystem oder infolge der Einsteinschen Relativitätstheorien, kommt es zu einer langsamen Drehung der Bahnachse (Periheldrehung). Wenn eine Abweichung so klein ist, dass ihr Quadrat vernachlässigt werden kann, so ist die Störung der Kepler-Bahn mit Hilfe des Runge-Lenz-Vektors elementar berechenbar. Es sei $\Phi (r)$ das Störpotential, das zum Kepler-Potential addiert wird. Für den Runge-Lenz-Vektor findet man (vgl. Beweis der Erhaltung)

${\frac {\mathrm {d} {\vec {A}}}{\mathrm {d} t}}=-\partial _{r}\Phi (r){\vec {e}}_{r}\times mr^{2}{\vec {\omega }}=mr^{2}\partial _{r}\Phi (r)~{\vec {e}}_{z}\times {\vec {e}}_{r}~{\frac {\mathrm {d} \varphi }{\mathrm {d} t}}.$

Die z-Richtung steht dabei senkrecht zur Bahnebene. Offenbar ist die Bewegung des Runge-Lenz-Vektors nicht zu jedem Zeitpunkt eine Drehung. Eine Drehung ergibt sich aber, wenn infinitesimale Änderungen über einen Umlauf integriert werden. Dafür findet man zunächst

$\left(\Delta {\vec {A}}\right)_{\mathrm {1\,Umlauf} }=\int _{0}^{T}{\frac {\mathrm {d} {\vec {A}}}{\mathrm {d} t}}\mathrm {d} t=m{\vec {e}}_{z}\times \int _{0}^{2\pi }r^{2}\partial _{r}\Phi (r)~{\vec {e}}_{r}~\mathrm {d} \varphi \qquad \qquad r=r(\varphi ).$

Da quadratische Effekte von $\Phi$ vernachlässigbar sein sollen, kann für $r(\varphi )$ die ungestörte Bahnkurve eingesetzt werden. Der radiale Einheitsvektor, zerlegt in Komponenten parallel und senkrecht zur Bahnachse, ist

${\vec {e}}_{r}(\varphi )={\vec {e}}_{A}\,\cos \varphi +{\vec {e}}_{\perp }\,\sin \varphi .$

Bei der Kepler-Ellipse ist $r(\varphi )$ eine Funktion von $\cos \varphi$ , daher ergibt das Integral über eine Periode mit dem Faktor $\sin \varphi$ für jedes Störpotential $\Phi (r)$ null. Es bleibt nur

$\left(\Delta {\vec {A}}\right)_{\mathrm {1\,Umlauf} }={\vec {e}}_{z}\Delta \varphi \times {\vec {A}},$

wobei ${\vec {A}}=mk\varepsilon {\vec {e}}_{A}$ eingesetzt wurde und der Drehwinkel $\Delta \varphi$ durch folgenden Ausdruck gegeben ist:

$\Delta \varphi ={\frac {1}{k\varepsilon }}\int _{0}^{2\pi }r(\varphi )^{2}\partial _{r}\Phi (r(\varphi ))\cos \varphi ~\mathrm {d} \varphi$

Bei der Störung einer Planetenbahn durch die Anwesenheit anderer Planeten ist das Störpotential nicht unmittelbar von der Form $\Phi (r)$ , erhält aber diese Form durch Mittelung über viele Umläufe von Planeten in einer gemeinsamen Bahnebene.

Quantenmechanik

In der Quantenmechanik kann im Wasserstoffproblem als Analogon zum Laplace-Runge-Lenz-Vektor der hermitesche Operator

${\hat {\vec {A}}}={\frac {1}{2m}}\left({\hat {\vec {p}}}\times {\hat {\vec {L}}}-{\hat {\vec {L}}}\times {\hat {\vec {p}}}\right)+Z\alpha \hbar c{\frac {\hat {\vec {x}}}{|{\hat {\vec {x}}}|}}$

definiert werden, wobei

${\hat {\vec {p}}}$ der Impulsoperator,
${\hat {{\vec L}}}$ der Drehimpulsoperator und
${\hat {\vec {x}}}$ der Ortsoperator sind, sowie
die Kernladungszahl,
$\alpha$ die Feinstrukturkonstante,
$\hbar$ das reduzierte Plancksche Wirkungsquantum,
die Lichtgeschwindigkeit und
die Masse des Elektrons sind.

Insbesondere ist in der Quantenmechanik $\textstyle {\hat {\vec {L}}}\times {\hat {\vec {p}}}\neq -{\hat {\vec {p}}}\times {\hat {\vec {L}}}$ , da der Kommutator zwischen Impuls- und Drehimpulsoperator nicht verschwindet. Der Hamilton-Operator des Coulomb-Problems ${\hat {H}}$ ist

${\hat {H}}={\frac {{\hat {\vec {p}}}^{2}}{2m}}-Z\alpha \hbar c{\frac {1}{|{\hat {\vec {x}}}|}}$

und aus der Definition des Drehimpulsoperators folgt die Kommutatorrelation

$[{\hat {H}},{\hat {A}}_{i}]=0$

für alle Komponenten des Laplace-Runge-Lenz-Operators. Da dieser selbst nicht zeitabhängig ist, folgt aus den Heisenbergschen Bewegungsgleichungen für quantenmechanische Operatoren

${\frac {\mathrm {d} }{\mathrm {d} t}}{\hat {\vec {A}}}={\frac {\partial }{\partial t}}{\hat {\vec {A}}}+{\frac {\mathrm {i} }{\hbar }}[{\hat {H}},{\hat {\vec {A}}}]=0$ .

Aus der Vertauschbarkeit des Hamilton-Operators und des Laplace-Runge-Lenz-Operators folgt, dass beide einen Satz gemeinsamer Eigenzustände besitzen und insbesondere ebenfalls der Hamilton-Operator und das Quadrat des Laplace-Runge-Lenz-Operators.

Die Kommutatorrelationen für die einzelnen Komponenten des Laplace-Runge-Lenz-Operators lauten

$[{\hat {A}}_{i},{\hat {A}}_{j}]=-2\mathrm {i} \hbar {\hat {H}}\varepsilon _{ijk}{\frac {{\hat {L}}_{k}}{m}}$

und für den Kommutator der Komponenten des Laplace-Runge-Lenz-Operator und des Drehimpulsoperators

$[{\hat {A}}_{i},{\hat {L}}_{j}]=\mathrm {i} \hbar c\varepsilon _{ijk}{\hat {L}}_{k}$

mit dem Levi-Civita-Symbol $\varepsilon$ . Insbesondere sind

$[{\hat {\vec {A}}}^{2},A_{j}]=[{\hat {\vec {A}}}^{2},L_{j}]=[{\hat {\vec {A}}}^{2},{\hat {\vec {L}}}^{2}]=[{\hat {\vec {L}}}^{2},{\hat {A}}_{j}]=[{\hat {\vec {L}}}^{2},{\hat {L}}_{j}]=0$ ,

also existiert ein Satz gemeinsamer Eigenzustände zu beiden Sätzen der Operatoren ${\hat {H}},{\hat {\vec {L}}}^{2},{\hat {L}}_{3},{\hat {\vec {A}}}^{2}$ und ${\hat {H}},{\hat {\vec {L}}}^{2},{\hat {\vec {A}}}^{2},{\hat {A}}_{3}$ .

Basierend auf einem Artikel in:

Wikipedia.de