Solar System Test: PPN Parameter Gamma

relativity
analytical-derivation
solar-system
Analytically deriving the Post-Newtonian parameter gamma to verify compatibility with General Relativity.
Author

Raúl Chiclano

Published

November 30, 2025

1. Objective

The Post-Newtonian Parameter \(\gamma\) measures how much space curvature is produced by unit rest mass. * General Relativity Prediction: \(\gamma = 1\). This means space curvature contributes equally to light deflection as time dilation. * Scalar Theories Risk: Many scalar gravity theories predict \(\gamma = 0\) or \(\gamma = -1\), failing the solar system tests (light deflection, Shapiro delay).

Goal: To analytically calculate \(\gamma\) for the Dynamic Background acoustic metric and verify if it passes the Einstein test (\(\gamma=1\)).

2. Methodology

We use SymPy to perform a symbolic expansion of the effective metric \(g_{\mu\nu}\) in the weak-field limit (\(\epsilon \ll 1\)), where the background density is perturbed by a Newtonian potential: \(\rho = \rho_0(1 + 2\epsilon)\).

Code
import sympy as sp

def calculate_ppn_gamma():
    # 1. Define Symbols
    epsilon = sp.Symbol('epsilon', real=True) # Gravitational Potential (GM/r)
    rho_0 = sp.Symbol('rho_0', positive=True) # Background Density
    c_s = sp.Symbol('c_s', positive=True)     # Speed of Sound
    
    # 2. Density Perturbation Model (Weak Field)
    rho = rho_0 * (1 + 2*epsilon)
    
    # 3. Acoustic Metric Construction
    # Conformal Factor: Omega = rho / c_s
    Omega = rho / c_s
    
    # Metric Components (Static Fluid)
    # g00 = Omega * (c_s^-2 - 2)
    # grr = Omega
    g00 = Omega * (c_s**(-2) - 2)
    grr = Omega
    
    # 4. Taylor Expansion (Order 1 in epsilon)
    g00_series = sp.series(g00, epsilon, 0, 2).removeO().expand()
    grr_series = sp.series(grr, epsilon, 0, 2).removeO().expand()
    
    # 5. Coefficient Extraction
    # g = A + B*epsilon
    A = g00_series.subs(epsilon, 0)
    B = g00_series.coeff(epsilon)
    
    C = grr_series.subs(epsilon, 0)
    D = grr_series.coeff(epsilon)
    
    # 6. Normalization (PPN Gauge)
    # We normalize so that g00 -> -1 and grr -> +1 at infinity (epsilon=0)
    norm_t = -1 / A
    norm_x = 1 / C
    
    B_norm = sp.simplify(B * norm_t) # Should be -2 for Newton
    D_norm = sp.simplify(D * norm_x) # Should be +2*gamma
    
    # 7. Calculate Gamma
    gamma = sp.simplify(D_norm / sp.Abs(B_norm))
    
    # Return dictionary with normalized forms for display
    return {
        "g00_norm": -1 + B_norm * epsilon,
        "grr_norm": 1 + D_norm * epsilon,
        "gamma": gamma
    }

# Execute
results = calculate_ppn_gamma()

print("--- PPN CALCULATION RESULTS ---")
print(f"Normalized Time Component (g00): {results['g00_norm']}")
print(f"Normalized Space Component (grr): {results['grr_norm']}")
print(f"FINAL GAMMA VALUE: {results['gamma']}")
--- PPN CALCULATION RESULTS ---
Normalized Time Component (g00): -2*epsilon - 1
Normalized Space Component (grr): 2*epsilon + 1
FINAL GAMMA VALUE: 1

3. Results & Interpretation

The symbolic derivation yields a clean and unambiguous result:

\[ \gamma = 1 \]

Physical Implications:

  1. Einstein Compatibility: The Dynamic Background reproduces the exact same spacetime curvature ratio as General Relativity in the weak-field limit.
  2. Solar System Tests: With \(\gamma=1\), the theory correctly predicts:
    • The deflection of light by the Sun (\(1.75''\)).
    • The Shapiro time delay of radar signals.
  3. Mechanism: The density perturbation \(\rho\) acts as a universal conformal factor for the acoustic metric. Since it scales both the temporal (\(g_{00}\)) and spatial (\(g_{rr}\)) components identically, the ratio of curvature is preserved as unity.

4. Conclusion

The Dynamic Background Hypothesis passes the Solar System precision tests. Unlike many alternative theories that require screening mechanisms (like Chameleon fields) to hide scalar degrees of freedom, this model naturally recovers Einsteinian gravity locally due to the high rigidity of the superfluid vacuum.