L33 LAPLACE RUNGE LENZ VECTOR LRL Vector Classical Mechanics B

In classical mechanics, the Laplace-Runge-Lenz vector (LRL vector) is a vector used chiefly to describe the shape and orientation of the orbit of one astronomical body around another, such as a binary star or a planet revolving around a star. For two bodies interacting by Newtonian gravity, the LRL vector.. Der Laplace-Runge-Lenz-Vektor (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 (z. B. Coulomb-Potential, Gravitationspotential), d. h. er ist auf jedem Punkt der Bahn gleich. Er zeigt vom Brennpunkt der Bahn (Kraftzentrum) zum nächstgelegenen Bahnpunkt (Perihel bei der.

Der Laplace-Runge-Lenz-Vektor (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 (z. B. Coulomb-Potential, Gravitationspotential), d. h. er ist auf jedem Punkt der Bahn gleich. 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
(
φ
)

{\displaystyle r(\varphi )}

eines Teilchens (z. B. Planet im Keplerproblem,

α

{\displaystyle \alpha }

-Teilchen gestreut an Atomkern) in diesem Kraftfeld, ohne eine einzige Bewegungsgleichung lösen zu müssen. Die Verwendung dieses Vektors geht mindestens bis auf Jakob Hermann (Mathematiker) zurück. Lenz selber verweist in seiner Originalarbeit auf ein Lehrbuch zur Vektoranalysis von Carl Runge, der seinerseits nie behauptet hat, es handele sich dabei um einen originalen Beitrag.

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

Sei

V
=
−

k
r

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

ein radialsymmetrisches anziehendes Potential, das mit einer Proportionalitätskonstante

k

{\displaystyle k}

umgekehrt proportional zum relativen Abstand zweier Objekte

r

{\displaystyle r}

ist. Dann ist der Laplace-Runge-Lenz-Vektor

A
→

{\displaystyle {\vec {A}}}

definiert als

A
→

=

p
→

×

L
→

−
m
k

e
→

r

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

,

wobei

• p
  →

  {\displaystyle {\vec {p}}}

  den Impuls des Körpers

• L
  →

  {\displaystyle {\vec {L}}}

  seinen Drehimpuls,

• m

  {\displaystyle m}

  seine Masse und

• e
  →

  r

  {\displaystyle {\vec {e}}_{r}}

  den radialen Einheitsvektor

bezeichnet.

Exzentrizitätsvektor

Der dimensionslose Vektor

e
→

=

A
→

m
k

=

p
→

×

L
→

m
k

−

e
→

r

{\displaystyle {\vec {e}}={\frac {\vec {A}}{mk}}={\frac {{\vec {p}}\times {\vec {L}}}{mk}}-{\vec {e}}_{r}}

heißt Exzentrizitätsvektor. Speziell im Keplerproblem mit

k
=
G
M
m

{\displaystyle k=GMm}

kann die Definitionsgleichung in eine Form überführt werden, in der die Masse

m

{\displaystyle m}

des bewegten Körpers (z. B. eines beobachteten Satelliten) eliminiert wurde:

e
→

=

v
→

×

h
→

G
M

−

e
→

r

{\displaystyle {\vec {e}}={\frac {{\vec {v}}\times {\vec {h}}}{GM}}-{\vec {e}}_{r}}

Dabei ist

• v
  →

  {\displaystyle {\vec {v}}}

  der Geschwindigkeitsvektor

• h
  →

  {\displaystyle {\vec {h}}}

  der spezifische Drehimpuls

• G

  {\displaystyle G}

  die Gravitationskonstante und

• M

  {\displaystyle M}

  die Masse im Gravitationszentrum, z. B. eines Sterns oder Planeten.

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:

L
→

=

r
→

×

p
→

=
m

r

2

ω
→

=

c
o
n
s
t

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

Die Winkelgeschwindigkeit bestimmt die Zeitableitung des zweiten Terms von

A
→

{\displaystyle {\vec {A}}}

, denn ein Einheitsvektor kann sich nur durch Drehung ändern:

d

d

t

m
k

e
→

r

=
m
k

ω
→

×

e
→

r

{\displaystyle {\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 (und wirkt Richtung Zentrum):

F
→

=
−

k

r

2

e
→

r

=

d

p
→

d

t

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

Für den ersten Term von

A
→

{\displaystyle {\vec {A}}}

gilt damit (beim Tauschen der Vektoren ändert sich ein Vorzeichen)

d

d

t

p
→

×

L
→

=

(

−

k

r

2

e
→

r

)

×

(

m

r

2

ω
→

)

=
m
k

ω
→

×

e
→

r

.

{\displaystyle {\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:

d

d

t

A
→

=

0
→

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

Herleitung mithilfe der bac-cab-Regel

A
=

p
→

×

L
→

−
m
k

e
→

r

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

Aus

L
→

˙

=
0

{\displaystyle {\dot {\vec {L}}}=0}

und der Produktregel bei Ableitungen folgt

d

d

t

A
→

=

p
→

˙

×

L
→

−
m
k

e
→

˙

r

{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\vec {A}}={\dot {\vec {p}}}\times {\vec {L}}-mk{\dot {\vec {e}}}_{r}}

mit

p
→

˙

=

F
→

=
−

k

r

2

e
→

r

{\displaystyle {\dot {\vec {p}}}={\vec {F}}=-{\frac {k}{r^{2}}}{\vec {e}}_{r}}

und

p
→

=

r
→

˙

m

{\displaystyle {\vec {p}}={\dot {\vec {r}}}m}

folgt:

d

A
→

d

t

=
−

k

e
^

r

r

2

×
(

r
→

×

r
→

˙

m

)
−

k

m

e
^

˙

r

=
−
m
k

r
→

r

3

×
(

r
→

×

r
→

˙

)
−
m
k

d

r
→

r

d

t

{\displaystyle {\frac {\mathrm {d} {\vec {A}}}{\mathrm {d} t}}=-{\frac {{k}{\hat {e}}_{r}}{r^{2}}}\times ({\vec {r}}\times {\dot {\vec {r}}}{m})-{k}{m}{\dot {\hat {e}}}_{r}=-mk{\frac {\vec {r}}{r^{3}}}\times ({\vec {r}}\times {\dot {\vec {r}}})-mk{\frac {\mathrm {d} {\frac {\vec {r}}{r}}}{\mathrm {d} t}}}

Nun wird

a
→

×
(

b
→

×

c
→

)
=

b
→

(
a
c
)
−

c
→

(
a
b
)

{\displaystyle {\vec {a}}\times ({\vec {b}}\times {\vec {c}})={\vec {b}}(ac)-{\vec {c}}(ab)}

angewendet. Beim zweiten Term wird die Quotientenregel und Kettenregel angewendet:

d

A
→

d

t

=
m
k

[

−

r
→

(

r
→

⋅

r
→

˙

r

3

)

+

r
→

˙

(

r
→

⋅

r
→

r

3

)

−

r
→

˙

1
r

+

r
→

r
→

⋅

r
→

˙

r

3

]

=
0

{\displaystyle {\frac {\mathrm {d} {\vec {A}}}{\mathrm {d} t}}=mk\left[-{\vec {r}}\left({\frac {{\vec {r}}\cdot {\dot {\vec {r}}}}{r^{3}}}\right)+{\dot {\vec {r}}}\left({\frac {{\vec {r}}\cdot {\vec {r}}}{r^{3}}}\right)-{\dot {\vec {r}}}{\frac {1}{r}}+{\vec {r}}{\frac {{\vec {r}}\cdot {\dot {\vec {r}}}}{r^{3}}}\right]=0}

Folgerung aus dem Noether-Theorem

In der Literatur findet sich teilweise die Aussage, es existiere zum Laplace-Runge-Lenz-Vektor keine zugehörige Symmetrietransformation der Lagrangefunktion. Dies ist nur dann korrekt, wenn man die Menge der Symmetrietransformationen auf orts- und zeitabhängige Transformationen einschränkt. Lässt man in einer allgemeineren Formulierung des Noether-Theorems auch geschwindigkeitsabhängige Transformationen zu, kann eine solche angegeben werden.

Die Lagrangefunktion eines attraktiven

1

/

r

{\displaystyle 1/r}

-Potentials lautet:

L

=

1
2

m

v
→

2

+

k
r

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

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

r

i

→

r

i

′

=

r

i

+
m

ε

k

(

v

i

r

k

−

1
2

r

i

v

k

−

1
2

r
→

⋅

v
→

δ

i
k

)

{\displaystyle 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

ε

k

{\displaystyle \varepsilon _{k}}

. Mithilfe der Bewegungsgleichungen kann die entsprechende Transformation der Geschwindigkeiten als

v

i

→

v

i

′

=

v

i

+

1
2

m

ε

k

(

v

i

v

k

−

v
→

2

δ

i
k

−

k
m

r

i

r

k

r

3

+

k
m

δ

i
k

r

)

{\displaystyle 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

O

(

ε

k

)

{\displaystyle {\mathcal {O}}(\varepsilon _{k})}

zeigt sich, dass sich diese wie

L

→

L

′

=

L

+
m
k

ε

k

(

v

k

r

−

r
→

⋅

v
→

r

3

r

k

)

=

L

+
m
k

ε

k

d

d

t

r

k

r

{\displaystyle {\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

(

v

i

r

k

−

1
2

r

i

v

k

−

1
2

r
→

⋅

v
→

δ

i
k

)

∂

L

∂

v

i

−
m
k

r

k

r

=

m

2

(

v
→

×
(

r
→

×

v
→

)

)

k

−
m
k

r

k

r

{\displaystyle 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

H

=

p
→

2

2
m

−

k
r

{\displaystyle {\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

∂

H

∂

r

i

=
k

r

i

r

3

∂

H

∂

p

i

=

p

i

m

∂

A
→

∂

r

i

=

p

2

e
→

i

−

p

i

p
→

−
m
k

e
→

i

r

+
m
k

r

i

r

3

r
→

∂

A
→

∂

p

i

=
2

p

i

r
→

−
(

r
→

⋅

p
→

)

e
→

i

−

r

i

p
→

{\displaystyle {\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

d

A
→

d

t

=

∂

A
→

∂
t

+
{

A
→

,

H

}
=

∂

A
→

∂
t

+

(

∂

A
→

∂

r

i

∂

H

∂

p

i

−

∂

A
→

∂

p

i

∂

H

∂

r

i

)

=
0

{\displaystyle {\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

r
→

{\displaystyle {\vec {r}}}

nun einfach nach der Kosinusbeziehung des Skalarprodukts (pfeillose Buchstaben kennzeichnen stets die Beträge des zugehörigen Vektors):

A
→

⋅

r
→

=
A
r
cos
⁡
φ
=

r
→

⋅

(

p
→

×

L
→

)

−
m
k
r
=

L
→

⋅

(

r
→

×

p
→

)

−
m
k
r
=

L

2

−
m
k
r

{\displaystyle {\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.

φ

{\displaystyle \varphi }

bezeichnet den Winkel zwischen Runge-Lenz- und Ortsvektor.

Durch Umschreiben entsteht eine Kegelschnittgleichung in Polarkoordinaten:

r
=

L

2

/

(
m
k
)

1
+

A

m
k

cos
⁡
φ

=

L

2

/

(
m
k
)

1
+

ε

k

cos
⁡
φ

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

,

wobei der Term

A

/

m
k

{\displaystyle A/mk}

als die numerische Exzentrizität des Kegelschnitts

ε

k

{\displaystyle \varepsilon _{k}}

, die die Bahnform Kreis (

ε

k

=
0

{\displaystyle \varepsilon _{k}=0}

), Ellipse (

0
<

ε

k

<
1

{\displaystyle 0<\varepsilon _{k}<1}

), Parabel (

ε

k

=
1

{\displaystyle \varepsilon _{k}=1}

) oder Hyperbel (

ε

k

>
1

{\displaystyle \varepsilon _{k}>1}

) bestimmt, identifiziert werden kann.

Weiterhin ist ebenfalls die Herleitung des Hodographen der Keplerbahn mithilfe des Laplace-Runge-Lenz-Vektors möglich. Da der Drehimpulsvektor senkrecht auf der Bewegungsebene steht und damit orthogonal zum Impuls ist, also

L
→

⋅

p
→

=
0

{\displaystyle {\vec {L}}\cdot {\vec {p}}=0}

, folgt nach

m
k

r

r
→

=

p
→

×

L
→

−

A
→

{\displaystyle {\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
(

A
→

×

p
→

)
⋅

L
→

+

A

2

{\displaystyle 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

z

{\displaystyle z}

-Richtung zeigt,

L
→

=
L

e
→

z

{\displaystyle {\vec {L}}=L{\vec {e}}_{z}}

, und der dazu orthogonale Laplace-Runge-Lenz-Vektor in

x

{\displaystyle x}

-Richtung,

A
→

=
A

e
→

x

{\displaystyle {\vec {A}}=A{\vec {e}}_{x}}

, folgt:

p

x

2

+

(

p

y

−

A
L

)

2

=

(

m
k

L

)

2

{\displaystyle 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

m
k

/

L

{\displaystyle mk/L}

, der um

A

/

L

{\displaystyle A/L}

vom Zentrum der Kraft verschoben ist. Für die Schnittpunkte des Hodographen mit der

x

{\displaystyle x}

-Achse

p

0

{\displaystyle p_{0}}

gilt:

p

0

=
±

m

2

k

2

L

2

−

A

2

L

2

=
±

2
m

|

E

|

{\displaystyle 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:

L
→

⋅

A
→

=

L
→

⋅

(

p
→

×

L
→

)

−
m
k

L
→

⋅

r
→

r

=

p
→

⋅

(

L
→

×

L
→

)

−
m
k

(

r
→

×

p
→

)

⋅

r
→

r

=
0

{\displaystyle {\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

  φ

  {\displaystyle \varphi }

  den Winkel zwischen Orts- und Runge-Lenz-Vektor darstellt und

  r

  {\displaystyle r}

  minimal ist für maximalen Nenner, d. h.

  cos
  ⁡
  φ
  =
  1
  ⇒
  φ
  =
  0

  {\displaystyle \cos \varphi =1\Rightarrow \varphi =0}

  .

• Der Runge-Lenz-Vektor hat als Betrag das

  m
  k

  {\displaystyle mk}

  -Fache der numerischen Exzentrizität der Bahnkurve. Dies wurde bereits bei der Herleitung derselben gezeigt. Da der Exzentrizitätsvektor der Runge-Lenz-Vektor geteilt durch

  m
  k

  {\displaystyle mk}

  ist, ist dessen Betrag gleich der numerischen Exzentrizität der Bahnkurve.

• 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, drei Komponenten des Drehimpulsvektors, Orientierung des Laplace-Runge-Lenz-Vektors) für sechs Anfangsbedingungen; es ist daher ein maximal superintegrables System.

Durch den Runge-Lenz-Vektor ist damit aus den Positions- und Geschwindigkeitvektoren eines beobachteten Objekts die vollständige Form und Orientierung seiner Bahn festgelegt, die in einer Ebene senkrecht zum Dreimpulsvektor liegt.

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

Φ
(
r
)

{\displaystyle \Phi (r)}

das Störpotential, das zum Kepler-Potential addiert wird. Für den Runge-Lenz-Vektor findet man (vgl. Beweis der Erhaltung)

d

A
→

d

t

=
−

∂

r

Φ
(
r
)

e
→

r

×
m

r

2

ω
→

=
m

r

2

∂

r

Φ
(
r
)
 

e
→

z

×

e
→

r

 

d

φ

d

t

.

{\displaystyle {\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

(

Δ

A
→

)

1

U
m
l
a
u
f

=

∫

0

T

d

A
→

d

t

d

t
=
m

e
→

z

×

∫

0

2
π

r

2

∂

r

Φ
(
r
)
 

e
→

r

 

d

φ

r
=
r
(
φ
)
.

{\displaystyle \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

Φ

{\displaystyle \Phi }

vernachlässigbar sein sollen, kann für

r
(
φ
)

{\displaystyle r(\varphi )}

die ungestörte Bahnkurve eingesetzt werden. Der radiale Einheitsvektor, zerlegt in Komponenten parallel und senkrecht zur Bahnachse, ist

e
→

r

(
φ
)
=

e
→

A

cos
⁡
φ
+

e
→

⊥

sin
⁡
φ
.

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

Bei der Kepler-Ellipse ist

r
(
φ
)

{\displaystyle r(\varphi )}

eine Funktion von

cos
⁡
φ

{\displaystyle \cos \varphi }

, daher ergibt das Integral über eine Periode mit dem Faktor

sin
⁡
φ

{\displaystyle \sin \varphi }

für jedes Störpotential

Φ
(
r
)

{\displaystyle \Phi (r)}

null. Es bleibt nur

(

Δ

A
→

)

1

U
m
l
a
u
f

=

e
→

z

Δ
φ
×

A
→

,

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

wobei

A
→

=
m
k
ε

e
→

A

{\displaystyle {\vec {A}}=mk\varepsilon {\vec {e}}_{A}}

eingesetzt wurde und der Drehwinkel

Δ
φ

{\displaystyle \Delta \varphi }

durch folgenden Ausdruck gegeben ist:

Δ
φ
=

1

k
ε

∫

0

2
π

r
(
φ

)

2

∂

r

Φ
(
r
(
φ
)
)
cos
⁡
φ
 

d

φ

{\displaystyle \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

Φ
(
r
)

{\displaystyle \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

A
→

^

=

1

2
m

(

p
→

^

×

L
→

^

−

L
→

^

×

p
→

^

)

+
Z
α
ℏ
c

x
→

^

|

x
→

^

|

{\displaystyle {\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

• p
  →

  ^

  {\displaystyle {\hat {\vec {p}}}}

  der Impulsoperator,

• L
  →

  ^

  {\displaystyle {\hat {\vec {L}}}}

  der Drehimpulsoperator und

• x
  →

  ^

  {\displaystyle {\hat {\vec {x}}}}

  der Ortsoperator sind, sowie

• Z

  {\displaystyle Z}

  die Kernladungszahl,

• α

  {\displaystyle \alpha }

  die Feinstrukturkonstante,

• ℏ

  {\displaystyle \hbar }

  die reduzierte Planck-Konstante,

• c

  {\displaystyle c}

  die Lichtgeschwindigkeit und

• m

  {\displaystyle m}

  die Masse des Elektrons sind.

Insbesondere ist in der Quantenmechanik

L
→

^

×

p
→

^

≠
−

p
→

^

×

L
→

^

{\displaystyle \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

H
^

{\displaystyle {\hat {H}}}

ist

H
^

=

p
→

^

2

2
m

−
Z
α
ℏ
c

1

|

x
→

^

|

{\displaystyle {\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

[

H
^

,

A
^

i

]
=
0

{\displaystyle [{\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

d

d

t

A
→

^

=

∂

∂
t

A
→

^

+

i

ℏ

[

H
^

,

A
→

^

]
=
0

{\displaystyle {\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

[

A
^

i

,

A
^

j

]
=
−
2

i

ℏ

H
^

ε

i
j
k

L
^

k

m

{\displaystyle [{\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

[

A
^

i

,

L
^

j

]
=

i

ℏ
c

ε

i
j
k

L
^

k

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

mit dem Levi-Civita-Symbol

ε

{\displaystyle \varepsilon }

. Insbesondere sind

[

A
→

^

2

,

A

j

]
=
[

A
→

^

2

,

L

j

]
=
[

A
→

^

2

,

L
→

^

2

]
=
[

L
→

^

2

,

A
^

j

]
=
[

L
→

^

2

,

L
^

j

]
=
0

{\displaystyle [{\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

H
^

,

L
→

^

2

,

L
^

3

,

A
→

^

2

{\displaystyle {\hat {H}},{\hat {\vec {L}}}^{2},{\hat {L}}_{3},{\hat {\vec {A}}}^{2}}

und

H
^

,

L
→

^

2

,

A
→

^

2

,

A
^

3

{\displaystyle {\hat {H}},{\hat {\vec {L}}}^{2},{\hat {\vec {A}}}^{2},{\hat {A}}_{3}}

.

Einzelnachweise

Figure 1.3 from Generalisations of the LaplaceRungeLenz vector in

Lecture 22 Laplace RungeLenz Vector YouTube

A free particle of energy E, momentum \vec{p} , angular momentum

The LaplaceRungeLenz Vector Short Notes YouTube

Evolution of the RungeLenz vector in an inspiralling binary. Initially

Solved 1) The RUNGE LENZ Vector operator is defined as 1

theorie106

Solved 8. (Bonus) LaplaceRungeLenz vector via Noether’s

Vector de RungeLenz

Solved The LaplaceRungeLenz (LRL) vector is an additional

L33 LAPLACE RUNGE LENZ VECTOR LRL Vector Classical Mechanics B

PPT Sect. 3.8 Motion in Time, Kepler Problem PowerPoint Presentation

The LaplaceRungeLenz Vector Classical Mechanics Homework DocsLib

Figure 4.1 from An Analogue of the LaplaceRungeLenz Vector for

Laplace Runge Lenz Vector YouTube

LaplaceRungeLenz vector Wikipedia

(PDF) LaplaceRungeLenz vector for arbitrary spin

LaplaceRungeLenz vector Context, History of rediscovery, Definition

[PhO]LaplaceRungeLenz vector YouTube

Practical Notes Of The Laplacerungelenz Vector 29. April 2013

Here Pauli refers to the vector as “previously utilized by Lenz” [1]. Since then, it has become popular to refer to the vector as the Runge-Lenz vector, as it is named in Shankar [7] and other popular quantum contexts.. The vector A lies in the plane of the orbit, points to the pericenter, and has magnitude A = m e, where e is the eccentricity of the orbit.