Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Constraining the Milky Way dark matter halo with LMC-induced reflex motion

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The outer slope of the Milky Way's dark matter halo shows up in the direction, contraction, and dipole/quadrupole balance of the LMC-induced reflex motion, while the dipole amplitude stays nearly unchanged.

desk verdict The abstract makes a plausible, useful claim about which LMC-reflex-motion observables constrain the outer halo, but the corrupt full text makes the numerics unverifiable; referee it only after a clean copy arrives. read the letter →

arxiv 2508.04781 v1 pith:4NDK6FAI submitted 2025-08-06 astro-ph.GA

Rashid Yaaqib , Michael Petersen , Jorge Peñarrubia This is my paper
classification astro-ph.GA
keywords LargeMagellanicCloudMilkyWaydarkmatterhaloreflexmotionvelocitydipoleN-bodysimulationsbasisfunctionexpansioninstabilityouterdensityprofile
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the Large Magellanic Cloud's gravitational pull moves the Milky Way's disc and halo relative to each other, and that this reflex motion can be used to measure the density profile of the Galaxy's outer dark matter halo. Using N-body simulations built on a basis-function expansion with halo profiles whose density falls as $\rho \propto r^{-\beta}$ beyond 50 kpc, the authors find that the amplitude of the reflex dipole is nearly independent of $\beta$, so it cannot by itself set the Milky Way–LMC mass ratio or the outer slope. The direction of the disc's motion, the halo contraction, and the balance between dipole and quadrupole distortions are all sensitive to $\beta$, giving alternative observables. The paper also reports a halo instability whose oscillation frequency grows with $\beta$ and would show up as a sinusoidal pattern in the mean radial velocity of halo stars. If right, future measurements should target these directional and oscillatory signatures rather than the dipole strength.

What carries the argument

The key machinery is a suite of basis-function expansion (BFE) N-body simulations, which represent the Milky Way halo as a family of truncated Navarro–Frenk–White profiles with $\rho \propto r^{-\beta}$ beyond $r = 50$ kpc. The BFE coefficients decompose the evolving gravitational potential into multipole moments, so the LMC-induced distortion is tracked as dipole and quadrupole terms; the direction and amplitude of the reflex motion, the contraction of the halo, and the oscillation frequency of the halo instability are all read off these coefficients. This setup lets the authors separate the $\beta$-dependent directional and quadrupole signals from the $\beta$-insensitive dipole amplitude.

What would settle it

Run the same simulation suite with a different truncation radius (e.g., 100 kpc instead of 50 kpc) and with an independent N-body method: if the instability frequency, the dipole-direction sensitivity, or the dipole/quadrupole balance changes materially, the central claim fails. Observationally, a long time series of mean radial velocities of halo stars should show whether the predicted sinusoidal pattern exists and whether its period tracks the outer-halo slope.

Watch

Extended reading notes

Core claim

The central claim is that the outer slope of the Milky Way's dark halo leaves a clear but selective imprint on the LMC-induced reflex motion. In the authors' BFE N-body models, the amplitude of the velocity dipole induced on halo-star kinematics varies with Galactocentric radius but barely changes as the outer halo slope $\beta$ is varied; the direction of the disc's reflex motion, by contrast, is strongly $\beta$-dependent, as is the contraction of the halo response. The same simulations show a halo instability whose oscillation frequency increases with $\beta$, which would appear as a sinusoidal variation of the mean radial velocity of halo stars, and steeper outer profiles suppress the di

Load-bearing premise

The load-bearing premise is that the simulated halos—a family with density $\rho \propto r^{-\beta}$ only beyond 50 kpc, evolved with a basis-function expansion—respond to the LMC the way the real Milky Way halo does, so the reported $\beta$-sensitivities and the instability are physical rather than artifacts of the truncation radius or the expansion method.

Editorial extensions

