Skip to main content
Sandbox Physics
Galaxy rotation curves and mass models

Galaxy rotation curves · L4

Reading feedback, halo and MOND comparisons

A halo–modified-dynamics comparison starts with its observations, fixed parameters and error treatment. A 2026 paper by Smith and colleagues puts these choices in one study that we can examine figure by figure.

Record the version and the 147-galaxy selection

The source is Amy Smith and colleagues’ “Dark Matter and MOND successes and degeneracies in reproducing SPARC rotation curves and weak lensing data”, arXiv:2610.33157v1. Its watermark is dated 23 October 2026. The abstract page says MNRAS Accepted; this guide follows the actual 26-page preprint. The PDF’s 2025 typesetting year is not used as a publication year.

Section 2 selects 147147 late-type SPARC galaxies, explaining quality 3, low-inclination and small-sample exclusions, plus choices where total speed squared becomes invalid. Those conditions define which population enters the comparison. Read them before interpreting a distribution of fit scores.

Our fixed NGC 2403 branch retains its own 73 points and does not import the paper’s multi-galaxy cuts. One named object’s NFW/core analysis and a population-level model comparison ask different questions.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Lelli, McGaugh & Schombert (2016), SPARC, AJ152,157

Figure 1 shows density; feedback enters through a simulation-based profile

Figure 1 plots radius in kiloparsecs against halo density on logarithmic axes. Identify curves by model name in the legend and note the illustrative halo mass/luminosity in the caption. A density core, cusp or transition is not already a speed-curve shape: enclosed mass must first be integrated.

Section 3.1 contrasts the NFW family fitted to collisionless simulations with DC14, calibrated on hydrodynamic simulations including baryonic feedback. Feedback involves star formation, explosions and gas flows changing the gravitational potential. Collisionless matter can change its orbit in that changing field; gas winds need not directly collide with it.

DC14 relates profile shape to stellar-to-halo mass ratio. The figure displays model predictions; the changing density profiles belong to the simulations used to calibrate that family.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

Read Figure 5’s inner residuals alongside Tables 1 and 3

Figure 5 compares NFW above and DC14 below for UGC01281, with speed and residual panels for each. Locate the departures at the same radii, then examine how stars, gas and halo share the force. The example shows how the feedback-motivated family alleviates this object’s inner mismatch.

Table 3 extends the comparison to sample medians. At fixed distance/inclination the NFW and DC14 working reduced chi-square medians are 1.3871.387 and 0.9340.934. The varying-geometry columns give 1.4851.485 and 1.1571.157, with parameter counts changing from 3 to 5. These are not just improvements of one unchanged fitted curve: conditions and degrees of freedom also change.

Return to Table 1 for parameter ranges and priors. A prior assigns weight to parameter choices before applying the likelihood. A width in dex for a lognormal mass-to-light distribution differs from a linear standard deviation. Recomputing the comparison requires the complete distribution definitions, width conventions and parameter ranges in the fitting code.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

Figure 6 uses different reference models in its panels

Section 4.1.6 uses BIC=−2ln⁡L+pln⁡N\mathrm{BIC}=-2\ln L+p\ln N, where pp is parameter count and NN the measurements in one galaxy. BIC approximately balances likelihood improvement against complexity under its sample/model conditions. It is not a probability for a dark-matter particle.

The horizontal axis is ΔBIC=BICX−BICref\Delta\mathrm{BIC}=\mathrm{BIC}_X-\mathrm{BIC}_{\rm ref}; the vertical axis is the fraction of galaxies whose comparison lies at or below a given horizontal value. Positive values mean the reference has smaller BIC. The upper reference is DC14. The lower comparisons use the appropriate fixed-acceleration Standard/Simple MOND counterpart. Subtracting curve positions across rows would mix reference models. This reading follows the separate reference for each function stated in section 4.1.6. The lower-right caption of Figure 6 names only Simple MOND; these two descriptions are not identical. Without the implementation, we leave that wording discrepancy unresolved.

Modified Newtonian Dynamics, MOND, changes the low-acceleration relation between motion and baryonic gravity. Here it is the authors’ modified-inertia realization, with functions connecting the high/low-acceleration behaviour and a one-dimensional approximation for environmental fields acting on internal motion. It is not a solution of every modified-gravity field equation. Letting the acceleration scale vary per galaxy is also used as a degeneracy experiment; that freedom should not silently become part of a universal-constant theory.

A small example makes the complexity term tangible. With the same 7373 data points and error coordinate, one extra parameter that reduces χ2\chi^2 by only 22 changes BIC by −2+ln⁡73≃2.29-2+\ln73\simeq2.29. The older model has smaller BIC. A real cross-theory comparison additionally needs common data definitions, defensible priors and model checks. Different physical parameter names alone do not mathematically prohibit statistical comparison.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Independent saved-observation rotation-curve calculation

