Symmetrien und Erhaltungsgrößen

Kontinuierliche Symmetrien und Erhaltungssätze

Betrachte kontinuierliche Transformationen, unter denen das physikalische System invariant ist.

In diesem Fall gibt es zu jeder kontinuierlichen Invarianz gegen infinitesimale Transformationen eine Erhaltungsgröße I ( Integral der Bewegung oder auch Konstante der Bewegung), das heißt, in diesem Fall gilt:

$\frac{dI}{dt}=0$ entlang der Bahn der angenommenen Bewegung ( längs der Bahn).

Dies ist die allgemeine Aussage des Theorems von Emmy Noether

Das Noether Theorem

Voraussetzung: Autonomes, das heißt, nicht explizit zeitabhängiges System mit f Freiheitsgraden und einer Lagrangefunktion

$L(q_1,...,{\dot{q}}_1,...,t)$

Theorem ( E.Noether, 1882-1935)

Die Lagrangefunktion $L(q_1,...,{\dot{q}}_1,...,t)$ eines autonomen Systems sei unter der Transformation

$\overline{q}\to h^s(\overline{q})$ invariant. Dabei ist s ein eindimensionaler Parameter und $h^{s=0}(\overline{q})=\overline{q}$ die Identität.

Dann gibt es ein Integral der Bewegung

$I(\overline{q},\dot{\overline{q}})=\sum _{i=1}^f\frac{\partial L}{\partial {\dot{q}}_i}{(\frac{d}{ds}{h}^s(q_i))}_{s=0}$

Beweis:

Sei $\overline{q}=\overline{q}(t)$ eine Lösung der Lagrangegleichung. Dann ist auch $\overline{q}(s,t):=h^s(\overline{q},t)$ Lösung, das heißt:

$\frac{d}{dt}\frac{\partial L(\overline{q}(s,t),\dot{\overline{q}}(s,t))}{\partial {\dot{q}}_i}=\frac{\partial L(\overline{q}(s,t),\dot{\overline{q}}(s,t))}{\partial q_i}$

Invarianz der Lagrangefunktion für beliebige s:

$\begin{array}{rl}& \frac{d}{ds}L(\overline{q}(s,t),\dot{\overline{q}}(s,t))=\sum _{i=1}^f(\frac{\partial L}{\partial q_i}\left(\frac{d{q}_i}{ds}\right)+\frac{\partial L}{\partial {\dot{q}}_i}{\left(\frac{d{\dot{q}}_i}{ds}\right)}_{})=0\\ & \Rightarrow \frac{d}{dt}I(\overline{q},\dot{\overline{q}})=\sum _{i=1}^f\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}_i}{(\frac{d}{ds}{h}^s(q_i))}_{s=0}\right)=\sum _{i=1}^f(\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}\left(\frac{d{q}_i}{ds}\right)+\frac{\partial L}{\partial {\dot{q}}_i}\frac{d}{dt}{\left(\frac{d{q}_i}{ds}\right)}_{})\\ \end{array}$

Mit

$\begin{array}{rl}& \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}=\frac{\partial L}{\partial q_i}\\ & \frac{d}{dt}\left(\frac{d{q}_i}{ds}\right)=\left(\frac{d{\dot{q}}_i}{ds}\right)\\ \end{array}$

und mit Hilfe von

$\frac{d}{ds}L(\overline{q}(s,t),\dot{\overline{q}}(s,t))=\sum _{i=1}^f(\frac{\partial L}{\partial q_i}\left(\frac{d{q}_i}{ds}\right)+\frac{\partial L}{\partial {\dot{q}}_i}{\left(\frac{d{\dot{q}}_i}{ds}\right)}_{})=0$

folgt dann:

$\frac{d}{dt}I(\overline{q},\dot{\overline{q}})=\frac{d}{ds}L=0$

Räumliche Translationsinvarianz

Seien die Kräfte konservativ und seien keine Randbedingungen:

$L=\frac{1}{2}\sum _{i=1}^N{m}_i{{\dot{\overline{r}}}_i}^2-V({\overline{r}}_1,...,{\overline{r}}_N)$

Eine Translation in Richtung x ist damit eine Translation der Form:

$h^s:{\overline{r}}_i{\to }_{}{\overline{r}}_i+s{\overline{e}}_{x}i=1,..,N$

Der Parameter s ist dabei beliebig.

Die Translationsinvarianz entlang der x- Achse bewirkt nun:

$\begin{array}{rl}& L(h^s({\overline{r}}_i),{\dot{\overline{r}}}_i)=\frac{1}{2}\sum _{i=1}^N{m}_i{{\dot{\overline{r}}}_i}^2-V({\overline{r}}_1+s{\overline{e}}_x,...,{\overline{r}}_N+s{\overline{e}}_x)=L({\overline{r}}_i,{\dot{\overline{r}}}_i)Forderung!\\ & \frac{dL}{ds}=-\sum _{i=1}^N({\mathrm{\nabla }}_{ri}\cdot {\overline{e}}_x)V=-\sum _{i=1}^N\frac{\partial }{\partial x_i}V=0Forderung!\\ \end{array}$

Das bedeutet aber: es darf keine äußere Kraft in x- Richtung geben !

Für die Transformation gilt:

$h^s({\overline{r}}_i)={\overline{r}}_i+s{\overline{e}}_{x}i=1,..,N$

$h^{s=0}({\overline{r}}_i)={\overline{r}}_i$ (Identität)

$\frac{d}{ds}{h}^s({\overline{r}}_i)={\overline{e}}_x$

Für unser Integral der Bewegung gilt jedoch:

$I=\sum _{i=1}^N{\mathrm{\nabla }}_{\dot{r}i}L\frac{d{h}^s}{ds}=\sum _i{m}_i{\dot{\overline{r}}}_i\cdot {\overline{e}}_x=\sum _i{m}_i{\dot{x}}_i=P_x$

Fazit: die Translationsinvarianz in x- Richtung bestimmt die Erhaltung der x-Komponente des Gesamtimpulses.

Dieser Zusammenhang ist leicht für die anderen Komponenten zu zeigen.

Dies kann auch umgekehrt betrachtet werden:

Wähle q1=s als verallgemeinerte Koordinate:

Nun gilt die Transformation:

${\overline{r}}_i={\overline{r}}_i(q_1,...,q_f,t)=q_1{\overline{e}}_x+\mathrm{\Delta }{\overline{r}}_i(q_1,...,q_f,t)$

mit

$q_1{\overline{e}}_x$

als Schwerpunktskoordinate und

$\mathrm{\Delta }{\overline{r}}_i(q_1,...,q_f,t)$

als Relativpositionen.

Es folgt:

$\frac{\partial }{\partial q_1}{\overline{r}}_i={\overline{e}}_x$

$\frac{\partial }{\partial {\dot{q}}_1}{\dot{\overline{r}}}_i=\frac{\partial }{\partial q_1}{\overline{r}}_i={\overline{e}}_x$

wegen  

${\dot{\overline{r}}}_i=\sum _k^{}\frac{\partial }{\partial q_k}{\overline{r}}_i{\dot{q}}_k+\frac{\partial }{\partial t}{\overline{r}}_i$

Invarianz Erhaltungssatz

${\frac{\partial L}{\partial q_1}}_{}=0\iff \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_1}=0$

 äquivalent zum Erhaltungssatz	

$\frac{\partial L}{\partial {\dot{q}}_1}=const$

Allgemein heißt $\frac{\partial L}{\partial {\dot{q}}_j}=p_j$ der zur Koordinate qj konjugierte verallgemeinerte Impuls.

Falls gilt dass ${\frac{\partial L}{\partial q_1}}_{}=0\iff \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_1}=0$ , wenn also die Lagrangefunktion invariant gegen q1- Änderungen ist, dann nennt man q1 eine zyklische Koordinate. der zu q1 konjugierte Impuls ist in diesem Fall eine Erhaltungsgröße .

Hier:

