Dunkle Energie

Konventionen: Für Gravitation und Energiedichte $\rho$ gilt $c=1$; Abschnitt 3 setzt $c$ in die Entfernungen ein. Wir wählen $a_0=1$, also $1+z=1/a$. Punkte bezeichnen Ableitungen nach der kosmischen Zeit. Die Entfernungsbeispiele setzen räumliche Flachheit voraus. Dunkle Energie bezeichnet die Komponente, die in der allgemeinen Relativitätstheorie (ART) die späte Beschleunigung erklärt; der Name identifiziert kein entdecktes Teilchen.

1. Das Phänomen: die beschleunigte Expansion des Universums

Warum überraschte die Beschleunigung? Gewöhnliche Materie zieht an und sollte eine bereits vorhandene Expansion bremsen. Die Supernovateams fanden ferne Supernovae vom Typ Ia nach der Standardisierung lichtschwächer und damit weiter entfernt als in verzögert expandierenden Vergleichsmodellen. Die Analyse von 1998 und die unabhängige Analyse von 1999 erschlossen die Beschleunigung aus der Entfernung-Rotverschiebungs-Beziehung. Sie beobachteten nicht in Echtzeit, wie eine einzelne Galaxie schneller wurde. Staub, Quellenentwicklung und Auswahl müssen geprüft werden; die Übereinstimmung verschiedener kosmologischer Verfahren stärkt die Deutung.

Der Skalenfaktor $a(t)$ beschreibt Abstände frei beweglicher Beobachter im mittleren kosmischen Fluss. Homogenität und Isotropie ergeben die FLRW-Metrik; $k$ beschreibt die räumliche Krümmung. Ein idealer Flüssigkeits-Energie-Impuls-Tensor in Einsteins Gleichungen liefert die beiden Friedmann-Gleichungen. Hier enthalten $\rho ,p$ die gesondert geschriebene kosmologische Konstante nicht:

$$\begin{array}{c}d{s}^2=-d{t}^2+a^2(t)\left[\frac{d{r}^2}{1-k{r}^2}+r^{2}d{\Omega }^2\right],H=\frac{\dot{a}}{a}.\\ H^2=\frac{8\pi G}{3}\rho -\frac{k}{a^2}+\frac{\Lambda }{3},\frac{\ddot{a}}{a}=-\frac{4\pi G}{3}(\rho +3p)+\frac{\Lambda }{3}.\end{array}$$

Druck wirkt in der ART gravitativ: In der Beschleunigungsgleichung steht $\rho +3p$, nicht nur die Dichte. Positive Energiedichte mit hinreichend negativem Druck kann deshalb beschleunigen, ohne negative Masse zu sein. Man kann $\Lambda$ auch als Flüssigkeit mit ${\rho }_{\Lambda }=\Lambda /(8\pi G)$ und $p_{\Lambda }=-{\rho }_{\Lambda }$ behandeln; dann entfallen die expliziten $\Lambda$-Terme. Beides zugleich wäre eine Doppelzählung. Bei positiver Gesamtdichte verlangt Beschleunigung $w_{\mathrm{total}}=p_{\mathrm{total}}/{\rho }_{\mathrm{total}}<-1/3$. Eine kleine Komponente mit $w<-1/3$ reicht bei dominierender Materie nicht aus.

Der Verzögerungsparameter $q$ quantifiziert diese Unterscheidung. Die Dichteanteile ${\Omega }_i$ beziehen sich jeweils auf $3H^2/(8\pi G)$ zur selben Epoche; dunkle Energie wird hier als Flüssigkeit behandelt. Beschleunigung bedeutet $q<0$, aber nicht zwingend ein wachsendes $H$. Abstände können immer schneller wachsen, während ihre relative Wachstumsrate sinkt.

$$q\equiv -\frac{\ddot{a}}{a{H}^2}=\frac{1}{2}{\Omega }_m+{\Omega }_r+\frac{1}{2}(1+3w){\Omega }_{\mathrm{DE}},\dot{H}=\frac{\ddot{a}}{a}-H^2.$$

