{"id":"aa899e19-d3cf-4133-b827-668e05a6bb8b","arxiv_id":"2508.04781","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"Basis-function-expansion N-body models show that the direction, not the amplitude, of the LMC-induced reflex motion of the Milky Way disc can constrain the outer dark matter halo, and reveal a slope-dependent halo instability.","lead":"Sizable N-body simulations of the Milky Way halo being pulled by the infalling Large Magellanic Cloud show that the amplitude of the resulting reflex motion barely depends on the outer halo's density slope, while its direction strongly does. The results tell astronomers which velocity measurements can actually constrain the outer dark matter halo and point to a new oscillatory signature in halo star motions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 50 kpc truncation of the BFE halo may produce the reported beta-sensitivity and 'instability' as numerical artifacts, since it coincides with the LMC's orbital radius; this must be tested before the constraints can be accepted.","rationale":"The reader identified the same load-bearing weakness: the BFE truncated-NFW models are assumed to faithfully represent the real Milky Way halo response, and the abstract offers no convergence or truncation tests. I agree, and I sharpen the concern by pointing out that the truncation radius (50 kpc) coincides with the LMC orbit, making it particularly prone to boundary artifacts. The claim that the instability frequency increases with beta is especially vulnerable, because a truncated halo's local oscillation frequency naturally scales with the enclosed mass profile; without a test at different truncation radii or with an independent code, the 'instability' could be a numerical breathing mode. The concrete test I propose—varying the truncation radius for a fixed beta—would directly determine whether the reported beta-dependencies survive changes in the boundary condition. If they do, the constraints become more credible; if they do not, claims (ii) and (iv) would be unsupported. Until such a test is performed, the paper remains unverified, matching the reader's verdict. I am not claiming the authors are wrong; I am identifying the single most decisive missing check that would settle whether the central claim is physical or an artifact.","tokens_in":13860,"tokens_out":4647,"duration_ms":56684,"concrete_test":"For a single beta value (e.g., beta=1), re-run the same BFE suite with the outer truncation radius moved from 50 kpc to 40 kpc and 80 kpc, holding all other parameters fixed. Compare the dipole direction angle and the period of the reported sinusoidal radial-velocity signal across these truncation radii. If either quantity shifts by an amount comparable to the variation seen across the fiducial beta grid, the claimed beta-sensitivity and instability frequency are dominated by the truncation artifact rather than by the outer halo slope.","verdict_should_be":"UNVERDICTED","load_bearing_attack":"The central claim—that the reflex dipole direction and a halo instability's frequency constrain beta—rests on the fidelity of BFE N-body models with NFW profiles truncated at r=50 kpc and power-law slope beta beyond. The choice of r=50 kpc is suspicious because the LMC's orbital radius is near 50 kpc; the truncated profile has a discontinuous density slope exactly where the perturber acts. The 'halo instability' in claim (iv) may be a breathing mode excited by the initial non-equilibrium from the truncation, and its frequency scaling with beta may simply follow the local epicyclic frequency of the truncated halo rather than a physical dark-matter response. Likewise, the reported sensitivity of the dipole direction to beta could be an artifact of fixing the truncation radius while changing the outer slope: since the mass distribution beyond 50 kpc is altered but the boundary condition at 50 kpc also changes, the comparison across beta mixes the effect of the outer slope with the truncation artifact. No convergence tests, truncation-radius variations, or independent-code comparisons are reported in the abstract, so these possibilities cannot be ruled out.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14099,"tokens_out":4557,"duration_ms":55959,"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":[{"comment":"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.","section":"Abstract (claim iv) / Methods: truncated NFW model"},{"comment":"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.","section":"Abstract (claim ii) and BFE-coefficient statements"},{"comment":"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.","section":"Full text / equations / figures"}],"minor_comments":[{"comment":"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.","section":"Abstract"},{"comment":"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.","section":"Abstract (claim iii)"},{"comment":"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.","section":"Abstract, final sentence"},{"comment":"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.","section":"Metadata"}],"recommendation":"major_revision","confidential_remarks":"I would ask the editor to obtain a clean, decodable version of the manuscript before sending it to a second referee; the supplied text is not readable. The scientific idea is interesting and the abstract states clear, falsifiable predictions, but the central numerical claims are unverifiable without the truncation/convergence tests requested in the major comments. If the clean version is identical to the arXiv PDF, the authors need to add those tests before this can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague —\n\nThe full text of this submission is mojibake; all I can assess is the abstract. On that basis, this is a serious and useful piece of work. The new content is the explicit mapping: the reflex dipole amplitude is largely insensitive to the outer halo slope beta, while the dipole direction, the induced halo contraction, and the frequency of a claimed halo instability are sensitive to beta. The instability's observable signature—a sinusoidal pattern in the mean radial velocity of halo stars—is genuinely interesting, if it holds up.\n\nCredit where due: the design is a parameter scan, not a fit to target, so the reported trends are emergent from the N-body integration rather than baked in by construction. The negative result on the dipole amplitude is valuable; it tells observers not to use that amplitude to constrain the outer profile.\n\nThe soft spots are numerical. The truncation at 50 kpc sits right at the LMC's orbital radius. The density-slope discontinuity there could excite breathing modes that masquerade as the claimed instability, and the beta-dependent frequency might simply track the local epicyclic frequency of the truncated model rather than a physical dark-matter response. The abstract reports no convergence tests, no truncation-radius variation, and no independent-code check. Those omissions are not damning by themselves, but they keep me from endorsing the instability claim without more.\n\nThere is also a production defect: the arXiv header says [cs.LG] while the paper is in astro-ph.GA. I attribute that to the same corruption, but it needs fixing.\n\nWho is this for? Anyone working on LMC-induced perturbations of the Milky Way halo. It deserves a serious referee if the equations and figures can actually be read. Right now the only honest path is to ask the authors for a clean manuscript before sending it out. If the corruption persists, desk-reject on technical grounds, not scientific merit.\n\nMy vote: conditional accept for peer review, pending a readable version.","headline":"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.","tokens_in":14621,"tokens_out":2626,"would_cite":false,"duration_ms":28850,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Large Magellanic Cloud","Milky Way dark matter halo","reflex motion","velocity dipole","N-body simulations","basis function expansion","dark halo instability","outer halo density profile"],"falsifier":"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.","tokens_in":13684,"feed_emoji":"🌌","tokens_out":9420,"duration_ms":97740,"temperature":0.7,"pith_summary":"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.","feed_headline":"Direction, not amplitude, of LMC reflex motion reveals halo slope","feed_subtitle":"Simulations find the outer halo slope controls reflex-motion direction and a halo oscillation; dipole amplitude barely moves","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[],"fun_headline_variants":["Reflex motion direction, not size, reveals Milky Way halo shape","LMC reflex direction probes dark halo slope, amplitude does not","New halo instability offers sine-wave test of outer dark matter","Steeper Milky Way halo: weaker dipole, stronger quadrupole","Amplitude blind to halo slope; direction and oscillation see it"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Reflex motion direction, not size, reveals Milky Way halo shape","LMC reflex direction probes dark halo slope, amplitude does not","New halo instability offers sine-wave test of outer dark matter","Steeper Milky Way halo: weaker dipole, stronger quadrupole","Amplitude blind to halo slope; direction and oscillation see it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1271,"prompt_tokens":823,"completion_tokens":448,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":360}},"tokens_in":567,"tokens_out":448,"duration_ms":5417,"temperature":1.0,"reasoning_tokens":360,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T23:47:37.607910+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}