Virialsatz

Der Virialsatz (lateinisch vis ‚Kraft‘) ist eine Beziehung zwischen dem zeitlichen arithmetischen Mittelwert der kinetischen Energie ${\overline {T}}$ und dem zeitlichen Mittel der potentiellen Energie ${\overline {U}}$ eines abgeschlossenen stationären physikalischen Systems. Der Virialsatz ist ursprünglich als Satz der klassischen Mechanik formuliert und ermöglicht allgemeine Abschätzungen der Anteile potentieller und kinetischer Energie auch in komplexen Systemen zum Beispiel in der Astrophysik. Es gibt auch einen quantenmechanischen Virialsatz.

Der Virialsatz wurde 1870 von Rudolf Clausius in dem Aufsatz Über einen auf die Wärme anwendbaren mechanischen Satz aufgestellt.

Skalare Formulierung

Das Virial G eines Systems aus N Teilchen ist die Summe der Skalarprodukte aus den Impulsen ${\vec {p_{i}}}$ und den Orten ${\vec {r_{i}}}$ dieser Teilchen, d. h.:

$G=\sum _{i=1}^{N}{\vec {p_{i}}}\cdot {\vec {r_{i}}}$.

Ist das Virial beschränkt, so gilt der Virialsatz:

${\overline {T}}=-{\frac {1}{2}}\sum _{i=1}^{N}{\overline {{\vec {F_{i}}}\cdot {\vec {r_{i}}}}}$

Dabei ist ${\vec {F_{i}}}$ die Resultierende der auf das i-te Teilchen einwirkenden Kräfte, die von anderen Teilchen des Systems ausgeübt werden. Da ein abgeschlossenes System betrachtet wird, existieren keine äußeren Kräfte. ${\overline {A}}$ bezeichnet den zeitlichen Mittelwert für Zeiten $\tau \rightarrow \infty$.

Ist die Kraft konservativ und besitzt sie ein Potential U, das homogen vom Grad k ist, d. h., für α > 0 gilt $U(\alpha {\vec {r_{i}}})=\alpha ^{k}\cdot U({\vec {r_{i}}})$, so vereinfacht sich die obige Form zu

${\overline {T}}={\frac {k}{2}}\,{\overline {U}}.$

Mit der Gesamtenergie ${\overline {E}}={\overline {T}}+{\overline {U}}=E$ folgt:

${\overline {T}}={\frac {k}{2}}\;{\overline {U}}={\frac {k}{k+2}}\;E$

bzw.

${\overline {U}}={\frac {2}{k+2}}\;E$

Für den bekanntesten Fall $k=-1$ (Gravitation, Coulombsche Kraft) ergibt sich z. B.:

${\overline {T}}=-\;{\frac {1}{2}}\;{\overline {U}}=-\;E.$

Insbesondere ergibt sich, dass die Gesamtenergie für die Anwendung des Virialtheorems im Fall k = -1 negativ sein muss (da ${\overline {T}}$ positiv oder Null ist), entsprechend der Grundvorraussetzung des Virialtheorems, dass das System räumlich beschränkt sein soll, also gravitativ bzw. durch Coulombkräfte gebunden.

Für den Fall harmonischer Schwingungen (k =2) ist:

${\overline {T}}={\overline {U}}$

Sonderfälle

Geschlossene Bahnen

In zwei Sonderfällen von harmonischen Potentialen, nämlich für das Potential des harmonischen Oszillators ($k=2$) und des Coulombpotentials ($k=-1$), erhält man geschlossene Bahnen. Dann kann das zeitliche Mittel auch über einen Bahndurchlauf erfolgen.

Vielteilchensystem

Befindet sich ein Vielteilchensystem im thermischen Gleichgewicht, so kann das System als ergodisch betrachtet werden, d. h., das Zeitmittel ist gleich dem Scharmittel für alle Beobachtungsgrößen. Da dies insbesondere für die kinetische und die potentielle Energie gilt und das Scharmittel der Energien aus der Summe der Einzelenergien, geteilt durch die Anzahl N der Objekte, gebildet wird, lässt sich das Scharmittel durch die Gesamtenergien ausdrücken. Wir erhalten daher für Gleichgewichtssysteme:

$T={\frac {k}{2}}\,U$

ohne Mittelung über die Zeit, denn die Werte sind zeitlich konstant.

Ableitung des Virialsatzes