If this is right

  • The reflex dipole amplitude should not be used alone to infer the outer halo slope or the Milky Way–LMC mass ratio, because it is nearly insensitive to $\beta$.
  • Measuring the direction of the disc's reflex motion, for example through proper motions of disc and halo tracers, could constrain the outer dark matter density slope.
  • A sinusoidal pattern in the mean radial velocity of halo stars would be a direct, potentially observable signature of the predicted halo instability, with its frequency set by $\beta$.
  • Steeper outer halos should show smaller dipole distortions but larger quadrupole distortions, so combining dipole and quadrupole measurements improves the constraint on the outer profile.
  • The LMC-induced contraction of the Milky Way halo depends strongly on the outer profile, so modeling that contraction is necessary when interpreting the reflex-motion signal.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the instability is physical, its frequency could act as a new dynamical clock for the outer halo: measuring the period of the sinusoidal radial-velocity pattern might pin down $\beta$ more precisely than the static dipole field could.
  • A direct numerical extension would be to rerun the same experiments with different truncation radii; if the instability frequency or the dipole-direction sensitivity shifts with truncation, the signal is an artifact of the 50 kpc cut rather than a real halo property.
  • The dipole/quadrupole trade-off suggests that future wide-area astrometric surveys should treat the reflex motion as a vector field, including direction and higher multipoles, rather than as a scalar amplitude, to maximize dark-matter constraining power.
  • The result implies that Milky Way–LMC mass estimates built on reflex-motion amplitude may have underestimated their uncertainty, because that amplitude is degenerate with the outer halo profile.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper uses basis function expansion (BFE) N-body simulations of the Milky Way halo with truncated NFW profiles (rho proportional to r^-beta beyond r=50 kpc) to study the LMC-induced reflex motion of the MW disc. The abstract reports four findings: (i) the reflex motion amplitude is largely insensitive to the outer DM slope beta; (ii) the direction of the disc motion is very sensitive to beta; (iii) the contraction of the MW halo depends strongly on beta; and (iv) a halo instability exists whose oscillation frequency increases with beta, producing a potentially observable sinusoidal pattern in the mean radial velocity of halo stars. The paper concludes that the reflex amplitude alone has limited constraining power for the outer MW profile and that the direction, dipole/quadrupole balance, and instability frequency are more promising diagnostics. The supplied body text is largely corrupted mojibake, so detailed verification of the methods, figures, and numerical results was not possible; this report is therefore based primarily on the readable abstract and the identifiable model description.

Significance. If the results hold, the paper would establish a useful new probe of the outer Milky Way dark matter halo: the direction of the LMC-induced reflex motion and the frequency of a halo oscillation could constrain the outer density slope beta, while the reflex amplitude would be shown to be a poor diagnostic. The predicted beta-dependent sinusoidal signal in halo-star radial velocities is falsifiable with current or near-future surveys, and the qualitative separation of amplitude insensitivity versus direction sensitivity is a clear, testable claim. The BFE methodology is standard in the field and the parameter scan over beta is a sensible design. However, the significance is currently contingent: the central claims rest on simulations with a fixed truncation radius at 50 kpc, and the manuscript as supplied provides no convergence tests, no truncation-radius variation, and no independent-code comparison. Until these are shown, the reported beta-sensitivities and the instability cannot be distinguished from numerical artifacts associated with the imposed density break.

major comments (3)
  1. [Abstract (claim iv) / Methods: truncated NFW model] The claimed halo instability is central to claim (iv), but the manuscript gives no evidence that it is a physical collective mode rather than a breathing/transient response to the imposed density break at r=50 kpc. The initial conditions must be in Jeans equilibrium for the truncated profile; if the outer power law is not initialized self-consistently, radial oscillations will have frequencies that scale with the local epicyclic frequency, which itself depends on beta. Please add (a) a test that the initial halo is in equilibrium in the absence of the LMC, (b) a variation of the truncation radius (e.g., 40 and 60 kpc) to show that the mode frequency and the dipole direction are stable, and (c) a comparison of the BFE result with a direct N-body or linear-response calculation for at least one beta.
  2. [Abstract (claim ii) and BFE-coefficient statements] The reported sensitivity of the reflex direction to beta is obtained while holding the truncation radius fixed at 50 kpc. Because the LMC's orbital radius is near 50 kpc, fixing the break radius while varying beta changes both the mass distribution beyond 50 kpc and the density slope at the perturbing satellite's location. These two effects are degenerate. The abstract reports no test with different truncation radii or with a smooth, non-broken outer profile. This is load-bearing because the main conclusion is that the reflex direction constrains beta; without separating the beta effect from the boundary effect, the claim is not established.
  3. [Full text / equations / figures] As supplied, the body text, equations, and figure captions are corrupted mojibake; I cannot verify the simulation parameters (particle number, force softening, basis order, time integration) or any quantitative results. This blocks a normal technical review. The authors should provide a readable manuscript and a reproducibility statement (code/data availability) so that the numerical claims can be checked directly.