Rechenbeispiel: Ein flaches Materie-$\Lambda$-Modell mit ${\Omega }_{m0}=0.30$, ${\Omega }_{\Lambda 0}=0.70$ und vernachlässigbarer Strahlung dient nur der Illustration. Heute ist $q_0=0.30/2-0.70=-0.55$. Wegen ${\rho }_m\propto a^{-3}$ beginnt die Beschleunigung bei ${\rho }_m=2{\rho }_{\Lambda }$: $a_{\mathrm{acc}}=(0.30/1.40{)}^{1/3}\simeq 0.598$, entsprechend $z_{\mathrm{acc}}\simeq 0.67$. Gleiche Materie- und Vakuumdichten treten erst bei $a\simeq 0.754$ auf. Dichtedominanz und Beschleunigungsbeginn sind also unterschiedliche Schwellen.

Beschleunigung hat eine eigene Schwelle. Im illustrativen flachen Modell $({\Omega }_{m0},{\Omega }_{\Lambda 0})=(0.30,0.70)$ erreicht $q(a)$ bei $a\simeq 0.598$ null (durchgezogene Hilfslinie); Dichtegleichheit folgt bei $a\simeq 0.754$ (gestrichelt). Die Kurve ist eine Modellrechnung, kein Datenfit.

2. Zustandsgleichung und Parametrisierungen

Wir definieren $w=p_{\mathrm{DE}}/{\rho }_{\mathrm{DE}}$. Für eine homogene Komponente ohne Energieaustausch mit anderen Spezies liefert die lokale Energie-Impuls-Erhaltung die Kontinuitätsgleichung. Division durch $H\rho$ und $d\ln a=Hdt$ führen zur integrierbaren Form:

$$\begin{array}{c}{\dot{\rho }}_{\mathrm{DE}}+3H({\rho }_{\mathrm{DE}}+p_{\mathrm{DE}})=0,\frac{d\ln {\rho }_{\mathrm{DE}}}{d\ln a}=-3(1+w).\\ {\rho }_{\mathrm{DE}}(a)={\rho }_{\mathrm{DE},0}\exp \left[-3{\int }_1^a(1+w(a^{\prime }))d\ln a^{\prime }\right].\end{array}$$

Für konstantes $w$ folgt $\rho \propto a^{-3(1+w)}$. Materie mit $w\simeq 0$ verdünnt sich mit dem Volumen; Strahlung mit $w=1/3$ verliert zusätzlich Energie durch Rotverschiebung. Bei $w=-1$ bleibt die Dichte konstant. Dasselbe folgt aus $d(\rho V)=-pdV$ für ein mitbewegtes Volumen $V\propto a^3$. Negativer Druck erlaubt wachsende Energie bei wachsendem Volumen. Das verletzt die lokale Erhaltung nicht: Eine allgemein expandierende Raumzeit besitzt nicht notwendigerweise die global erhaltene Gesamtenergie einer statischen Raumzeit.

Die Chevallier–Polarski–Linder-Form (CPL) beschreibt mögliche Zeitabhängigkeit durch den heutigen Wert $w_0$ und den Änderungsparameter $w_a$. Einsetzen in das Integral liefert die folgende Dichte; der Exponentialfaktor darf nicht fehlen. Siehe die Originalarbeiten von Chevallier–Polarski und Linder.

$$\begin{array}{c}w(a)=w_0+w_a(1-a)=w_0+\frac{w_{a}z}{1+z},\\ \frac{{\rho }_{\mathrm{DE}}(z)}{{\rho }_{\mathrm{DE},0}}=(1+z{)}^{3(1+w_0+w_a)}\exp \left[-\frac{3w_{a}z}{1+z}\right].\end{array}$$

CPL ist eine Anpassungsfamilie, keine mikroskopische Theorie. $(-1,0)$ entspricht $\Lambda$. Die Form kann $w=-1$ kreuzen, was ein einzelnes minimal gekoppeltes kanonisches Skalarfeld positiver Energie nicht kann. Bei $w<-1$ verletzt eine ideale Flüssigkeit positiver Dichte mit $\rho +p<0$ die Nullenergiebedingung; eine effektive Anpassung bedeutet aber nicht zwangsläufig ein fundamentales Geistfeld. Für $w_a\ne 0$ divergiert CPL bei $a\to \infty$ und rechtfertigt keine Vorhersage der fernen Zukunft.