Hier wird der Darstellung im Lehrbuch von Landau und Lifschitz gefolgt, wo der Virialsatz in Zusammenhang mit Skalierungsverhalten mechanischer Größen (mechanische Ähnlichkeit) diskutiert wird. Dabei wird nur ausgenutzt, dass die kinetische Energie $T$ quadratisch in den Geschwindigkeiten ${\vec {v}}_{i}$ ist und die Impulse werden formal über ${\vec {p}}_{i}={\frac {\partial T}{\partial {\vec {v}}_{i}}}$ eingeführt. Dann ist nach dem Satz von Euler über homogene Funktionen:

$\sum _{i}{\frac {\partial T}{\partial {\vec {v}}_{i}}}\cdot {\vec {v}}_{i}=2T$

woraus folgt:

$2T=\sum _{i}{\vec {p}}_{i}\cdot {\vec {v}}_{i}={\frac {d}{dt}}(\sum _{i}{\vec {p}}_{i}\cdot {\vec {v}}_{i})-\sum _{i}{\vec {r}}_{i}\cdot {\frac {d}{dt}}{\vec {p}}_{i}$

Nun bildet man den asymptotischen Grenzwert des zeitlichen Mittelwerts:

${\overline {f}}=\lim _{\tau \to \infty }{\frac {1}{\tau }}\int _{0}^{\tau }f(t)dt$

Insbesondere gilt für den zeitlichen Mittelwert der Zeitableitung des Virials G:

${\overline {{\frac {d}{dt}}(\sum _{i}{\vec {p}}_{i}\cdot {\vec {v}}_{i})}}={\overline {({\frac {dG}{dt}})}}=\lim _{\tau \to \infty }{\frac {G(\tau )-G(0)}{\tau }}$

Hat man es mit einem System zu tun, in dem Geschwindigkeiten und Orte der Teilchen beschränkt sind (z.B. periodische Bahnen), so folgt

${\overline {({\frac {dG}{dt}})}}=0$

und es folgt mit ${\vec {F}}_{i}={\frac {d}{dt}}{\vec {p}}_{i}=-{\frac {\partial U}{\partial {\vec {r}}_{i}}}$ der Virialsatz:

$2{\overline {T}}={\overline {(\sum _{i}{\vec {r}}_{i}\cdot {\frac {\partial U}{\partial {\vec {r}}_{i}}})}}=k{\overline {U}}$

wenn man annimmt, dass das Potential U eine homogene Funktion der Ortskoordinaten vom Grad k ist.

Anwendungsbeispiel: Massenbestimmung astronomischer Haufen

Anwendung findet der Virialsatz beispielsweise in der Astrophysik und der Himmelsmechanik. Dort benutzt man das Newton’sche Gravitationspotential, das homogen vom Grad −1 ist. Dann gilt:

$2T=-U$

Der Virialsatz erlaubt es, recht gute Abschätzungen für die Gesamtmassen dynamisch gebundener Systeme wie Sternhaufen oder Galaxienhaufen zu finden. Die Gesamtmasse eines solchen Haufens kann dann vollständig durch Beobachtungsgrößen wie Radialgeschwindigkeiten, Winkelabstände und scheinbare Helligkeiten der Einzelobjekte ausgedrückt werden. Die einzige Voraussetzung für die Anwendung des Virialsatzes ist die Kenntnis des Abstandes des Haufens. Wir wollen das Vorgehen bei einer Massenbestimmung eines solchen Haufens hier skizzieren:

Die kinetische Gesamtenergie eines Stern- oder Galaxienhaufens ist durch

$T={\frac {1}{2}}\sum _{i}m_{i}\|{\vec {v}}_{i}\|^{2}$

gegeben. Aber weder die Einzelmassen m_i noch die Geschwindigkeitsbeträge |v_i| sind Beobachtungsgrößen. Die Einführung der Gesamtmasse $M=\sum m_{i}$ erlaubt die Umformung

$T={\frac {M}{2}}\sum _{i}{\frac {m_{i}}{M}}\|{\vec {v}}_{i}\|^{2}.$

Nun trifft man zwei Annahmen:

a) Die Einzelmassen m_i sind proportional zu den Einzelleuchtkräften l_i und daher gilt

$T={\frac {M}{2}}\sum _{i}{\frac {l_{i}}{L}}\|{\vec {v}}_{i}\|^{2}={\frac {M}{2}}\langle \|{\vec {v}}\|^{2}\rangle _{L},$

wobei der letzte Term das leuchtkraftgewichtete Mittel über die Geschwindigkeiten bezeichnet.

b) Das System ist sphärisch symmetrisch und befindet sich im Gleichgewicht (man sagt auch es ist virialisiert). Daher sind die Geschwindigkeiten über die Raumrichtungen gleichverteilt (Gleichverteilungssatz). Dann gilt

$\langle \|{\vec {v}}\|^{2}\rangle _{L}=3\langle v_{R}^{2}\rangle _{L},$

