← Home | Kann ein „schiefes“ Universum wie ein symmetrisches aussehen?
Testprojekt · Lokales KI-Ökosystem · Berlin · Juli 2026

Kann ein „schiefes“ Universum wie ein symmetrisches aussehen?

Eine offene Frage aus einem Dezember-2025-Paper von Lavinia Heisenberg (Universität Heidelberg): Kann eine Raumzeit, die selbst nicht überall gleich aussieht, trotzdem Einsteins Gravitationsgleichungen so erfüllen, dass sie sich für die Gravitation vollkommen gleichförmig verhält? Symbolische Rechnung mit SageMath/SageManifolds auf dem Mac Studio — Stand Juli 2026: der Fall A(r)=1 ist vollständig klassifiziert (siehe unten), die Untersuchung des allgemeinen Falls A(r) läuft.

Angestoßen durch den ZEIT-Podcast „Alles gesagt?“ mit Lavinia Heisenberg (Juli 2026) · SageMath 10.9 / SageManifolds · Mac Studio M3 Ultra (256 GB) · Juli 2026

Worum es geht — ohne Formeln

Die moderne Kosmologie beruht auf einer Annahme, die man das Kosmologische Prinzip nennt: Auf großen Skalen sieht das Universum überall gleich aus (homogen) und in jede Richtung gleich aus (isotrop). Das ist keine bewiesene Tatsache, sondern eine gut bestätigte Annahme, die die Mathematik der Kosmologie erst handhabbar macht.

Lavinia Heisenberg hat in einem Paper vom Dezember 2025 (arXiv:2512.16934) die Grundlagen dieser Annahme neu durchdacht und dabei eine Frage aufgeworfen, die zunächst harmlos klingt, aber nicht trivial ist. Einsteins Feldgleichungen verbinden zwei Größen: die Metrik (das vollständige geometrische „Maßsystem“ der Raumzeit — sie legt fest, wie lang jede Strecke in jede Richtung an jedem Ort ist) und den Einsteintensor (der Teil dieser Geometrie, der tatsächlich mit Materie und Energie wechselwirkt). Wenn die Metrik selbst homogen und isotrop ist, folgt zwingend, dass auch der Einsteintensor es ist — das ist ein bewiesener, unstrittiger Sachverhalt.

Heisenbergs Frage geht in die andere Richtung, und die ist nicht automatisch beantwortet: Kann der Einsteintensor homogen und isotrop sein, obwohl die zugrundeliegende Metrik es nicht ist? Anders gesagt: Könnte es eine Raumzeit geben, die selbst schief, ungleichmäßig oder mit einer Vorzugsrichtung ausgestattet ist — die sich aber, gemessen an dem, was tatsächlich auf die Materie wirkt, exakt so verhält wie unser gewöhnliches, symmetrisches Universum? Heisenberg schreibt selbst, dass sie wegen der algebraischen Komplexität kein konkretes Beispiel konstruieren konnte, und nennt die Frage ein offenes Problem.

Warum das mehr als eine Spielerei ist: Wäre die Antwort „ja“, gäbe es prinzipiell verborgene geometrische Struktur, die sich hinter dem beobachtbaren, perfekt symmetrischen kosmischen Hintergrund verstecken könnte — unsichtbar für jede Messung, die nur über die Feldgleichungen läuft. Wäre die Antwort „nein“, würde das bestätigen, dass die beobachtete Symmetrie tatsächlich auf eine symmetrische Geometrie schließen lässt, nicht nur auf ein symmetrisches Nebenprodukt.

Praezisierungen (14. Juli 2026, lokales Komitee qwen3:235b-a22b):

Geschaerfte Motivation: Die ursprüngliche Fragestellung, ob verborgene geometrische Strukturen hinter dem beobachtbaren, perfekt symmetrischen kosmischen Hintergrund existieren könnten, erfordert eine präzisere Einordnung: Solche Strukturen wären nicht unsichtbar für kinematische Beobachtungen (Geodäten von Testkörpern und Licht), die stets die volle Metrik g erfassen, sondern lediglich vor der Schlussfolgerung „symmetrische Materie impliziert symmetrische Geometrie“ verborgen. Im Bianchi-I-Staubmodell würde etwaige Scherung unmittelbar in der Anisotropie der Rotverschiebung des kosmischen Mikrowellenhintergrunds sichtbar – ein Phänomen, das durch den klassischen Satz von Ehlers, Geren & Sachs (1968) erklärt wird: Isotropie der Hintergrundstrahlung für alle mitbewegten Beobachter erzwingt bereits die FLRW-Symmetrie. Die ursprüngliche Frage bleibt somit bestehen, ob eine geometrische Degeneriertheit existieren könnte, die durch die Feldgleichungen nicht erfasst wird, jedoch durch geodätische Effekte nachweisbar wäre.

Invariante Formulierung: Im Kern vergleicht das Projekt zwei Lie-Algebren von Vektorfeldern: die Killing-Algebra KV(g) = {ξ : Lξ g = 0} und die Materie-Kollineationsalgebra MC(g) = {ξ : Lξ G = 0}. Stets gilt KV(g) ⊆ MC(g), doch die entscheidende Frage lautet, wann diese Inklusion strikt ist und wie groß der Überschuss MC(g)  ∖  KV(g) werden kann. Elemente dieser Differenzmenge werden in der Literatur als „proper matter collineations“ bezeichnet. Ist der Einstein-Tensor G als Bilinearform nichtdegeneriert, so ist MC(g) isomorph zur Killing-Algebra der zugehörigen Pseudometrik G – endlichdimensional mit einer maximalen Dimension von 10. Degeneriert G hingegen, wie im projektgefundenen Fall pt = 0 mit Gθθ = 0, kann MC(g) unendlichdimensional werden, was die zentrale Signatur unserer Analyse darstellt.

Warum die CAS-Identitaet strukturell gilt: Für statische, sphärisch symmetrische Metriken g mit diagonalem Einstein-Tensor G (nur r-abhängige Komponenten) und rein winkelartigem, statischem Vektorfeld ξ (ξt = ξr = 0) verschwindet die Lie-Ableitung Lξ Gμν komponentenweise für alle Indizes außer θθ und φφ. Diese beiden Komponenten tragen einen gemeinsamen Faktor Gθθ = Gφφ, sodass Lξ G = 0 äquivalent ist zur Bedingung Gθθ = 0 (d.h. pt = 0). Die Invarianz der verbleibenden Komponenten (Energiedichte ρ, radialer Druck pr) folgt automatisch aus der Diagonalstruktur und der Symmetrie von g, ohne zusätzliche Feldgleichungsbedingungen. Dies erklärt strukturell die CAS-Identität Gθθ = -rW/(4A² C²), die somit nicht als Zufall, sondern als direkte Konsequenz der geometrischen Rahmenbedingungen auftritt.

Definition "stark nichttrivial": Bisher wurde nichttrivial definiert durch Lξ g ≠ 0 , Lξ G = 0 und G ≠ 0 . Die Verschärfung "stark nichttrivial" erfordert zusätzlich, dass die Einschränkung von G auf die von ξ deformierten Richtungen nicht verschwindet – es liegt kein "Nullblock-Mechanismus" vor, bei dem die Deformation exakt dort wirkt, wo G nichts trägt. Daraus folgt der präzisierte Klassifikationssatz: Im statischen A=1 -System existiert keine stark nichttriviale Konfiguration; die pt=0 -Klasse ist der vollständige Katalog der blockweise-trivialen Fälle.