minor comments (4)
  1. [Abstract] The range of beta explored and the assumed MW-LMC mass ratio should be stated explicitly; claim (i) says the amplitude is insensitive to beta, but the reader needs to know the parameter range over which this null result is asserted.
  2. [Abstract (claim iii)] The phrase 'contraction of the MW halo' should be defined quantitatively (e.g., change in monopole coefficient, radial density profile, or potential depth) and distinguished from the instantaneous reflex motion of the disc.
  3. [Abstract, final sentence] The statement that 'steeper truncations produce smaller dipole distortions, while amplifying the quadrupole distortion' is qualitative. Specify the metric (e.g., normalized BFE coefficients), the radial range, and the uncertainty on the distortion amplitudes.
  4. [Metadata] The arXiv header displayed in the supplied text reads 'cs.LG' while the paper is submitted as astro-ph.GA; verify the correct archive category. Also, ensure all equations and special characters are encoded so the text is readable.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: beta is a scanned input parameter and the dipole/quadrupole response, direction sensitivity, and instability frequency are emergent N-body outputs, not re-expressions of the input profile.

full rationale

The paper's load-bearing claims (i)-(iv) are derived from N-body simulations in which the outer halo slope beta is a varied input parameter, not a fitted quantity. The reflex-motion amplitude, the disc-motion direction, the halo contraction, and the claimed sinusoidal mean radial velocity signature are all measured from the time-evolved phase-space distribution after the LMC perturber is introduced. Nothing in the abstract or visible text indicates that any of these outputs is defined in terms of beta or is recovered from the input profile by construction. The statement 'using BFE coefficients we find that steeper truncations produce smaller dipole distortions, while amplifying the quadrupole distortion' describes a diagnostic decomposition of the simulated response; BFE coefficients are a basis expansion tool, not a fitting target. Concerns that the r=50 kpc truncation or the authors' own BFE machinery could introduce numerical artifacts are potential correctness or robustness issues, not circularity: they do not show that a prediction is equivalent to its inputs by construction. Self-citation of the BFE method is normal and is not load-bearing in the sense of importing an unverified uniqueness theorem or ansatz; the core results rest on the integration itself. Therefore the derivation chain is self-contained with respect to the paper's stated inputs, and the circularity score is 0.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The central claims rest on the simulation setup: the halo model family (single power-law beyond 50 kpc), the BFE truncation, and the mapping from simulated particles to observable halo-star kinematics. These are all inputs or domain assumptions rather than derived results, and none could be audited because the methods text is unreadable. No parameters appear to be fitted to the target signal; beta is scanned as an input. No new physical entities are introduced; the reported halo instability is a dynamical phenomenon, not an entity.

free parameters (3)
  • beta, outer halo power-law slope = input, varied across runs (values not given in abstract)
    Central claims are sensitivity statements with respect to beta; it is a design choice, not fitted, but the conclusions' scope is defined by it.
  • truncation radius / power-law start radius = 50 kpc
    All halos are truncated NFW with the beta power-law beyond r = 50 kpc; no robustness test against this chosen radius is visible in the abstract.
  • MW-LMC mass ratio = not stated in abstract
    The paper argues this ratio alone does not set dipole amplitude; the adopted ratio is an input that scales all simulated signals.
assumptions (3)
  • domain assumption BFE simulations with the adopted expansion accurately capture the self-gravitating response of the MW halo to LMC infall
    The whole study rests on basis function expansion N-body models; if basis truncation or the 50 kpc boundary dominates the dipole response, findings (i)-(v) would be artifacts. Methods section unreadable so this cannot be audited.
  • ad hoc to paper Truncated NFW halos with a single power law beyond 50 kpc span the relevant outer-halo profile space
    Beta is the only outer-halo shape freedom; realistic MW outer halos could have more complex shapes (flattening, baryonic feedback, streams), excluded by construction.
  • domain assumption Simulated particles map onto observable halo-star kinematics
    The observable claims (velocity dipole of halo stars, sinusoidal mean radial velocity) require the simulation tracer population to match the real stellar halo, which is not specified in the abstract.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Constraining the Milky Way dark matter halo with LMC-induced reflex motion." pith.science (2026). https://pith.science/paper/4NDK6FAI

@misc{pith2026250804781,
  author       = {Pith},
  title        = {Pith review of: Constraining the Milky Way dark matter halo with LMC-induced reflex motion},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4NDK6FAI}},
  note         = {Machine review of arXiv:2508.04781}
}
abstract

