A Velocity Coupled Radial Acceleration Ansatz for Disk-Galaxy Rotation Curves: Fits to SPARC, Bayesian Inference, and Parameter Identifiability

This paper introduces and evaluates a phenomenological "velocity-coupled radial acceleration" model for disk-galaxy rotation curves, demonstrating that while it fits SPARC data comparably to standard NFW and Burkert halo models, its parameters suffer from significant degeneracy that limits their identifiability for the majority of galaxies.

The Gist

The Big Picture: The "Missing Weight" Mystery

Imagine you are watching a carousel. If you spin it slowly, the horses on the outside stay put. But if you spin it fast, the horses on the outside should fly off because there isn't enough "grip" (gravity) holding them in.

In real galaxies, stars on the outer edges are spinning incredibly fast. According to our current understanding of gravity (based on the visible stars and gas), they should fly off into space. But they don't. They stay in orbit.

For decades, astronomers have solved this by saying, "There must be invisible 'Dark Matter' holding them together." They imagine a giant, invisible cloud of weight surrounding the galaxy.

This paper asks a different question: What if we don't need to invent a new invisible cloud? What if the rules of how things move are just slightly different than we thought?

The New Idea: The "Speed-Sensitive Tether"

The author, Nalin Dhiman, proposes a new way to describe the extra force holding these stars in.

• The Old Way (Dark Matter Halos): Imagine the galaxy is wrapped in a giant, invisible fishing net (the Dark Matter halo). The net has a specific shape and weight. You have to guess the shape of the net to make the math work.
• The New Way (VCA - Velocity-Coupled Acceleration): Imagine the galaxy is attached to the center by a special, magical rubber band. The tighter the band pulls, the faster the star is moving.
  • If the star moves slowly, the band is loose.
  • If the star moves fast, the band pulls harder to keep it from flying away.

The author calls this a "Velocity-Coupled" force. It's a simple mathematical rule: The faster you go, the harder you are pulled in.

The Experiment: Testing the Theory

The author took 171 real galaxies from a famous database (called SPARC) and ran a massive simulation.

1. The Setup: He took the visible stuff (stars, gas, dust) and calculated how fast the stars should move based on normal gravity.
2. The Test: He added his new "Speed-Sensitive Tether" rule to the math.
3. The Comparison: He compared his new rule against the two most popular "Dark Matter" shapes (called NFW and Burkert).

Think of it like a car race.

• Car A (NFW): A standard sports car with a known engine.
• Car B (Burkert): A modified truck with a known suspension.
• Car C (VCA): A new prototype car with a weird, speed-sensitive engine.

The author drove all three cars on the same 171 tracks (galaxies) to see which one stayed on the road best.

The Results: A Surprising Tie

Here is what happened:

1. It Works: The new "Speed-Sensitive Tether" model fits the data almost as well as the standard Dark Matter models. It successfully predicts how fast the stars are moving without needing to invent a complex invisible cloud.
2. The Winner: The "Burkert" model (the modified truck) was slightly better at fitting the data overall. However, the new "VCA" model was very competitive, often beating the "NFW" model (the sports car).
3. The Catch (The "Fuzzy" Parameters): While the model works, it's hard to pin down the exact settings.
   • Imagine trying to tune a radio. You have two knobs: one for "Volume" and one for "Tone."
   • In this model, the knobs are linked. If you turn one up, you have to turn the other down to get the same sound.
   • For many galaxies, the data wasn't clear enough to tell the difference between the two knobs. The author calls this a "degeneracy." We know the result works, but we aren't 100% sure what the specific settings are for every single galaxy.

The "So What?" Conclusion

The author is very honest about what this means. He says: "This isn't a complete theory of the universe yet."