wobei v_R die radialen Pekuliargeschwindigkeiten bezeichnet, d. h. die Abweichungen der Radialgeschwindigkeit vom Haufenmittelwert. Damit erhält man:

$T={\frac {3M}{2}}\langle v_{R}^{2}\rangle _{L}$

Andererseits gilt für die potentielle Gesamtenergie unter der Bedingung der sphärischen Symmetrie

$U=-{\frac {GM^{2}}{\alpha R}}$

mit

• der Gravitationskonstanten G,
• dem Gesamtradius R des Systems,
• und dem morphologischen Faktor α, der von der radialen Verteilungsfunktion, also der Geometrie des Haufens, abhängt. Für eine (allerdings unrealistische) Gleichverteilung innerhalb des Radius R ist beispielsweise α = 5/3. Im Allgemeinen ist der Faktor aus den beobachteten Winkelabständen der Einzelsysteme zum Haufenzentrum zu bestimmen.

Durch Anwendung des Virialsatzes für die Gravitation erhalten wir die Gesamtmasse des Haufens zu

$M={\frac {3\alpha R}{G}}\langle v_{R}^{2}\rangle _{L}.$

Eine weitere astrophysikalische Anwendung ist die Abschätzung der Jeans-Masse.

Verallgemeinerung auf Tensoren

Im Rahmen der Kontinuumsmechanik wird der tensorielle Virialsatz aus der stoßfreien Boltzmann-Gleichung und den daraus abgeleiteten Jeans-Gleichungen bewiesen und in der Astrophysik verwendet. Wenn als Wechselwirkung wiederum die Gravitation angenommen wird, hat der Satz die Form

${\frac {1}{2}}{\frac {d^{2}}{dt^{2}}}I_{ij}=2T_{ij}+\Pi _{ij}+U_{ij}$

mit

• dem Trägheitstensor $I_{ij},$
• dem Tensor $T_{ij}$ der kinetischen Energie,
• dem Spannungstensor $\Pi _{ij}$ und
• dem Tensor $U_{ij}$ der potentiellen Energie.

Im statischen Fall fällt die Zeitableitung auf der linken Seite der Gleichung weg, und da der Spannungstensor spurfrei ist, ergibt die Spur der Gleichung wieder den skalaren Virialsatz.

Quantenmechanik

Für die Quantenmechanik behält der Virialsatz seine Gültigkeit, wie von Fock gezeigt wurde.

Der Hamiltonoperator des Systems aus Punktteilchen sei

$H=V(\{X_{i}\})+\sum _{n}P_{n}^{2}/2m.$

Man bilde den Kommutator von $H$ mit $X_{n}P_{n}$, gebildet aus dem Ortsoperator $X_{n}$ und dem Impulsoperator $P_{n}=-i\hbar d/dX_{n}$ des n-ten Teilchens:

$[H,X_{n}P_{n}]=X_{n}[H,P_{n}]+[H,X_{n}]P_{n}=i\hbar X_{n}{\frac {dV}{dX_{n}}}-i\hbar {\frac {P_{n}^{2}}{m}}$

Bildet man das Virial durch Summierung über die Teilchen $Q=\sum _{n}X_{n}P_{n}$, so folgt

${\frac {i}{\hbar }}[H,Q]=2T-\sum _{n}X_{n}{\frac {dV}{dX_{n}}}$

mit der kinetischen Energie $T=\sum _{n}P_{n}^{2}/2m$. Nach den Heisenbergschen Bewegungsgleichungen ist die linke Seite geich $-dQ/dt$. Der Erwartungswert $\langle dQ/dt\rangle$ verschwindet in einem stationären Zustand, sodass mit

$2\langle T\rangle =\sum _{n}\langle X_{n}dV/dX_{n}\rangle$

die Quantenversion des Virialsatzes folgt, wobei die spitzen Klammern für quantenmechanische Erwartungswerte der jeweiligen Operatoren für einen stationären Zustand stehen.

Literatur

• Lew Dawidowitsch Landau, Jewgeni Michailowitsch Lifschitz: Mechanik. Lehrbuch der theoretischen Physik, Bd. 1. Deutsch, Frankfurt/M. 2004, ISBN 3-8171-1326-9.

Gibt eine einfache Herleitung des skalaren Virialsatzes.

• James Binney, Scott Tremaine: Galactic Dynamics. Princeton Series in Astrophysics. Princeton University Press, Princeton, N.J. 1988, ISBN 0-691-08445-9.

Hier findet man die tensorielle Verallgemeinerung und Anwendungen.

Weblinks

• Eine weitere Herleitung der skalaren Variante (englisch).