$\begin{array}{rl}& p_1=\frac{\partial L}{\partial {\dot{q}}_1}=\frac{\partial }{\partial {\dot{q}}_1}(T-V)=\frac{\partial T}{\partial {\dot{q}}_1}=\frac{\partial }{\partial {\dot{q}}_1}\left(\sum _i\frac{1}{2}{m}_i{{\dot{\overline{r}}}_i}^2\right)=\sum _i{m}_i{\dot{\overline{r}}}_i\frac{\partial }{\partial {\dot{q}}_1}{\dot{\overline{r}}}_i\\ & mit\frac{\partial }{\partial {\dot{q}}_1}{\dot{\overline{r}}}_i={\overline{e}}_x\\ & p_1=\sum _i{m}_i{\dot{\overline{r}}}_i{\overline{e}}_x=P_x\\ \end{array}$

Verallgemeinerung auf Nichtkonservative Kräfte

$\frac{d}{dt}\frac{\partial T}{\partial {\dot{q}}_1}-\frac{\partial T}{\partial q_1}=Q_1=\sum _i{\overline{X}}_i\frac{\partial }{\partial q_1}{\overline{r}}_i={\overline{e}}_x\sum _i{\overline{X}}_i$

Xi kennzeichnet dabei die Kraft. Nun steht rechts also die resultierende Kraft in x- Richtung. Existiert keine resultierende Kraft in x- Richtung ( Translationsinvarianz in x- Richtung), so gilt:

$\frac{d}{dt}\frac{\partial T}{\partial {\dot{q}}_1}-\frac{\partial T}{\partial q_1}=Q_1=\sum _i{\overline{X}}_i\frac{\partial }{\partial q_1}{\overline{r}}_i={\overline{e}}_x\sum _i{\overline{X}}_i=0$

Invarianz sagt

$\frac{\partial T}{\partial q_1}=Q_1=0\Rightarrow \frac{d}{dt}\frac{\partial T}{\partial {\dot{q}}_1}=0\iff P_x=\frac{\partial T}{\partial {\dot{q}}_1}=const$

Nebenbedingung für das fehlen konservativer Kräfte ( Falls Q1 konservative Kraft ist)

$Q_1=0\Rightarrow \frac{\partial }{\partial q_1}V({\overline{r}}_1+q_1{\overline{e}}_x,...,{\overline{r}}_N+q_1{\overline{e}}_x)=\sum _i{\mathrm{\nabla }}_{ri}V\frac{\partial }{\partial q_1}\left(q_1{\overline{e}}_x\right)={\overline{e}}_x\sum _i{\mathrm{\nabla }}_{ri}V=-{\overline{e}}_x\sum _i{\overline{X}}_i=0$

Beispiel: ein Teilchen im Potenzial V=V(y,z)

Das Potenzial hänge nicht von x ab: ${\frac{\partial L}{\partial x}}_{}=0$

Daraus folgt: ${\frac{\partial L}{\partial \dot{x}}}_{}=m\dot{x}=P_x=const$

In diesem Fall existiert ein Integral der Bewegung:

$I(\overline{r},\dot{\overline{r}})=\frac{\partial L}{\partial \dot{\overline{r}}}\cdot {\frac{d{h}^s}{ds}}_{}=\frac{\partial L}{\partial \dot{x}}=P_x=const$

wegen

$\begin{array}{rl}& \frac{\partial L}{\partial \dot{\overline{r}}}={\mathrm{\nabla }}_{\dot{r}}L\\ & {\frac{d{h}^s}{ds}}_{}={\overline{e}}_x\\ \end{array}$

Beispiel: 2 Teilchen mit innerer Paarwechselwirkung

$V({\overline{r}}_1,{\overline{r}}_2)=V({\overline{r}}_1-{\overline{r}}_2)$

 Das Potenzial kann auch anisotrop sein.

Es sollen keine äußeren Kräfte wirken, so dass das Potenzial unabhängig von den Schwerpunktskoordinaten wird.

Gleichzeitig soll Translationsinvarianz entlang x-, - und z- Richtung vorliegen:

$\begin{array}{rl}& L({\overline{r}}_1,{\overline{r}}_2,{\dot{\overline{r}}}_1,{\dot{\overline{r}}}_2)=\frac{m_1}{2}{{\dot{\overline{r}}}_1}^2+\frac{m_2}{2}{{\dot{\overline{r}}}_2}^2-V({\overline{r}}_1-{\overline{r}}_2)\\ & L(h^s\left({\overline{r}}_1\right),h^s\left({\overline{r}}_2\right),{\dot{\overline{r}}}_1,{\dot{\overline{r}}}_2)=\frac{m_1}{2}{{\dot{\overline{r}}}_1}^2+\frac{m_2}{2}{{\dot{\overline{r}}}_2}^2-V(({\overline{r}}_1-s{\overline{e}}_i)-({\overline{r}}_2-s{\overline{e}}_i))=L({\overline{r}}_1,{\overline{r}}_2,{\dot{\overline{r}}}_1,{\dot{\overline{r}}}_2)\\ \end{array}$ für alle i = x,y,z

Somit existieren gleich drei Integrale der Bewegung:

$\begin{array}{rl}& I_x=\frac{\partial L}{\partial {\dot{\overline{r}}}_1}{\overline{e}}_x+\frac{\partial L}{\partial {\dot{\overline{r}}}_2}{\overline{e}}_x=\frac{\partial L}{\partial {\dot{x}}_1}+\frac{\partial L}{\partial {\dot{x}}_2}=m_1{\dot{x}}_1+m_2{\dot{x}}_2=P_x=const\\ & I_y=\frac{\partial L}{\partial {\dot{\overline{r}}}_1}{\overline{e}}_y+\frac{\partial L}{\partial {\dot{\overline{r}}}_2}{\overline{e}}_y=\frac{\partial L}{\partial {\dot{y}}_1}+\frac{\partial L}{\partial {\dot{y}}_2}=m_1{\dot{y}}_1+m_2{\dot{y}}_2=P_y=const\\ & I_z=\frac{\partial L}{\partial {\dot{\overline{r}}}_1}{\overline{e}}_z+\frac{\partial L}{\partial {\dot{\overline{r}}}_2}{\overline{e}}_z=\frac{\partial L}{\partial {\dot{z}}_1}+\frac{\partial L}{\partial {\dot{z}}_2}=m_1{\dot{z}}_1+m_2{\dot{z}}_2=P_z=const\\ \end{array}$

Dies ist, aufgrund des Fehlens äußerer Kräfte, gerade der Schwerpunkts- Erhaltungssatz:

$M\dot{\overline{R}}={\overline{P}}_{}=const$

Mit den Schwerpunktskoordinaten

$\overline{R}:=\frac{1}{M}\sum _{i=1}^2{m}_i{\overline{r}}_i$

Und der Gesamtmasse

$M:=\sum _{i=1}^2{m}_i$

Räumliche Isotropie

Nebenbedingung: konservative Kräfte, keine Zwangsbedingungen

Es erfolgt eine Drehung des Bezugssystems um den Winkel $\varphi =s$ um die z- Achse.

An einer Skizze kann man sich schnell verdeutlichen:

$h^s:{\overline{r}}_i=(x_i,y_i,z_i)\to {\overline{r}}_i\stackrel{´}{}=(x{\stackrel{´}{}}_i,y{\stackrel{´}{}}_i,z{\stackrel{´}{}}_i)$

Dabei gilt:

$\begin{array}{rl}& x_i\stackrel{´}{}=x_i\cos s+y_i\sin s\\ & y_i\stackrel{´}{}=y_i\cos s-x_i\sin s\\ & z_i\stackrel{´}{}=z_i\\ \end{array}$

Rotationsinvarianz für die Drehung um die z- Achse:

Betrachten wir infinitesimale Transformationen ( Drehungen um die z- Achse mit kleinen Winkeln $\delta \varphi =\delta s$

$\left(\begin{array}{c}{x}_i\stackrel{´}{}\\ y_i\stackrel{´}{}\\ z_i\stackrel{´}{}\\ \end{array}\right)=\left(\begin{array}{ccc}\cos s& \sin s& 0\\ -\sin s& \cos s& 0\\ 0& 0& 1\\ \end{array}\right)\left(\begin{array}{c}{x}_i\\ y_i\\ z_i\\ \end{array}\right)\approx [\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\\ \end{array}\right)+\left(\begin{array}{ccc}0& s& 0\\ -s& 0& 0\\ 0& 0& 0\\ \end{array}\right)]\left(\begin{array}{c}{x}_i\\ y_i\\ z_i\\ \end{array}\right)$