Aktueller Stand (Juli 2026, nach Schritten 69–74): Das statische Programm ist vollständig klassifiziert — mit einer wesentlichen Korrektur früherer Ergebnisse. (1) Eich-Lemma: Der Kreuzterm B(r) ist im statischen Fall immer reine Eichung (T=t−f(r), f'=B/A; Konvention: ds² enthält 2B dt dr, also g01=B). (2) Hauptresultat: Die Lösungsklasse der Projektfrage im statischen Diagonalsystem mit rein winkelartigem Deformationsfluss ist exakt die Klasse aller Metriken mit verschwindendem Tangentialdruck pt=0 — sie enthält die Defizitwinkel-Familie und die volle Barriola–Vilenkin-Lösung mit Massenterm; der frühere Literatur-Dissens ist damit aufgelöst (vom Komitee selbst nachgerechnet). (3) Klassifikationssatz (A=1): Für radiale Flussanteile (Fr≠0) existiert keine nichttriviale Konfiguration; die einzige Kandidatenfamilie C=1/(1−kr²) hat einen Killing-Fluss (LT3g=0) und ist trivial. (4) Korrektur: Die frühere Shooting-Trajektorie (Abschnitte unten, Schritte 27–58 samt Nachträgen) erfüllt die 12er-Bedingung des vollen Systems nicht und ist keine gültige LT3G=0-Konfiguration. Details: Schritte 69–74; die betroffenen Abschnitte bleiben transparent dokumentiert, sind aber überholt.
Update (10. Juli 2026, Schritt 75): Im offenen Fr≠0-Zweig bei allgemeinem A(r) (siehe Punkt 3 oben) wurde eine erste konkrete, vollständig verifizierte nichttriviale Konfiguration gefunden: A(r)=r², C(r)=3 (konstant), mit radialem Fluss Fr(r)=λr², Fth=0. CAS-Kette: E0=0 und die volle 12er-Restbedingung R=0 sind beide erfüllt (unabhängig reproduziert), LT3g≠0 (kein Killing-Fluss) und ein nichtverschwindender Riemann-Tensor (echt gekrümmt, keine flache Verkleidung) wurden direkt auf der konkreten Metrik nachgerechnet, ebenso LT3G=0 in allen zehn Komponenten unabhängig von der allgemeinen Herleitung. Komitee-Review mit dokumentiertem Dissens: ein Gutachten (qwen3) lehnte formal ab, bestätigte in der eigenen Rechnung aber den Kernbefund und identifizierte einen einzigen tragenden Einwand (fehlender unabhängiger Vier-Komponenten-Nachweis), der durch die genannte direkte Nachrechnung inzwischen geklärt ist; das zweite Gutachten (deepseek-r1:70b) blieb methodisch zu schwach für nennenswertes Gewicht. Offen: ob A=r² ein isolierter Punkt oder Teil einer größeren Lösungsfamilie ist (die zugrundeliegende Bedingung ist eine ODE 3. Ordnung, generisch eine Drei-Parameter-Familie). Details: Schritt 75.
Update (11. Juli 2026, Schritt 76): Der Punkt A=r² ist kein Einzelfall: Die W2=0-Bedingung reduziert sich (Skalensymmetrie A→μA, dann Übergang zu x=ln r) auf eine autonome ODE 2. Ordnung mit drei Fixpunkten: h=0 (A=1, Sattel), h=−1 (entartet, C=0), h=2 (A=r², stabile Spirale). Ein CAS-Test bestätigt: die gesamte 3-Parameter-Lösungsfamilie erfüllt tatsächlich E0=0 UND die volle 12er-Bedingung R=0 — keine Scheinfamilie. Numerische Integration (breiter Startbereich) zeigt: A~r² ist ein robuster, aber nicht universeller Attraktor für r→∞; die Familie erfüllt dort die Energiebedingungen (p_r≡0 identisch, ρ>0, p_t>0). Bei r→0 dagegen entwickelt sich eine echte Krümmungssingularität — kein sanfter Grenzwert (korrigiert in Schritt 79: nicht der Ricci-Skalar divergiert, wie hier zunächst stand, sondern der Kretschmann-Skalar, K~r⁻⁶; der Ricci-Skalar verschwindet tatsächlich). Komitee-Review (zwei Runden, dokumentierter Dissens): ein Gutachten identifizierte dabei einen echten Fehler in einer Zwischenversion (fälschliche Gleichsetzung von Energiebedingungsverletzung und numerischer Instabilität, korrigiert). Details: Schritt 76.
Update (11. Juli 2026, Schritt 77): Der parallel untersuchte AC=1-Zweig (Schwarzschild-Klasse-Ansatz, Schritt 75b) erwies sich als überbestimmt (zwei unabhängige Bedingungen für eine Funktion A(r)). Die einzige gefundene geschlossene Lösungsfamilie, A=1+K/r, ist — nach Elimination und Invarianztest — exakt die Schwarzschild-Metrik: reines Vakuum (G≡0), damit trivial im Projektsinn, obwohl der gefundene Fluss nicht-Killing ist. Ein informatives Negativresultat, sauber identifiziert. Auch der dritte Faktor (kubisch in A') erwies sich bei genauerer Prüfung als Scheinlösung (nicht invariant unter dem E0-Fluss) — der gesamte AC=1-Zweig ist damit vollständig und negativ klassifiziert. Details: Schritt 77.

Die Idee, symbolisch anzugreifen

Diese Frage lässt sich nicht durch Nachdenken allein beantworten — sie erfordert, tatsächlich eine Metrik hinzuschreiben und mit Computeralgebra zu prüfen, ob und wie ihr Einsteintensor symmetrisch werden kann, ohne dass die Metrik selbst es ist. Genau das haben wir versucht, mit SageMath (genauer: dessen Erweiterung SageManifolds, die symbolische Tensorrechnung auf gekrümmten Räumen beherrscht), vollständig lokal auf dem Mac Studio, ohne Cloud-Zugriff.

Metrik g vollständige Geometrie Einsteintensor G wirkt auf Materie g symmetrisch ⇒ G symmetrisch (bewiesen, unstrittig) G symmetrisch ⇒ g symmetrisch? (Heisenbergs offene Frage — nicht bewiesen)
Die Hinrichtung ist gesichert. Die Rückrichtung ist die offene Frage, die Heisenberg stellt.

Als konkreten Testfall haben wir eine kugelsymmetrische Metrik gewählt, mit einer Materiequelle, die symmetrisch sein soll, aber nicht notwendig verschwindet — das ist die einfachste Klasse, in der die Frage überhaupt sinnvoll gestellt werden kann. (Der reine Vakuumfall, dazu unten mehr, stellt sich als eine andere, trivialere Frage heraus.) Sie ist groß genug, um nicht automatisch auf die triviale (bereits als symmetrisch bekannte) Lösung zu kollabieren.

Vier weitere Einschränkungen kamen hinzu — drei davon sind echte Koordinatenwahlen ohne Informationsverlust, die vierte ist etwas anderes und wurde später zum eigenständigen Forschungsgegenstand. Zuerst der statische Fall (keine Zeitabhängigkeit) — die in der Literatur übliche Reihenfolge. Dann die diagonale Form (kein t-r-Kreuzterm): inzwischen kein Ansatz mehr, sondern ein bewiesenes Lemma (Schritt 69) — im gesamten statischen Fall entfernt eine reine Zeitverschiebung T=t−f(r) den Kreuzterm B(r) exakt, ohne jede Einschränkung. Dann die Areal-Radius-Eichung (Winkelkomponente = r²) — ebenfalls eine echte, verlustfreie Koordinatenwahl. Und schließlich A(r)=1: anders als die drei vorigen war das nie eine bewiesene Eichung, sondern eine plausible, aber ungeprüfte Vereinfachung. Der gesamte im Folgenden geschilderte Projektverlauf — von der Potenzgesetz-Analyse bis zur numerischen Shooting-Trajektorie — fand innerhalb dieser vierten Einschränkung statt, und die spätere vollständige Klassifikation (Schritte 69–74) zeigt: Für A=1 ist der dort untersuchte nichttriviale Zweig tatsächlich leer. Ob das an A=1 liegt oder eine echte Eigenschaft der Fragestellung ist, klärt erst die Lockerung genau dieser einen Einschränkung — die laufende Untersuchung des allgemeinen Falls A(r), seit Juli 2026: erster verifizierter nichttrivialer Befund: W₂=0-Familie (Schritte 75–79); AC=1-Zweig negativ klassifiziert (Schritt 77).

Allgemeine kugelsymm. Metrik B=0 (kein Kreuzterm) Diagonale Form bewiesen für alle statischen Metriken (Schritt 69) statisch (kein t) Zeitunabhängiger Fall D=r² (Eichung, verlustfrei) Areal-Radius-Eichung A(r)=1 (Ansatz, ungeprüft) Eine freie Funktion C(r) übrig unter A(r)=1 Zweig 1: C=1 = flaches Minkowski (trivial) Zweig 2: Fr≠0 radialer Deformationsfluss Schritte 9–58 Potenzgesetz, Phasenraum, Shooting-Trajektorie 12er-Bedingung nie geprüft Widerrufen Schritte 27–58 (Retro-Befund, Schritt 74): keine gültige Lösung exakt für A=1 (Schritte 72–73) A=1: Zweig ist leer einzige Lösung C=1/(1−kr²) ist Killing-trivial jetzt: A(r) frei Laufende Untersuchung: A(r) allgemein W₂=0-Familie verifiziert nichttrivial (Schritte 75–79)
Der Weg der Vereinfachung — und der weitere Verlauf des rechten Zweigs: ausführlich untersucht (Schritte 9–58), widerrufen (Retro-Befund, Schritt 74), für A=1 abschließend als leer klassifiziert (Schritte 72–73). Offen bleibt, ob sich das bei allgemeinem A(r) ändert.

Was tatsächlich herauskam

Kernbefund
Eine saubere Verzweigung: trivial versus möglich nichttrivial

Nach mehreren Vereinfachungsschritten (keine Kreuzterme, keine Zeitabhängigkeit, eine naheliegende Koordinatenwahl) reduzierte sich die gesamte Fragestellung auf eine einzige unbekannte Funktion. Die Bedingung, dass der Einsteintensor symmetrisch wird, zerfällt dabei sauber in zwei Fälle: Entweder die Metrikfunktion nimmt den Wert 1 an — das ist exakt der flache, leere Raum (Minkowski-Raumzeit), also keine echte Abweichung von einem symmetrischen Universum. Oder eine zweite, algebraische Bedingung ist erfüllt, die eine tatsächlich neue, vom trivialen Fall verschiedene Möglichkeit eröffnet.

Test
Ein naheliegender Versuch scheitert — kontrolliert

Der naheliegendste Kandidat für den nichttrivialen Zweig (eine Zusatzfunktion, die einfach proportional zum Radius wächst) führt auf eine entartete, physikalisch unzulässige Metrik — eine Konfiguration, bei der eine Komponente der Metrik null wird, was geometrisch nicht erlaubt ist. Das ist ein negatives, aber informatives Ergebnis: Es zeigt, dass der nichttriviale Zweig kein Freibrief ist, sondern selektiv, und dass ein Quadrieren-Schritt in der Rechnung (technisch notwendig, um eine störende Wurzel zu beseitigen) genau diese Art von Scheinlösung einführen kann — wir haben sie identifiziert und verworfen, statt sie unkommentiert stehen zu lassen.

Stand: Die ursprüngliche Frage wurde auf eine einzige, explizite gewöhnliche Differentialgleichung dritter Ordnung für eine Funktion reduziert. Ob reguläre (nicht-singuläre) Lösungen dieser Gleichung existieren, wurde im Rahmen dieses Projekts nicht entschieden. Symbolisch (mit Computeralgebra) ist damit ein sauberer Endpunkt erreicht — die Fortsetzung erfordert numerische Methoden (etwa ein Shooting-Verfahren), um das zu prüfen.
Was dieses Projekt nicht ist: Es ist keine Lösung von Heisenbergs offener Frage. Es ist eine präzise Reduktion eines Spezialfalls (kugelsymmetrisch, statisch, eine bestimmte Koordinatenwahl, symmetrische aber nicht notwendig verschwindende Materiequelle) auf ein einziges, explizites mathematisches Objekt — mit dem ehrlichen Befund, dass ein erster Kandidat für ein Gegenbeispiel sich als Artefakt der Rechnung entpuppt hat, nicht als echte Lösung.

Die ODE-Analyse: Potenzgesetz-Lösungen

Bevor ein numerisches Verfahren sinnvoll ist, muss man wissen, wonach man überhaupt sucht. Deshalb wurde zuerst geprüft, ob die gefundene Gleichung eine Skalensymmetrie besitzt (Invarianz unter r → λr) — wenn ja, existieren einfache Potenzgesetz-Lösungen Fr(r)=r^p, die als natürlicher Ausgangspunkt für ein Randwertproblem dienen könnten. Die Gleichung besitzt diese Symmetrie tatsächlich. Einsetzen von Fr(r)=r^p liefert eine algebraische Bedingung, unabhängig von r:

2(4p²−12p+25)(2p−3)(p−1)³p = 0

Reelle Wurzeln: p=0, p=1 (dreifach — das ist die bereits verworfene entartete Lösung), und p=3⁄2. Beide neuen Kandidaten wurden einzeln geprüft, mit einer Selbstkorrektur unterwegs: eine erste Rechnung ohne die Annahme r>0 lieferte ein irreführendes Ergebnis, weil √(r⁶) nicht automatisch zu r³ vereinfacht wurde; mit der korrekten Annahme ergab sich ein anderes, verlässliches Bild.

Vollständiges Negativresultat für diese Ansatzklasse: p=0 liefert nur die triviale Lösung C=1 (Minkowski) — die einzige zweite algebraische Möglichkeit (C=−4) verletzt die Zusatzbedingung und scheidet aus. p=3⁄2 liefert ebenfalls nur C=1 als zulässige Lösung; die zweite algebraische Möglichkeit (C=−1) verletzt die Zusatzbedingung nachweislich. Da die Skalensymmetrie eine vollständige Suche innerhalb der Potenzgesetz-Klasse garantiert, ist damit diese gesamte Ansatzklasse erschöpft — ohne ein einziges nichttriviales Ergebnis.

Konsequenz: Ein numerisches Shooting-Verfahren lässt sich nicht auf ein Potenzgesetz-Verhalten als Randbedingung stützen, weil es keins gibt. Falls überhaupt eine nichttriviale Lösung existiert, muss ihr Verhalten komplizierter sein als ein einfaches Potenzgesetz. Deshalb wurde das Shooting-Verfahren zurückgestellt, bis diese Asymptotik verstanden ist — ein numerischer Versuch ohne begründetes Startintervall wäre geratenes Herumprobieren, kein begründeter nächster Schritt.

Gegenprüfung durch das lokale Modellkomitee: Der konkrete Rechenschritt für p=3⁄2 (Diskriminante, beide Lösungszweige, Prüfung der Zusatzbedingung) wurde unabhängig an qwen3:235b-a22b und deepseek-r1:70b (lokal, Ollama) gegeben, ohne Werkzeugzugriff, nur eigene Rechnung. qwen3:235b-a22b bestätigte alle Zwischenwerte exakt. deepseek-r1:70b unterlief ein eigener Rechenfehler bei einem Koeffizienten, mit dennoch formaler Zustimmung zum Gesamtergebnis — ein dokumentiertes Beispiel dafür, dass auch das lokale Komitee nicht unfehlbar ist, und ein Beleg dafür, warum unabhängige Gegenrechnung mehr wert ist als eine einzelne Bestätigung.

Nachtrag: logarithmische Korrektur, ebenfalls negativ

Da p=1 eine dreifache Wurzel der Potenzgesetz-Gleichung ist, schlug qwen3:235b-a22b — nach eigener Rechnung, ohne Vorgabe — den in der ODE-Theorie üblichen nächsten Schritt für mehrfache charakteristische Wurzeln vor: eine logarithmische Korrektur Fr(r)=r·(1+a·ln r+b·(ln r)²). Mathematisch ein sauberer, korrekt begründeter Vorschlag. Eingesetzt in die tatsächliche (nichtlineare) Gleichung ergibt sich aber: Da 1, ln r, (ln r)², (ln r)³ als Funktionen von r linear unabhängig sind, muss jeder ln(r)-Koeffizient einzeln verschwinden — und das erzwingt a=b=0. Auch dieser Ansatz kollabiert also auf den trivialen Fall. Kein Fehler des Komitees: Der Trick funktioniert bei linearen ODEs mit mehrfachen Wurzeln zuverlässig, hier verhindert die Nichtlinearität der Gleichung genau das. Informativ, weil es zeigt, dass der naheliegendste Nächste-Schritt-Kandidat ebenfalls ausscheidet.

Wichtiger Vorbehalt, inzwischen geprüft: Eichabhängigkeit der bisherigen Struktur

Alle bisherigen Ergebnisse (die saubere Verzweigung, die Potenzgesetz-Faktorisierung mit ganzzahligen Wurzeln, das perfekte-Quadrat-Diskriminanten-Muster) wurden zunächst nur in der Areal-Radius-Eichung (A=1, D=r²) erzielt. Zur Kontrolle wurde dieselbe Herleitung in einer anderen Standardeichung wiederholt: isotrope Koordinaten (A=1, gρρ=gθθ/ρ²=Ψ(ρ)⁴, ebenfalls nur eine freie Funktion). Die Algebra sieht dort tatsächlich anders aus: Die analoge Bedingung liefert zunächst keine saubere quadratische, sondern eine kubische Gleichung (Cardano-Struktur), und die Potenzgesetz-Analyse (Fr(ρ)=ρ^p, Ψ(ρ)=ρ^q) ergab zunächst eine ganze Kurve möglicher (p,q)-Paare statt isolierter Punkte — das sah zwischenzeitlich nach einer echten, in der Areal-Eichung unsichtbaren nichttrivialen Lösung aus.

Die vollständige Parallel-Herleitung (alle vier Bedingungen (LT3G)00,11,12,22=0 gemeinsam, nicht nur die ersten beiden) löst das auf: Die vermeintliche Kurve erfüllt die dritte Bedingung ((LT3G)12=0) an keiner Stelle — diese Bedingung faktorisiert unabhängig zu −4(q+1)q·ρ^(p−2)·sin(θ), null nur für q=0 oder q=−1. Beide überlebenden Fälle wurden explizit auf Flachheit geprüft: q=0 (Ψ=const) ist trivial per Konstruktion; q=−1 (Ψ=1/ρ) erwies sich bei direkter Berechnung des Riemann-Tensors als identisch verschwindend — flacher Raum, nur in invertierten isotropen Koordinaten geschrieben (die bekannte Kelvin-Inversion ρ→1/ρ, unter der flacher Raum mit genau diesem Ψ⁴-Konformfaktor auf sich selbst abgebildet wird).

Ergebnis der Eichkontrolle: Innerhalb der Potenzgesetz-Ansatzklasse ist das vollständige Negativresultat jetzt in zwei unabhängigen Eichungen bestätigt — mit völlig unterschiedlicher Algebra unterwegs (quadratisch vs. kubisch, ein zusätzlicher, versteckter Inversions-Zweig, der in der Areal-Eichung gar nicht auftaucht), aber demselben Endergebnis: keine nichttriviale Lösung. Das ist ein deutlich stärkerer Beleg für Eichunabhängigkeit als eine einzelne Rechnung — nicht bewiesen für jede denkbare Eichung, aber nicht mehr nur eine möglicherweise zufällig günstige Parametrisierung.
Korrekturvermerk (Schritte 72–73, Juli 2026): Die folgenden Abschnitte bis einschließlich „Schritt 48–55“ (samt der Nachträge zu Schritt 64/66) beruhen auf der Shooting-Trajektorie, die konstruktionsbedingt nur aus der 00er- und der 11er-Bedingung gebaut wurde; die 12er-Bedingung wurde entlang der Trajektorie nie geprüft. Die spätere vollständige Klassifikation (Schritte 69–74) zeigt: Diese Trajektorie verletzt die 12er-Bedingung und ist keine LT3G=0-Konfiguration. Die ODE- und Krümmungsanalysen unten sind als Mathematik der reduzierten Gleichung korrekt, beantworten aber nicht die Projektfrage. Die Abschnitte bleiben dem Fehler-Transparenz-Prinzip dieses Projekts folgend unverändert stehen. Komitee-geprüft (qwen3:235b-a22b mit eigener Gegenrechnung der Asymptotik, deepseek-r1:70b).

Phasenraumanalyse: zwei Sättel, keine Spiralen

Zurück zur Areal-Eichung. Nächster Schritt: die ODE dritter Ordnung via x=ln(r) und einer logarithmischen Hilfsvariable u=d(ln Fr)/dx in ein autonomes System umschreiben. Potenzgesetz-Lösungen entsprechen darin Fixpunkten u=p=const — ein Kontrolltest (u=p, alle Ableitungen null einsetzen) reproduziert exakt die bekannte Quintik 2(4p²−12p+25)(2p−3)(p−1)³p, was die Transformation bestätigt.

Linearisierung um jeden Fixpunkt liefert eine charakteristische Gleichung mit den Eigenwerten (Stabilitätsexponenten):

Befund: Beide nichtentarteten Fixpunkte sind Sättel mit reellen Eigenwerten, keine Spiralen. Das bedeutet: es gibt keinen lokalen Hinweis auf versteckte, oszillierende Nicht-Potenzgesetz-Lösungen in der Nähe dieser selbstähnlichen Punkte. Eine echte, nichttriviale Lösung — falls sie existiert — wäre keine kleine Störung um p=0 oder p=3/2, sondern müsste eine globale Bahn sein (etwa eine heterokline Verbindung zwischen den Sätteln); das kann lokale Linearisierung prinzipiell nicht zeigen. Sie liefert aber die korrekten lokalen Eigenvektoren, mit denen ein Shooting-Verfahren bei r→0 oder r→∞ jetzt tatsächlich begründet ansetzen könnte.

Numerisches Shooting: eine echte, verifizierte nichttriviale Bahn — und eine wichtige Korrektur

Hinweis: überholt — siehe Korrekturvermerk oben und Schritte 69–74.

Ausgehend von der instabilen Richtung des Sattels bei p=0 (Julia, DifferentialEquations.jl) wurde die Trajektorie numerisch integriert. Eine technische Schwierigkeit musste dabei sauber gelöst werden: Die autonome Gleichung ist quadratisch in der zweiten Ableitung u_xx — geerbt aus dem Quadrieren, das schon in Schritt 9 nötig war —, sodass bei jedem Integrationsschritt zwei algebraisch mögliche Werte existieren, von denen einer eine Scheinlösung sein kann. Die Zweigwahl erfolgte per Kontinuität (stets der zum vorherigen Wert nächstgelegene).

Für eine Störung in Richtung des instabilen Eigenvektors ergab sich eine stabile, bis x=8 laufende Trajektorie, die sich u→1 annähert (dem entarteten Fixpunkt, nicht p=3/2). Bevor das als Ergebnis zählen konnte, wurde es unabhängig gegen die ungequadrierte Bedingung aus Schritt 8 geprüft — diese lässt sich elegant rein in (u,u_x) ausdrücken als dC/dx=2(C−1)(1−u). Bei grober Auflösung zeigte der Vergleich kleine Restdifferenzen; bei feinerer Auflösung (Δx=0,02 statt 0,5) schrumpften diese auf 10&supmin;⁵–10&supmin;⁶ — numerisches Rauschen, keine echte Verletzung. Die Trajektorie ist damit als echte Lösung bestätigt, kein Artefakt des Quadrierens.

Die notwendige Korrektur: Bei dieser Verifikation fiel eine grundlegendere Ungenauigkeit auf, die die gesamte bisherige Darstellung betrifft. Die Bezeichnung „Vakuumfall“ oben ist unpräzise: Wir haben nie Gμν=0 gefordert, sondern nur, dass Gμν symmetrisch ist. In der Areal-Radius-Eichung gilt exakt Grr=−(C−1)/r² — das ist nur für C=1 (flach) null. Jede nichttriviale Lösung, einschließlich der eben verifizierten Trajektorie, hat also zwangsläufig Gμν≠0 — es ist Materie im Spiel, kein Vakuum.

Das wirft die Frage auf, ob der echte Vakuumfall (G≡0 für alle Komponenten) eigenständig untersucht werden sollte. Die Antwort ist aufschlussreich: Mit A(r) und C(r) beide frei (statt A=1 fest) liefert Gtt=Grr=0 exakt die Schwarzschild-Metrik (Birkhoffs Theorem, hier direkt bestätigt: A=1−2M/r, C=1/(1−2M/r) ergibt G≡0 identisch). Aber der Vakuumfall von Heisenbergs Frage ist dabei trivial — aus einem definitorischen, nicht physikalischen Grund: Der Nulltensor ist gegenüber jeder Transformation symmetrisch, unabhängig davon, ob diese Transformation überhaupt eine Symmetrie der Metrik ist. Schwarzschild selbst ist bereits ein (uninteressantes) Beispiel: gμν ist nicht translationssymmetrisch, aber Gμν=0 ist trivial „symmetrisch“.

Auflösung: Nicht die Rechnung war falsch gerahmt, sondern nur das Wort „Vakuum“. Die einzige Version der Frage mit echtem Inhalt verlangt eine symmetrische, aber nicht notwendig verschwindende Materiequelle — genau das, was dieses Projekt von Anfang an tatsächlich untersucht hat. A=1 war dabei eine legitime Eichwahl für diesen Fall (eine Zeitkoordinatenwahl bei vorgegebener Materie), kein Fehler. Die heutige, verifizierte Trajektorie ist damit der erste echte, nichttriviale Kandidat des gesamten Projekts — mit Materie, nicht im Vakuum. Alle früheren Formulierungen auf dieser Seite wurden entsprechend korrigiert.

Geometrische Charakterisierung: negative Energiedichte, Zustandsgleichung w=−1/3

Hinweis: überholt — siehe Korrekturvermerk oben und Schritte 69–74.

Aus Gμν=8πTμν lassen sich Energiedichte ρ und Drücke direkt ablesen: ρ·8π=(C²+rC′−C)/(r²C²), Radialdruck pr·8π=−(C−1)/(r²C), Tangentialdruck pt·8π=−C′/(2rC²). Ausgewertet entlang der verifizierten Trajektorie (r von 1,65 bis über 1000): ρ·8π ist durchgehend negativ (≈−0,003), ebenso ρ+pr und ρ+pt (≈−0,002). Die schwache Energiebedingung ist verletzt — die Materie, die diese Geometrie trägt, ist exotisch, keine gewöhnliche Substanz.

Auffällig: ρ·8π bleibt über einen r-Bereich von 2,8 bis 1141 (Faktor >400) nahezu konstant, mit dem Verhältnis w=pr/ρ≈−1/3 — die Zustandsgleichung kosmischer Strings/Defektnetzwerke, zugleich die kritische Schwelle zwischen Bremsung und Beschleunigung in der Friedmann-Gleichung. Analytisch bestätigt: im Grenzfall C→0 reduziert sich w auf C/(rC′−C); w=−1/3 ist dazu äquivalent zu C′=−2C/r, also C∝1/r² — und genau dieses Verhalten nähert sich die numerische Trajektorie an (die Konstante K=C·r² wächst von 7,9 über 962 auf über 1000, noch nicht konvergiert, aber sichtbar einem Grenzwert zulaufend). Die Konsistenz wurde zusätzlich gegen die ursprüngliche, ungequadrierte Bedingung geprüft: C∝1/r² impliziert dort u≈1−C, also u→1 im selben Grenzfall — Fixpunktanalyse, Energiebedingung und Integrabilität passen strukturell zusammen, nicht nur numerisch.

Konvergenzbeweis: Der Fixpunkt u=1 ist nicht nur linear (Schritt 26), sondern auch quadratisch entartet — erst die kubische Ordnung der Störungsentwicklung (u=1−ε) liefert eine nichttriviale Gleichung. Der Ansatz ε(x)=A·e−γx führt zu einer sauber faktorisierenden Bedingung: −(γ²+2γ+17)(γ+1)(γ−2)=0. Reelle Wurzeln: γ=2 und γ=−1. γ=2 ist exakt die Rate, die C∝1/r² erfordert — jetzt algebraisch bewiesen, nicht nur numerisch beobachtet. Die Amplitude A bleibt frei (der erwartete Integrationskonstanten-Freiheitsgrad, entspricht dem numerisch beobachteten K).

Literaturvergleich: w=−1/3 als bekannte, aber strukturell andersartige Signatur

Bronnikov & Zaslavskii (Phys. Rev. D 78, 021501(R), 2008, arXiv:0801.0889; erweitert in arXiv:0904.4904; populär zusammengefasst als „black holes can have curly hair“, arXiv:0905.4425) zeigen: Materie im Gleichgewicht mit einem Killing-Horizont ist nur für diskrete w=−n/(n+2k) möglich (n=Horizont-Ordnung, k=Dichteabfall-Rate). Der generische Fall (n=1, k=1) ergibt exakt w=−1/3 — von ihnen explizit als „Gas ungeordneter kosmischer Strings“ identifiziert.

Der Unterschied zu unserem Fall ist wichtig: Ihr Ergebnis gilt an einem Killing-Horizont (dort verschwindet die Zeitkomponente der Metrik). Unsere Trajektorie hat A=1 fest, nie null — kein Horizont in diesem Sinn. Bei uns verschwindet stattdessen die Radialkomponente C(r) bei r→∞, ein strukturell anderer entarteter Punkt (räumliche statt zeitliche Entartung, Unendlichkeit statt Horizont). Die Übereinstimmung liegt nicht im geometrischen Ort, sondern in der Methode: Beide Male erzwingt eine lokale Entwicklung der Einsteingleichungen um einen entarteten Punkt mit der einfachsten Ordnung denselben Wert w=−1/3 — eine echte strukturelle Parallele, keine Identität der Aufbauten.

Die vollständige globale Lösung: bestätigt bei r→∞, ungeklärt bei r→0

Hinweis: überholt — siehe Korrekturvermerk oben und Schritte 69–74.

Die Trajektorie wurde über einen deutlich größeren Bereich integriert (x=0 bis x=20 statt x=8). Vorwärts (r→∞): K=C·r² konvergiert sauber — 0,999 bei r=1, dann 51,8, 748,0, 992,5, 998,5 bei r bis knapp 3000. Das bestätigt numerisch exakt, was zuvor analytisch bewiesen wurde (γ=2): K läuft klar auf einen Grenzwert nahe 1000 zu. Die Integration bricht bei x=8 nur aus numerischen Gründen ab (die Dynamik wird am entarteten Fixpunkt extrem langsam) — kein physikalisches Problem.

Rückwärts (r→0), erstmals untersucht: Hier ist der Befund mehrdeutig, und das soll er auch bleiben, bis er geklärt ist. C(r) bleibt bis r≈0,0067 fast exakt bei 1 (0,999998 bis 1,0) — die Metrik selbst scheint dort praktisch flach zu bleiben. Die Hilfsgröße u (die logarithmische Ableitung von Fr) divergiert dagegen dramatisch, von −0,003 bei r=0,05 auf −10,27 bei r=0,0067, dann bricht die Integration instabil ab.

Was das bedeuten könnte, ohne dass ich es entscheiden kann: Da Grr=−(C−1)/r² und C praktisch bei 1 bleibt, könnten Energiedichte und Druck dort tatsächlich gegen null gehen — die Geometrie selbst bliebe dann nahe flach, während nur die Hilfskonstruktion (der Translationsgenerator T3, ausgedrückt über u) an dieser Stelle zusammenbricht: kein physikalischer Krümmungsdefekt, sondern ein Versagen der Koordinatenbeschreibung. Es könnte aber auch sein, dass die eigentliche Physik hier noch nicht sichtbar ist, weil bisher nur C und u geprüft wurden, nicht die tatsächlichen Krümmungsinvarianten. Das bleibt offen.

Krümmungsinvarianten bei r→0: Korrektur eines eigenen Fehlers

Bei der Vorbereitung der feinen r→0-Integration wurde der zuvor verwendete Kretschmann-Skalar gegengeprüft — und als fehlerhaft entlarvt. Eine echte Krümmungsinvariante darf nicht von θ abhängen; die Formel aus der ersten Rechnung enthielt aber einen nicht herausfallenden cos²θ-Term, ein Indexkontraktionsfehler bei der automatisierten Index-Hebung in SageManifolds. Die saubere Neuberechnung über vollständige, explizite Indexkontraktion ergibt die tatsächlich θ-unabhängige Formel:

K = 2(2C⁴+r²C′²−4C³+2C²)/(r⁴C⁴)

Mit dieser korrigierten Formel an denselben r→0-Punkten neu ausgewertet:

r≈0,252: Ricci=−0,00647, Kretschmann=1,41·10&supmin;⁵
r≈0,0925: Ricci=−0,00638, Kretschmann=1,36·10&supmin;⁵
r≈0,0341: Ricci=−0,00546, Kretschmann=0,99·10&supmin;⁵

Korrigierter Befund, das Gegenteil des vorherigen: Mit der richtigen Formel bleibt auch der Kretschmann-Skalar klein und nahezu konstant (leicht fallend statt wachsend). Die zuvor berichtete „echte Krümmungszunahme“ existierte nicht — sie war ein Rechenfehler, keine Physik. Damit ist die ursprünglich vorsichtige Hypothese (nur ein Versagen der T3/u-Hilfskonstruktion, keine echte geometrische Auffälligkeit in diesem Bereich) wieder die naheliegendere Lesart, nicht widerlegt.
Warum das hier offen dokumentiert wird: Ein publizierter Zwischenbefund erwies sich als falsch, nicht nur als vorläufig. Das wird nicht stillschweigend ersetzt, sondern ausdrücklich als Korrektur festgehalten — derselbe Standard, der in diesem Projekt durchgängig gelten soll.

Feine, verifizierte Integration: das eigentliche Bild

Die grobe Dreipunkt-Nachprüfung oben war selbst irreführend — zu wenige, zu weit auseinanderliegende Punkte mit grobem Differenzenquotienten für C′. Die feine (effektiv Δx=0,02), gegen die ungequadrierte Bedingung verifizierte Integration (Restdifferenzen 5·10&supmin;⁵ bis 2,6·10&supmin;¹⁰ — echte Lösung bestätigt) zeigt: Der Ricci-Skalar bleibt über den gesamten untersuchten Bereich (r=0,368 bis r=0,0067) praktisch exakt konstant bei −0,0060095. Der Kretschmann-Skalar bleibt bis r≈0,05 ebenfalls praktisch exakt konstant (1,2038·10&supmin;⁵), zeigt dann aber einen echten, wenn auch bescheidenen Anstieg auf 1,659·10&supmin;⁵ bei r=0,0067 — etwa 38 % Zuwachs über die letzte Dekade in r.

Ehrliches Fazit nach drei verschiedenen Zwischenständen zu dieser Frage: Die Krümmung bleibt über den größten Teil des untersuchten Bereichs bemerkenswert flach, zeigt aber am unteren Rand des numerisch erreichbaren Bereichs einen echten, moderaten Anstieg — weder das dramatische Wachstum der ersten (fehlerhaften) Rechnung noch die leichte Abnahme der zweiten (zu groben) Nachprüfung. Ob das bei noch kleinerem r zu einer echten Singularität anwächst oder sich stabilisiert, lässt sich mit den vorliegenden Daten nicht entscheiden — die Integration wird bei r≈0,0067 instabil (Divergenz der u/T3-Hilfskonstruktion), bevor sich das klärt.

Schritt 40–47: Kretschmann-Vorhersage bestätigt, dann grundlegende Korrektur — r* ist eine Koordinatensingularität

Der oben offen gelassene Punkt (milder Kretschmann-Anstieg nahe r≈0,0067) wurde in einer Folgesitzung weiterverfolgt. Der Weg dorthin ging über zwei Zwischenkorrekturen, die hier vollständig dokumentiert werden, nicht nur das Endergebnis.

Erste Bestätigung: Eine Fixpunkt-Linearisierung um C=1, u=0 (derselbe Sattel wie in der Phasenraumanalyse oben) sagt K→12A² für r→0 voraus, mit A aus dem Verhältnis ε=C−1 zu r² gefittet. Numerisch bestätigt: A₀≈−9,955·10⁻⁴, 12·A₀²≈1,189·10⁻⁵ — trifft den beobachteten Plateauwert (1,2038·10⁻⁵) auf besser als 1 % genau.

Numerische Nachprüfung mit höherer Präzision: Um zu klären, ob der Anstieg bei noch kleinerem r ein Rechenpräzisions-Artefakt ist, wurde dieselbe Integration mit 256-Bit-BigFloat-Arithmetik wiederholt (statt gewöhnlicher Fließkomma-Genauigkeit), mit zwei unterschiedlichen Integrationsverfahren (feste Schrittweite und ein etablierter adaptiver Löser). Beide brechen praktisch am selben Ort zusammen wie die ursprüngliche Rechnung: bei x≈−5,045 (r*≈0,00644). Drei methodisch unabhängige numerische Ansätze, die exakt an derselben Stelle scheitern, legten zunächst nahe, dass dort tatsächlich eine echte Singularität der Lösung vorliegt — kein Präzisionsproblem.

Analytische Klärung (Frobenius-Ansatz): Um das zu entscheiden, wurde ein Polstellenansatz u(x)~A·(x−x*)⁻ᵖ direkt in die algebraische Bestimmungsgleichung für die zweite Ableitung eingesetzt. Die führende Ordnungsbalance ist konsistent bei p=1; die resultierende Koeffizientengleichung (computeralgebraisch gelöst, nicht per Hand — eine erste Handrechnung enthielt einen Fehler, der bei der Gegenprüfung auffiel) liefert 2A³(A−1)(A−½)²(A+½)=0, mit A=−½ als einzig physikalisch relevanter Wurzel. Das ergibt u(x)~−1/(2(x−x*)) — ein exakter, rationaler Koeffizient, keine numerische Schätzung, und mit den numerischen u-Werten auf unter 1 % Abweichung bestätigt.

Das eigentliche Ergebnis: Dieses u(x)-Verhalten wurde durch die Formel für C(r) und weiter durch Ricci-Skalar R und Kretschmann-Skalar K propagiert (Reihenentwicklung um die Singularität). Ergebnis: C(r) divergiert dort tatsächlich (C~−1/(2ε²)), aber sowohl R als auch K bleiben endlich — beide sind reguläre Taylor-Reihen ohne negative Potenzen, mit K→9/(2r*⁴) und R→3/r*² im Grenzwert.

Die korrekte Einordnung, nach zwei vorherigen Zwischenständen ("Artefakt vermutet", dann "vermutlich echte Krümmungssingularität"): r*≈0,0064 ist eine Koordinatensingularität der Arealradius-Eichung — die Metrikkomponente C(r) divergiert, während die tatsächliche Raumzeitkrümmung (Ricci, Kretschmann) endlich bleibt. Die Analogie ist die divergierende grr-Komponente von Schwarzschild in Schwarzschild-Koordinaten am Horizont, wo Krümmungsinvarianten ebenfalls endlich sind. Die drei numerischen Integrationsverfahren brachen zu Recht zusammen — nicht wegen echter Krümmungsdivergenz, sondern weil die Koordinate selbst dort unbrauchbar wird.

Offen bleibt: ob sich die Lösung in besseren Koordinaten (analog zu Eddington-Finkelstein bei Schwarzschild) glatt über r* hinaus fortsetzen lässt, und wie sie sich für r<r* verhält.

Schritt 48–55: Verifikation, Fortsetzung über r* hinaus, und eine echte Singularität bei r=0

Hinweis: überholt — siehe Korrekturvermerk oben und Schritte 69–74.

Zwei Lücken aus den Schritten 45b/46 wurden zuerst geschlossen. Erstens: Die Wahl A=−½ für die Frobenius-Wurzel war nie tatsächlich belegt worden — das Sage-Skript in Schritt 45b brach mit einem internen Fehler exakt an der Stelle ab, an der die Auswahl hätte begründet werden müssen, und Schritt 46 übernahm A=−½ ohne Beleg. Beide Wurzeln (A=−½ und A=+½) wurden jetzt vollständig durchgerechnet: A=+½ liefert nur einen einfachen Pol mit irrationalen (√41-behafteten) Koeffizienten, A=−½ eine rein rationale Reihe. Entscheidend ist das Vorzeichen von u(ε) bei Annäherung an x* — die tatsächliche BigFloat-Trajektorie (Schritt 41–43) fällt auf u=−10,27, was nur A=−½ reproduziert. A=+½ ist damit widerlegt, nicht nur verworfen.

Wichtiger Nebenbefund: Beide Wurzeln liefern in erster Näherung endliches K — "Krümmung bleibt endlich" unterscheidet die Zweige NICHT. Nur das u-Vorzeichen und die Rationalität der Reihenkoeffizienten tun das.

Regularisierung durch v=1/u, und ein dokumentierter eigener Fehler

Da u selbst divergiert (nicht nur C), wurde die Bewegungsgleichung direkt (nicht nur ihre lokale Reihe) auf v=1/u umgeschrieben. Bei v=0 (=x*) fällt die zweite Ableitung vxx aus der Gleichung heraus; übrig bleibt eine reine Bedingung an die Steigung: −2(vx−2)²(vx−1)(vx+2)=0, mit den Wurzeln vx∈{2 (doppelt), 1, −2}.

Eine erste Handrechnung wählte hier den falschen Zweig (vx=+2, aus einem Vorzeichenfehler bei der Ableitung dv/dx=1/A). Die daraus berechneten Reihenkoeffizienten widersprachen der numerischen Trajektorie im Vorzeichen komplett (v ist durchgehend negativ, die Vorhersagen waren positiv) — der Fehler wurde dadurch aufgedeckt und korrigiert, nicht stillschweigend behoben. Der korrekte Zweig (vx=−2) liefert zwei algebraische Kandidaten für die nächsten Reihenkoeffizienten; gegen die BigFloat-Trajektorie (mit dem präzisen, nicht gerundeten x* aus der Solver-Warnung) getestet, bestätigt sich genau einer auf besser als 0,1 % Übereinstimmung — und zwar mit irrationalen (√17-behafteten) Koeffizienten. Die "rationale Koeffizienten"-Heuristik aus dem A=±½-Fall war hier also kein verlässliches Kriterium, nur ein unterstützendes Indiz.

Ergebnis: v(ε) ist eine gewöhnliche, konvergente Potenzreihe durch ε=0 hindurch, numerisch bestätigt. Die Lösung setzt sich analytisch über r* hinaus fort — nicht nur C(r) war regularisierbar, die gesamte lokale Lösung ist es.

Numerische Fortsetzung nach r<r*: eine echte Singularität bei r=0

Mit der exakt aus der Bewegungsgleichung hergeleiteten ODE für v (analog zu den früheren Acoef/Bcoef/Ccoef-Koeffizienten für u) wurde die Integration über den Punkt hinaus fortgesetzt, ausgehend von der oben bestätigten Taylor-Reihe als Anfangsbedingung. Über eine Dekade in r verläuft alles glatt und monoton — C(r) fällt von ≈0,995 auf ≈0,558 — aber der Kretschmann-Skalar wächst dabei bereits deutlich, von ≈1,8·10⁷ auf ≈7,5·10¹³.

Eine deutlich längere Integration (r schrumpft dabei um 26 Größenordnungen, bis 8,9·10&supmin;¹⁴, durchgehend stabil) zeigt: Kretschmann folgt einem sauberen Potenzgesetz K∼r⁻⁶ — der Exponent konvergiert auf −6,000000000168 im am weitesten von r* entfernten Bereich. Das ist keine Tendenz, das ist eine harte Konvergenz auf einen ganzzahligen Exponenten.

Analytische Erklärung, nicht nur Fit: C(r) selbst folgt auf demselben Ast C∼a·r (Exponent im am weitesten von r* entfernten Bereich: 0,9999999), während u(r) gegen exakt 1 konvergiert. Mit C∼a·r, C′∼a (konstant) dominieren im Kretschmann-Zähler bei kleinem r die Terme r²C′² und 2C² (beide Ordnung r²) über die höherordnungen C³- und C⁴-Terme — das ergibt Zähler∼6a²r², Nenner r⁴C⁴∼a⁴r⁸, also K∼(6/a²)·r⁻⁶. Algebraisch hergeleitet, nicht nur numerisch beobachtet — und deckungsgleich mit dem Fit.

Physikalisch: grr=C(r)→0 bei r→0 ist eine echte, nicht behebbare Entartung der Metrik, begleitet von echt divergierender Krümmung — eine genuine Krümmungssingularität, strukturell das Gegenteil von r* (dort divergierte C bei endlichem u; hier verschwindet C bei u→1). Bemerkenswert: u→1 ist exakt derselbe entartete Fixpunkt (dreifache Wurzel p=1), der seit den Schritten 25/26 das ganze Projekt durchzieht.

Was das bedeutet: Der r<r*-Ast führt nicht in ein ruhiges Gebiet, sondern von einer harmlosen Koordinatensingularität (bei r*) direkt auf eine echte Krümmungssingularität (bei r=0) zu. Offen bleibt die physikalische Klassifikation dieser Singularität (Vergleich mit bekannten Klassen wie Kasner-artigen Anisotropie-Singularitäten — C∼r könnte darauf hindeuten, ist aber ungeprüft).

Schritt 56: Physikalische Klassifikation über den Eigenabstand

Der Eigenabstand ρ=∫√C(r)dr wurde entlang des r<r*-Astes numerisch berechnet. Zentrales Ergebnis: ρ konvergiert gegen einen endlichen Grenzwert (ρend≈0,006009), obwohl r selbst gegen 0 geht — die Singularität liegt in endlichem Eigenabstand, ist also eine echte, erreichbare Grenze der Raumzeit (wie Schwarzschild r=0), keine bloß asymptotische Eigenschaft. Der Skalierungsexponent r∼ρ^q bestätigt dabei die aus C∼a·r (Schritt 55) erwartete Beziehung q=2/3.

Bild nahe der Singularität: ds²≈−dt²+dρ²+const·ρ^(4/3)dΩ² — beide Winkelrichtungen schrumpfen gemeinsam wie ρ^(2/3). Kein Kasner im strengen kosmologischen Sinn (dort drei unabhängige räumliche Exponenten mit Summenbedingungen aus der Vakuumgleichung), sondern eine anisotrope statische Punktsingularität mit nichttrivialem Exponenten 2/3.

Bronnikov-Zaslavskii-Vergleich: ein struktureller, nicht nur negativer Befund

Die einschlägigen Arbeiten (Bronnikov & Zaslavskii, PRD 78, 021501, arXiv:0801.0889; Bronnikov, Elizalde, Odintsov & Zaslavskii, PRD 78, 064049, arXiv:0805.1095) wurden vollständig gelesen, nicht nur die Abstracts. Ihre Klassifikation ("usual"/"naked"/"truly naked" Horizonte) setzt voraus, dass die Zeitkomponente der Metrik an einem Kandidatenhorizont gegen null geht — bei uns ist A=1 exakt fest, nie null. Unser r=0 ist also gar kein Horizont-Kandidat in ihrem Sinn.

Der eigentliche Grund, warum diese Literatur hier nicht greift: Die Autoren grenzen ihre Analyse selbst ausdrücklich von gewöhnlichen Krümmungssingularitäten (divergierende Invarianten) ab — genau das ist unser Fall (K∼r⁻⁶ divergiert). Ihre "truly naked horizons" sind gerade dadurch interessant, dass dort die Krümmungsinvarianten endlich bleiben. Kein negativer Befund aus unzureichender Suche, sondern ein struktureller Grund: diese Arbeiten klassifizieren bewusst ein anderes, subtileres Phänomen. Der früher gefundene w=−1/3-Treffer betrifft weiterhin nur den r→∞-Ast, nicht diesen r=0-Punkt.

Schritt 57–58: Ricci-Tensor-Komponenten und das Tipler-Kriterium

Die vollen Ricci-Tensor-Komponenten (nicht nur der Skalar) wurden für die Metrik hergeleitet: Rtt=0 (exakt, algebraische Folge von A=1, unabhängig von C(r)), Rrr=C′/(rC). Damit lässt sich die Tipler-Frage entscheiden: Ist die Singularität "stark" (zerstört jeden ausgedehnten Körper) oder "schwach"?

Für eine radiale zeitartige Geodäte (Erhaltungsgröße E, (ur)²=(E²−1)/C) wurde die Größe τ²·Rabuaub (τ = verbleibende Eigenzeit bis zur Singularität) analytisch und numerisch bestimmt. Sie konvergiert auf einen endlichen, von Null verschiedenen Wert (≈0,4444, konstant über rund 15 Größenordnungen in r).

Nachtrag (Schritt 64, lokales Modellkomitee + CAS-Gegenprüfung): Der Grenzwert ist keine numerische Eigenart des gewählten Referenzfalls E²−1=1, sondern eine universelle Konstante 4/9 — exakt unabhängig von E und dem Koeffizienten a in C∼a·r nahe r=0 (Sage-Verifikation: die symbolische Ableitung des Ergebnisses nach E und nach a ergibt beide Male 0). Ein unabhängiges lokales Modell (qwen3:235b-a22b) widersprach zunächst einer im internen Rechenkommentar genannten Formel (E²−1)/a; ein zweites Modell (DeepSeek-R1-671B) stimmte dieser — wie sich zeigte, falschen — Formel zu. Die CAS-Nachrechnung bestätigte den Einwand des ersten Modells. Nachtrag (Schritt 66): Qwen3 hatte zusätzlich eingewandt, dass das Clarke-Królak-Kriterium strenggenommen nur für Geodäten gilt, die tatsächlich an der Singularität enden — nicht-radiale Geodäten (Drehimpuls L≠0) erreichen r=0 aber nie (Umkehrpunkt bei r_turn=L/√(E²−1)>0, CAS-bestätigt). Auf gezielte Nachfrage bestätigten beide lokalen Modelle unabhängig voneinander (qwen3:235b-a22b und DeepSeek-R1-671B, kein Dissens): Da nicht-radiale Geodäten nicht an der Singularität enden, gehören sie nicht zur relevanten Menge des Kriteriums — die Prüfung der radialen Geodäten (die das Kriterium mit 4/9>0 erfüllen) ist somit vollständig und hinreichend.
Ergebnis: Das erfüllt das hinreichende Kriterium von Clarke & Królak (1985) für eine starke Krümmungssingularität im Sinne von Tipler (1977). Ein ausgedehnter, radial einfallender Körper wird nicht nur unendlich gezeitenverzerrt — sein durch Jacobi-Felder definiertes Volumen geht tatsächlich gegen null.

Damit sind die physikalische Klassifikation der r=0-Singularität und der Vergleich mit der Bronnikov-Zaslavskii-Literatur abgeschlossen.

Schritt 59–62: Erste Lockerung — Kreuzterm B(r)≠0 (statisch)

Begonnen mit der kleineren der beiden verbleibenden Vereinfachungen: der Kreuzterm B(r) in ds²=−A dt²+2B dt dr+C dr²+r²d&Omega² wird wieder zugelassen (A=1, D=r² bleiben). Kontrolltest bestanden: bei B=0 werden die früheren Ergebnisse reproduziert. Zentraler neuer Befund: Es taucht eine vollständig neue Bedingung (LT3G)02=0 auf, die bei B=0 identisch verschwindet — bei B≠0 aber echt ist.

Diese Bedingung liefert, nach Herausdividieren des trivialen Faktors B, eine saubere gewöhnliche Differentialgleichung erster Ordnung für B(r) in Abhängigkeit von C(r). B=0 ist nur die triviale Lösung dieser Gleichung — für jedes vorgegebene C(r) existiert ein echtes, nichttriviales B(r). Der Kreuzterm-Zweig ist also real, nicht nur formal möglich.

Struktureller Fund: Nach Einsetzen dieser B'-Beziehung in die übrigen Gleichungen zeigt sich (LT3G)01 = −B·(LT3G)00 exakt — diese beiden sind nicht unabhängig. Das System reduziert sich von sechs auf effektiv vier Bedingungen.

Die (LT3G)12-Gleichung entkoppelt vollständig: Sie faktorisiert (von Sage bestätigt) zu Q·[2(B²+C)Fr(r) + Fth(r) − r·Fth'(r)] = 0, mit Q:=B²+C−1. Im generischen Fall Q≠0 ist das eine lineare ODE erster Ordnung, die sich direkt durch Quadratur lösen lässt: Fth(r) folgt aus einem expliziten Integral über B, C und Fr — kein eigenständiges Problem mehr. Der Sonderfall Q=0 entspricht bemerkenswerterweise exakt der Bedingung −det(g)=r⁴sin²θ.

Ehrlicher Stand: Das eigentlich harte System — drei gekoppelte Funktionen B(r), C(r), Fr(r) in den verbleibenden drei Gleichungen — ist damit noch nicht gelöst. Die Entkopplung von Fth ist der Fortschritt dieser Sitzung, das Kernsystem bleibt offen. Das ist kein kurzer nächster Schritt mehr, sondern ein eigenständiges Vorhaben von ähnlichem Umfang wie das bisherige Projekt.

Schritt 67: Der Sonderzweig Q=0 ist reine Eichung — abgeschlossen

Der in Schritt 62 identifizierte Sonderfall Q:=B²+C−1=0 (äquivalent zur dort gefundenen Determinantenbedingung) wurde als erster Zugriff auf das Kernsystem untersucht: B=±√(1−C) direkt in die Metrik eingesetzt (A=1, D=r²). Ergebnis (SageManifolds, symbolisch): Nicht nur alle (LT3G)-Bedingungen, sondern der gesamte Einstein-Tensor und sogar der volle Riemann-Tensor verschwinden identisch — für beliebiges C(r) mit 0<C≤1. Die Metrik ist flach.

Analytische Erklärung (CAS-verifiziert): Die Koordinatentransformation T=t−f(r), R=r mit f'(r)=B(r) überführt Minkowski exakt in diese Metrik (g01=f'=B, g11=1−f'²=C); beide Vorzeichen von B sind durch das Vorzeichen von f' abgedeckt, C=1 (B=0, Standard-Minkowski) gehört zum Zweig. Der Q=0-Zweig ist die reine Eichbahn von Minkowski unter radiusabhängigen Zeitverschiebungen — dieselbe Konstruktion, die bei Schwarzschild die Gullstrand-Painlevé-Koordinaten liefert, hier auf den flachen Raum angewandt. Physikalisch ist der Zweig leer: LT3g≠0, aber G=0 ist der triviale Vakuumfall. Die Vermutung aus Schritt 62 („möglicherweise eine natürliche Eichbedingung“) ist damit im stärksten Sinn bestätigt: Q=0 ist eine Eichung. Die B≠0-Untersuchung muss im generischen Fall Q≠0 fortgesetzt werden; Q fungiert dabei als Diskriminator flach/nichtflach, nicht als Krümmungsmaß — die Krümmung hängt von Ableitungen von Q ab (Präzisierung aus dem lokalen Komitee-Review).
Lokales Komitee (Schritt 67): Die Interpretation wurde vor Publikation von zwei lokalen Modellen unabhängig geprüft (qwen3:235b-a22b, deepseek-r1:70b, beide via Ollama auf dem Mac Studio): einstimmige Zustimmung, kein Dissens. qwen3 deckte dabei eine Vorzeichen-Inkonsistenz in der ursprünglichen Formulierung der Behauptung auf (f'=−√(1−C) gepaart mit B=+√(1−C); korrekt: f'=B) — die Substanz des Befunds war davon nicht betroffen, die Korrektur ist dennoch dokumentiert. Ein erster Review-Lauf scheiterte technisch (Thinking-Budget verbraucht, done_reason=length, leere Antwort) und wurde mit größerem Token-Budget wiederholt.

Schritte 68–68d: Das generische Kernsystem kollabiert — und eine neue Lösungsfamilie

Mit S:=B²+C und Q:=S−1 wurde das verbliebene Kernsystem strukturell reduziert (jeder Einzelschritt CAS-verifiziert): Nach Einsetzen der B'-Beziehung gilt exakt (LT3G)01=−B·(LT3G)00; die geeignet normierte Differenz der 00- und 22-Gleichung ist rein algebraisch und faktorisiert zu −2B²SQ·(Fr+2Fth). Für Q≠0 folgt Fth=−Fr/2, aus der 12er-Quadratur dann r·Fr'=(1−4S)·Fr — und damit ist die 00-Gleichung identisch erfüllt (Bianchi-konsistent), während die 11-Gleichung auf 9·S·Q·Fr/r³=0 kollabiert. Da S>0 aus der Signatur: Der Zweig Fr≠0, Q≠0 ist leer. (Hinweis zum Geltungsbereich: die Schritt-62-Faktorisierung setzt die B'-Beziehung und damit Fr≠0 voraus. Stetigkeits-Zusatz: bei C¹-Metrik ist Q stetig, stückweise gemischte Konfigurationen liefern daher nichts Neues.)

Der Fall Fr=0 (rein winkelartiger Deformationsfluss) liefert dagegen etwas: Alle Bedingungen tragen den gemeinsamen Faktor (B²+C)' — entweder Fth=0 (trivial) oder B²+C=s0=konstant, und dann sind alle (LT3G)-Bedingungen für beliebiges Fth(r) erfüllt, bei LT3g≠0. Dieselbe Zeitverschiebung wie in Schritt 67 (f'=B) bildet die Familie auf die diagonale Metrik −dT²+s0dr²+r²dΩ² ab. Deren Geometrie: G00=G11=−(s0−1)/(s0r²), Tangentialkomponenten null, Kretschmann K=4(s0−1)²/(s0²r4) — gekrümmt für s0≠1, mit konstantem Raumwinkeldefizit (Reskalierung R=√s0·r ergibt den Defizitfaktor 1/√s0).

Physikalischer Kern: Die Materiequelle der s0-Familie ist 8πρ=(s0−1)/(s0r²), pr=−ρ, pt=0 (Zustandsgleichung wr=−1). Für s0>1 sind alle klassischen Energiebedingungen erfüllt: NEC radial saturiert, WEC, SEC (marginal), DEC — die erste Konfiguration dieses Projekts, deren Materie keine Energiebedingung verletzt. Geometrie und Quelle entsprechen der asymptotischen Form des globalen Monopols (Barriola–Vilenkin) bzw. einer Letelier-String-Cloud; zu dieser Literatur-Zuordnung siehe den offenen Dissens unten. Mit der Definition „nichttrivial“ := (LT3g≠0 ∧ LT3G=0 ∧ G≠0) ist die s0≠1-Familie die einzige nichttriviale Konfiguration des gesamten Kreuzterm-Systems; für s0≠1 muss der Fluss rein winkelartig sein (Fr=0 erzwungen). Ehrliche Einschränkung: G ist T3-invariant, aber räumlich inhomogen (∼1/r²) — „homogen“ gilt nur im LT3-operationalisierten Sinn dieses Projekts. Ob der Spezialfall B=0, C=konstant im früheren B=0-Programm je betrachtet wurde, ist ungeprüft und offen markiert.
Lokales Komitee (Schritte 68–68d) und ein offener Dissens: Vier Gutachten (qwen3:235b-a22b und deepseek-r1:70b, je zwei Claims, sequentiell orchestriert, ein Claim pro Aufruf): Rechenkette und CAS-Resultate einstimmig bestätigt; drei Präzisierungen der Gutachter wurden übernommen (Geltungsbereich der Schritt-62-Faktorisierung, Stetigkeits-Zusatz, Energiebedingungs-Formulierung statt des Etiketts „nicht-exotisch“). Offener Dissens zur Literatur-Benennung: qwen3 hält die Identifikation mit Barriola–Vilenkin/String-Cloud für ungenau (volle BV-Lösung habe nichttriviales gtt), deepseek-r1 stimmt ihr zu; die hiesige Einordnung bezieht sich auf die asymptotische Defizitform bzw. die sphärisch-symmetrische Letelier-Quelle. Die Zuordnung wird mit dokumentiertem Dissens publiziert (Entscheidung des Projektträgers) und kann durch Prüfung der Originalarbeiten nachgeschärft werden.

Schritte 69–74: Eich-Lemma, pt=0-Klasse, Klassifikationssatz — und die Retro-Korrektur

Eich-Lemma (Schritt 69): Im gesamten statischen Fall (ds²=−A dt²+2B dt dr+C dr²+r²dΩ², Konvention g01=B, A>0) eliminiert T=t−f(r) mit f'=B/A den Kreuzterm exakt (Diagonalform A, C+B²/A; CAS-verifiziert). B ist statisch immer reine Eichung; die Kreuzterm-Resultate der Schritte 59–68d sind rückblickend der mühsam bestätigte Spezialfall A=1. Der statische Rest reduziert sich modulo Eichung auf die Diagonalmetrik mit zwei Funktionen A(r), C(r).

Die pt=0-Klasse (Schritte 70–71): Für rein winkelartigen Fluss (Fr=0) tragen alle Bedingungen den gemeinsamen Faktor W, und es gilt die exakte CAS-Identität Gθθ=−rW/(4A²C²): Die Lösungsklasse ist genau die Klasse aller statischen sphärisch-symmetrischen Metriken mit verschwindendem Tangentialdruck, mit beliebigem Fth(r) bei LT3g≠0. Sie enthält die Defizitwinkel-Familie (Schritte 68c/d) und die volle Barriola–Vilenkin-Monopol-Lösung mit Massenterm (A=1−α−2M/r, C=1/A) — das Komitee (qwen3) hat Letzteres in seinem Gutachten selbst nachgerechnet und damit seinen früheren Literatur-Einwand in der Substanz zurückgezogen. Für radiale Flussanteile (Fr≠0) wird das Diagonalsystem von einer Auflösungshierarchie regiert (00er-Gleichung: A''-freie Kernbedingung E0; 22er: Fth algebraisch; 11er: Fr' linear; 12er: Restbedingung); bei allgemeinem A liegt der Rest nicht im Differentialideal von E0 — ob dort, wie im Fall A=1, versteckte gemeinsame Faktoren die scheinbare Zusatzbedingung kollabieren lassen, ist offen und ausdrücklich nicht behauptet.

Klassifikationssatz A=1 (Schritte 72–73): Für Fr≠0, C'≠0 hinterlässt die Eliminationskette eine Fr-freie Restbedingung P(r,C,C'), kubisch in C' (Fr-Unabhängigkeit doppelt CAS-geprüft und formal aus der Homogenität vom Grad 1 aller vier Bedingungen in (Fr,Fth) belegt — Schritt 74, auf Komitee-Forderung nachgereicht; die Herleitung wurde zusätzlich über einen zweiten, unabhängigen Rechenweg kreuzvalidiert). P faktorisiert als P=g·qP: Der Zweig g=0 (C'=2C(C−1)/r) ist unter dem E0-Fluss invariant und geschlossen lösbar, C(r)=1/(1−kr²) — Raumschnitte konstanter Krümmung mit konstantem gemischtem Einstein-Tensor diag(−3k,−k,−k,−k), perfektes Fluid w=−1/3. Aber die 11er/22er-Bedingungen erzwingen (bis auf Gesamtskala) Fr=Fth=√(1−kr²), und dieser Fluss ist ein Killing-Vektor der Konstantkrümmungs-Geometrie: LT3g=0 identisch — die Familie ist im Projektsinn trivial. Der Zweig qP=0 ist nicht invariant; die Resultante mit seinem Fluss-Rest, 192(5C²−4C+4)(C−1)5C³r², hat keine reellen nichttrivialen Nullstellen. Satz: Die einzige nichttriviale Konfiguration des A=1-Systems ist die Defizitwinkel-Familie (Fr=0); für Fr≠0 existiert keine.

Retro-Korrektur: Die frühere Shooting-Trajektorie wurde nur aus der 00er- und 11er-Bedingung konstruiert (Integrabilitäts-ODE, Schritte 6–9); die einzigen 12er-Prüfungen des Altprogramms betrafen Potenzansätze in isotroper Eichung, nie die Trajektorie selbst. Nach dem Klassifikationssatz kann sie das volle System nicht lösen (sie ist nicht von der Form 1/(1−kr²): C∼a·r bei r→0, r*-Struktur); das Komitee bestätigte dies einstimmig, qwen3 mit eigener Asymptotik-Gegenrechnung (Rest ∼−2akr³≠0) und explizitem Ausschluss von Rettungspfaden (Normierung, additive Killing-Anteile, stückweise Konstruktionen). Damit sind die Ergebnisse der Schritte 27–58 samt Nachträgen 64/66 als Antworten auf die Projektfrage widerrufen; als Analyse der reduzierten {00,11}-ODE bleiben sie mathematisch korrekt und dokumentiert.
Lokales Komitee (Schritte 69–74) und übernommene Präzisierungen: Vier orchestrierte Gutachten (qwen3:235b-a22b, deepseek-r1:70b; ein Claim pro Aufruf; ein done_reason=length-Ereignis wurde vom Auto-Retry mit verdoppeltem Token-Budget aufgefangen). Übernommen: expliziter Beleg der Fr-Unabhängigkeit (Schritt 74), explizite Konventionsangabe g01=B, vorzeichenpräzise W-Identität. Nicht übernommen (dokumentierte Gutachter-Überkonfidenz): qwen3s Behauptung, im allgemeinen A-Fall gäbe es keine versteckten gemeinsamen Faktoren — das ist ungeprüft und bleibt offen. Eigene Fehler dieser Etappe, transparent: ein konzeptionell falscher Vorab-Anker (widerlegt durch die eigene Rechnung), ein Sage-API-Fehlgebrauch (.has() auf Funktionsableitungen) und ein NameError; alle dokumentiert in den Sicherungsdateien.

Schritt 75: Der W2-Zweig — ein erster nichttrivialer Fall bei A≠1

Der A''''-Koeffizient des reduzierten 12er-Rests R (Fr≠0-Zweig, allgemeines A(r), aus Schritt 71/74) faktorisiert exakt als −4r⁵A³C·W·W2, mit W wie aus Schritt 70/71 (W=0 ⇔ pt=0) und der neuen Größe W2=AC−rA'−A. Auf {W=W2=0} lässt sich R=0 nicht mehr direkt nach A'''' auflösen — genau dort, wo im A=1-Fall (Schritt 72) eine scheinbare Zusatzbedingung auf eine Killing-Familie kollabierte, war die naheliegende nächste Frage.

Spot-Test AC=1 (Schwarzschild-Klasse): hier faktorisiert der 12er-Rest als [Polynom in A und Ableitungen bis A''''] · Fr(r) — Fr fällt vollständig heraus, die Nichttrivialität reduziert sich auf eine reine, von Fr entkoppelte ODE 4. Ordnung für A allein. Ob diese Lösungen hat, ist noch nicht untersucht.

W2=0 legt C algebraisch durch A fest: C=(A+rA')/A. Anders als der g-Zweig im A=1-Fall löst dies E0 nicht automatisch — es entsteht eine echte ODE 3. Ordnung für A. Der Potenzgesetz-Ansatz A=rp faktorisiert diese Bedingung exakt zu (p+1)²(p−2)p·r3p−1=0: p=−1 ist entartet (C=0, unzulässig), p=0 der triviale A=1-Fall, p=2 ist neu.

A=r², C=3 — vollständig verifiziert: E0=0 und der volle 12er-Rest R=0 (beide unabhängig reproduziert). Aus 11er/22er: Fr'=2Fr/r ⇒ Fr(r)=λr², Fth=0. Direkt auf der konkreten Metrik ds²=−r²dt²+3dr²+r²dΩ² berechnet: LT3g≠0 (kein Killing-Vektor; explizite Komponenten ungleich 0), Ricci-Skalar=0 aber Ricci-Tensor≠0 (kein Vakuum) und Riemann-Tensor≠0 (echt gekrümmt), und — unabhängig von der allgemeinen Herleitungskette direkt nachgerechnet — LT3G=0 in allen zehn Komponenten. Erste im Projekt verifizierte nichttriviale Konfiguration mit A≠1.
Lokales Komitee (Claim K5): qwen3:235b-a22b und deepseek-r1:70b, sequentiell, adversarialer Prüfauftrag. qwen3 lehnte formal ab; die eigene Rechnung im Gutachten (u.a. ein selbst berechneter Kretschmann-Skalar ≠0) bestätigte jedoch den Kernbefund, und von fünf Prüfpunkten war nur einer substantiell — der fehlende unabhängige Nachweis, dass E0=0 und R=0 alle vier LT3G-Komponenten garantieren. Dieser wurde durch die oben genannte direkte Neuberechnung (unabhängig vom allgemeinen A(r),C(r)-Formalismus) geschlossen. deepseek-r1:70b lieferte eine methodisch schwache, wenig adversariale Antwort mit geringem Evidenzwert. Ausdrücklich offen, nicht behauptet: ob A=r² ein isolierter Punkt oder Teil einer größeren Lösungsfamilie ist (die W2=0-Bedingung ist eine ODE 3. Ordnung, generisch eine Drei-Parameter-Familie); Verhalten bei r→0 und r→∞ sowie Energiebedingungen sind nicht untersucht.

Schritt 76: Von einem Punkt zu einer Familie — und ihre Grenzen

Die W2=0-Bedingung (Schritt 75, C=(A+rA')/A) reduziert sich, ausgenutzt über die Skalensymmetrie A→μA (Substitution A=exp(f), dann g=f', dann g=h(r)/r, x=ln r), auf eine autonome Gleichung 2. Ordnung: (h+1)h''−2h'²+(h²−1)h'+h⁴−3h²−2h=0. Bei h'=h''=0 faktorisiert das exakt zu h(h+1)²(h−2)=0 — die drei bekannten Potenzgesetz-Wurzeln p∈{−1,0,2}. Das Phasenraumsystem (h,w=dh/dx) hat drei Fixpunkte: h=0 (A=1) ist ein hyperbolischer Sattel, h=2 (A=r²) eine stabile Spirale (Eigenwerte −½±i√23/2), h=−1 (C=0) ist singulär.

Keine Scheinfamilie: Der A''''-unabhängige Rest der vollen 12er-Bedingung R verschwindet identisch für JEDES A(r) auf W2=0 — R=0 folgt allein aus E0=0. Die gesamte 3-Parameter-Familie (zwei Integrationskonstanten aus dem Phasenraum, eine aus der A-Reskalierung) besteht vollständig aus echten Lösungen, nicht nur der isolierte Punkt A=r².

Numerisches Phasenporträt (scipy, breiter Startbereich): Für r→∞ konvergiert ein breiter Bereich von Anfangsbedingungen auf h=2 — außer stark negative Startwerte, die stattdessen in einen Blow-up nahe h=−1 laufen (endliches r!). Für r→0 laufen alle zu h=2 spiralenden Trajektorien auf h→−1 zu. Kein universeller, aber ein robuster Attraktor.

Energie-Impuls-Tensor: p_r≡0 identisch auf der ganzen Familie (strukturell, da Grr∝−W2). Entlang der Attraktor-Familie (h zwischen 1,5 und 2,5): ρ>0, p_t>0 — Standard-Energiebedingungen (WEC, NEC) durchgehend erfüllt.

Bei r→0 (korrigiert, siehe Schritt 79): Der Ricci-Skalar verschwindet entlang der Attraktor-Familie (~r, an zwei unabhängigen Trajektorien numerisch bestätigt) — über diesen Invarianten zeigt sich keine Singularität. Der Kretschmann-Skalar dagegen divergiert echt, K~r⁻⁶ — dieselbe Potenz, die das Projekt bereits einmal bei seiner früheren, inzwischen widerrufenen r=0-Singularität fand. Es gibt also weiterhin eine echte Krümmungssingularität bei r→0, strukturell ähnlich wie bei Schwarzschild (Ricci verschwindet, Kretschmann divergiert) — ohne dass die Geometrie damit selbst Schwarzschild wäre.

Lokales Komitee (Claims K6, K7): Zwei Runden, qwen3:235b-a22b und deepseek-r1:70b, adversarialer Prüfauftrag. In K6 identifizierte qwen3 zwei eigene Rechenfehler beim Nachrechnen der Skalen-Homogenität und der Autonomie-Behauptung (bei genauem Nachrechnen widerlegt), aber einen echten offenen Punkt: fehlender unabhängiger Nachweis von LT3G=0 entlang der ganzen Familie (nicht nur am Einzelpunkt) — geschlossen durch eine allgemeine Formel LT3g[t,t]=−A'·Fr·cosθ, die direkt zeigt, dass die Metrik überall auf der Familie (außer am bekannten Trivialpunkt) nicht-Killing ist. In K7 identifizierte qwen3 einen echten Fehler in einer Zwischenversion: die Behauptung, Energiebedingungsverletzung falle mit numerischer Instabilität zusammen, war zu pauschal (Gegenbeispiel: h<−1 ist energiebedingungskonform, aber ebenso instabil) — korrigiert. qwen3s Vorschlag eines Krümmungstests führte direkt zum Befund der echten Singularität bei r→0. deepseek-r1:70b blieb in beiden Runden methodisch schwach, ohne erkennbare eigenständige Nachrechnung. Offen: eine Frobenius-Analyse der Singularität bei r→0, der AC=1-Zweig (Schritt 75b) und eine systematischere (gitterbasierte statt punktuelle) Phasenraum-Kartierung.

Schritt 77: Der AC=1-Zweig — Schwarzschild in Verkleidung

Der in Schritt 75b als billiger Spot-Test eingeführte AC=1-Ansatz (C=1/A, Schwarzschild-Klasse) wurde weiterverfolgt. Erste Lücke geschlossen: Die 00er-Bedingung (E0) war für diesen Zweig nie geprüft worden. Ergebnis: Der Zweig ist überbestimmt — die "reine ODE 4. Ordnung für A" aus Schritt 75b (Q=0) muss gemeinsam mit E0=0 (2. Ordnung) gelten.

Elimination von A'',A''',A'''' über E0 reduziert Q auf einen Rest, der sauber faktorisiert: [kubisches Polynom in A'] · (rA'+A−1)³ · (rA'−A). Ein Invarianztest (dieselbe Methodik wie beim A=1-Fall, Schritt 72) zeigt: (rA'−A)=0 ist nicht invariant unter dem E0-Fluss — Scheinlösung. (rA'+A−1)=0 dagegen ist invariant und löst sich geschlossen: A=1+K/r.

Vollständig verifiziert, dann als trivial erkannt: E0=0 und Q=0 sind für A=1+K/r beide exakt erfüllt (nicht nur das Invarianzargument). Fth=Fr/2, Fr(r)=λ√(r(r+K)) explizit hergeleitet. Aber: A=1+K/r, C=1/A ist exakt die Schwarzschild-Metrik (K=−2M) — und Schwarzschild ist reines Vakuum (G≡0 CAS-bestätigt: Ricci-Tensor, Ricci-Skalar und voller Einstein-Tensor identisch 0). Bei G≡0 gilt LT3G=0 automatisch für jedes Vektorfeld, nicht wegen echter Struktur — dieselbe Vakuum-Falle, die früh im Projekt bei der ursprünglichen Schwarzschild-Erkenntnis dokumentiert wurde. Der Kandidat ist damit trivial im Projektsinn, trotz LT3g≠0.

Auch der dritte Faktor verworfen: Ein Invarianztest (Resultante von F und dF/dr bezüglich A', mit A'' aus E0 substituiert) zeigt: Die Resultante verschwindet nicht identisch, sondern nur bei isolierten A-Werten (A=0 entartet, A=1 trivial, oder Nullstellen einer Quartik) — der Faktor ist nicht invariant unter dem E0-Fluss, ebenfalls eine Scheinlösung. Damit ist der gesamte AC=1-Zweig vollständig und negativ klassifiziert: keiner der drei Faktoren liefert eine echte, nichttriviale, nicht-Vakuum-Lösung — methodisch analog zum Klassifikationssatz für A=1 (Schritte 72–73).

Schritt 79: Korrektur — Ricci verschwindet, aber Kretschmann divergiert bei r→0

Die in Schritt 76 berichtete Aussage "der Ricci-Skalar divergiert bei r→0" beruhte auf einer fehlerhaften Methodik: ausgewertet wurde bei festem r=1 für verschiedene h-Werte, mit jeweils künstlich bei r=1 auf A=1 normiert — das entspricht nicht dem tatsächlichen r→0-Verhalten EINER echten Trajektorie, bei der A, A', A'', h, w gemeinsam mit r variieren.

Korrigierte Methodik: Zunächst wurde die Annäherung h→−1 (x=ln r→−∞) analytisch charakterisiert — Ansatz ε=h+1~B·eγx in die autonome Gleichung eingesetzt, charakteristische Gleichung γ²+2γ−3=0 (CAS-bestätigt), relevante Wurzel γ=1, also h+1~B·r. Danach zwei unabhängige Trajektorien numerisch integriert (h, w UND f=ln A gemeinsam, A',A'' konsistent aus den etablierten Formeln rekonstruiert, nicht unabhängig normiert) — beide zeigen: die Steigung von ln|Ricci| gegen x konvergiert sauber auf ≈1,0 (Ricci~r, verschwindet), keine Divergenz.

Ein zweiter, unabhängiger eigener Fehler bei der Gegenprüfung gefunden: Die Kretschmann-Berechnung (hier und bereits beim Komitee-Gutachten K5 in Schritt 75) beruhte auf einer falschen Annahme über den von SageManifolds gelieferten Tensortyp (angenommen (0,4), tatsächlich (1,3)) — die alte Formel enthielt einen nicht herausfallenden sin(θ)-Term, unmöglich für eine echte Invariante. Korrigiert über Riem.down(g,0) (macht (1,3) zu (0,4)), dann vierfaches .up(), dann Kontraktion; die neue Formel ist nachweislich θ-unabhängig und dreifach verifiziert (allgemeine Formel, direkte Rechnung an der konkreten Metrik A=r²,C=3, Handrechnung aus den einzelnen Riemann-Komponenten — alle drei ergeben K=8/3 bei r=1, nicht den zuvor im Komitee-Gutachten K5 genannten Wert 80/9, der ohne CAS-Werkzeuge "aus eigener Rechnung" entstand und sich als fehlerhaft erweist). Mit der korrigierten Formel an derselben Trajektorie ausgewertet: die Steigung von ln(Kretschmann) gegen x konvergiert auf −6,000 — K~r⁻⁶, dieselbe Potenz, die das Projekt bereits einmal bei seiner früheren, widerrufenen r=0-Singularität fand (Schritt 55).

Korrigiertes Ergebnis: Ricci-Skalar → 0 bei r→0 (verschwindet, keine Singularität über diesen Invarianten), aber Kretschmann-Skalar divergiert echt, K~r⁻⁶ — es gibt weiterhin eine echte Krümmungssingularität, nur über den richtigen Invarianten nachgewiesen. Strukturell ähnlich wie bei Schwarzschild (Ricci≡0 überall, Kretschmann divergiert bei r=0) — passend zum unabhängigen Befund, dass der AC=1-Zweig (Schritt 77) tatsächlich Schwarzschild ist, aber keine Identität der beiden Geometrien (hier ist die Gesamtgeometrie nicht Vakuum).

Lokales Komitee (Claim K8): qwen3:235b-a22b empfahl zunächst den vollständigen Rückzug der Korrektur. Bei genauer Prüfung erwiesen sich die zentralen technischen Einwände selbst als nicht haltbar: Der Vorwurf, die charakteristische Gleichung sei durch nichtlineare Terme ungültig, verkennt, dass alle drei führenden Terme konsistent von derselben Ordnung sind — das ist die korrekte Dominant-Balance-Methode an einem entarteten Punkt, keine Verletzung davon. Der Vorwurf eines Fehlers in der Kretschmann-Korrektur enthielt einen sich selbst widersprechenden Satz und schlug einen im Kern identischen Code als "Korrektur" vor; die dreifache unabhängige Verifikation (s.o.) widerlegt den Einwand mit einer harten Zahl. Übernommen wurden die berechtigten Vorsichtspunkte: mehr als zwei Trajektorien wären für eine vollständig generische Aussage wünschenswert, und die Schwarzschild-Analogie ist strukturell, keine Identität. deepseek-r1:70b stimmte pauschal und ohne erkennbare eigene Rechnung zu, wie in fast allen bisherigen Runden.

Nächste Schritte

Das Projekt ist damit nicht abgeschlossen, sondern der erste Abschnitt eines längeren Untersuchungsprogramms. Konkret anstehend:

Einordnung

Lavinia Heisenberg ist Professorin am Institut für Theoretische Physik der Universität Heidelberg und Mitbegründerin mehrerer aktueller Forschungsrichtungen in der Gravitationsphysik (unter anderem der sogenannten „Generalisierten Proca-Theorie“ und der „Geometrischen Trinität der Gravitation“, die zeigt, dass sich Einsteins Gravitation gleichwertig über Krümmung, Torsion oder Nichtmetrizität beschreiben lässt). Das hier behandelte Paper (arXiv:2512.16934) ist eher konzeptioneller, grundlagenorientierter Natur als ihre übrige, meist modellbauende Arbeit — und genau deshalb ein besonders geeigneter Ausgangspunkt für ein Testprojekt: Die offene Frage stammt direkt und ausdrücklich von der Autorin selbst, nicht aus einer eigenen, unabgesicherten Vermutung.

Werkzeuge und Ablage

T. Riepe, Berlin. Juli 2026. Quelle der offenen Frage: L. Heisenberg, „A Re-Examination Of Foundational Elements Of Cosmology“, arXiv:2512.16934 (Dezember 2025), Abschnitt 5 / Conclusio. Angeregt durch den ZEIT-Podcast „Alles gesagt?“ mit Lavinia Heisenberg (Juli 2026); Podcast-Inhalt selbst wurde für dieses Projekt nicht ausgewertet, nur das Paper direkt. Berechnung: SageMath/SageManifolds (lokal, Mac Studio). Koordination / Interpretation (Schritte 59–62): Claude Sonnet 5 (Anthropic), auf Anweisung und mit Gegenprüfung durch T. Riepe. Interpretation der Kernbehauptungen zu Schritt 48–58 (Wurzelwahl, Fortsetzung über die Koordinatensingularität, Krümmungssingularität bei r=0, Tipler-Kriterium, Bronnikov-Zaslavskii-Vergleich): nachträglich vollständig durch ein lokales Modell-Komitee (qwen3:235b-a22b, DeepSeek-R1-0528-671B, beide lokal auf dem Mac Studio) geprüft und per CAS (Sage) gegengerechnet — Claude fungierte dabei nur noch als Koordinator zwischen den lokalen Werkzeugen, nicht als interpretierende Instanz (siehe Schritt 63–64, inklusive eines dabei aufgedeckten und korrigierten Fehlers bei Behauptung 4). Interpretation der Schritte 59–62 (Kreuzterm B≠0, Fth-Quadratur) erfolgte weiterhin durch Claude direkt, ohne nachträgliche lokale Komitee-Prüfung. Unabhängige Gegenrechnung eines konkreten Zwischenschritts (Vorzeichen von C für Fr=r³⁄²) durch das lokale Modellkomitee: qwen3:235b-a22b bestätigte das Sage-Ergebnis exakt; deepseek-r1:70b unterlief dabei ein eigener Rechenfehler (c=−5r³ statt −r³) bei formal unverändertem „Stimme zu“ — dokumentiert als Beleg dafür, dass auch das lokale Komitee nicht unfehlbar ist. Schritte 48–55 (Verifikation der Wurzelwahl, Fortsetzung über r* hinaus, Krümmungssingularität bei r=0): zusätzlich Julia/DifferentialEquations.jl (BigFloat, 256 Bit); diese Sitzung enthält einen dokumentierten eigenen Vorzeichenfehler (falscher Zweig in der v=1/u-Regularisierung), der durch Vorzeichen-Widerspruch mit der numerischen Trajektorie aufgedeckt und korrigiert wurde — explizit als Korrektur festgehalten, nicht stillschweigend ersetzt. Schritte 56–58 (Eigenabstand, Literaturvergleich, Ricci-Tensor-Komponenten, Tipler-Kriterium): Literatur vollständig gelesen (arXiv:0801.0889, arXiv:0805.1095), nicht nur Abstracts. Schritte 59–62 (Kreuzterm B≠0, Fth-Quadratur): Kernsystem noch offen, ehrlich als solches gekennzeichnet. Schritt 67 (Q=0-Sonderzweig): Interpretation vor Publikation durch das lokale Komitee geprüft (qwen3:235b-a22b, deepseek-r1:70b, einstimmig), CAS-Verifikation per Sage (Riemann-Tensor, Pullback-Identität); eine von qwen3 aufgedeckte Vorzeichen-Inkonsistenz in der Claim-Formulierung ist dokumentiert. Schritte 68–68d: vier orchestrierte Komitee-Gutachten (qwen3:235b-a22b, deepseek-r1:70b), Kern einstimmig, drei Präzisierungen übernommen, ein offener Dissens zur Literatur-Benennung transparent dokumentiert; Publikationsentscheidung (Option C) durch Thomas. Schritte 69–74: vier orchestrierte Gutachten (qwen3:235b-a22b, deepseek-r1:70b), Retro-Korrektur einstimmig bestätigt, drei Präzisierungen übernommen, eine Gutachter-Überkonfidenz dokumentiert zurückgewiesen; Korrektur-Rahmung (Leseführung auf aktuelle Resultate) nach Entscheidung von Thomas.