3. Beobachtungsgrößen des kosmologischen Hintergrunds

Teleskope messen Flüsse, Winkel, Spektren und Korrelationen. Ein Modell übersetzt sie in Entfernungen und Expansionsraten; kein Instrument liest die Dichte dunkler Energie direkt ab. Jedes Verfahren liefert einen anderen Teil dieser Übersetzung.

3.1 Hubble-Expansion

Mit ${\rho }_{c0}=3H_0^2/(8\pi G)$ definieren wir ${\Omega }_{i0}={\rho }_{i0}/{\rho }_{c0}$. Bei räumlicher Flachheit ist ihre Summe eins. Separat erhaltene Materie, Strahlung und dunkle Energie ergeben:

$$\frac{H^2(z)}{H_0^2}={\Omega }_{r0}(1+z{)}^4+{\Omega }_{m0}(1+z{)}^3+{\Omega }_{\mathrm{DE},0}\frac{{\rho }_{\mathrm{DE}}(z)}{{\rho }_{\mathrm{DE},0}}.$$

Bei $z=1$ liefert das Beispiel $H/H_0=\sqrt{0.30\times 8+0.70}\simeq 1.761$. Ersetzt man $w=-1$ durch konstantes $w=-0.8$ bei gleichen heutigen Anteilen und gleichem $H_0$, wird der Beitrag dunkler Energie $0.70\times 2^{0.6}$ und $H/H_0\simeq 1.860$. Ihre Dichte ändert sich dort um rund 52%, die gesamte Expansionsrate aber nur um 5.7%. Materie schwächt den Effekt ab, und integrierte Entfernungen glätten ihn zusätzlich. Zeitabhängigkeit ist daher schwer zu messen.

3.2 Entfernungen (SNe Ia, BAO, CMB)

Für einen radialen Lichtstrahl verbindet $d{s}^2=0$ die Entfernung mit dem Integral über $1/H$. Im flachen Raum hängen die mitbewegte Entfernung $\chi$, die Leuchtkraftentfernung $d_L$ und die Winkeldurchmesserentfernung $D_A$ wie folgt zusammen. Bei Krümmung ersetzt eine krümmungsabhängige transversale Entfernung $\chi$ in den letzten beiden Beziehungen. Die übliche Entfernungsdualität setzt außerdem Nullgeodäten für Photonen und Erhaltung ihrer Anzahl voraus.

$$\begin{array}{c}\chi (z)=c{\int }_0^z\frac{d{z}^{\prime }}{H(z^{\prime })},d_L=(1+z)\chi ,D_A=\frac{\chi }{1+z},\\ \mu =5{\log }_{10}\left(\frac{d_L}{10\mathrm{pc}}\right),D_V={\left[(1+z{)}^2{D}_A^2\frac{cz}{H(z)}\right]}^{1/3}.\end{array}$$

Kalibrierte Leuchtkraft $L$ und Fluss $F$ erfüllen $F=L/(4\pi d_L^2)$. SNe Ia sind standardisierbare, nicht identische Kerzen: Lichtkurvenform, Farbe und Wirtsgalaxie erfordern Korrekturen. Der Modul $\mu$ drückt die Entfernung logarithmisch aus. Ohne absolute Kalibrierung liefert eine Supernovastichprobe vorwiegend relative Entfernungen; absolute Leuchtkraft und $H_0$ sind entartet. Eine Verschiebung von 0.10 Magnituden verändert die Entfernung um $10^{0.10/5}-1\simeq 4.7\%$: Kalibrierung zählt.

Baryonische akustische Oszillationen (BAO) nutzen einen statistisch bevorzugten Abstand, den Schallwellen im frühen Plasma hinterließen. Der mitbewegte Schallhorizont $r_d$ wird zur Epoche des baryonischen Strahlungswiderstands bestimmt. Quer zur Sichtlinie messen BAO $D_M/r_d$, längs davon $D_H/r_d$ mit $D_H=c/H$. Für kleine Trennungen gilt schematisch:

$${\theta }_{\mathrm{BAO}}\simeq \frac{r_d}{D_M(z)},\Delta z_{\mathrm{BAO}}\simeq \frac{H(z)r_d}{c},D_M=(1+z)D_A.$$

