FINAL_PHYSICS_FIX_SUMMARY.md

FINAL Physics Plots Fix Summary

Probleme identifiziert:

1. g₁/g₂ Plot

✅ KORREKT - Verwendet G79 Daten in pc, piecewise fit funktioniert ⚠️ Überlappungen - Fixed mit spacing 0.18/0.13

2. Time Dilation

❌ FALSCH - Crossover bei 6.0 statt 1.387 Ursache: gamma_seg verwendet r/(r_c*r_s) was für METER ist, aber r_ratio ist dimensionslos

3. Radial Stretch

❌ ERROR - Funktion existiert und ist korrekt definiert Ursache: Gleiche gamma_seg Problem

4. Combined Analysis

❌ LEER - Gleiche gamma_seg Problem

Root Cause:

Die gamma_seg Funktion in ssz_physics_plots.py Zeile 37-39:

def gamma_seg(r, r_s, alpha=ALPHA, r_c=R_C):
    return 1 - alpha * np.exp(-(r/(r_c*r_s))**2)

Verwendet r/(r_c*r_s) was bedeutet:

• r in METER
• r_c dimensionslos (in Einheiten von r_s)
• r_s in METER

ABER: In PAPER-RESTORED g79_temperature_plots.py wird verwendet:

def gamma_seg(r, alpha=ALPHA, r_c=R_C):
    return 1.0 - alpha * np.exp(-(r / r_c)**2)

Mit:

• r in PARSEC
• r_c in PARSEC

Lösung:

Zwei separate Funktionen:

1. gamma_seg_meter(r, r_s, alpha, r_c) - für r in Meter, r_c in r_s
2. gamma_seg_pc(r, r_c, alpha) - für r in pc, r_c in pc

G79 Plot verwendet pc-Version, andere Plots verwenden meter-Version.

Die RICHTIGE Normalisierung für dimensionslosen r/r_s Plot:

• gamma_seg sollte (r/r_s) / r_c verwenden
• NICHT r / (r_c * r_s)

Mit r_c = 1.9 und r/r_s = 1.387:

• Normalisiert: 1.387 / 1.9 = 0.73
• Xi = alpha * exp(-(0.73)²) = 0.12 * exp(-0.53) = 0.07
• D = 1/(1+0.07) = 0.93 -D_time_dilation = sqrt(0.93 * (1-1/1.387)) = sqrt(0.93 * 0.28) = sqrt(0.26) = 0.51

Das ist NAH an 0.528!

Ich muss die Formel korrigieren.