Virialsatz

Der Virialsatz ({{#invoke:Vorlage:lang|full|CODE=la |SCRIPTING=Latn |SERVICE=lateinisch}} ‚Kraft‘) ist eine Beziehung zwischen den zeitlichen arithmetischen Mittelwerten der kinetischen Energie <math>\overline{T}</math> und der potentiellen Energie <math>\overline{U}</math> eines abgeschlossenen physikalischen Systems. Der Virialsatz wurde 1870 von Rudolf Clausius aufgestellt in dem Aufsatz Über einen auf die Wärme anwendbaren mechanischen Satz.

Das Virial ist dabei nach Clausius der Ausdruck<ref>R. Clausius: Über einen auf die Wärme anwendbaren mechanischen Satz. Annalen der Physik, Band 217, 1870, S. 124–130.</ref><ref name=Goldstein>H. Goldstein: Klassische Mechanik. Akademische Verlagsgesellschaft, 1978, S. 76 f.</ref><ref>Die Definitionen des Virials variieren etwas, z. B. lassen sowohl Wolfgang Pauli in seinen Vorlesungen über Thermodynamik (ETH Zürich 1958) als auch das unten zitierte Buch von Honerkamp den Vorfaktor −1/2 in der Definition des Virials weg und Pauli lässt auch die Mittelbildung weg.</ref>

<math>-\frac{1}{2} \sum_{i=1}^N \overline{\vec{F_i} \cdot \vec{r_i}}.</math>

Hierbei bezeichnet

• <math>\vec{F_i}</math> die auf das <math>i</math>-te Teilchen wirkende Kraft
• <math>\vec r_i</math> den Ortsvektor des <math>i</math>-ten Teilchens
• der Querstrich einen unten näher erläuterten Mittelwert, z. B. ein Zeit- oder Scharmittel.

Der Virialsatz wurde von Clausius ursprünglich als Satz der klassischen Mechanik formuliert (als Gleichheit von Virial und mittlerer kinetischer Energie). Er ermöglicht allgemeine Abschätzungen der Anteile potentieller und kinetischer Energie auch in komplexen Systemen, z. B. in Mehrkörperproblemen der Astrophysik. Es gibt auch einen quantenmechanischen Virialsatz, einen Virialsatz der statistischen Mechanik, aus dem u. a. das ideale Gasgesetz und Korrekturen für reale Gase abgeleitet wurden, sowie einen relativistischen Virialsatz.

Der Virialsatz gilt nur unter gewissen Voraussetzungen, etwa im Fall des Virialsatzes der Mechanik, dass mit zeitlicher Mittelwertbildung Orte und Geschwindigkeiten der Teilchen beschränkt sind, oder dass ein thermisches Gleichgewicht herrscht.

Virialsatz der Mechanik

Teilchen in einem konservativen Kraftfeld

Einen einfachen Fall stellen <math>N</math> untereinander nicht wechselwirkende Teilchen in einem äußeren Kraftfeld <math>\vec{F_E}</math> dar, das konservativ, also von einem Potential <math>\Phi(\vec r)</math> abgeleitet, ist (die dazugehörende Ladung sei mit <math>q</math> bezeichnet, sie ist für den Fall der Gravitation gerade die Masse):

<math>-\vec{F_E} (\vec r)= q \, \nabla\Phi( \vec r)</math>

Darin ist <math>\nabla \Phi(\vec{r})</math> der Gradient des Feldes bzw. des Potentials.

Der Virialsatz gilt, wie unten dargelegt wird, falls die Bewegung im Endlichen bleibt, also Ort und Impuls für alle Zeiten beschränkt sind, und lautet

<math>\begin{alignat}{2}

\overline T &= -&&\frac{1}{2} \sum_{i=1}^N \overline{\vec{F_i} \cdot \vec{r_i}}\\

           &=  &&\frac{q}{2} \sum_{i=1}^N \overline{\nabla \Phi(\vec{r_i}) \cdot \vec{r_i}},

\end{alignat}</math>

wobei

• <math>T</math> die kinetische Energie des Teilchens ist
• der Querstrich den zeitlichen Mittelwert für Zeiten <math>\tau \to \infty</math> bezeichnet.

Nimmt man zusätzlich ein in der Ortsvariablen homogenes Potential vom Grad <math>k</math> an, d. h. es gilt <math>\Phi(\alpha \, \vec {r}) = \alpha^k \cdot \Phi(\vec {r})</math> für <math>\alpha > 0</math> (Werte für k finden sich weiter unten in Folgerungen und Beispiele), dann vereinfacht sich obige Gleichung mit der Eulerschen Gleichung für homogene Funktionen:<ref name=Honerkamp>{{#invoke:Vorlage:Literatur|f}}</ref>

<math>\nabla \Phi(\vec{r}) \cdot \vec{r} = k \, \Phi(\vec{r})</math>

zu

<math>\overline{T} = \frac{k}{2} \, \overline U,</math>

wobei <math>\textstyle U = \sum q_i \Phi(\vec{r_i})</math> die gesamte potentielle Energie der Teilchen ist. Der Virialsatz ist daher eine Beziehung zwischen mittlerer kinetischer und mittlerer potentieller Energie.

Untereinander wechselwirkende Teilchen

Für die Ableitung der Gasgesetze und die Anwendung in der Astrophysik ist der Fall eines abgeschlossenen Systems von <math>N</math> miteinander wechselwirkenden Teilchen von besonderem Interesse. Wie oben ergibt sich unter der Voraussetzung einer im Endlichen ablaufenden Bewegung der Virialsatz:

<math>\overline T = -\frac{1}{2} \sum_{i=1}^N \overline{\vec{F_i} \cdot \vec{r_i}}</math>

Dabei ist <math>\vec {F_i}</math> die Resultierende der auf das <math>i</math>-te Teilchen einwirkenden Kräfte, die von anderen Teilchen des Systems ausgeübt werden. Da ein abgeschlossenes System betrachtet wird, existieren diesmal keine äußeren Kräfte. Da <math>\textstyle \sum_i \vec {F_i}=0</math> gilt, ist die Wahl des Ursprungs für die Ortsvektoren <math>\vec r_i</math> im Virial beliebig. Auf den ersten Blick sieht der Ausdruck im Virial kompliziert aus, lässt sich aber unter der Annahme, dass die paarweise zwischen den Teilchen wirkenden Kräfte jeweils von homogenen Potentialen vom Grad <math>k</math> abgeleitet werden können, wie oben auf die Form

<math>\overline{T} = \frac{k}{2}\,\overline{U}</math>

bringen.

Folgerungen und Beispiele

Mit der Gesamtenergie <math>\overline{E} = \overline{T} + \overline{U} = E</math> folgt aus dem Virialsatz:

<math>\overline{T} = \frac{k}{2}\,\overline {U} = \frac{k}{k+2}\,E</math>

<math>\overline{U} = \frac{2}{k+2}\,E</math>

Für den bekannten Fall <math>k = -1</math> (Gravitation, Coulombsche Kraft) ergibt sich z. B.:

<math>\overline{T} = -\frac{1}{2}\,\overline{U} = -E</math>

Insbesondere ergibt sich, dass die Gesamtenergie für die Anwendung des Virialtheorems im Fall <math>k = -1</math> negativ sein muss (da <math>\overline {T}</math> positiv ist).

Für harmonische Schwingungen (<math>k = 2</math>) gilt:

<math>\overline{T} = \overline{U} = \frac{1}{2}\,E</math>

Ableitung

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

<math>\sum_i \frac{\partial T}{\partial \vec v_i} \cdot \vec v_i =2 T,</math>

woraus

<math>2 T = \sum_i \vec p_i \cdot \vec v_i = \frac{d}{dt} ( \sum_i \vec p_i \cdot \vec r_i ) - \sum_i \vec r_i \cdot \frac{d}{dt} \vec p_i = \frac{d G}{dt} - \sum_i \vec r_i \vec F_i</math>

folgt, wobei <math>G</math> die Summe der Skalarprodukte aus den Impulsen <math>\vec{p}_i</math> und den Orten <math>\vec{r}_i</math> aller Teilchen ist:

<math>G= \sum_{i=1}^N \vec{p}_i \cdot \vec{r}_i</math>

Nun bildet man den asymptotischen Grenzwert des zeitlichen Mittelwerts:

<math>\overline {f}=\lim_{\tau \to\infty} \frac{1}{\tau} \int_0^{\tau} f(t) dt</math>

Insbesondere gilt für den zeitlichen Mittelwert der Zeitableitung von <math>G</math>:

<math>\overline {\frac{d}{dt} (\sum_i \vec p_i \cdot \vec r_i)} = \overline {\left(\frac{d G}{dt}\right)} = \lim_{\tau \to\infty} \frac{G(\tau) - G(0)}{\tau}</math>

Hat man es mit einem System zu tun, in dem die Geschwindigkeiten und Orte der Teilchen beschränkt sind (z. B. bei periodischen Bahnen),<ref name=Honerkamp /> so folgt

<math>\overline {\left(\frac{dG}{dt} \right)} = 0</math>

und mit <math>\vec F_i = \frac{d}{dt} \vec p_i = -\frac{\partial U}{\partial \vec r_i}</math> weiter der Virialsatz

<math>2 \overline {T} =\overline {\left( \sum_i \vec r_i \cdot \frac{\partial U}{\partial \vec r_i} \right)}= k \overline {U},</math>

wenn man annimmt, dass das Potential <math>U</math> eine homogene Funktion der Ortskoordinaten vom Grad <math>k</math> ist. In dieser Sicht drückt der Satz eine Gleichheit von Mittelwerten von kinetischer und potentieller Energie aus mit Vorfaktoren, die sich aus dem Skalierungsverhalten ergeben: 2 bei der kinetischen Energie, da die Geschwindigkeiten oder Impulse quadratisch eingehen, <math>k</math> beim Potential, da die Ortsvariablen mit Potenz <math>k</math> eingehen.

Eine ähnliche Ableitung findet sich schon bei Clausius und in dem Lehrbuch der klassischen Mechanik von Herbert Goldstein.<ref name=Goldstein /> Goldstein weist auch darauf hin, dass der Virialsatz mit Potentialterm auch dann gilt, wenn zusätzlich zu den Potentialkräften Reibungskräfte vorhanden sind, die proportional zur Geschwindigkeit sind, da diese keinen Beitrag zum Virialsatz liefern. Das gilt aber nur, falls sich ein Fließgleichgewicht einstellt, also Energie zugeführt wird, sodass die Bewegung nicht vollständig zum Erliegen kommt, denn dann würden alle Zeitmittelwerte verschwinden.

Sonderfälle der Mittelwertbildung

Gewöhnlich bezeichnet der Querstrich wie schon bei Clausius den zeitlichen Mittelwert für Zeiten <math>\tau \to \infty</math>. In bestimmten Sonderfällen kann das aber auch vereinfacht werden.

Geschlossene Bahnen

Liegen geschlossene Bahnen vor, so kann das Zeitmittel durch die Mittelung über eine Periode ersetzt werden. Der Virialsatz folgt hier unmittelbar aus der Periodizität der Bewegung.

In zwei Sonderfällen homogener Potentiale, nämlich für das Potential des harmonischen Oszillators (<math>k = 2</math>) und für das Coulombpotential (<math>k = -1</math>), erhält man für finite (d. h. nicht ins Unendliche gehende) Bewegungen im Ein- oder Zweikörperproblem immer geschlossene Bahnen.<ref>{{#invoke:Vorlage:Literatur|f}}</ref>

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 gebildet wird aus der Summe der Einzelenergien, geteilt durch die Anzahl <math>N</math> der Objekte, lässt sich das Scharmittel durch die Gesamtenergien ausdrücken. Wir erhalten daher für Gleichgewichtssysteme

<math>T = \frac{k}{2}\,U</math>

ohne Mittelung über die Zeit, denn die Werte sind zeitlich konstant (siehe auch unten die Behandlung des Virialsatzes im Rahmen der statistischen Mechanik).

Astrophysik

Für das gravitative <math>N</math>-Teilchensystem in der Astrophysik (z. B. als Modell von Galaxien- und Sternhaufen) ist die o. g. Grundvoraussetzung in der Ableitung des Virialsatzes, nämlich dass das System räumlich beschränkt bleibt, für große Zeiträume nicht gegeben. All diese Haufen lösen sich irgendwann auf, da immer wieder Teilchen durch die gegenseitige Wechselwirkung (Störung) mit den anderen genug Energie aufsammeln, um zu entkommen.

Allerdings sind die Zeiträume, in denen das geschieht, sehr lang: In der Astrophysik definiert die Relaxationszeit <math>T_\text{relax}</math> eines Sternhaufens oder einer Galaxie die Zeit, in der sich eine Gleichgewichtsverteilung einstellt.<ref name="Voigt">H. Voigt: Abriss der Astronomie. BI Verlag, 1980, S. 367 ff., S. 487.</ref> Sie beträgt bei der Milchstraße <math>T_\text{relax} \approx 7 \cdot 10^{13}</math> Jahre (bei einem Alter von <math>13{,}6 \cdot 10^9</math> Jahren) und für typische Kugelsternhaufen <math>10^{10}</math> Jahre. Innerhalb des Zeitraums <math>T_\text{relax}</math> erreichen 0,74 Prozent der Sterne nach der Maxwellschen Geschwindigkeitsverteilung die Fluchtgeschwindigkeit und entweichen.

Numerische Rechnungen zeigten, dass der Anteil sogar noch etwas höher liegt,<ref>Sebastian von Hoerner: Zeitschrift für Astrophysik. Band 50, 1960, 184. Danach etwa fünfmal höher.</ref> und dass der Virialsatz in den Haufen aufgrund des sich einstellenden Gleichgewichts (mit einer Anlaufzeit von zwei bis drei Relaxationszeiten) gut erfüllt ist. Nach dem Ablauf von <math>42 \cdot T_\text{relax}</math> sind 90 Prozent der Sterne abgewandert.

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:

<math>2 T = -U</math>

Der Virialsatz erlaubt es, recht gute Abschätzungen für die Gesamtmassen dynamisch gebundener Systeme wie Sternhaufen, Galaxien 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 anhand der Massenbestimmung eines solchen Haufens hier skizzieren:

Die kinetische Gesamtenergie eines Stern- oder Galaxienhaufens ist durch

<math>T = \frac{1}{2} \sum_i m_i v_i^2</math>

gegeben. Aber weder die Einzelmassen <math>m_i</math> noch die Geschwindigkeitsbeträge <math>v_i</math> sind Beobachtungsgrößen. Um diese einzuführen, müssen die Beiträge der einzelnen Objekte durch die Gesamtmasse <math>\textstyle M = \sum m_i</math> und geeignete Mittelwerte ausgedrückt werden. Zum Beispiel kann man annehmen, dass die Einzelmassen <math>m_i</math> proportional zu den Einzelleuchtkräften <math>l_i</math> sind und ein leuchtkraftgewichtetes Mittel bilden (durch den Index <math>L</math> angedeutet):

<math>T = \frac{M}{2} \sum_i \left(\frac{m_i}{M} \cdot v_i^2 \right) = \frac{M}{2} \sum_i \left( \frac{l_i}{L} \cdot v_i^2 \right) = \frac{M}{2} \langle v^2 \rangle_L</math>

Nimmt man an, dass das System sphärisch symmetrisch ist und sich im Gleichgewicht befindet (man sagt dann auch, es ist virialisiert), dann sind die Geschwindigkeiten über die Raumrichtungen gleichverteilt und es gilt

<math>\langle v^2 \rangle = 3 \langle v_{R}^2 \rangle,</math>

wobei <math>\textstyle \sqrt{\langle v^2 \rangle}</math> bzw. <math>\textstyle \sqrt{\langle v_R^2 \rangle}</math> die Streuungen (Abweichungen vom Mittelwert) der Geschwindigkeiten sind, das heißt die räumlichen bzw. Radialgeschwindigkeiten relativ zum Schwerpunkt des Haufens.<ref name="Voigt" /> Beispielsweise haben die Galaxien des Coma-Haufens eine Gaußverteilung der Geschwindigkeiten mit einer Streuung <math>v</math> von 1000 km/s. Damit erhält man:

<math>T = \frac{3 M}{2} \langle v_{R}^2 \rangle</math>

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

<math>U= -\frac{\alpha G M^2}{R}</math>

mit

• der Gravitationskonstanten <math>G</math>,
• dem Gesamtradius <math>R</math> des Systems und
• einem Faktor <math>\alpha</math>, der von der Größenordnung 1 ist und von der radialen Verteilungsfunktion, also der Geometrie des Haufens, abhängt. Für eine (allerdings unrealistische) Gleichverteilung innerhalb des Radius <math>R</math> ist beispielsweise <math>\alpha = 3/5</math>. 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:<ref>Roger Tayler: Galaxien. Aufbau und Entwicklung. Vieweg, 1986, S. 120.</ref><ref name="Voigt" />

<math>M = \frac{3 R}{\alpha G} \langle v_{R}^2 \rangle</math>

Die sich aus der Beobachtung ergebende Masse heißt Virialmasse. Da <math>\alpha</math> von der Größenordnung 1 ist, sieht man außerdem, dass die mittlere Geschwindigkeit <math>\langle v \rangle</math> etwa der Fluchtgeschwindigkeit entspricht (mit genauer Übereinstimmung für <math>\alpha = 2</math>).

Obwohl diese Methode der Massenbestimmung mit Unsicherheiten behaftet ist, merkte mit ihr bei der Messung von stark abweichenden Fluchtgeschwindigkeiten von Galaxienhaufen und der Deutung der Rotverschiebung Fritz Zwicky schon 1933 an, dass ein Großteil der Masse sehr dicht in Form Dunkler Materie vorliegen könne: Die Summe der Massen der sichtbaren Galaxien des Haufens lag eine Größenordnung niedriger.<ref>A. Unsöld, B. Baschek: Der neue Kosmos. Springer, 1988, S. 346.</ref> Denn zur Erklärung der Rotverschiebung sei eine 400-mal größere Massendichte erforderlich, als die aus den Massen der leuchtenden Materie abgeleitete Dichte. „Falls sich dies bewahrheiten sollte, würde sich also das überraschende Resultat ergeben, dass dunkle Materie in sehr viel größerer Dichte vorhanden ist als leuchtende Materie.“<ref>F. Zwicky, Die Rotverschiebung von extragalaktischen Nebeln. Helvetica Physica Acta, Band 6, 1933, S. 125. Online</ref> Auch bei elliptischen Galaxien ergab sich, dass die Virialmasse um Faktoren 10 bis 100 größer als die leuchtende Masse ist. Im Gegensatz zu Spiralgalaxien, wo man die Masse aus der Rotationskurve bestimmen kann, ist die Virialmethode bei elliptischen Galaxien häufig die einzige Methode der Massenbestimmung.

Eine weitere astrophysikalische Anwendung ist die Abschätzung der Jeans-Masse und der Satz findet auch Anwendung in Untersuchungen zur Stabilität von Gaskugelmodellen für Sterne.<ref name="Chandrasekhar">S. Chandrasekhar: An introduction to the study of stellar structure. Chicago 1939, S. 51 ff.</ref> Für ein durch Gravitation zusammengehaltenes ideales Gas als Sternmodell lässt sich mit dem Virialsatz zeigen, dass der Stern in der Endphase (wenn alle Fusionsprozesse zum Erliegen gekommen sind) nicht abkühlen kann. Erhöht sich der Betrag der gravitativen Bindungsenergie <math>U</math> durch die Kontraktion des Sterns, geht die Hälfte des Zuwachses in die kinetische Energie der als ideales Gas aufgefassten Sternmaterie und erhöht somit die Temperatur, der Rest wird abgestrahlt.<ref>Wolfgang Hillebrandt, Ewald Müller: Einführung in die Theoretische Astrophysik. Skript der TU München, 2008, Kapitel 2 (PDF).</ref> Wird der Druck im Innern zu hoch, bricht die Beschreibung als klassisches ideales Gas allerdings zusammen, da sich ein entartetes Fermigas bildet (Weißer Zwerg).

Tensor-Form

Im Rahmen der Kontinuumsmechanik wird der tensorielle Virialsatz aus der stoßfreien Boltzmann-Gleichung bewiesen und in der Astrophysik verwendet.

Wenn als Wechselwirkung wiederum die Gravitation angenommen wird, hat der Satz die Form

<math>\frac{1}{2} \frac{d^2}{dt^2} I_{ij} = 2 T_{ij} + \Pi_{ij} + U_{ij}</math>

mit

• dem Trägheitstensor <math>I_{ij},</math>
• dem Tensor <math>T_{ij}</math> der kinetischen Energie,
• dem Spannungstensor <math>\Pi_{ij}</math> und
• dem Tensor <math>U_{ij}</math> 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.

Das Auftreten der zweiten Zeitableitung des Trägheitstensors kann aus folgender Umformulierung von <math>G</math> im skalaren Fall motiviert werden:

<math>G = \sum_{k=1}^N \vec p_k \cdot \vec x_k = \sum_{k=1}^N m_{k} \, \frac{d \vec x_k}{dt} \cdot \vec x_k = \frac{1}{2} \frac{d}{dt} \sum_{k=1}^N m_{k} \, \vec x_k \cdot \vec x_k = \frac{1}{2} \frac{dI}{dt}</math>

mit dem skalaren Trägheitsmoment <math>I = \sum_{k=1}^{N} m_{k} {\vec x_k}^{2}.</math>

Varianten in der Astrophysik

Für Anwendungen in der Astrophysik wurde folgende Form des Virialsatzes

<math>\frac{1}{2} \frac{d^2 I}{d t^2} = 2 T + \Omega</math>

zuerst von Henri Poincaré<ref>H. Poincaré: Leçons sur les hypothèses cosmogoniques. Paris 1911.</ref> und Arthur Eddington<ref>A. Eddington: Monthly Notices Roy. Astron. Soc. 76, 1916, 528.</ref> abgeleitet.<ref name="Chandrasekhar" />

Für stationäre Systeme verschwindet die linke Seite, und in der betrachteten Anwendung war <math>\Omega</math> die potentielle gravitative Energie der Teilchen einer Gaswolke oder der Sterne in Galaxien:

<math>\Omega= -\sum_{i \neq j} \frac{G m_i m_j}{r_{ij}}</math>

In der Himmelsmechanik war diese Form des Virialsatzes schon Joseph-Louis Lagrange bekannt (1772, in einer Abhandlung zum Dreikörperproblem) und von Carl Gustav Jacobi verallgemeinert worden (Vorlesungen über Dynamik).<ref>S. Chandrasekhar: Hydrodynamic and hydromagnetic stability. Oxford University Press, 1961, S. 596.</ref>

Eine Aufteilung der kinetischen Energie <math>T</math> in

• einen Anteil <math>E_\mathrm{kin}</math> der hydrodynamischen Flüsse und
• einen Anteil <math>E_W</math> der zufälligen Wärmebewegung

sowie bei der potentiellen Energie eine zusätzliche Betrachtung

• eines Anteil <math>E_M</math> von Magnetfeldern

liefert den Virialsatz in folgender skalarer Form:<ref>Henrik Beuther: Sternentstehung. Skript, 2009 (PDF; 2,8 MB).</ref>

<math>\frac{1}{2} \frac{d^2 I}{d t^2} = 2 E_\mathrm{kin} + 2 E_W + \Omega + E_M</math>

Eine Tensorform dieses Virialsatzes für astrophysikalische Anwendungen in Anwesenheit magnetischer Felder wurde 1954 von Eugene N. Parker gegeben<ref>E. Parker: Tensor Virial Equations. Physical Review 96, 1954, 1686–1689.</ref> sowie 1953 von Subramanyan Chandrasekhar und Enrico Fermi.<ref>S. Chandrasekhar, E. Fermi: Problems of Gravitational Stability in the Presence of a Magnetic Field. Astrophysical Journal, 118, 1953, 116.</ref> Chandrasekhar entwickelte auch spezialisierte Virialsätze für seine Diskussion der Gleichgewichtsfiguren rotierender Flüssigkeiten.<ref>S. Chandrasekhar: Ellipsoidal figures of equilibrium. Yale University Press, 2009.</ref>

In der Plasmaphysik lässt sich als Anwendung des Virialsatzes zeigen, dass es keine stationären endlichen, durch die eigenen Magnetfelder eingeschlossenen Plasmakonfigurationen (Plasmoide) gibt.<ref>George Schmidt: Physics of High Temperature Plasmas. Academic Press, 1979, S. 72.</ref> Stattdessen sind für den Einschluss des Plasmas z. B. äußere Wände oder äußere Magnetfelder erforderlich.

Der Virialsatz der Quantenmechanik

Für die Quantenmechanik behält der Virialsatz seine Gültigkeit, wie bereits von Heisenberg, Jordan und Born in der Dreimännerarbeit 1926 gezeigt<ref>{{#invoke:Vorlage:Literatur|f}}{{#if:

       | {{#if: Vorlage:Cite book/ParamBool
               | Vorlage:Toter Link/archivebot
               | Vorlage:Webarchiv/archiv-bot
         }}
  }}{{#invoke:TemplatePar|check
   |all    = title=
   |opt    = vauthors= author= author1= authorlink= author-link= author-link1= author1-link= author2= author3= author4= author5= author6= author7= author8= author9= editor= last= first= last1= first1= last2= first2= last3= first3= last4= first4= last5= first5= last6= first6= last7= first7= last8= first8= last9= first9= last10= first10= last11= first11= last12= first12= last13= first13= last14= first14= last15= first15= others= script-title= trans-title= date= year= volume= issue= number= series= page= pages= at= issn= arxiv= bibcode= doi= pmid= pmc= jstor= oclc= id= url= url-status= format= access-date= archive-date= archive-url= archivebot= offline= location= publisher= language= quote= work= journal= newspaper= magazine= periodical=  name-list-style= url-access= doi-access= display-authors= via= s2cid= mr= type= citeseerx=  accessdate= archivedate= archiveurl= coauthors= month= day= last16= first16= last17= first17= last18= first18= last19= first19= last20= first20= last21= first21= last22= first22= last23= first23= last24= first24= last25= first25= last26= first26= last27= first27= last28= first28= last29= first29= last30= first30= last31= first31=
   |cat      = Wikipedia:Vorlagenfehler/Vorlage:Cite journal
   |errNS    = 0
   |template = Vorlage:Cite journal
   |format   = 
   |preview  = 1
  }}Vorlage:Cite book/URL{{#if:  | Vorlage:Cite book/Meldung }}{{#if:        | Vorlage:Cite book/Meldung }}{{#if: Zeitschrift für Physik
     || Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
        | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
       | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}Vorlage:Cite book/Meldung2{{#ifexpr: 0{{#ifeq:^^|^^||+1}}{{#ifeq:Max Born|^^||+1}}{{#ifeq:^^|^^||+1}}{{#ifeq:^^|^^||+1}} > 1
    | Vorlage:Cite book/Meldung
  }}</ref> und später unter anderem von Fock verallgemeinert wurde.<ref>{{#invoke:Vorlage:Literatur|f}}{{#if: 
       | {{#if: Vorlage:Cite book/ParamBool
               | Vorlage:Toter Link/archivebot
               | Vorlage:Webarchiv/archiv-bot
         }}
  }}{{#invoke:TemplatePar|check
   |all    = title=
   |opt    = vauthors= author= author1= authorlink= author-link= author-link1= author1-link= author2= author3= author4= author5= author6= author7= author8= author9= editor= last= first= last1= first1= last2= first2= last3= first3= last4= first4= last5= first5= last6= first6= last7= first7= last8= first8= last9= first9= last10= first10= last11= first11= last12= first12= last13= first13= last14= first14= last15= first15= others= script-title= trans-title= date= year= volume= issue= number= series= page= pages= at= issn= arxiv= bibcode= doi= pmid= pmc= jstor= oclc= id= url= url-status= format= access-date= archive-date= archive-url= archivebot= offline= location= publisher= language= quote= work= journal= newspaper= magazine= periodical=  name-list-style= url-access= doi-access= display-authors= via= s2cid= mr= type= citeseerx=  accessdate= archivedate= archiveurl= coauthors= month= day= last16= first16= last17= first17= last18= first18= last19= first19= last20= first20= last21= first21= last22= first22= last23= first23= last24= first24= last25= first25= last26= first26= last27= first27= last28= first28= last29= first29= last30= first30= last31= first31=
   |cat      = Wikipedia:Vorlagenfehler/Vorlage:Cite journal
   |errNS    = 0
   |template = Vorlage:Cite journal
   |format   = 
   |preview  = 1
  }}Vorlage:Cite book/URL{{#if:  | Vorlage:Cite book/Meldung }}{{#if:        | Vorlage:Cite book/Meldung }}{{#if: Zeitschrift für Physik
     || Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
        | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
       | Vorlage:Cite book/Meldung
  }}{{#if: Vorlage:Cite book/ParamBool
     | Vorlage:Cite book/Meldung
  }}Vorlage:Cite book/Meldung2{{#ifexpr: 0{{#ifeq:W. A. Fock|^^||+1}}{{#ifeq:^^|^^||+1}}{{#ifeq:^^|^^||+1}}{{#ifeq:^^|^^||+1}} > 1
    | Vorlage:Cite book/Meldung
  }}</ref>

Der Hamiltonoperator des Systems aus Punktteilchen sei

<math>H = V(\{X_i\})+\sum_n P_n^2/2m.</math>

Man bilde den Kommutator von <math>H</math> mit <math>X_n P_n</math>, gebildet aus dem Ortsoperator <math>X_n</math> und dem Impulsoperator <math>P_n = -i\hbar d/dX_n</math> des <math>n</math>-ten Teilchens:

<math>[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}</math>

Bildet man durch Summierung über die Teilchen <math>\textstyle Q = \sum_n X_n P_n</math>, so folgt

<math>\frac{i}{\hbar} [H,Q] = 2 T-\sum_n X_n \frac{dV}{dX_n}</math>

mit der kinetischen Energie <math>\textstyle T = \sum_n P_n^2/2m</math>.

Nach den Heisenbergschen Bewegungsgleichungen ist die linke Seite gleich <math>-dQ/dt</math>. Der Erwartungswert <math>\langle dQ/dt\rangle</math> verschwindet in einem stationären Zustand, sodass mit

<math>2\langle T\rangle = \sum_n\langle X_n dV/dX_n\rangle</math>

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

Der Virialsatz der statistischen Mechanik

Wie der Gleichverteilungssatz gehört auch eine Version des Virialsatzes zu den allgemeinen Aussagen der klassischen statistischen Mechanik.

Als Mittelbildung mit Hilfe des kanonischen Ensembles erhält man (vgl. den Gleichverteilungssatz):

<math>\left\langle x_i \frac{\partial H}{\partial x_i} \right\rangle = k_\mathrm{B} T</math>

<math>\left\langle p_i \frac{\partial H}{\partial p_i} \right\rangle = k_\mathrm{B} T</math>

mit <math>H = H_\mathrm{kin} + U (x)</math>.

Die untere Gleichung liefert:

<math>\frac{1}{2} \left\langle p_i \frac{\partial H}{\partial p_i} \right\rangle = \left\langle \frac{p_i^2}{2m} \right\rangle = \frac{1}{2} k_\mathrm{B} T</math>,

also einen Beitrag <math>\frac{1}{2} k_\mathrm{B} T</math> pro Freiheitsgrad für die mittlere kinetische Energie (Gleichverteilungssatz).

Die untere und obere Gleichung zusammen liefern den Virialsatz der statistischen Mechanik:

<math>\left\langle H_\mathrm{kin} \right\rangle = \left\langle \sum_i \frac{p_i^2}{2m} \right\rangle = \frac{1}{2} \sum_i \left\langle \vec x_i \frac{\partial U}{\partial \vec x_i} \right\rangle\, ,</math>

der auch in der Quantenstatistik gilt.

Es ist nach Clausius üblich, den Beitrag des Potentials aufzuteilen in

• das innere Virial, d. h. den Beitrag <math>V_\mathrm{int} (\vec x_i)</math> des Potentials der inneren Kräfte (Wechselwirkung der Teilchen untereinander) und
• das äußere Virial, d. h. den Beitrag <math>\textstyle W = \sum_i W (\vec x_i)</math> des Wandpotentials bzw. der Kräfte auf die Wand.

Das äußere Virial liefert:

<math>\sum_i \left\langle \vec x_i \frac{\partial W}{\partial \vec x_i} \right\rangle = p \int d \vec f \cdot \vec x = p \int d V (\operatorname{div} \vec x ) = 3 p V</math>

mit

• dem Druck <math>p</math> und
• dem Volumen <math>V</math>.

Dabei wurde über die Oberfläche (Wand) integriert und der Gaußsche Integralsatz angewandt.

Damit erhält man die Virialform der thermischen Zustandsgleichung:

<math>3 p V = 2 \left\langle H_\mathrm{kin} \right\rangle - \sum_i \left\langle \vec x_i \frac{\partial V_{int}}{\partial \vec x_i} \right\rangle\,</math>,

also für <math>N</math> Teilchen mit dem Gleichverteilungssatz:

<math>p V = N k_\mathrm{B} T - \frac{1}{3} \sum_i \left\langle \vec x_i \frac{\partial V_{int}}{\partial \vec x_i} \right\rangle</math>

Das ist die ideale Gasgleichung mit dem Virial der inneren Kräfte als Zusatzterm. Das Virial kann nach Potenzen der Teilchendichte <math>N/V</math> entwickelt werden (siehe: Virialentwicklung) für die Entwicklung von Zustandsgleichungen für reale Gase.

Die Ableitung der Gasgleichung war das Hauptziel der ursprünglichen Arbeit von Clausius, wobei er den Virialsatz der Mechanik als Grundlage benutzte.

Der Virialsatz der Relativitätstheorie

Es gibt auch einen relativistischen Virialsatz. Für Teilchen in Wechselwirkung mit elektromagnetischen Feldern findet er sich im Lehrbuch der theoretischen Physik von Landau und Lifschitz,<ref>Landau, Lifschitz, Klassische Feldtheorie. Band 2, Akademie Verlag, 1977, S. 99 f., § 34.</ref> er lässt sich aber auch für andere Wechselwirkungen formulieren.<ref>J. Gaite: The relativistic virial theorem and scale invariance. Physics Uspekhi, Band 56, 2013, S. 919.</ref>

Da die Spur des Energie-Impuls-Tensors des elektromagnetischen Feldes verschwindet, kann man – unter Verwendung des vierdimensionalen Energieerhaltungssatzes für Systeme mit beschränkter Bewegung (Impulse, Koordinaten u. a. variieren zwischen endlichen Schranken, die elektromagnetischen Felder verschwinden im Unendlichen) – ähnlich wie beim klassischen Virialsatz durch Mittelung über die Zeit zeigen:

<math>E = \langle \sum_i m_i c^2 \frac{1} {\sqrt {1 - \left( \frac{v_i}{c} \right)^2}}\rangle </math>

mit

• der Gesamtenergie <math>\textstyle E = \int \overline {T^{0}_{0}} \mathrm dV = \int \overline {T^{\alpha}_{\alpha}} \mathrm dV</math> des Systems
  • dem Energie-Impuls-Tensor <math>T^{\alpha \, \beta}</math> des Gesamtsystems aus Teilchen und Feldern
  • dem vierdimensionalen Index <math>\alpha = 0, 1, 2, 3</math>
  • der Spur <math>T^{\alpha}_{\alpha}</math>, wobei die Einsteinsche Summationskonvention verwendet wird.

Für kleine Geschwindigkeiten <math>v \ll c</math> ergibt sich die klassische Form des Virialsatzes für das Coulombpotential:

<math>E - \sum_i m_i c^2 = -\overline E_\mathrm{kin}</math>

wobei die Ruheenergien der Teilchen von der Gesamtenergie abgezogen werden.

Relativistische Versionen des Virialsatzes wurden schon von Chandrasekhar angewandt auf Weiße Zwerge. Er untersuchte auch Versionen in der allgemeinen Relativitätstheorie im Rahmen der Post-Newton-Näherung.<ref>S. Chandrasekhar: The Post-Newtonian Equations of Hydrodynamics in General Relativity. Astrophysical Journal, Band 142, 1965, S. 1488–1512, {{#invoke:Vorlage:bibcode|f|errHide=1|errNS=0|errClasses=editoronly error|errCat=Wikipedia:Vorlagenfehler/Parameter:bibcode}}.</ref><ref>George W. Collins: The Virial Theorem in Stellar Astrophysics. Pachart Press, 1978, Kapitel 2.</ref>

Literatur

• L. D. Landau, E. M. Lifschitz: Lehrbuch der theoretischen Physik, Bd. 1. Mechanik. 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.

• Wilhelm Brenig: Statistische Theorie der Wärme. 3. Auflage, Springer 1992, S. 144 f. (Virialsatz in statistischer Mechanik).
• George W. Collins: The Virial Theorem in Stellar Astrophysics. Pachart Press, 1978, Online.
• R. Becker: Theorie der Wärme. 1961, S. 85 (zum äußeren Virial).
• Albrecht Unsöld: Der neue Kosmos. Springer, 2. Aufl., 1974, S. 283, Ableitung und Bedeutung für die Berechnung des Aufbaus von Sternen. (Nicht im 1966er B.I.-Taschenbuch.)

Weblinks

• John Baez: The Virial Theorem Made Easy.

Einzelnachweise und Anmerkungen

<references />

{{#ifeq: s | p | | {{#if: 4218966-4sh85143793 | |

}} }}{{#ifeq:||{{#if: | [[Kategorie:Wikipedia:GND fehlt {{#invoke:Str|left|{{{GNDCheck}}}|7}}]] }}{{#if: | {{#if: | | }} }} }}{{#if: | {{#ifeq: 0 | 2 | | }} }}{{#if: | {{#ifeq: 0 | 2 | | }} }}{{#ifeq: s | p | {{#if: 4218966-4 | | {{#if: {{#statements:P227}} | | }} }} }}{{#ifeq: s | p | {{#if: 4218966-4 | {{#if: {{#invoke:Wikidata|pageId}} | {{#if: {{#statements:P227}} | | }} }} }} }}{{#ifeq: s | p | {{#if: sh85143793 | | {{#if: {{#statements:P244}} | | }} }} }}{{#ifeq: s | p | {{#if: sh85143793 | {{#if: {{#invoke:Wikidata|pageId}} | {{#if: {{#statements:P244}} | | }} }} }} }}{{#ifeq: s | p | {{#if: | | {{#if: {{#statements:P214}} | | }} }} }}{{#ifeq: s | p | {{#if: | {{#if: {{#invoke:Wikidata|pageId}} | {{#if: {{#statements:P214}} | | }} }} }} }}Vorlage:Wikidata-Registrierung