Ein statistischer Maßstab, zwei Richtungen. Der orange Querabstand verbindet $r_d$ mit $\theta$ und $D_M$; der blaue Längsabstand verbindet $r_d$ mit $\Delta z$ und $H$. Die Geometrie ist schematisch, der Maßstab mitbewegt.

$D_V$ fasst eine isotrope BAO-Auswertung zusammen; anisotrope Analysen behalten radiale und transversale Information. BAO sind ein Korrelationsmerkmal, keine starren Stäbe zwischen einzelnen Galaxien. Die kosmische Mikrowellenhintergrundstrahlung (CMB) beschränkt frühe Dichten und akustische Skalen und kalibriert $r_d$. Späte Parameter folgen weiterhin nur unter Annahmen über die frühe Physik: Andere frühe Expansion kann den Maßstab ändern. Die Planck-Parameteranalyse zeigt solche bedingten Schlussfolgerungen.

4. Wachstum kosmischer Strukturen

Expansion ist nur die Hälfte des Tests. Materieüberdichten $\delta =\delta {\rho }_m/{\rho }_m$ wachsen durch Anziehung. Für lineare, unterhalb des Horizonts liegende Störungen druckloser Materie in der ART, bei glatter dunkler Energie und vernachlässigbarer Strahlung, ergeben Kontinuitäts-, Euler- und Poisson-Gleichung:

$$\ddot{\delta }+2H\dot{\delta }-4\pi G{\rho }_m\delta =0,f\equiv \frac{d\ln D}{d\ln a}.$$

$\delta \propto D(a)$ beschreibt die wachsende Lösung. $2H\dot{\delta }$ dämpft das Wachstum, der Dichteterm treibt es an. In einer flachen reinen Materieära gilt $D\propto a$; glatte dunkle Energie verlangsamt es relativ zu $a$. $f\simeq {\Omega }_m(a{)}^{0.55}$ ist eine Näherung nahe $\Lambda$CDM, keine allgemeine Identität. Rotverschiebungsraum-Verzerrungen beschränken etwa $f{\sigma }_8$. ${\sigma }_8$ ist die quadratische mittlere lineare Dichtefluktuation in Kugeln mit Radius $8h^{-1}$ Mpc, wobei $h=H_0/(100\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1})$.

Galaxienbewegungen reagieren auf das Newtonsche Potential $\Psi$, während Linseneffekte $\Phi +\Psi$ prüfen. Mit der Fourierkonvention ${\nabla }^2\to -k^2$ parametrisiert man Abweichungen von einem Referenzmodell oft durch:

$$k^2\Psi =-4\pi G_{\mathrm{eff}}(k,a)a^2{\rho }_m\Delta ,\eta =\frac{\Phi }{\Psi }.$$

$\Delta$ ist die mitbewegte Materiedichtestörung; $k$ bezeichnet hier die Wellenzahl, nicht die Hintergrundkrümmung. In der ART ohne nennenswerte anisotrope Spannungen gilt $G_{\mathrm{eff}}=G$ und $\eta =1$. Neutrinospannungen oder andere nichtideale Komponenten können auch in der ART unterschiedliche Potentiale erzeugen. Vor einer Deutung als modifizierte Gravitation müssen Materiephysik, Galaxienbias und systematische Messfehler geprüft werden.

5. Mikroskopische Modelle im Überblick

Ein Hintergrundfit für $w$ identifiziert weder ein Feld noch eine Gravitationstheorie eindeutig. Modelle mit gleichem $H(z)$ können verschiedene Störungen, Stabilitätsbedingungen und lokale Gravitationseffekte haben.

5.1 Kosmologische Konstante $\Lambda$

Ein Lorentz-invariantes Vakuum hat einen zum metrischen Tensor proportionalen Energie-Impuls-Tensor. Sein Druck ist negativ, obwohl seine Energiedichte positiv sein kann:

$$T_{\mu \nu }^{(\Lambda )}=-{\rho }_{\Lambda }{g}_{\mu \nu },p_{\Lambda }=-{\rho }_{\Lambda },{\rho }_{\Lambda }=\frac{\Lambda }{8\pi G}.$$

