{"id":"63a29528-65f2-46d3-b433-a407412c65fd","arxiv_id":"2608.10781","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Dynamic dark matter substructure in Milky Way-like halos enhances neutron star dark matter accretion by at most a factor of about two, so it cannot explain the twelve-order-of-magnitude gap between equation-of-state and accretion estimates.","lead":"This paper uses a cosmological N-body simulation to track how dark matter density around neutron stars changes as Milky Way-like halos evolve, and finds that substructure boosts accreted dark matter by at most a factor of about two. A smart generalist might read it because it closes off one proposed explanation for why neutron stars seem to hold far more dark matter than simple accretion estimates allow.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative factor-2 claim is a density-exposure ratio, not a capture-rate calculation; Eq. (2.1) assumes velocity-independent linear scaling that is untested for substructure encounters.","rationale":"The reader's weakest_assumption points to precisely the gap I identify: Eq. (2.1) converts time-integrated density into accreted mass with a fixed prefactor that absorbs the velocity distribution, and the paper never verifies that this prefactor is the same for smooth-halo and substructure encounters. My independent reading of Sections 2.3 and 3 confirms that all quoted ratios are density ratios. I therefore agree with the reader's concern. I do not think this warrants changing the CONDITIONAL verdict: the qualitative conclusion that environmental effects cannot explain the orders-of-magnitude discrepancy is robust, because even if velocity corrections changed the factor from 2 to 10 or 0.2, the gap to 10^-2 M_sun remains enormous. The paper's ownFootnote 3 cross-check at 10 h^-1 kpc also supports the robustness of the order-of-magnitude statement. However, the quantitative headline 'mass accretion can be enhanced by a factor ~2' should be presented as a density-exposure result unless the velocity-dependent check is run, which is exactly the conditionality the reader already assigned.","tokens_in":8701,"tokens_out":11411,"duration_ms":122565,"concrete_test":"Using the same Sahyadri snapshots and NS trajectories, recompute the ratios in Fig. 2 with a velocity-dependent capture coefficient instead of Eq. (2.1). For each snapshot, estimate the local DM velocity dispersion from particles in the NS's Voronoi cell or a few-kpc kernel, adopt C ∝ rho_chi / sqrt(sigma_v^2 + v_NS^2) (or the full Gould integral including gravitational focusing), integrate over the NS age, and compare the 95th percentile of M_Vor/M_NFW with the density-only ratio. If the percentile stays within a factor of a few, the linear-density approximation is validated for the paper's purpose; if it shifts by an order of magnitude, the factor-2 claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.3 computes M_acc with Eq. (2.1), a linear, velocity-averaged capture formula. The ratios in Fig. 2 are therefore ratios of time-integrated local DM densities (Voronoi versus evolving NFW), with the same fixed prefactor for every NS and every snapshot; the DM velocity distribution and the NS orbital velocity never enter the calculation. The headline statement that 'mass accretion can be enhanced by a factor ~2 at 95% confidence' is thus strictly a statement about density exposure, not about accreted mass. Standard capture physics (e.g., Gould 1987) has a capture rate C that depends on the relative DM--NS velocity, roughly C ∝ rho_chi / v_eff with gravitational focusing at low velocities. A subhalo encounter with a lower velocity dispersion than the smooth halo could give a mass-enhancement factor larger than the density ratio, while optically thick or saturated capture in a very dense clump would make it smaller. The paper does not test this mapping, so the quantitative factor-2 claim is under-supported. This does not threaten the qualitative conclusion that environmental effects cannot close a roughly 12-order-of-magnitude gap, but it does mean the central quantitative result is not yet established as a statement about DM accretion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper asks whether time-varying local dark-matter densities in the Galactic halo, caused by substructure evolution, can substantially increase the DM mass accreted by a neutron star (NS) relative to the standard smooth-NFW expectation. Using the Sahyadri N-body simulation, the authors select 255 Milky Way analog halos at z=0, follow their most massive progenitors back in time, place a NS at 20 h^-1 kpc from the halo center in either a stationary (Case 1) or circular-orbit (Case 2) configuration, and compute the local DM density at the NS position at every snapshot using both an evolving NFW profile and a Voronoi-tessellation density field. The time-integrated densities are converted to an accreted DM mass via Eq. (2.1), and the ratio M_Vor^acc/M_NFW^acc is presented as the environmental enhancement factor. The median ratios are 1.14 (Case 1) and 0.81 (Case 2), with a long tail toward higher values; the 95th percentile is about 2. The paper concludes that environmental effects enhance accreted mass by at most a factor of about two and therefore cannot explain the orders-of-magnitude discrepancy between equation-of-state-based DM fractions in NSs and smooth-accretion estimates.","tokens_in":8945,"tokens_out":7138,"duration_ms":72110,"significance":"If the result is robust, this is a useful negative result: it shows that dynamical halo substructure, even in a realistic cosmological simulation, does not reconcile the ~10^-14 M_sun accretion-derived DM mass with the ~10^-2 M_sun values inferred from TOV-based analyses of DM-admixed NSs. The paper's use of a high-resolution cosmological simulation and a well-defined Voronoi density estimator is a strength, as is the transparent construction of the ratio, which cancels the common prefactor in Eq. (2.1). The consistency check at 10 h^-1 kpc reported in the footnote strengthens the qualitative conclusion. However, the central quantitative claim of a factor-of-two enhancement is currently presented as a statement about accreted mass, whereas the calculation actually yields a ratio of time-integrated densities under a deliberately linear capture formula; the sensitivity of this mapping to the DM velocity distribution is not tested. This does not threaten the paper's main negative conclusion, which is robust to order-of-magnitude caveats, but it does mean the quantitative factor-of-two claim needs either reframing or additional analysis.","major_comments":[{"comment":"The ratio M_Vor^acc/M_NFW^acc is computed using Eq. (2.1), which assumes that the accreted mass is proportional to the time-averaged DM density times the NS age with a fixed prefactor that is identical for numerator and denominator. Consequently, the quantity plotted in Figure 2 is strictly a ratio of time-integrated local DM densities, not a ratio of actual accreted masses. Real capture rates depend on the DM–NS relative-velocity distribution (e.g., gravitational focusing at low velocities) and can in principle saturate in very dense environments. The paper does not test how such velocity-dependent effects would alter the factor-of-two estimate. The headline statement in the Abstract and Conclusion that 'mass accretion can be enhanced by a factor ~2' is therefore not yet established as a statement about accreted DM mass. The authors should either rephrase the claim as a density-exposure enhancement, or augment the analysis with a simple, physically motivated velocity-dependent capture model for substructure encounters to show that the factor of about two remains a valid proxy for the accreted-mass enhancement.","section":"Section 3, Figure 2"},{"comment":"The quoted factor of about two at the 95th percentile is a point estimate from the distribution of ratios across 255 halos, with no propagated uncertainty. The Voronoi densities themselves are subject to shot noise, because the simulation particle mass is 8.1e7 M_sun/h and the NS is placed at 20 h^-1 kpc, only about six times the force-softening length. A portion of the width of the ratio distribution could therefore be numerical rather than physical. The authors should provide a bootstrap or jackknife error on the 95th percentile and ideally a convergence check varying the Voronoi tracer density or the simulation resolution. Without this, the quantitative claim that the enhancement is 'about a factor 2' at 95% confidence is under-supported, even though the qualitative negative conclusion would likely survive such uncertainties.","section":"Section 2.2 and Abstract"},{"comment":"The Voronoi density field measures the full local density including both subhalo encounters and the aspherical, triaxial shape of the smooth halo, while the baseline is a spherical NFW profile. The ratio shown in Figure 2 is therefore a measure of total deviation from spherical symmetry, not specifically 'the dynamics of substructure' as stated in the Abstract. This conflation does not affect the paper's main conclusion that environmental effects are insufficient to close the gap, but it does affect the physical interpretation. The authors should either soften the substructure-specific language or perform a comparison between the Voronoi densities and the azimuthally averaged density of the simulated halo to isolate the substructure contribution.","section":"Section 3"}],"minor_comments":[{"comment":"The phrase '95% confidence' is a sample percentile of the halo distribution, not a confidence interval from a statistical inference; it should be rephrased as '95th percentile of the distribution' to avoid implying a formal confidence statement.","section":"Section 3"},{"comment":"The sentence 'at most factor ~ 2 enhancement' is too strong, since the 95th percentile leaves a 5% tail of the distribution above that value; a safer phrasing would be 'the 95th percentile of the enhancement is a factor of about two.'","section":"Section 2.3"},{"comment":"The fixed 20 h^-1 kpc radius and the circular-orbit model are simplified treatments of NS trajectories; the footnote checking 10 h^-1 kpc is reassuring, but the paper would benefit from a brief justification of why 20 h^-1 kpc is representative for NS populations, particularly because many NSs receive natal kicks and are born in the disk.","section":"Section 2.3"},{"comment":"The capture formula should cite the original velocity-dependent capture-rate formalism (e.g., Gould 1987) in addition to the review reference [41], so that the limitations of the linear, density-only scaling are clear to the reader.","section":"Section 2.3"},{"comment":"There are several presentation issues: the shaded 68% and 90% regions in Figure 2 are described in the text but not visible in the manuscript as provided; the reference list contains formatting anomalies (e.g., [12] shows 'JCAP 2023 (\"2023\") 073'); and the figure caption for Figure 1 should define the color scale and the meaning of the arrow more explicitly.","section":"Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JCAP and the core negative result is likely robust. The main issue is that the quantitative factor-of-two claim is presented as a statement about accreted mass when the calculation actually yields a density-exposure ratio; this needs either reframing or a velocity-dependent sensitivity test. The lack of error bars on the 95th percentile is a secondary but important issue. The self-citation of the Sahyadri simulation and the Voronoi method is transparent and appropriate, and the authors appear to have performed a reasonable consistency check at 10 h^-1 kpc. With careful revisions, the paper could become a solid contribution; I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a clean negative result that deserves a serious referee. The authors use time-resolved Voronoi densities from the Sahyadri cosmological simulation to ask whether Milky Way substructure can boost the DM mass accreted by a neutron star. They find the boost is at most a factor ~2 relative to an evolving NFW baseline, so environment cannot explain the 12-order-of-magnitude gap between TOV-based DM fractions and simple accretion estimates. I think the qualitative conclusion holds up.\n\nWhat's genuinely new: this is the first time the time-dependent, spatially resolved density field from an MW-analog simulation is used for this question. The comparison to an evolving NFW profile is the right baseline—better than a static halo. The two NS placement cases (stationary and circular orbit) and the 10 h^-1 kpc check cover the reasonable simplifications. The writing is transparent about the pipeline and the Sahyadri method is cited properly.\n\nThe soft spots are real but not fatal. The headline factor ~2 is a ratio of time-integrated DM densities, not a calculation of accreted mass. Eq. (2.1) uses a fixed linear prefactor that cancels in the ratio, so the velocity dependence of capture never enters. Standard capture physics has the rate scaling roughly as rho/v_eff; a dense, cold subhalo encounter could give an enhancement larger than the density ratio, while saturation could make it smaller. The paper doesn't test this mapping. I'd like the authors to either compute the capture integral along the trajectories or state clearly that the result is about density exposure only. That said, even a factor of 10 or 100 in capture efficiency would not close a 12-order-of-magnitude gap, so the central null conclusion is safe.\n\nMinor: no error bars on the factor ~2, and the fixed 20 h^-1 kpc placement is a simplification; the 10 h^-1 kpc check mitigates it. The paper also overstates slightly by calling the ratio 'mass accretion' in the abstract; it's really density-weighted exposure.\n\nWho is it for? Anyone working on DM in compact objects, and it's a useful methodological template for using cosmological simulations in this context. I'd send it to peer review—a referee can push for the capture-physics caveat and a sharper abstract. I'd likely cite it as the reference for environmental enhancement bounds.","headline":"Clean null result: substructure boosts NS dark-matter exposure by at most ~2, but the factor is density-only and capture physics is left untested.","tokens_in":9491,"tokens_out":2302,"would_cite":true,"duration_ms":22846,"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":"This paper shows that the clumpy, evolving dark matter halo of the Milky Way enhances the dark matter mass a neutron star accretes by at most a factor of about two, so environmental effects cannot explain the large gap between…","keywords":["dark matter accretion","neutron stars","Galactic halo substructure","N-body simulations","Voronoi tessellation","Milky Way analogues","NFW profile"],"falsifier":"A calculation that replaces the linear $M_{\\rm acc}\\propto\\langle\\rho_\\chi\\rangle$ scaling in Eq. (2.1) with a velocity-dependent capture formalism applied to the simulated phase-space distribution of dark matter encounters, or an observation of a neutron star whose dark matter fraction exceeds the simulation's 95th-percentile bound for its environment, would falsify the paper's central claim.","tokens_in":8467,"feed_emoji":"⭐","tokens_out":7881,"duration_ms":73259,"temperature":0.7,"pith_summary":"The paper asks whether the time-varying, clumpy dark matter environment of the Milky Way's halo can raise the dark matter mass a neutron star accretes enough to explain the gap between equation-of-state predictions and simple accretion estimates. It uses a high-resolution cosmological N-body simulation of Milky Way-like halos, tracks the local dark matter density around model neutron stars via Voronoi tessellation, and integrates that density over each star's age. The 95th-percentile enhancement over the smooth Navarro-Frenk-White baseline is only a factor of about two for both stationary and orbiting neutron stars. The authors conclude that environmental effects cannot account for the orders-of-magnitude discrepancy between ~$10^{-2}\\,M_\\odot$ and ~$10^{-14}\\,M_\\odot$, so the explanation must lie elsewhere.","feed_headline":"Dark matter substructure can't close the neutron star gap","feed_subtitle":"Even with realistic evolving halo densities, accretion stays within a factor of two of smooth-halo predictions.","key_machinery":"The central object is the time-averaged local dark matter density $\\langle\\rho_\\chi\\rangle = (1/t)\\int_0^t dt'\\,\\rho_\\chi(\\mathbf{x}(t'),t')$ sampled along the neutron star's trajectory. The paper estimates $\\rho_\\chi$ in two ways: from an evolving Navarro-Frenk-White profile whose mass and scale radius come from the simulation's halo catalog, and from a Voronoi tessellation of simulation particles in a 100 $h^{-1}\\,$kpc sub-box, a method that assigns each particle a cell volume and thereby yields the local density field. These densities enter Eq. (2.1), $M_{\\rm acc}\\approx 10^{-14}(\\langle\\rho_\\chi\\rangle/0.3\\,{\\rm GeV\\,cm^{-3}})(\\sigma_{\\chi n}/10^{-45}\\,{\\rm cm^2})(t/{\\rm Gyr})\\,M_\\odot$, so the entire comparison reduces to how the Voronoi density integrated over time differs from the NFW density integrated over time. Two neutron star placements are treated: stationary at 20 $h^{-1}\\,$kpc, and on a circular orbit of that radius.","core_discovery":"The central discovery is that the dynamically evolving, spatially structured dark matter distribution in a Milky Way-like halo does not significantly change the dark matter mass accreted by a neutron star. Measured as the ratio $M_{\\rm acc}^{\\rm Vor}/M_{\\rm acc}^{\\rm NFW}$, the median is 1.14 for a stationary neutron star and 0.81 for one on a circular orbit at 20 $h^{-1}\\,$kpc, with a 95th-percentile value of about 2 in both cases. A footnote reports that moving the stationary star to 10 $h^{-1}\\,$kpc still yields at most a factor of about 2 enhancement. The authors therefore state that environmental effects cannot explain the discrepancy between equation-of-state estimates of ~$10^{-2}\\,M_\\odot$ of dark matter inside a neutron star and accretion-based estimates of ~$10^{-14}\\,M_\\odot$.","pith_inferences":["The paper's linear scaling assumption in Eq. (2.1) is the main lever: if capture efficiency depends on the velocity distribution of dark matter particles in subhalo encounters, a factor of about two in time-integrated density could translate to a different factor in accreted mass; testing this with a phase-space-aware capture calculation is a natural next step.","Because the simulation resolves only halos above roughly $3.2\\times10^9\\,M_\\odot$, the densest small subhalos are absent; a higher-resolution run could produce a longer tail of rare high-density encounters, though the paper's box-size convergence check suggests such a tail would not overturn the main conclusion.","An implication the authors leave implicit is that neutron-star dark matter searches should prioritize mechanisms that convert baryonic matter into dark matter inside the star, or accumulation during earlier stellar phases, rather than the galactic environment."],"forward_implications":["If the factor of about two ceiling holds, then smooth Navarro-Frenk-White accretion estimates are order-of-magnitude reliable, and the missing dark matter mass in neutron stars must come from earlier evolutionary stages or from microphysical channels like neutron-to-dark-matter conversion.","The distribution of $M_{\\rm acc}^{\\rm Vor}/M_{\\rm acc}^{\\rm NFW}$ is skewed above 1 for stationary stars and has a tail toward 2.5 for orbiting stars, so rare substructure encounters do add mass, but not enough to change the overall picture.","Because the conclusion is phrased as a 95th-percentile bound, it provides a quantitative target: any proposed environmental mechanism must produce more than a factor of about two to matter.","The analysis at 10 $h^{-1}\\,$kpc extends the conclusion inward to higher densities, reinforcing that the halo environment is not the decisive factor.","The result redirects attention from the ambient dark matter density to the capture physics and to dark matter accumulation during the main-sequence and supernova phases of the neutron star's progenitor."],"supporting_citations":[{"why":"Defines the Navarro-Frenk-White density profile that serves as the smooth baseline accretion environment.","marker":"[31]"},{"why":"Supplies Eq. (2.1), the linear accretion-mass scaling that converts time-averaged dark matter density into accreted mass.","marker":"[41]"},{"why":"Provides the high-resolution N-body simulation suite and the Voronoi-based density estimation method used throughout the analysis.","marker":"[32]"},{"why":"Supplies the phase-space halo finder used to identify halos and subhalos in the simulation.","marker":"[34]"},{"why":"Provides the merger trees used to trace Milky Way analogs and their progenitors back in time.","marker":"[35]"},{"why":"Calculates dark matter accumulation in neutron stars using local densities and NFW profiles, establishing the comparison the paper builds on.","marker":"[27, 28]"},{"why":"Gives earlier estimates of dark matter accumulation using static NFW densities, defining the discrepancy the paper addresses.","marker":"[14, 25]"},{"why":"Estimates mass accretion from dark matter clumps in a Milky Way simulation, the closest prior treatment of dynamical clumpy environments.","marker":"[26]"}],"fun_headline_variants":["Halo substructure barely boosts neutron star dark matter","Galactic halo dynamics fail to explain neutron star dark matter","Dark matter accretion on neutron stars unaffected by halo clumps","Neutron star dark matter gap remains despite halo structure","Milky Way-like halo substructure won't solve neutron star enigma"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's central estimate rests on the assumption that accreted dark matter mass is directly proportional to the time-averaged local dark matter density, so that a factor-of-two density enhancement translates one-to-one into a factor-of-two mass enhancement, with no saturation or velocity dependence in capture.","fun_headline_variants_meta":{"raw":{"variants":["Halo substructure barely boosts neutron star dark matter","Galactic halo dynamics fail to explain neutron star dark matter","Dark matter accretion on neutron stars unaffected by halo clumps","Neutron star dark matter gap remains despite halo structure","Milky Way-like halo substructure won't solve neutron star enigma"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000518,"raw_usage":{"total_tokens":2522,"prompt_tokens":968,"completion_tokens":1554,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":584,"completion_tokens_details":{"reasoning_tokens":1471}},"tokens_in":584,"tokens_out":1554,"duration_ms":10313,"temperature":1.0,"reasoning_tokens":1471,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:26:44.554231+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation that replaces the linear $M_{\\rm acc}\\propto\\langle\\rho_\\chi\\rangle$ scaling in Eq. (2.1) with a velocity-dependent capture formalism applied to the simulated phase-space distribution of dark matter encounters, or an observation of a neutron star whose dark matter fraction exceeds the simulation's 95th-percentile bound for its environment, would falsify the paper's central claim.","supporting_citations":[{"cited_title":"The Dark Side of Neutron Stars","cited_arxiv_id":"1308.3222","evidence_quote":"Supplies Eq. (2.1), the linear accretion-mass scaling that converts time-averaged dark matter density into accreted mass."},{"cited_title":"Behroozi, R.H","cited_arxiv_id":null,"evidence_quote":"Supplies the phase-space halo finder used to identify halos and subhalos in the simulation."},{"cited_title":"Behroozi, R.H","cited_arxiv_id":null,"evidence_quote":"Provides the merger trees used to trace Milky Way analogs and their progenitors back in time."}],"review_version":1}