Dabei gilt die rechtsseitige Taylorentwicklung für kleine Winkel. Wir schreiben

$\left(\begin{array}{ccc}0& s& 0\\ -s& 0& 0\\ 0& 0& 0\\ \end{array}\right)=-s{\overline{\overline{J}}}_z$

Mit ${\overline{\overline{J}}}_z$ als Erzeugende für infinitesimale Drehungen um die z- Achse.

Somit folgt:

$\left(\begin{array}{c}{x}_i\stackrel{´}{}\\ y_i\stackrel{´}{}\\ z_i\stackrel{´}{}\\ \end{array}\right)=\left(\begin{array}{c}{x}_i\\ y_i\\ z_i\\ \end{array}\right)+s\left(\begin{array}{c}{y}_i\\ -x_i\\ 0\\ \end{array}\right)=\left(\begin{array}{c}{x}_i\\ y_i\\ z_i\\ \end{array}\right)+s({\overline{r}}_i\times {\overline{e}}_z)$

Formal schreibt man:

${\overline{r}}_i\stackrel{´}{}=h^s({\overline{r}}_i)|{s=0}_{}+s{(\frac{d}{ds}{h}^s({\overline{r}}_i))}_{s=0}+O(s^2)$

mit ${(\frac{d}{ds}{h}^s({\overline{r}}_i))}_{s=0}={\overline{r}}_i\times {\overline{e}}_z$

Rotationsinvarianz der Lagrange-Funktion

$T=\frac{1}{2}\sum _i{m}_i{{\dot{\overline{r}}}_i}^2$ ist rotationsinvariant, da nur von $\left|{\dot{\overline{r}}}_i\right|$ abhängig und die Drehmatrix ändert die Abstände nicht.

( Drehungen sind orthogonale Transformationen).

${\left(\frac{\partial L}{\partial s}\right)}_{s=0}=-{\left(\frac{\partial V}{\partial s}\right)}_{s=0}=-\sum _{i=1}^N\left({\mathrm{\nabla }}_{ri\stackrel{´}{}}V\right){\left(\frac{d{\overline{r}}_i\stackrel{´}{}}{ds}\right)}_{s=0}=\sum _{i=1}^N{\overline{F}}_i({\overline{r}}_i\times {\overline{e}}_z)$

wegen:

$\begin{array}{rl}& \left({\mathrm{\nabla }}_{ri\stackrel{´}{}}V\right)=-{\overline{F}}_i\\ & {\left(\frac{d{\overline{r}}_i\stackrel{´}{}}{ds}\right)}_{s=0}={\left(\frac{d{h}^s}{ds}\right)}_{s=0}\\ \end{array}$

Als zyklische Permutation gilt dann jedoch:

${\left(\frac{\partial L}{\partial s}\right)}_{s=0}=-{\left(\frac{\partial V}{\partial s}\right)}_{s=0}={\overline{e}}_z\cdot \sum _i({\overline{F}}_i\times {\overline{r}}_i)=-{\overline{e}}_z\cdot \sum _i({\overline{r}}_i\times {\overline{F}}_i)$

Mit $\sum _i({\overline{r}}_i\times {\overline{F}}_i)$ als gesamtes Drehmoment und der Tatsache, dass die z-Komponente des äußeren resultierenden Drehmomentes verschwindet:

$-{\overline{e}}_z\cdot \sum _i({\overline{r}}_i\times {\overline{F}}_i)={\left(\frac{\partial L}{\partial s}\right)}_{s=0}=-{\left(\frac{\partial V}{\partial s}\right)}_{s=0}=0$

Interpretation nach dem Noetherschen Theorem

$I(\overline{r},\dot{\overline{r}})=\sum _{i=1}^N\frac{\partial L}{\partial {\dot{\overline{r}}}_i}\cdot {\left(\frac{d{h}^s}{ds}\right)}_{s=0}=\sum _i{m}_i{\dot{\overline{r}}}_i\cdot ({\overline{r}}_i\times {\overline{e}}_z)=-{\overline{e}}_z\sum _i({\overline{r}}_i\times m_i{\dot{\overline{r}}}_i)=-{\overline{e}}_z\overline{l}=-l_z$

Also: Rotationsinvarianz entspricht Drehimpulserhaltung

Andere Betrachtungsweise

Wähle $q_1=\varphi =s$ als verallgemeinerte Koordinate

Trafo: ${\overline{r}}_i={\overline{r}}_i(\varphi ,q_2,...,q_f,t)$

mit $\frac{\partial }{\partial q_1}{\overline{r}}_i={(\frac{d}{ds}{h}^s({\overline{r}}_i))}_{s=0}={\overline{r}}_i\times {\overline{e}}_z$

Für infinitesimale Drehung um z-Achse.

Invarianz Erhaltungssätze

${\frac{\partial L}{\partial q_1}}_{}=0\iff \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_1}=0$

 äquivalent zum Erhaltungssatz	

$p_1=\frac{\partial L}{\partial {\dot{q}}_1}=const$

Der Winkel ist also eine zyklische Variable.

Berechnet man den verallgemeinerten konjugierten Impuls zu $q_1=\varphi =s$ , so ergibt sich:

$p_1=\frac{\partial L}{\partial {\dot{q}}_1}=\frac{\partial L}{\partial {\dot{q}}_1}=\frac{1}{2}\sum _i{m}_i\frac{\partial }{\partial {\dot{q}}_1}{{\dot{\overline{r}}}_i}^2=\sum _i{m}_i{\dot{\overline{r}}}_i\frac{\partial }{\partial {\dot{q}}_1}{\dot{\overline{r}}}_i=\sum _i{m}_i{\dot{\overline{r}}}_i({\overline{r}}_i\times {\overline{e}}_z)=-{\overline{e}}_z\sum _i({\overline{r}}_i\times m_i{\dot{\overline{r}}}_i)=-l_z$

wegen

$\frac{\partial }{\partial {\dot{q}}_1}{\dot{\overline{r}}}_i=\frac{\partial }{\partial q_1}{\overline{r}}_{i}da{}_{}{\dot{\overline{r}}}_i=\sum _k\frac{\partial {\overline{r}}_i}{\partial q_k}{\dot{q}}_k+\frac{\partial {\overline{r}}_i}{\partial t}$

Es ergibt sich also wieder die z-Komponente des Drehimpulses als verallgemeinerter Impuls.

Nebenbedingung:

Wir betrachteten hier eine passive Drehung des Korodinatensystems. Die Aktive Drehung des Koordinatensystems ist jedoch äquivalent. Das bedeutet, wir drehen aktiv alle Massenpunkte mit $\tilde{\varphi }=-\varphi$ .

Dazu gehören dann die konjugierten Impulse +lz

Beispiel:

N Teilchen mit einer inneren Paarwechselwirkung, die nur vom Abstand abhängt:

$V({\overline{r}}_1,...,{\overline{r}}_N)=V(r_{12},...,r_{ij},...)$ mit $r_{ij}=|{\overline{r}}_i-{\overline{r}}_j|$

Rotationsinvarianz gegen Drehung um alle Achsen:

$\frac{\partial V(r_{12},...,r_{ij},...)}{\partial \varphi }=\sum _{i,j}\frac{\partial V}{\partial r_{ij}}\cdot \frac{\partial }{\partial \varphi }{r}_{ij}=0$ für beliebige Achsen, da