Figure 10’s megaparsec curve comes from stacked lensing

The horizontal unit now becomes Mpc\mathrm{Mpc}. The inner blue curve averages HI rotation curves in baryonic-mass bins. Outer points come from stacked weak-lensing surface-density signals, deprojected to acceleration and expressed as equivalent circular speed. They are not gas tracers directly followed at those outer radii.

Section 4.3 and the caption describe mass bins, at least three contributors to each HI radial mean, exclusion of the inner two lensing points and the isolation limit beyond one megaparsec. The halo curves also use different mass–concentration prescriptions. Each choice belongs to this particular comparison.

Separate the full-range and outer-only scores. In Table 6 a small full reduced score can accompany a large outer score; numerous inner points can conceal an outskirts mismatch. Figure 14 then checks the fitted halo parameters against the selected scaling relations.

The authors find that their MOND cases sustain the extended flat curves more naturally, and discuss further effects such as material outside bound haloes. We read this fixed version in full, but did not acquire its privately supplied full lensing covariance or rerun its 147-galaxy MCMC/deprojection. Read the reported comparison together with those conditions.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

Worked example: changing baryonic gravity changes a collisionless orbit

A small model explains why feedback can influence dark matter. Start a test particle on a circular orbit about a central mass, then remove one tenth of that mass instantaneously. Its velocity and angular momentum are unchanged at that instant. Its potential becomes shallower, and the orbit becomes more elongated.

For initial central mass MM and circular radius r0r_0, v02=GM/r0v_0^2=GM/r_0. Removing fraction ff changes specific energy to E′=(f−12)GM/r0E'=(f-\tfrac12)GM/r_0. The specific angular momentum is h=r0v0h=r_0v_0. At an apocentre, radial speed is zero, so E′=h2/(2ra2)−G(1−f)M/raE'=h^2/(2r_a^2)-G(1-f)M/r_a. Substitution gives the original turning point r0r_0 and, for 0<f<1/20<f<1/2, the outer one ra=r0/(1−2f)r_a=r_0/(1-2f). Thus f=0.1f=0.1 gives ra=1.25r0r_a=1.25r_0.

The independent material checks three removed fractions by solving the turning-point energy equation numerically. This illustrates an orbit responding to a time-dependent potential. Real feedback also requires the gas-loss/return history, repeated cycles and self-gravity. This point-mass example does not produce a DC14 halo.

Reading: Independent saved-observation rotation-curve calculation · Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

Keep the conditions beside the conclusion

After reading, state whether the evidence concerns individual rotation or stacked lensing; how cusp/core structure enters; whether geometry is fixed; which reference defines each score; and how outer material and covariance are handled. Those details preserve the question the study actually answers.

They also identify useful next measurements. Better spatial resolution probes the inner curve, independent geometry and stellar-mass information narrows compensating choices, and lensing extends the radial scale. Combined evidence can help separate physical explanations that give similar local curves.

Reading: Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Independent saved-observation rotation-curve calculation

Exercises and explained answers

  1. An upper-panel NFW comparison in Figure 6 has ΔBIC=+8\Delta\mathrm{BIC}=+8. Which model does the sign favour, and does it give a MOND probability?

    Read the explanation

    The upper reference is DC14. Positive delta means NFW has larger BIC, favouring DC14 in that specific comparison. Lower panels have separate MOND references. This number cannot be converted into a cross-theory probability.

  2. Why distinguish the outer “speed” in Figure 10 from gas speed in an individual galaxy?

    Read the explanation

    The outer values come from stacked weak lensing, deprojection and a dynamical conversion, involving populations, selection and covariance. They are not directly observed gas orbits at a megaparsec in the same galaxy.

  3. For an instantaneous central-mass loss f=0.2f=0.2, find the apocentre. Why does this not calculate a real galactic core?

    Read the explanation

    ra/r0=1/(1−0.4)=5/3r_a/r_0=1/(1-0.4)=5/3. It describes one test orbit in a prescribed changing central potential. A density profile requires many orbits, self-gravity, gas histories and repeated feedback to evolve together.

Try the prediction

Compare the stellar disk and conditional halo

  1. Inspect the manual cusp and core forward models and describe their inner component curves.
  2. Change the common geometry and observe the coherent residual changes.
  3. Use the Lab as a finite Newtonian mass-model example, then list the additional models, observations and statistical conditions in the paper’s Figures 6 and 10.

Standalone scientific materials

Download the scientific ZIP · Read setup and calculation scope