• It's not a replacement for physics: The new rule doesn't come from a deep, fundamental law of nature (like Einstein's gravity). It's a "phenomenological" rule—meaning it's a mathematical trick that works well to describe what we see, even if we don't know why it works physically yet.
• It's a useful tool: It shows that we can describe galaxy rotation curves with a very simple, compact formula that doesn't require a complex invisible cloud.
• The Mystery Remains: The fact that this simple "speed-dependent" rule works so well suggests that there is a deep, hidden connection between how much visible matter a galaxy has and how fast it spins. Whether that connection is caused by Dark Matter or a modification of gravity is still an open question.

The Takeaway Analogy

Imagine you are trying to explain why a kite stays in the air.

• Standard View: "There is an invisible wind current (Dark Matter) pushing it up."
• This Paper's View: "Maybe the kite is just designed so that the faster it flies, the more lift it generates automatically."

The paper proves that the "automatic lift" idea fits the data just as well as the "invisible wind" idea for many kites. It doesn't prove the wind doesn't exist, but it shows that we might be able to explain the flight with a simpler rule we haven't fully understood yet.

Technical Summary

1. Problem Statement

Disk galaxy rotation curves are a primary empirical probe of the relationship between baryonic matter and galactic dynamics. While standard models attribute the discrepancy between observed rotation speeds and baryonic predictions to Dark Matter (DM) halos (e.g., NFW or Burkert profiles), phenomenological regularities like the Radial Acceleration Relation (RAR) suggest a tight coupling between baryonic and total acceleration.

The paper addresses a specific question: Can a minimal, two-parameter phenomenological ansatz, distinct from a density-profile halo, provide a self-consistent description of circular dynamics that competes with standard DM models in fitting SPARC rotation curves? The authors aim to test a model where the "missing" acceleration is not derived from a mass distribution but is directly coupled to the local tangential velocity.

2. Methodology

2.1 The VCA Model

The authors propose a Velocity-Coupled Acceleration (VCA) model. Instead of adding a mass density profile, they introduce an additional inward radial acceleration term ($a_{vca}$) proportional to the local tangential speed ($v$):
$$a_{vca}(r) = \gamma(r) v(r)$$
where the coupling function $\gamma(r)$ saturates with radius:
$$\gamma(r) = \frac{v_\infty}{r + r_0}$$

• Parameters: $v_\infty$ (asymptotic velocity scale) and $r_0$ (saturation length scale).
• Dynamics: By combining this with the circular motion condition ($v^2/r = g_{bar} + a_{vca}$), the authors derive a quadratic equation for $v(r)$. This yields a closed-form analytical solution for the rotation curve:
  $$v(r) = \frac{1}{2} \left[ A(r) + \sqrt{A(r)^2 + 4v_{bar}^2(r)} \right]$$
  where $A(r) = v_\infty r / (r + r_0)$. This algebraic self-consistency is a key technical feature, avoiding the need for iterative numerical integration of differential equations.

2.2 Data and Pipeline

• Dataset: The SPARC sample (Lelli et al. 2016), comprising $N_{gal} = 171$ galaxies with high-quality rotation curves and homogeneous baryonic decompositions (gas, disk, bulge).
• Baryonic Input: Fixed stellar mass-to-light ratios ($\Upsilon_{disk}=0.5, \Upsilon_{bulge}=0.7$).
• Comparison Models: The VCA model is fitted alongside two standard two-parameter DM halos:
  1. NFW (Navarro-Frenk-White): Cuspy halo motivated by $\Lambda$CDM.
  2. Burkert: Empirically cored halo.
• Optimization: All models use an identical pipeline involving bounded nonlinear least squares (SciPy) and a systematic error floor ($\sigma_0$) added to observational uncertainties.
• Inference: For VCA, the authors perform MCMC inference (affine-invariant ensemble sampler) to map posterior distributions, quantify parameter degeneracies, and assess identifiability.
• Validation:
  • Information Criteria: AIC and BIC for model selection.
  • Cross-Validation: A radial holdout test (fitting inner 70%, predicting outer 30%) to assess generalization.
  • RAR Analysis: Evaluating the model's ability to reproduce the Radial Acceleration Relation ($g_{obs}$ vs. $g_{bar}$).

3. Key Contributions

1. Closed-Form Phenomenological Ansatz: The derivation of an exact, closed-form solution for rotation curves where the "dark" component is a velocity-dependent force rather than a mass density. This allows for rapid, stable fitting without numerical integration.
2. Rigorous Identifiability Analysis: Unlike many phenomenological studies, this paper explicitly quantifies parameter identifiability using MCMC. It reveals that for the majority of galaxies, the parameters $v_\infty$ and $r_0$ are strongly degenerate.
3. Apples-to-Apples Benchmarking: The study applies a unified optimization and error model to VCA, NFW, and Burkert, ensuring a fair comparison of fit quality and predictive power.
4. Physical Interpretation Limits: The authors clearly delineate the model as a kinematic fit rather than a dynamical theory, acknowledging that the velocity-dependent force is non-conservative and does not derive from a scalar potential.

4. Key Results

4.1 Fit Quality and Model Comparison

• Performance: With a fiducial systematic error floor of $\sigma_0 = 5$ km/s, the VCA model is competitive with the NFW halo and performs comparably to the Burkert halo in terms of Information Criteria (AIC/BIC).
• Win Rates: VCA outperforms NFW in ~61% of galaxies (at $\sigma_0=5$) but is less frequently preferred over Burkert. However, a significant fraction of galaxies show $\Delta AIC \le 2$, indicating statistical equivalence between VCA and Burkert.
• Predictive Power: Radial cross-validation (outer 30% holdout) shows that VCA, NFW, and Burkert have comparable median RMSE on unseen outer radii, suggesting VCA does not overfit the inner data significantly.

4.2 Parameter Identifiability

• Degeneracy: MCMC results show a strong $v_\infty$–$r_0$ degeneracy (banana-shaped posteriors) for most galaxies.
• Constrained Subset: Only 47 out of 171 galaxies (27.5%) yield "well-constrained" posteriors (defined by narrow width and sufficient radial lever arm). For the remaining 72.5%, the data cannot uniquely separate the velocity scale from the saturation radius.
• Robust Quantities: While individual parameters are often unconstrained, derived quantities like the effective amplitude $A(R_{max})$ or the velocity at the outermost radius are well-constrained.

4.3 Radial Acceleration Relation (RAR)

• The VCA model naturally reproduces the overall locus of the RAR.
• The RMS scatter of VCA predictions relative to observed acceleration is comparable to NFW (0.196 dex vs. 0.199 dex) but slightly larger than Burkert (0.177 dex).

4.4 Physical Implications

• Asymptotic Behavior: The model naturally transitions to flat rotation curves ($v \to v_\infty$) at large radii without explicitly invoking a halo density profile.
• Effective Mass: In the outer regions, the VCA term mimics an isothermal sphere ($\rho \propto r^{-2}$), while in the inner regions, it mimics a constant density core.
• Limitations: The model is not derived from a conservative potential. It implies a "drag-like" force that does no work on circular orbits but would exchange energy on non-circular trajectories, limiting its applicability to full dynamical simulations or lensing.

5. Significance and Conclusion

The paper establishes that a velocity-coupled acceleration is a viable, compact descriptive tool for disk galaxy rotation curves. It demonstrates that the tight correlation between baryons and dynamics (the RAR) can be captured by a simple kinematic rule without assuming a specific Dark Matter density profile.

Key Takeaways:

• Competitiveness: VCA is a strong empirical competitor to standard DM halos (NFW/Burkert) for fitting rotation curves.
• Phenomenology vs. Theory: The authors emphasize that VCA is a phenomenological ansatz, not a complete dynamical theory. It serves as a diagnostic for baryon-kinematics couplings rather than a replacement for physical halo models.
• Identifiability Warning: The strong parameter degeneracy suggests that interpreting $v_\infty$ and $r_0$ as fundamental physical constants for individual galaxies is often unjustified; the model effectively constrains the function of the acceleration rather than its specific parameters.

The study concludes that while VCA does not yet replace physical halo models, it provides a valuable, mathematically tractable framework for exploring the empirical regularities of galaxy dynamics and hints at underlying mechanisms connecting baryonic matter to kinematics. Future work would require a covariant formulation to address energy conservation and applicability to non-circular motions.