$\begin{array}{rl}& \frac{\partial }{\partial \varphi }{r}_{ij}=\frac{\partial }{\partial \varphi }{\left[({\overline{r}}_i-{\overline{r}}_j)({\overline{r}}_i-{\overline{r}}_j)\right]}^{1/2}=\frac{1}{r_{ij}}({\overline{r}}_i-{\overline{r}}_j)\frac{\partial }{\partial \varphi }({\overline{r}}_i-{\overline{r}}_j)=\frac{{\overline{r}}_i-{\overline{r}}_j}{r_{ij}}(\frac{\partial }{\partial \varphi }{\overline{r}}_i-\frac{\partial }{\partial \varphi }{\overline{r}}_j)\\ & \frac{\partial }{\partial \varphi }{\overline{r}}_i={\overline{r}}_i\times {\overline{e}}_k\\ & \Rightarrow \frac{{\overline{r}}_i-{\overline{r}}_j}{r_{ij}}(\frac{\partial }{\partial \varphi }{\overline{r}}_i-\frac{\partial }{\partial \varphi }{\overline{r}}_j)=\frac{{\overline{r}}_i-{\overline{r}}_j}{r_{ij}}[({\overline{r}}_i-{\overline{r}}_j)\times {\overline{e}}_k]=\frac{1}{r_{ij}}{\overline{e}}_k[({\overline{r}}_i-{\overline{r}}_j)\times ({\overline{r}}_i-{\overline{r}}_j)]=0\\ \end{array}$

Also ist der resultierende Drehimpuls $\overline{l}$ eine Erhaltungsgröße

Erzeugende der infinitesimalen Drehung um z-Achse

Die infinitesimale Drehung läßt sich schreiben als:

${\overline{r}}_i\stackrel{´}{}=h^s({\overline{r}}_i)=(\overline{\overline{1}}-s{\overline{\overline{J}}}_z){\overline{r}}_i$

Mit der Erzeugenden ${\overline{\overline{J}}}_z=\left(\begin{array}{ccc}0& -1& 0\\ 1& 0& 0\\ 0& 0& 0\\ \end{array}\right)$

Bei einer Drehung um den endlichen Winkel $\varphi$ gilt:

${\overline{r}}_i\stackrel{´}{}={\overline{\overline{R}}}_z(\varphi ){\overline{r}}_i=\left(\begin{array}{ccc}\cos \varphi & \sin \varphi & 0\\ -\sin \varphi & \cos \varphi & 0\\ 0& 0& 1\\ \end{array}\right){\overline{r}}_i$

Es gilt:

${\overline{\overline{R}}}_z(\varphi )=\exp (-{\overline{\overline{J}}}_z\varphi )$

mit Definition

$\exp (-{\overline{\overline{J}}}_z\varphi ):=\overline{\overline{1}}+(-{\overline{\overline{J}}}_z\varphi )+\frac{1}{2}{(-{\overline{\overline{J}}}_z\varphi )}^2+...+\frac{1}{k!}{(-{\overline{\overline{J}}}_z\varphi )}^k$

Beweis:

Für

$\begin{array}{rl}& \overline{\overline{M}}=\left(\begin{array}{cc}0& -1\\ 1& 0\\ \end{array}\right)\Rightarrow {\overline{\overline{M}}}^2=-\overline{\overline{1}},{\overline{\overline{M}}}^3=-\overline{\overline{M}},{\overline{\overline{M}}}^4=\overline{\overline{1}}\\ & {\overline{\overline{M}}}^{2n}={(-1)}^n\overline{\overline{1}}\\ & {\overline{\overline{M}}}^{(2n+1)}={(-1)}^n\overline{\overline{M}}\\ \end{array}$

Mit Hilfe der Taylorreihen für Sinus und Cosinus folgt dann:

$\begin{array}{rl}& \left(\begin{array}{cc}\cos \varphi & \sin \varphi \\ -\sin \varphi & \cos \varphi \\ \end{array}\right)=\overline{\overline{1}}\sum _{n=0}^{\mathrm{\infty }}\frac{{(-1)}^n}{\left(2n\right)!}{\varphi }^{2n}-\overline{\overline{M}}\sum _{n=0}^{\mathrm{\infty }}\frac{{(-1)}^n}{(2n+1)!}{\varphi }^{2n+1}\\ & =\sum _{n=0}^{\mathrm{\infty }}\frac{1}{\left(2n\right)!}{\overline{\overline{M}}}^{2n}{\varphi }^{2n}-\overline{\overline{M}}\sum _{n=0}^{\mathrm{\infty }}\frac{1}{(2n+1)!}{\overline{\overline{M}}}^{2n+1}{\varphi }^{2n+1}\\ & =\exp (-\overline{\overline{M}}\varphi )\\ \end{array}$

Analog behandelbar ist die Drehung um die x-Achse

Erzeugende:

${\overline{\overline{J}}}_x=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -1\\ 0& 1& 0\\ \end{array}\right)$

Hier gewinnen wir die Drehmatrix:

${\overline{\overline{R}}}_x(\varphi )=\exp (-{\overline{\overline{J}}}_x\varphi )$

Bei der y- Achse gilt:

Erzeugende:

${\overline{\overline{J}}}_y=\left(\begin{array}{ccc}0& 0& 1\\ 0& 0& 0\\ -1& 0& 0\\ \end{array}\right)$

Hier gewinnen wir die Drehmatrix:

${\overline{\overline{R}}}_y(\varphi )=\exp (-{\overline{\overline{J}}}_y\varphi )$

Beliebige Drehungen um den Winkel $\varphi$ mit der Drehachse $\overline{n}$

${\overline{\overline{R}}}_{}(\overline{\varphi })=\exp (-\varphi \sum _{i=1}^3{n}_i{\overline{\overline{J}}}_i)$

mit $\overline{\varphi }:=\varphi \overline{n}$

Die Drehmatrizen ${\overline{\overline{R}}}_{}(\overline{\varphi })=\exp (-\varphi \sum _{i=1}^3{n}_i{\overline{\overline{J}}}_i)$ bilden nun eine 3- parametrige $({\varphi }_1,{\varphi }_2,{\varphi }_3)$ , stetige, diffbare $\left(in\varphi \right)$

und orthogonale Gruppe.

Eine solche Gruppe heißt Lie- Gruppe oder kontinuierliche Gruppe in drei reellen Dimensionen

SO(3)

$SO\left(3\right)=\{\overline{\overline{R}}:R^3\to R^{3}linear|{\overline{\overline{R}}}^t\overline{\overline{R}}=1|\det \overline{\overline{R}}=1\}$

Mit ${\overline{\overline{R}}}^t\overline{\overline{R}}=1$ als Orthogonalitätsbedingung, so dass $\left|\overline{r}\stackrel{´}{}\right|=\left|\overline{r}\right|$

und 

$\det \overline{\overline{R}}=1$ zum Ausschluß von Raumspiegelungen.

Die Erzeugenden ${\overline{\overline{J}}}_i$ der Drehgruppe bilden eine Lie- Algebra mit dem Lieschen Produkt (=Kommutator):

$[{\overline{\overline{J}}}_i,{\overline{\overline{J}}}_k]={\overline{\overline{J}}}_i{\overline{\overline{J}}}_k-{\overline{\overline{J}}}_k{\overline{\overline{J}}}_i$ i,k=x,y,z

Dabei vertauschen 2 Drehungen um unterschiedliche Achsen nicht. Das bedeutet, das Ergebnis hängt von der Reihenfolge ab !:

$\begin{array}{rl}& [{\overline{\overline{J}}}_x,{\overline{\overline{J}}}_y]={\overline{\overline{J}}}_z\\ & [{\overline{\overline{J}}}_z,{\overline{\overline{J}}}_x]={\overline{\overline{J}}}_y\\ & [{\overline{\overline{J}}}_y,{\overline{\overline{J}}}_z]={\overline{\overline{J}}}_x\\ \end{array}$

-> zyklische Permutation des Lieschen Produktes

Zeitliche Translationsinvarianz

Die Zeit spielt in der klassischen Mechanik im Ggstz zur relativistischen Mechanik gegenüber dem Ort eine Sonderrolle.

Deshalb ist eine direkte Anwendung des Noether- Theorems nicht moeglich.

Zeitliche Translationsinvarianz ist erfüllt, falls:

1. die Zwangsbedingungen die Zeit t nicht explizit enthalten:

$\begin{array}{rl}& {\overline{r}}_i={\overline{r}}_i(q_1,...,q_f)\\ & \frac{\partial }{\partial t}{\overline{r}}_i=0\Rightarrow {\dot{\overline{r}}}_i=\sum _j^{}\frac{\partial }{\partial q_j}{\overline{r}}_i{{\dot{q}}_j}_{}\\ \end{array}$

Dabei ist $\frac{\partial }{\partial q_j}{\overline{r}}_i$ Funktion von q1...qf

$\frac{\partial }{\partial t}L=0$

1. Nebenbedingung: Aus der Existenz eines Potenzials der eingeprägten Kräfte folgt NICHT automatisch die Erhaltung der Energie, da die Zwangsbedingungen die Zeit enthalten könnten.

Wenn die Zwangsbedingungen die Zeit enthalten, so ist die Energie nicht enthalten.

${\overline{r}}_i={\overline{r}}_i(q_1,...,q_f,t)$

Kinetische Energie:

$T=\frac{1}{2}\sum _i^{}{m}_i{{\dot{\overline{r}}}_i}^2=\frac{1}{2}\sum _{j,k}^{}{T}_{jk}{\dot{q}}_j{\dot{q}}_k$

Mit

$T_{jk}=\sum _{i=1}^N{m}_i\left(\frac{\partial {\overline{r}}_i}{\partial q_j}\right)\left(\frac{\partial {\overline{r}}_i}{\partial q_k}\right)$ ist abhängig von den q1...qf im Gegensatz zum Fall der kleinen Schwingungen, der eingangs behandelt wurde.

T ist eine homogene quadratische Funktion der ${\dot{q}}_1...{\dot{q}}_f$

Also $T(\lambda {\dot{q}}_1,...,\lambda {\dot{q}}_f)={\lambda }^{2}T({\dot{q}}_1,...,{\dot{q}}_f)$

Nach $\lambda$ wird partiell abgelitten, dann wird $\lambda =1$ gesetzt.

$\begin{array}{rl}& \sum _{k=1}^N\left(\frac{\partial T}{\partial \left(\lambda {\dot{q}}_k\right)}\right)\left(\frac{\partial \left(\lambda {\dot{q}}_k\right)}{\partial \lambda }\right)|{\lambda =1}_{}=2\lambda T|{\lambda =1}_{}\iff \sum _{k=1}^N\left(\frac{\partial T}{\partial \left({\dot{q}}_k\right)}\right){\dot{q}}_k=2T\\ & \left(\frac{\partial \left(\lambda {\dot{q}}_k\right)}{\partial \lambda }\right)={\dot{q}}_k\\ \end{array}$

Obere Äquivalenz ist der sogenannte Eulersche Satz

Da V unabhängig von ${\dot{q}}_1...{\dot{q}}_f$ gilt auch:

$\sum _{k=1}^N\left(\frac{\partial L}{\partial \left({\dot{q}}_k\right)}\right){\dot{q}}_k=2T$

Zur totalen Zeitableitung von L:

$\begin{array}{rl}& \frac{dL}{dt}=\sum _k^{}(\frac{\partial L}{\partial {\dot{q}}_k}{\ddot{q}}_k+\frac{\partial L}{\partial q_k}{\dot{q}}_k)+\frac{\partial L}{\partial t}\\ & \frac{\partial L}{\partial q_k}=\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_k}und\frac{\partial L}{\partial t}=0wegen2.(oben)\\ \end{array}$

Somit:

$\frac{dL}{dt}=\sum _k^{}(\frac{\partial L}{\partial {\dot{q}}_k}{\ddot{q}}_k+\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_k}{\dot{q}}_k)=\frac{d}{dt}\sum _k^{}\frac{\partial L}{\partial {\dot{q}}_k}{\dot{q}}_k=2\frac{dT}{dt}$ wegen $\sum _{k=1}^N\left(\frac{\partial L}{\partial \left({\dot{q}}_k\right)}\right){\dot{q}}_k=2T$

Somit:

$0=\frac{d}{dt}(2T-L)=\frac{d}{dt}(T+V)\Rightarrow T+V=konst$

Zeitranslationsinvarianz bedingt also Energieerhaltung !

Oder: Skleronome Zwangsbedingungen: $\frac{\partial L}{\partial t}=0$ bedingen: E=T+V=constant

Nebenbemerkung

Die Aussage folgt auch aus dem Noether-Theorem, wenn man noch den folgenden Trick anwendet: (Scheck, Aufgabe 2.17)

Mache t zu einer q-artigen Variablen durch eine parametrisierte Darstellung: $q_k=q_k(\tau ),t=t(\tau )$

Als Lagrangefunktion muss man sich definieren:

$\overline{L}(q_k,t,\frac{d{q}_k}{d\tau },\frac{d{t}_{}}{d\tau }):=L(q_k,\frac{1}{\left({}^{dt}╱{}_{d\tau }\right)}\frac{d{q}_k}{d\tau },t,\frac{dt}{d\tau })$

soll invariant unter Zeittranslationen sein:

$h^s(\overline{q},t)=(\overline{q},t+s)$

Dann gilt:

1. Hamiltonsches Prinzip auf

$\overline{L}$

angewandt:

$0=\delta \underset{\tau 1}{\overset{\tau 2}{\int }}\overline{L}d\tau =\delta \underset{t1}{\overset{t2}{\int }}Ldt\iff \frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_k}-\frac{\partial L}{\partial q_k}=0$

2. Noethersches Theorem für $\overline{L}$

Integral der Bewegung I:

$\begin{array}{rl}& I=\sum _{i=1}^{f+1}\frac{\partial L}{\partial {\dot{q}}_i}{(\frac{d}{ds}{h}^s(q_1,...,q_{f+1}))}_{s=0}=\frac{\partial \overline{L}}{\partial {\dot{q}}_{f+1}}\\ & mit(\frac{d}{ds}{h}^s(q_1,...,q_{f+1}))=(0,...,0,1)fNullen,1anStellef+1mit{q}_{f+1}=t\\ \end{array}$

$\begin{array}{rl}& I=\frac{\partial \overline{L}}{\partial {\dot{q}}_{f+1}}=\frac{\partial \overline{L}}{\partial \left(\frac{dt}{d\tau }\right)}=L+\sum _{k=1}^f\frac{\partial L}{\partial {\dot{q}}_k}(-\frac{1}{{\left(\frac{dt}{d\tau }\right)}^2})\frac{d{q}_k}{d\tau }\frac{dt}{d\tau }\\ & =L-\sum _{k=1}^f\left(\frac{\partial L}{\partial \left({\dot{q}}_k\right)}\right){\dot{q}}_k=T-V-2T=-(T-V)\\ \end{array}$

Also Erhaltung der Energie durch zeitliche Translationsinvarianz

Das Zweikörperproblem

Hier werden die Erhaltungssätze zur Lösung der Bewegungsgleichung verwendet.

Idee:

f Freiheitsgrade -> f Differenzialgleichungen 2. Ordnung

• 2f Integrationskonstanten nötig ! ( jeweils zweifaches Integrieren). ( Anfangsbedingungen).
• Also existieren auch 2f Integrale der Bewegung

Falls alle 2f Integrale der Bewegung bekannt wären:

$I_r(q_1,...,q_f,{\dot{q}}_1,...{\dot{q}}_f)=c_{r}r=1,...2f$

So wäre das Problem vollständig gelöst:

$q_k=q_k(c_1,...,c_{2f},t)k=1,...,f$

Also ist es das Ziel, möglichst viele Integrale der Bewegung zu finden.

Beispiel: Zweikörperproblem

2 Massen, m1 und m2 unter dem Einfluss Ihrer inneren Wechselwirkung: V(|r1-r2|) ( Zentralpotenzial).

Beispiel: Sonne / Erde unter Gravitationswechselwirkung

Zahl der Freiheitsgrade: f=6

Also: es muessten 12 Integrale der Bewegung existieren

Erhaltungssätze