Modelling perturbations of the Milky Way (MW) halo induced by the infall of the Large Magellanic Cloud (LMC) offers new avenues to constrain the dark matter (DM) distribution in our Galaxy. A key observable is the reflex motion of the MW disc with respect to the halo induced by the LMC's infall, which imprints a velocity dipole on kinematics of halo stars. Here we investigate how the dipole varies with Galactocentric radius, and study the sensitivity of the reflex motion signal to different DM outer-halo profiles. Using a suite of basis function expansion (BFE) simulations with truncated NFW profiles ($\rho \propto r^{-\beta}$ beyond $r=50$ kpc), our $N$-body models show that (i) The reflex motion amplitude varies with Galactocentric radius but is largely insensitive to the outer DM slope, implying that the MW-LMC mass ratio alone does not set the dipole strength. (ii) In contrast, the direction of the disc motion is very sensitive to the density distribution of the outer DM halo. (iii) The contraction of the MW halo induced by the LMC's gravitational pull also depends strongly on the outer DM halo profile. (iv) We find a halo instability whose oscillation frequency increases with $\beta$ producing a potentially observable signature - a sinusoidal pattern of the mean radial velocity of halo stars. Finally, using BFE coefficients we find that steeper truncations produce smaller dipole distortions, while amplifying the quadrupole distortion. These results highlight the limited constraining power of the reflex motion amplitude alone for outer MW profile parameters.

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. The Milky Way - Large Magellanic Cloud Interaction with Simulation Based Inference

    astro-ph.GA 2025-10 conditional novelty 5.0 of 10

    Simulation-based inference on outer-halo star velocities gives a Milky Way reflex speed of 26.4 km/s and an LMC enclosed mass of 9.2×10^10 solar masses within 50 kpc.

Reference graph

Works this paper leans on

30 extracted references · 2 canonical work pages · cited by 1 Pith paper

  1. [1]

    Bland-Hawthorn J., Gerhard O., 2016, @doi [Annual Review of Astronomy and Astrophysics] 10.1146/annurev-astro-081915-023441 , 54, 529

    doi:10.1146/annurev-astro-081915-023441 2016

  2. [2]

    Brooks R. A. N., Sanders J. L., Dillamore A. M., Garavito-Camargo N., Price-Whelan A. M., 2025, @doi [arXiv e-prints] 10.48550/arXiv.2507.10667 , https://ui.adsabs.harvard.edu/abs/2025arXiv250710667B p. arXiv:2507.10667

    doi:10.48550/arxiv.2507.10667 2025

  3. [3]

    arXiv:2410.09149

    Bystr \"o m A., et al., 2024, @doi [arXiv e-prints] 10.48550/arXiv.2410.09149 , p. arXiv:2410.09149

    arXiv doi:10.48550/arxiv.2410.09149 2024

  4. [4]

    Chandra V., et al., 2025, @doi [The Astrophysical Journal] 10.3847/1538-4357/addab6 , 988, 156

    doi:10.3847/1538-4357/addab6 2025

  5. [5]

    B., 2021, Jupyter Book , Zenodo, @doi 10.5281/zenodo.4539666

    Community E. B., 2021, Jupyter Book , Zenodo, @doi 10.5281/zenodo.4539666

    doi:10.5281/zenodo.4539666 2021

  6. [6]

    Correa Magnus L., Vasiliev E., 2021, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stab3726 , 511, 2610

    doi:10.1093/mnras/stab3726 2021

  7. [7]

    The Dipole Instability in Gravitational $N$-body Systems: A Natural Explanation for Lopsidedness and Off-Centered Nuclei in Galaxies

    Dattathri S., van den Bosch F. C., Weinberg M. D., Banik U., 2025, @doi [arXiv e-prints] 10.48550/arXiv.2505.23905 , p. arXiv:2505.23905

    work page Pith review arXiv doi:10.48550/arxiv.2505.23905 2025

  8. [8]

    Erkal D., et al., 2019, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stz1371 , 487, 2685

    doi:10.1093/mnras/stz1371 2019