Das Problem der kosmologischen Konstante ist die Kleinheit des renormierten gesamten Vakuumbeitrags einschließlich des gravitativen Konstantenterms und der Quantenkorrekturen. Eine Abschätzung mit Abschneideskala ist keine gemessene Vorhersage; die Abweichung hängt von Energieskala und Regularisierung ab. Ein freies $\Lambda$ passt den Wert an, erklärt aber nicht seine Stabilität gegenüber großen Korrekturen. Das Problem bleibt auch bei künftig bestätigtem $w=-1$; siehe die theoretische Diskussion von Peebles und Ratra.

Davon verschieden ist die Koinzidenzfrage: Warum beobachten wir eine Epoche mit ähnlichen Materie- und Vakuumdichten, obwohl ihr Verhältnis wie $a^3$ variiert? Sie motiviert dynamische Modelle, beweist sie aber nicht.

5.2 Quintessenz

Quintessenz führt ein kanonisches, minimal gekoppeltes Skalarfeld ein. Für ein homogenes Feld bei Signatur $(-,+,+,+)$ ergeben kinetische Energie $K$ und Potential $V$:

$$\begin{array}{c}\mathcal{L}=-\frac{1}{2}{g}^{\mu \nu }{\partial }_{\mu }\varphi {\partial }_{\nu }\varphi -V(\varphi ),\\ K=\frac{1}{2}{\dot{\varphi }}^2,{\rho }_{\varphi }=K+V,p_{\varphi }=K-V,w_{\varphi }=\frac{K-V}{K+V}.\end{array}$$

Variation der Wirkung liefert die Feldgleichung; für $\dot{\varphi }\ne 0$ folgt sie auch aus der Kontinuitätsgleichung mit ${\rho }_{\varphi },p_{\varphi }$. Der Term $3H\dot{\varphi }$ wirkt als kosmologische Reibung.

$$\ddot{\varphi }+3H\dot{\varphi }+V_{,\varphi }=0,{\rho }_{\varphi }+3p_{\varphi }=4K-2V.$$

Dominiert das Feld, verlangt Beschleunigung $K<V/2$. Bei $K\ll V$ ähnelt es einer kosmologischen Konstante; bei $V=0$ gilt $w=1$, also beschleunigt ein Skalar nicht automatisch. Für $K,V\ge 0$ und positive Dichte ist $-1\le w\le 1$. Inverse Potenz- oder Exponentialpotentiale können für passende Parameter nachführende Lösungen besitzen, die einige Anfangsbedingungen weniger wichtig machen, aber nicht jede Feinabstimmung beseitigen. Potential, Kopplungen und Störungsstabilität bleiben physikalische Annahmen, die ein CPL-Fit nicht herleitet.

5.3 Modifizierte Gravitation

Modifizierte Gravitation ändert die Beziehung zwischen Materie und Krümmung. Metrische $f(R)$-Theorien führen im Allgemeinen einen zusätzlichen skalaren Freiheitsgrad ein. Chamäleon-Abschirmung kann lokale Effekte unterdrücken; Vainshtein-Abschirmung tritt in anderen Theorien mit Ableitungswechselwirkungen auf. Keine gilt universell. Horndeski- und degenerierte Skalar-Tensor-Theorien höherer Ordnung strukturieren weitere Möglichkeiten. GW170817 und sein Gammastrahlen-Gegenstück beschränken im beobachteten Bereich den Unterschied zwischen Gravitations- und Lichtgeschwindigkeit stark. Unter ihren Voraussetzungen schließen sie viele Wechselwirkungen aus, aber nicht jede modifizierte Gravitation.

6. Spannungen und Konsistenztests

Die Hubble-Spannung vergleicht modellabhängige frühkosmische $H_0$-Schlüsse mit lokalen Entfernungsleitern. Die Planck-Analyse ist ein Beispiel für Erstere, keine direkte lokale Expansionsmessung. Späte dunkle Energie zu ändern löst das Problem nicht automatisch, da BAO und Supernovaentfernungen ebenfalls passen müssen. Frühe dunkle Energie betrifft eine frühere Epoche und ist kein anderer Name für die heute beschleunigende Komponente.