1. V(|r1-r2|) ist translationsinvariant.

Somit ist der Impuls: $\overline{P}={\overline{p}}_1+{\overline{p}}_2$ =konstant

Der Schwerpunkt: $\overline{R}=\frac{1}{M}\overline{P}t+{\overline{R}}_0$ bewegt sich gleichförmig und geradlinig.

Dies folgt aus: $M\dot{\overline{R}}=\overline{P}=const$

M:=m1 + m2

Somit sind 6 Integrationskonstanten gefunden: $\overline{P},\overline{R}$

1. V(|r1-r2|) ist rotationsinvariant:

Damit ist der Drehimpuls $\overline{l}=m_1{\overline{r}}_1\times {\overline{v}}_1+m_2{\overline{r}}_2\times {\overline{v}}_2=const$

Es sind drei weitere Integrationskonstanten $\overline{l}$ gefunden.

1. Die zeitliche Translationsinvarianz bei konservativer Kraft:

$E=\frac{1}{2}{m}_1{{\overline{v}}_1}^2+\frac{1}{2}{m}_2{{\overline{v}}_2}^2+V\left(|{\overline{r}}_1-{\overline{r}}_2|\right)=const$

Eine Integrationskonstante E

Insgesamt sind 10 Integrale der Bewegung gefunden. Es bleiben nur 2 Integrationskonstanten, nämlich der Nullpunkt der Zeit- und Winkelskala. Diese ergeben sich aus den ANfangsbedingungen.

Impuls- und Drehimpulserhaltung

Lagrange- Formulierung:

$L=T-V=\frac{1}{2}{m}_1{{\overline{v}}_1}^2+\frac{1}{2}{m}_2{{\overline{v}}_2}^2-V\left(|{\overline{r}}_1-{\overline{r}}_2|\right)$

Verallgeminerte Koordinaten: Schwerpunktskoordinaten:

$\left(\begin{array}{c}{q}_1\\ q_2\\ q_3\\ \end{array}\right):=\overline{R}=\frac{1}{M}(m_1{\overline{r}}_1+m_2{\overline{r}}_2)$ Schwerpunktskoordinate

$\left(\begin{array}{c}{q}_4\\ q_5\\ q_6\\ \end{array}\right):=\overline{r}={\overline{r}}_1-{\overline{r}}_2$ Relativkoordinate

Die Umkehrung liefert dann die gesuchten Größen:

$\begin{array}{rl}& {\overline{r}}_1=\overline{R}+\frac{m_2}{M}{\overline{r}}_{}{\overline{r}}_2=\overline{R}-\frac{m_1}{M}\overline{r}\\ & {\dot{\overline{r}}}_1=\dot{\overline{R}}+\frac{m_2}{M}{\dot{\overline{r}}}_{}{\dot{\overline{r}}}_2=\dot{\overline{R}}-\frac{m_1}{M}{\dot{\overline{r}}}_{}\\ & L=\frac{M}{2}{\dot{\overline{R}}}^2+\frac{1}{2}m{\dot{\overline{r}}}^2-V(r)\\ \end{array}$

Dabei bezeichnet

$r:=\left|\overline{r}\right|$

den Abstand und

$m=\frac{m_1{m}_2}{m_1+m_2}$ die relative Masse

$L=\frac{M}{2}{\dot{\overline{R}}}^2+\frac{1}{2}m{\dot{\overline{r}}}^2-V(r)$

$\overline{R}$ ist zyklische Koordinate: $\frac{\partial L}{\partial R_k}=0\Rightarrow \frac{\partial L}{\partial {\dot{R}}_k}=M{\dot{R}}_k=P_k=const$ mit k= x,y,z

$\Rightarrow \overline{R}=\frac{1}{M}\overline{P}t+{\overline{R}}_0$

Verwende das Schwerpunktsystem als Inertialsystem:

o.B.d.A: $\overline{R}=\dot{\overline{R}}=0$

Damit ergibt sich die vereinfachte Lagrangegleichung

$L=\frac{1}{2}m{\dot{\overline{r}}}^2-V(r)$

mit:

$\begin{array}{rl}& {\overline{r}}_1=+\frac{m_2}{M}{\overline{r}}_{}{\overline{r}}_2=-\frac{m_1}{M}\overline{r}\\ & {\dot{\overline{r}}}_1=+\frac{m_2}{M}{\dot{\overline{r}}}_{}{\dot{\overline{r}}}_2=-\frac{m_1}{M}{\dot{\overline{r}}}_{}\\ \end{array}$

Der Drehimpuls berechnet sich gemäß:

$\overline{l}=m_1{\overline{r}}_1\times {\overline{v}}_1+m_2{\overline{r}}_2\times {\overline{v}}_2=(\frac{m_1{m_2}^2}{M^2}+\frac{m_2{m_1}^2}{M^2})\overline{r}\times \dot{\overline{r}}=m\overline{r}\times \dot{\overline{r}}=const$ (Rotationsinvarianz)

Somit folgt aber auch (zyklische Vertauschbarkeit):

$\overline{l}\cdot \overline{r}=\overline{l}\cdot \dot{\overline{r}}=0$

Beide, Radiusvektor und Geschwindigkeitsvektor $\overline{r},\dot{\overline{r}}$ liegen in der Ebene senkrecht zu $\overline{l}$ ( Im Schwerpunktsystem).

Übergang zu Polarkoordinaten. Wir legen das Koordinatensystem so, dass der Drehimpuls parallel zur z- Achse zeigt:

$\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 \\ \end{array}$

Somit:

${\dot{\overline{r}}}^2={\dot{x}}^2+{\dot{y}}^2=...={\dot{r}}^2+r^2{\dot{\varphi }}^2$

Nun wählen wir neue verallgemeinerte Koordinaten statt x,y : $(r,\varphi )$

$L=\frac{1}{2}m({\dot{r}}^2+r^2{\dot{\varphi }}^2)-V(r)$

$\varphi$ ist zyklische Koordinate: $\frac{\partial L}{\partial \varphi }=0\Rightarrow \frac{\partial L}{\partial \dot{\varphi }}=m{r}^2\dot{\varphi }=l=const$

Hier: l = lz, da lx = ly =0

Also: $m{r}^2\dot{\varphi }=l_z=m(x\dot{y}-y\dot{x})=const$

Flächensatz: 2. keplersches Gesetz

Geometrische Interpretation von $m{r}^2\dot{\varphi }=l_z=m(x\dot{y}-y\dot{x})=const$

Radiusvektor überstreicht in gleichen Zeiten gleiche Flächen.

Das heißt: Die Flächengeschwindigkei ist konstant:

Für die Fläche gilt:

$\delta F=\frac{1}{2}\left|\overline{r}\right|\cdot |\overline{r}+\delta \overline{r}|\sin \delta \varphi \approx \frac{1}{2}{r}^2\delta \varphi$

Dabei gilt die rechtsseitige Näherung für sehr kleine Änderungen in Radiusvektor und Winkel. Bleibt richtig für infinitesimale Betrachtung:

$\frac{d}{dt}F=\frac{1}{2}{r}^2\frac{d\varphi }{dt}=\frac{l}{2m}=const$

Energieerhaltung und Bahngleichung

Bestimmen wir die Lagranggleichung 2. Art für den radius r:

$\frac{d}{dt}\frac{\partial L}{\partial \dot{r}}-\frac{\partial L}{\partial r}=0$

$\begin{array}{rl}& \frac{\partial L}{\partial \dot{r}}=m\dot{r}\\ & \frac{\partial L}{\partial r}=mr{\dot{\varphi }}^2-V\stackrel{´}{}(r)\\ \end{array}$

Somit gilt:

$m\ddot{r}-mr{\dot{\varphi }}^2+V\stackrel{´}{}(r)=0$

Mit der Zentrifugalkraft $mr{\dot{\varphi }}^2$

Die Zeitableitung des Winkels können wir eliminieren durch die Bewegungskonstante l:

$\dot{\varphi }=\frac{l}{m{r}^2}$