Show all 30 references
  1. [9]

    Erkal D., et al., 2021, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stab1828 , 506, 2677

    2021 doi

  2. [10]

    Garavito-Camargo N., Besla G., Laporte C. F. P., Johnston K. V., G \'o mez F. A., Watkins L. L., 2019, @doi [The Astrophysical Journal] 10.3847/1538-4357/ab32eb , 884, 51

    2019 doi

  3. [11]

    A., Besla G., Carpintero D

    G \'o mez F. A., Besla G., Carpintero D. D., Villalobos \'A ., O'Shea B. W., Bell E. F., 2015, @doi [The Astrophysical Journal] 10.1088/0004-637x/802/2/128 , 802, 128

    2015 doi

  4. [12]

    R., et al., 2020, @doi [Nature] 10.1038/s41586-020-2649-2 , 585, 357

    Harris C. R., et al., 2020, @doi [Nature] 10.1038/s41586-020-2649-2 , 585, 357

    2020 doi

  5. [13]

    D., 2007, @doi [Computing in Science & Engineering] 10.1109/MCSE.2007.55 , 9, 90

    Hunter J. D., 2007, @doi [Computing in Science & Engineering] 10.1109/MCSE.2007.55 , 9, 90

    2007 doi

  6. [14]

    Iorio G., Belokurov V., 2021, @doi [ ] 10.1093/mnras/stab005 , https://ui.adsabs.harvard.edu/abs/2021MNRAS.502.5686I 502, 5686

    2021 doi

  7. [15]

    E., et al., 2023, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad551 , 521, 4936

    Koposov S. E., et al., 2023, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad551 , 521, 4936

    2023 doi

  8. [16]

    Lilleengen S., et al., 2022, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stac3108 , 518, 774

    2022 doi

  9. [17]

    J., 2016, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stw2759 , 465, 76

    McMillan P. J., 2016, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stw2759 , 465, 76

    2016 doi

  10. [18]

    A., Besla G., Erkal D., Ma Y.-Z., 2015, @doi [Monthly Notices of the Royal Astronomical Society: Letters] 10.1093/mnrasl/slv160 , 456, L54

    Pe \ n arrubia J., G \'o mez F. A., Besla G., Erkal D., Ma Y.-Z., 2015, @doi [Monthly Notices of the Royal Astronomical Society: Letters] 10.1093/mnrasl/slv160 , 456, L54

    2015 doi

  11. [19]

    E., 2007, @doi [Computing in Science and Engineering] 10.1109/MCSE.2007.53 , 9, 21

    P \'e rez F., Granger B. E., 2007, @doi [Computing in Science and Engineering] 10.1109/MCSE.2007.53 , 9, 21

    2007 doi

  12. [20]

    S., Pe \ n arrubia J., 2020, @doi [Monthly Notices of the Royal Astronomical Society: Letters] 10.1093/mnrasl/slaa029 , 494, L11

    Petersen M. S., Pe \ n arrubia J., 2020, @doi [Monthly Notices of the Royal Astronomical Society: Letters] 10.1093/mnrasl/slaa029 , 494, L11

    2020 doi

  13. [21]

    S., Pe \ n arrubia J., 2021, @doi [Nature Astronomy] 10.1038/s41550-020-01254-3 , 5, 251

    Petersen M. S., Pe \ n arrubia J., 2021, @doi [Nature Astronomy] 10.1038/s41550-020-01254-3 , 5, 251

    2021 doi

  14. [22]

    S., Weinberg M

    Petersen M. S., Weinberg M. D., 2025, @doi [Journal of Open Source Software] 10.21105/joss.07302 , 10, 7302

    2025 doi

  15. [23]

    Shipp N., et al., 2021, @doi [The Astrophysical Journal] 10.3847/1538-4357/ac2e93 , 923, 149

    2021 doi

  16. [24]

    Vasiliev E., 2023a, @doi [Galaxies] 10.3390/galaxies11020059 , 11

    doi

  17. [25]

    Vasiliev E., 2023b, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad2612 , 527, 437

    doi

  18. [26]

    Vasiliev E., Belokurov V., Erkal D., 2020, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/staa3673 , 501, 2279

    2020 doi

  19. [27]

    Virtanen P., et al., 2020, @doi [Nature Methods] 10.1038/s41592-019-0686-2 , 17, 261

    2020 doi

  20. [28]

    D., 2023, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad2591 , 525, 4962

    Weinberg M. D., 2023, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad2591 , 525, 4962

    2023 doi

  21. [29]

    S., Pe \ n arrubia J., 2024, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stae1363 , 531, 3524

    Yaaqib R., Petersen M. S., Pe \ n arrubia J., 2024, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stae1363 , 531, 3524

    2024 doi

  22. [30]

    Zhao H., 1996, @doi [ ] 10.1093/mnras/278.2.488 , https://ui.adsabs.harvard.edu/abs/1996MNRAS.278..488Z 278, 488

    1996 doi

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.