$S_8={\sigma }_8\sqrt{{\Omega }_m/0.3}$ beschreibt eine wichtige Entartung von Linsenmessungen, aber nicht ihren gesamten Informationsgehalt. Unterschiede gegenüber CMB-Schlüssen hängen von Datensatz, nichtlinearer Modellierung, baryonischer Rückwirkung, Rotverschiebungskalibrierung und intrinsischen Galaxienausrichtungen ab. Vergleiche vollständiger Likelihoods und gemeinsamer Störparameter sind aussagekräftiger, als jede Spannung als unabhängigen Beleg neuer Physik zu zählen.

7. Ablauf der Parameterschätzung: praktische Hinweise

1. Modell und Priorverteilungen festlegen: Krümmung, Neutrinomassen, frühe Physik und mögliche Klumpung dunkler Energie; $w_0,w_a$ allein sind keine vollständige Theorie.
2. $H$, Entfernungen, frühe akustische Skalen und lineare Störungen konsistent berechnen. Präzisionsanalysen benötigen einen Boltzmann-Löser statt nur eines Wachstumsindex oder komprimierter CMB-Entfernungsinformation.
3. Messstatistiken mit Auswahl, Kalibrierung und Störparametern vorhersagen. Kovarianzen überlappender Daten berücksichtigen; scheinbar unabhängige Multiplikation korrelierter Likelihoods zählt Information doppelt.
4. Posteriorverteilungen bestimmen und Anpassungsgüte, Priorabhängigkeit sowie Stabilität beim Weglassen einzelner Verfahren testen. Zusätzliche Parameter verbessern einen Fit, ohne allein neue Physik nachzuweisen.

8. Stand der Forschung

Befunde geprüft bis 13. Oktober 2026. Die DESI-DR2-BAO-Arbeit von 2025 beschreibt BAO allein gut durch flaches $\Lambda$CDM. In der CPL-Erweiterung bevorzugten BAO+CMB zeitabhängige dunkle Energie mit $3.1\sigma$; je nach zusätzlicher Supernovastichprobe ergaben sich $2.8$–$4.2\sigma$. Das sind bestimmte Vergleiche mit $(w_0,w_a)=(-1,0)$, keine direkte Entdeckung eines Skalarfelds.

Die DES-Analyse mehrerer Verfahren von 2026 kombiniert Geometrie und Wachstum. Innerhalb von CPL ergaben sich $2.2\sigma$ Abweichung für DES allein, $2.3\sigma$ mit DESI DR2 BAO und $3.0\sigma$ mit primärer CMB-Information. Die Kombinationen überschneiden sich mit anderen Analysen; ihre Signifikanzen sind keine unabhängig multiplizierbaren Bestätigungen.

Die DESI-Lyα-Analyse der vollständigen Korrelationsform vom Juli 2026, im August überarbeitet, ergänzt eine Alcock–Paczyński-Geometriemessung bei effektiver Rotverschiebung $2.33$. Der neue Anker rückt in Richtung $\Lambda$CDM. Die Arbeit nennt $2.7\sigma$ CPL-Präferenz mit DESI+CMB und $3.2\sigma$ mit ihrem zusätzlichen Supernovadatensatz. Der Befund entwickelt sich somit mit der verwendeten Information.

Eine Zahl von Standardabweichungen hängt von Modell, Likelihood und systematischer Fehlerbehandlung ab; sie ist nicht die Wahrscheinlichkeit einer zeitabhängigen dunklen Energie. Offen bleibt, ob unabhängige Kalibrierungen und Wachstumstests das Muster erhalten, welche Theorie es stabil erzeugen könnte und ob einfachere Erklärungen bestehen bleiben. Die Hinweise verdienen Untersuchung; der mikroskopische Ursprung der Beschleunigung ist ungelöst.

Häufige Missverständnisse: Dunkle Energie ist nicht die in Halos verklumpende dunkle Materie. Negativer Druck bedeutet weder negative Energiedichte noch einen Wind aus einem kosmischen Zentrum. Gebundene Atome, Sonnensysteme und Galaxien folgen nicht einfach dem Skalenfaktor. Weder heutige Beschleunigung noch ein CPL-Fit bei endlicher Rotverschiebung bestimmt das endgültige Schicksal des Universums: Konstantes Vakuum, veränderliches Potential und künftige Instabilität erlauben unterschiedliche Fortsetzungen.