$m\ddot{r}-\frac{l^2}{m{r}^3}+V\stackrel{´}{}(r)=0$

1. Integral: Trick: Wir müssen die Gleichung auf zeitliche Änderung bringen. Zu diesem zweck multiplizieren wir alles mit

$\dot{r}$

$\begin{array}{rl}& m\ddot{r}\dot{r}-\frac{l^2}{m{r}^3}\dot{r}+\dot{r}V\stackrel{´}{}(r)=0\\ & m\ddot{r}\dot{r}=\frac{d}{dt}\left(\frac{m}{2}{\dot{r}}^2\right)\\ & \frac{l^2}{m{r}^3}\dot{r}=\frac{d}{dt}(-\frac{l^2}{2m{r}^2})\\ & \dot{r}V\stackrel{´}{}(r)=\frac{d}{dt}V(r)\\ \end{array}$

Somit können wir Integration über die zeit ausführen und es ergibt sich:

$\frac{m}{2}{\dot{r}}^2+\frac{l^2}{2m{r}^2}+V(r)=const=E$ Energieerhaltung mit $T=\frac{m}{2}({\dot{r}}^2+\frac{l^2}{m^2{r}^2})=\frac{m}{2}({\dot{r}}^2+r^2{\dot{\varphi }}^2)$

Andere Interpretation

Die Bewegung der beiden Körper ist ebenfalls als eindimensionale Bewegung in einem effektiven

Radialpotenzial

$\tilde{V}(r):=\frac{l^2}{2m{r}^2}+V(r)$

Dabei wird $\frac{l^2}{2m{r}^2}$ als Zentrifugalbarriere bezeichnet.

Es ergibt sich: $\frac{m}{2}{\dot{r}}^2+\tilde{V}(r)=const=E$

Somit:

$\dot{r}=\sqrt{\frac{2}{m}(E-\tilde{V}(r))}=\frac{dr}{dt}$

Integration liefert:

$\underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{\sqrt{\frac{2}{m}(E-\tilde{V}(r\stackrel{´}{}))}}=\underset{t_o}{\overset{t}{\int }}dt\stackrel{´}{}$

Es sind somit t( r) und r( t) berechenbar.

Der Winkel folgt dann aus:

$\dot{\varphi }=\frac{d\varphi }{dt}=\frac{l}{mr{(t)}^2}$ durch Einsetzen:

$\underset{{\varphi }_o}{\overset{\varphi }{\int }}d\varphi \stackrel{´}{}=\underset{t_o}{\overset{t}{\int }}\frac{l}{m{r}^2(t\stackrel{´}{})}dt\stackrel{´}{}$

Es ergibt sich also: $\varphi (t)$ .

Die Bahngleichung wird gewonnen gemäß:

$\frac{dr}{d\varphi }=\frac{\dot{r}}{\dot{\varphi }}=\frac{m{r}^2\sqrt{\frac{2}{m}(E-\tilde{V}(r))}}{l}=r^2\sqrt{\frac{2m}{l^2}(E-\tilde{V}(r))}$

Es folgt:

$\underset{{\varphi }_o}{\overset{\varphi }{\int }}d\varphi \stackrel{´}{}=\underset{r_o}{\overset{r}{\int }}dr\stackrel{´}{}\frac{1}{r{\stackrel{´}{}}^2\sqrt{\frac{2m}{l^2}(E-\tilde{V}(r\stackrel{´}{}))}}$

Daraus erhält man als Bahngleichung $\varphi (r)$ bzw. $r(\varphi )$ .

Die Bahngleichung.

Planetenbewegung und Keplersche Gesetze

Betrachten wir speziell das Gravitationspotenzial als Wechselwirkung:

$V(r)=-\frac{\gamma m_1{m}_2}{r^{}}$ mit $r=|{\overline{r}}_1-{\overline{r}}_2|$

Somit ergibt sich ein effektives Radialpotenzial gemäß

$\tilde{V}(r)=-\frac{k}{r^{}}+\frac{l^2}{2m{r}^2}k:=\gamma m_{1}m>0$

ALs Grenzwert folgt:

$\begin{array}{rl}& r\to 0:\tilde{V}(r)=\frac{l^2}{2m{r}^2}\to \mathrm{\infty }\\ & r\to \mathrm{\infty }:\tilde{V}(r)=-\frac{k}{r}\to 0\\ \end{array}$

Differenziation findet ein Minimum:

$\begin{array}{rl}& \frac{d\tilde{V}(r)}{dr}=\frac{k}{r^2}-\frac{l^2}{m{r}^3}=0\to r_o=\frac{l^2}{mk}\\ & \tilde{V}(r_o)=\frac{-m{k}^2}{2l^2}\\ \end{array}$

Wegen

$\frac{m}{2}{\dot{r}}^2=E-\tilde{V}(r)$ ist eine Bewegung nur für $E-\tilde{V}(r)\ge 0$ möglich. Also muss $E\ge \tilde{V}(r)$

Es gilt:

$0>E\ge \tilde{V}(r_o)\Rightarrow 0>E\ge \frac{-m{k}^2}{2l^2}$

Bahnen sind geschlossen ( Ellipse, Spezialfall: Kreis)

$E>0$ Bahnen sind offen. ( Hyperbeln)

Wir werden sehen, dass für E=0 eine Parabelbahn folgt.

Das Potenzial hat die folgende Gestalt:

Für $0>E\ge \tilde{V}(r_o)\Rightarrow 0>E\ge \frac{-m{k}^2}{2l^2}$

Sind die Umkehrpunkte durch $\tilde{V}(r)=-\frac{k}{r^{}}+\frac{l^2}{2m{r}^2}=E$

bestimmt ( quadratisch Gleichung in r mit zwei Lösungen):

$r_{\min /\max }=\frac{1}{2\left|E\right|}(k\mp \sqrt{k^2-\frac{2l^2\left|E\right|}{m}})$

Für E>0 gibt es nur noch eine Lösung für r, die positiv und damit physikalisch sinnvoll ist.

Aus

gewinnt man den inneren Umkehrpunkt:  

Die Bahngleichung kann nun explizit berechnet werden:

$\underset{{\varphi }_o}{\overset{\varphi }{\int }}d\varphi \stackrel{´}{}=\varphi -{\varphi }_o=\underset{r_o}{\overset{r}{\int }}dr\stackrel{´}{}\frac{1}{r{\stackrel{´}{}}^2\sqrt{\frac{2m}{l^2}(E-\tilde{V}(r\stackrel{´}{}))}}=\underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\frac{1}{\sqrt{\frac{2mE}{l^2}+\frac{2mk}{l^{2}r\stackrel{´}{}}-\frac{1}{r{\stackrel{´}{}}^2}}}$

Dieses Integral ist nicht leicht zu berechnen, jedoch lediglich ein mathematisches Problem. Es gelingt mit einer geschickten Substitution:

Zunächst soll der Ausdruck unter der Wurzel quadratisch ergänzt werden:

$\frac{2mE}{l^2}+\frac{2mk}{l^{2}r\stackrel{´}{}}-\frac{1}{r{\stackrel{´}{}}^2}=-{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^2+\frac{m^2{k}^2}{l^4}+\frac{2mE}{l^2}$

mit

$\begin{array}{rl}& -{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^2+\frac{m^2{k}^2}{l^4}+\frac{2mE}{l^2}:=D[1-\frac{1}{D}{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^2]\\ & D:=\frac{2m}{l^2}(\frac{m{k}^2}{2l^2}+E)\\ \end{array}$

Dabei gilt:

$\begin{array}{rl}& \frac{m{k}^2}{2l^2}=\tilde{V}(r_o)\\ & \Rightarrow D:=\frac{2m}{l^2}(\frac{m{k}^2}{2l^2}+E)\ge 0\\ \end{array}$

Substitution:

$\begin{array}{rl}& \cos \vartheta \stackrel{´}{}:=\frac{1}{\sqrt{D}}(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})\Rightarrow \frac{d\cos \vartheta }{dr\stackrel{´}{}}=-\frac{1}{\sqrt{D}}\left(\frac{1}{r{\stackrel{´}{}}^2}\right)\\ & \frac{d\cos \vartheta }{d\vartheta }=-\sin \vartheta \stackrel{´}{}\Rightarrow -\sin \vartheta d\vartheta =d\cos \vartheta \\ & -\sin \vartheta \stackrel{´}{}d\vartheta =-\frac{1}{\sqrt{D}}\left(\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\right)\\ \end{array}$

Somit folgt:

$\begin{array}{rl}& \varphi -{\varphi }_o=\underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\frac{1}{\sqrt{\frac{2mE}{l^2}+\frac{2mk}{l^{2}r\stackrel{´}{}}-\frac{1}{r{\stackrel{´}{}}^2}}}=\underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\frac{1}{\sqrt{D[1-\frac{1}{D}{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^2]}}=\underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\frac{1}{\sqrt{D}[1-\frac{1}{D}{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^{}]}\\ & \underset{r_o}{\overset{r}{\int }}\frac{dr\stackrel{´}{}}{r{\stackrel{´}{}}^2}\frac{1}{\sqrt{D}[1-\frac{1}{D}{(\frac{1}{r{\stackrel{´}{}}^{}}-\frac{mk}{l^2})}^{}]}=\underset{{\vartheta }_0}{\overset{\vartheta }{\int }}d\vartheta \stackrel{´}{}\sin \vartheta \stackrel{´}{}\frac{1}{\sqrt{1-{\cos }^2\vartheta }\stackrel{´}{}}=\underset{{\vartheta }_0}{\overset{\vartheta }{\int }}d\vartheta \stackrel{´}{}=\vartheta -{\vartheta }_0\\ & \vartheta -{\vartheta }_0=\arccos \frac{1}{\sqrt{D}}(\frac{1}{r^{}}-\frac{mk}{l^2})-\arccos \frac{1}{\sqrt{D}}(\frac{1}{{r_o}^{}}-\frac{mk}{l^2})\\ \end{array}$

Also in Summary:

$\varphi -{\varphi }_o=\arccos \frac{1}{\sqrt{D}}(\frac{1}{r^{}}-\frac{mk}{l^2})-\arccos \frac{1}{\sqrt{D}}(\frac{1}{{r_o}^{}}-\frac{mk}{l^2})$

Eine der Integrationskonstanten,

${\varphi }_o$ oder $r_o$ kann frei eingesetzt werden.

Wir wählen den Winkel willkürlich:

Mit der vereinfachenden Wahl von

${\varphi }_o=\arccos \frac{1}{\sqrt{D}}(\frac{1}{{r_o}^{}}-\frac{mk}{l^2})$

ergibt sich:

$\begin{array}{rl}& \varphi (r)=\arccos \frac{1}{\sqrt{D}}(\frac{1}{r^{}}-\frac{mk}{l^2})\\ & \Rightarrow \frac{1}{r(\varphi )}=\frac{mk}{l^2}+\sqrt{D}\cos \varphi =\frac{mk}{l^2}(1+\epsilon \cos \varphi )\\ & mit\epsilon :=\sqrt{D}\frac{l^2}{mk}=\sqrt{1+\frac{2E{l}^2}{m{k}^2}}\\ \end{array}$

Wesentlich ist unsere Bahngleichung:

$\frac{1}{r(\varphi )}=\frac{mk}{l^2}+\sqrt{D}\cos \varphi =\frac{mk}{l^2}(1+\epsilon \cos \varphi )$

Dies ist nämlich, wie jedem Mathematiker bekannt ist, die Gleichung eines Kegelschnitts in Polarkoordinaten:

$\begin{array}{rl}& \epsilon >1\cong E>0Hyperbel(offeneBahn)\\ & \epsilon =1\cong E=0Parabel(offeneBahn)\\ & \epsilon <1\cong -\frac{m{k}^2}{2l^2}<E<0Ellipse(geschlosseneBahn)\\ \end{array}$

Für da zweidimensionale Problem ist die Umrechnung auf kartesische Korodinaten sehr einfach:

Dies ist nämlich, wie jedem Mathematiker bekannt ist, die Gleichung eines Kegelschnitts in Polarkoordinaten:

$\begin{array}{rl}& \cos \varphi =\frac{x}{r}\\ & \sin \varphi =\frac{y}{r}\\ & r=\sqrt{(x^2+y^2)}\\ \end{array}$

Für $\epsilon <1$ folgt:

$\begin{array}{rl}& {(\frac{mk}{l^2}(1-{\epsilon }^2)x+\epsilon )}^2+\frac{m^2{k}^2}{l^4}(1-{\epsilon }^2)y^2=1\\ & \frac{m^2{k}^2}{l^4}{(1-{\epsilon }^2)}^2{(x+\frac{l^2}{mk}\frac{\epsilon }{(1-{\epsilon }^2)})}^2+\frac{m^2{k}^2}{l^4}(1-{\epsilon }^2)y^2=1\\ \end{array}$

Dies kann vereinfacht werden zu:

$\frac{{(x+e)}^2}{a^2}+\frac{y^2}{b^2}=1$

mit der Exzentrizität

$e=\sqrt{a^2-b^2}$

Dies ist die Gleichung einer Ellipse mit einem Brennpunkt im Ursprung.

Die Hauptachsen lauten:

$\begin{array}{rl}& a=\frac{l^2}{mk(1-{\epsilon }^2)}=\frac{k}{2\left|E\right|}\\ & b=\frac{l^2}{mk\sqrt{1-{\epsilon }^2}}=\frac{l}{\sqrt{2m\left|E\right|}}\\ \end{array}$

Die relative Exzentrizität:

$\epsilon =\frac{e}{a}=\sqrt{1-\frac{b^2}{a^2}}$

e, die absolute Exzentrizität ist der absolute Abstand zwischen Mittelpunkt der Ellipse und einem Brennpunkt.

Keplersches Gesetz

Folgt also aus der Bewegungsgleichung mit Gravitationspotenzial bei negativen Energien:

Die Planetenbahnen sind Ellipsen, in deren einen Brennpunkt die Sonne steht.

Keplersches Gesetz

T²~a³

Beweis:

Für die Fläche einer Ellipse gilt:

$F=\pi ab$

Wenn wir das zweite Keplersche Gesetz verwenden ( Flächensatz), so gilt:

$\frac{dF}{dt}=\frac{l}{2m}=const$

Es ergibt sich der folgende Zusammenhang mit der Umlaufzeit:

$\underset{0}{\overset{T}{\int }}dt\frac{dF}{dt}=\frac{l}{2m}T=F=\pi ab$

Aus der Herleitung des ersten Keplerschen Gesetzes ist bekannt:

$\begin{array}{rl}& \frac{b^2}{a}=\frac{l^2}{mk}\Rightarrow b=\frac{l}{\sqrt{mk}}\sqrt{a}\\ & T=\frac{2m\pi ab}{l}=\frac{2\sqrt{m}\pi a^{\frac{3}{2}}}{\sqrt{k}}\\ & \frac{T^2}{a^3}=\frac{4{\pi }^{2}m}{k}\\ \end{array}$

Die zweiten Potenzen der Umlaufdauer sind somit nicht exakt proportional zur dritten Potenz der großen Halbachsen, da auch die Masse des Planeten noch eingeht:

$\begin{array}{rl}& k=\gamma m_1{m}_2\\ & m=\frac{m_1{m}_2}{m_1+m_2}\\ & \frac{m}{k}=\frac{1}{\gamma (m_1+m_2)}\\ \end{array}$

Falls die Planeten jedoch deutlich leichter sind als die Zentralgestirne, so gilt:

$\begin{array}{rl}& \frac{m}{k}\approx \frac{1}{\gamma m_2}\\ & \frac{T^2}{a^3}\approx \frac{4{\pi }^2}{\gamma m_2}\\ \end{array}$

Leitet man dies aus dem Kraftansatz ab, so steckt der Fehler der Vernachlässigung der Planetenmasse in der Annahme einer kreisförmigen Bewegung um das Zentralgestirn. Das Ergebnis ist ebenso fehlerbelastet.