Is Space "Lumpy"? Probing Quantum Gravity with Python

This Python script explores a frontier of modern physics known as Lorentz Invariance Violation (LIV), which investigates whether the "smooth" fabric of Einstein's spacetime becomes "lumpy" or "foamy" at incredibly small scales. By calculating the arrival time differences of gamma rays from distant cosmic sources like quasars, the code determines a characteristic length scale for this spatial graininess. While classical physics assumes all light travels at exactly $c$ in a vacuum, quantum gravity theories suggest that higher-energy photons might "bump" into the quantum foam of space, causing a measurable lag over billions of light-years.

The Python Implementation

To ensure Google and other search engines can "read" this code, I have embedded it as text within a preformatted block. This allows the logic to be searchable while remaining easy to read.

# Constants
C = 29979245800.0          # Speed of light in cm/sec
LY_TO_CM = 9.46073e17      # 1 Light Year in cm
PLANCK_LENGTH = 1.616e-33  # Theoretical "quantum" grain size (cm)
NUCLEUS_SIZE = 1e-13       # Average size of an atomic nucleus (cm)

def analyze_space_grain(delta_t, distance, unit, wl_diff_cm):
    """
    Calculates spatial lumpiness (L) based on the observed time delay 
    of photons with different wavelengths.
    """
    dist_cm = distance * LY_TO_CM if unit.lower() in ["ly", "light year"] else distance
    lumpiness = (delta_t * C * wl_diff_cm) / dist_cm
    return lumpiness

# Dataset of Cosmic Observations
rays = {
    "Quasar": (1.0, 1e-9, 5e9, "ly"),           
    "Distant Object": (0.7, 4e-12, 10e9, "ly")
}

print(f"{'Source':<18} | {'Lumpiness (cm)':<15} | {'Classification'}")
print("-" * 65)

for name, (td, wld, d, u) in rays.items():
    lump = analyze_space_grain(td, d, u, wld)
    
    if lump > NUCLEUS_SIZE:
        status = "Large (Nuclear Scale)"
    elif lump > 1e-20:
        status = "Sub-Nuclear Scale"
    elif lump > PLANCK_LENGTH:
        status = "Quantum Scale"
    else:
        status = "Smooth (Planck Limit)"
        
    print(f"{name:<18} | {lump:.2e} cm    | {status}")

Physics Significance

By comparing our results to the size of an atomic nucleus or the theoretical Planck Length, we can set experimental "upper limits" on the smoothness of the universe. If our observations show almost no time delay, it proves that spacetime remains perfectly smooth even at sub-particle dimensions.