The offline calculation needs Python, NumPy and SciPy. Default synthetic data are fully generated and fitted. Optional named SPARC files are fetched only on your computer; the ZIP has no original rows or browser Engine.

Individual scripts, sizes and SHA-256
  • README.md (6785 B)
    5fbf0d29e006fc06f35d2290fd8fee7d5abc2301b601814165495331f9c37a84
  • LICENSE.txt (1325 B)
    4f15facc82fd0218a74907ae83f1716569277d9c40f55b249451dd8d89191b2f
  • requirements.txt (185 B)
    359f4c650ca6dd5a2705535e32377aa86db27d8d6ebd9b47333070850b9c5980
  • source-identities.json (724 B)
    409cf4931c19dc97049f943eb005165a59f88d1842a3dd2596eaa1dabcfbd458
  • own-observation.json (9009 B)
    48cd32150a8eff02b4c3e870cb82f8df008d90c2b5aed3f9bca743506852a5e1
  • analyze.py (2750 B)
    5baed281117d3980af0f0e19e428c4163f79a97a84d1db71955ecbee7993ac3a
  • conditional_fit.py (3618 B)
    37700d28edc852a10d09f949a7270a652b61f1c9d30cc2d8d0e99d4347e92d70
  • fetch_named_source.py (5014 B)
    240b29a1c6e929895b5ead31269b495b1223be0615da62c703e2962028e6de5a
  • generate_observation.py (849 B)
    c5442424f9cecdcd81022ec31cae0b81bdc21400f06eea6dde4108c7062de743
  • observation_record.py (2693 B)
    2d791013b5c8b289e6bfdc5ba9ed408e18b9fc4334d44ba32091546ed297e752
  • physical_units.py (255 B)
    2db6448e5f4e657ef48ffd9d2a58c7487c6f7655a8033dc001327fe43192507c
  • plot.py (5734 B)
    c30e8e160f1f73f6451ebb903402a316ecba858814e6b1d5059ad31926abb4e1
  • rotation_model.py (2003 B)
    97c2c5ca56373f20996b26ef179267549ad3a341817b29666ce368029cfa96a1
  • save_native.py (655 B)
    bbf0881e28cf1dd4d13105af656cc4072e3a42a62a5f7013b9f48075b764c082
  • spherical_halos.py (1459 B)
    4678e06bf172167f3a8a3a5d199447658eb51860507405b80110807165c5ef83
  • synthetic_record.py (1538 B)
    6285cda07f2b9324806d25ecc89712eb6b94479e3caaa1f24a88c004d449a063
  • thin_disk_reference.py (1298 B)
    67ba0512e47c10d10c5b8d17f435d808812eec25f0edce296b6c80eaf72e0bc0
  • validate.py (5709 B)
    7cfbf116c5e51e9a3e1aa0ecd113a24f7a069706336ae01221de67dbf331cb6d
  • verify_analytic_limits.py (3073 B)
    68bca45f2d3d8a1f8cda804bb45c3e165025effa0221c2f45d7c1416009bae6a
  • verify_nuisance_mapping.py (4184 B)
    04b8587eb21a257b809fe930c33d43ddb2e008e4ebf2026134ee0a1d35b4015e
  • SHA256SUMS.json (3306 B)
    2fb52c937f2026d872f053e17343e55595fb3568862c20573d59123f3d897b49
  • galaxy-rotation-science-v1.zip (29081 B)
    f389268ed651be3a66ad7d87ee6760ebe39d87f6bd2498b2fb8904c9f10893a3

Sources and calculation scope

  1. Lelli, McGaugh & Schombert (2016), SPARC, AJ152,157

    Official fixed Table1 and Rotmod_LTG.zip; named NGC2403 member only. The master paper sections3.1–3.3/5.1–5.3 and original Table2 notes establish signed gas, unit stellar mass-to-light, finite thickness and error meaning. Fresh source acquisition is separate from cached replay. Original rows and measured-derived plots stay local.

  2. Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

    Full fixed26-page preprint read; Figures1/5/6/10/14 and Table3 visually inspected. Watermark23October2026; abstract says MNRAS Accepted, while the PDF retains placeholder acceptance and2025 style year. Guide follows sections2–6, fixed/varying geometry, internal BIC references and stacked-lensing conditions. No147-object MCMC or privately supplied covariance is reproduced. Original PDF/figures are excluded.

  3. Independent saved-observation rotation-curve calculation

    Own MIT Python produces synthetic observations, actual six-start fits, profiles, commonD/i scenarios, positive-density/ring-potential/Jacobian references and SVGs. Optional exact named-source acquisition and73-row conversion/fits were actually run. Default is offline; ZIP contains no App Engine, original SPARC rows, paper or measured-derived graph.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy