Galactic Dynamics: Deriving the MOND Regime

Investigating if the superfluid vacuum naturally generates non-Newtonian gravity at low accelerations.

1. Objective

The Missing Mass Problem in galaxies is usually solved by adding Dark Matter. An alternative is Modified Newtonian Dynamics (MOND), which changes gravity at low accelerations (\(a < a_0\)). * AQUAL Theory: A relativistic version of MOND where the Lagrangian depends non-linearly on the kinetic term \(X = (\nabla \phi)^2\). * Goal: To check if the effective Lagrangian of the Dynamic Background contains non-linear terms in \(X\) that could mimic Dark Matter.

2. Methodology

We perform a hydrodynamic reduction of the scalar field Lagrangian. By integrating out the density degrees of freedom (\(\rho\)), we derive an effective action \(\mathcal{L}_{eff}(X)\) for the gravitational potential \(\theta\).

import sympy as sp

def derive_mond_lagrangian():
    # 1. Define Fields
    rho = sp.Symbol('rho', positive=True) # Fluid Density
    X = sp.Symbol('X', positive=True)     # Kinetic Term ~ (Gravity)^2
    alpha = sp.Symbol('alpha', real=True) # Mass Parameter
    beta = sp.Symbol('beta', positive=True) # Self-interaction
    
    # 2. Scalar Field Lagrangian (Hydrodynamic Form)
    # L = Kinetic - Potential
    V = alpha * rho + beta * rho**2
    L = -0.5 * rho * X - V
    
    # 3. Equation of State (dL/drho = 0)
    # The fluid adjusts its density to minimize action
    dL_drho = sp.diff(L, rho)
    rho_sol = sp.solve(dL_drho, rho)[0]# 4. Effective Lagrangian L_eff(X)
    # Substitute rho back into L
    L_eff = L.subs(rho, rho_sol)
    L_eff = sp.expand(L_eff) # Expandir para ver el término X^2
    
    # 5. Response Function mu(X)
    # In Newton, mu = constant. In MOND, mu depends on X.
    mu_function = sp.diff(L_eff, X) # Derivar respecto a X
    
    return {
        "L_original": L,
        "rho_equilibrium": rho_sol,
        "L_effective": L_eff,
        "mu_response": mu_function
    }

# Execute
results = derive_mond_lagrangian()

print("--- EFFECTIVE LAGRANGIAN RESULTS ---")
print(f"L_eff(X): {results['L_effective']}")
print(f"Response Function mu(X): {results['mu_response']}")

--- EFFECTIVE LAGRANGIAN RESULTS ---
L_eff(X): 0.0625*X**2/beta + 0.25*X*alpha/beta + 0.25*alpha**2/beta
Response Function mu(X): 0.125*X/beta + 0.25*alpha/beta

3. Results & Interpretation

The symbolic derivation reveals a non-trivial structure:

\[ \mathcal{L}_{eff}(X) \propto X^2 + \dots \]

Physical Implications:

1. Non-Newtonian Behavior: The effective Lagrangian is quadratic in X, not linear. This means the theory belongs to the AQUAL (A Quadratic Lagrangian) class of modified gravity theories.
2. Variable Stiffness: The response function \(\mu(X)\) is not constant. The “stiffness” of spacetime depends on the local gravitational field strength.
3. Dark Matter Mimicry: This non-linearity provides a natural mechanism to explain the flat rotation curves of galaxies and the anomalies in Wide Binaries (Chae et al.) without invoking particle Dark Matter. The vacuum itself acts as a “phantom dark matter” halo.

4. Conclusion

The Dynamic Background naturally transitions into a Modified Gravity regime at low accelerations. This unifies the success of General Relativity in the Solar System (high acceleration) with the MOND-like phenomenology observed in galaxies (low acceleration).