{"id":"9704d25e-9e8f-4ad3-905c-79ab377fd0d3","arxiv_id":"2504.12393","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A new Monte Carlo scheme embeds dark matter-baryon scattering into SPH and MFM simulations and reproduces analytic heat and momentum transfer, including a 1:1000 physical mass ratio.","lead":"This paper presents the first numerical scheme for simulating dark matter particles scattering off ordinary gas particles inside cosmological N-body and hydrodynamics simulations. The method passes analytic test problems and opens the way to new astrophysical constraints on dark matter-baryon interactions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The r=1000 validation in §3.3 appears to violate the scheme's own mass-ratio constraint Eq. (12), so the headline mass-ratio capability may rest on an invalid test configuration.","rationale":"Good-faith reading: the paper is a numerical-methods paper whose central claim is that this is the first scheme that simulates DM-baryon scattering with correct angular and velocity dependence and mass ratio. The scheme is described in detail, validated against analytic heat-conduction and momentum-transfer solutions, and the code is built on an existing SIDM implementation; the SPH tests agree with the analytic solutions, and the r=1000 test is specifically the evidence for the mass-ratio claim. The weakest point is not only the thermalization assumption already flagged by the reader, although that is real and acknowledged, but the internal consistency of the r=1000 validation: the stated particle numbers imply mbary/mDM≈2.14, which is incompatible with the derivation of the virtual-particle construction (Eqs. (10)–(12)) for r=1000. If this is a typo in the paper, it is a serious one because the test is the sole demonstration of the mass-ratio capability; if it is not a typo, the simulation must be relying on rejection and artificial heat conduction, which the paper itself says can underestimate interactions by preferentially discarding high-velocity scatterings. The proposed check settles which case holds. The reader's conditional verdict is reasonable; the additional condition should be a corrected and transparent r=1000 benchmark (or an explicit exemption for the frequent-scattering branch), plus reporting of rejection rates. The velocity-dependent cross-section claim is also untested, but that is secondary to the mass-ratio inconsistency because all present tests use velocity-independent cross-sections.","tokens_in":33792,"tokens_out":20508,"duration_ms":213601,"concrete_test":"Rerun the §3.3 unequal-mass heat-conduction test with numerical masses that satisfy Eq. (12): for r=1000 and ζ=5 take NDM=10^8 and Nbary=10^5 (mbary/mDM≈10^4), keeping σT/mχ=1000 cm²/g and the same physical setup, and compare the DM kinetic-energy curve with Eqs. (19)–(20). If the result shifts by more than the few-percent agreement shown in Fig. 11, the published r=1000 validation is not a clean test of unequal-mass scattering. In addition, report the fraction of rejected scatterings in the published setup; if it is not negligible, the analytic match is not evidence that the physical scattering is being simulated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim includes accurate treatment of the DM/baryon mass ratio, and §2.1 requires the virtual-particle mass to be mvirt = r mDM with the numerical mass ratio matching the physical one. Section 2.2 then derives Eq. (12): for ζ=5, positive internal energy after a large-angle event requires mbary/mDM > r(ζ²/3+1) ≈ 9.33r. The r=1000 heat-conduction test in §3.3 uses Nbary=46656, NDM=10^5, and equal total masses (10^10 M☉ each), giving mbary/mDM ≈ 2.14 and mvirt = 1000 mDM ≈ 467 mbary. This is a factor of roughly 4400 below the Eq. (12) bound and, more fundamentally, makes Eq. (11) have no real solution for |vrand|, so the virtual particle cannot be created from the baryonic particle under the paper's own energy bookkeeping. The paper's fallback (rejecting scatterings that produce non-positive internal energy and relying on artificial heat conduction, §2.2, §2.3, §3.3) removes preferentially the high-relative-velocity events that carry the largest energy and momentum transfer. Thus the agreement in Fig. 11 may reflect the rejection/heat-conduction correction rather than the physical unequal-mass scattering kernel, leaving the mass-ratio part of the central claim unsupported. If the small-angle/frequent branch is meant to be exempt from the mass-splitting requirement, that exemption is not stated for the creation step, which §2.1 describes for every pair.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a numerical scheme for simulating dark matter-baryon (DM-baryon) scattering in N-body and hydrodynamics codes. The method creates a 'virtual' particle from each baryonic numerical particle, draws its velocity from the baryon's Maxwell-Boltzmann distribution, and scatters it with a DM particle using an SIDM scattering routine; after scattering, the virtual particle is destroyed and its momentum and energy changes are folded back into the baryonic particle. The authors validate the scheme against analytic solutions for heat conduction in both directions, momentum transfer, unequal-mass scattering with r=1000, isotropic scattering, and a comoving test, and then apply it to the collapse of an overdensity to study halo formation. They report good agreement for the SPH implementation, discuss the role of viscosity and heat conduction, and identify the negative-internal-energy problem and energy non-conservation in MFM as key limitations.","tokens_in":34065,"tokens_out":12721,"duration_ms":120878,"significance":"The manuscript presents a conceptually new Monte Carlo coupling of DM-baryon scattering to SPH/MFM codes via a virtual-particle construction, with a physically motivated probability and drag-force formulation. The heat-conduction tests (Figs. 2-4) and the equal-mass small-angle tests are checked against standard analytic solutions (Dvorkin et al. 2014; Muñoz et al. 2015) and reproduce them well, especially with SPH. The paper is commendably transparent about the thermalisation assumption, the role of viscosity and heat conduction, the negative-internal-energy problem, and the energy non-conservation of the MFM implementation. The halo-formation application, while idealized, demonstrates a concrete use case and yields interesting non-monotonic baryon-density behavior with cross-section. If the mass-ratio capability can be established in a valid parameter regime, the scheme would be a valuable tool for astrophysical DM probes. However, the current r=1000 validation does not respect the scheme's own constraints, so the headline claim of accurate mass-ratio handling is not yet supported.","major_comments":[{"comment":"The r=1000 test configuration violates the scheme's own mass-ratio constraint. With mbary/mDM about 2.14 and mvirt = 1000 mDM about 467 mbary, the right-hand side of Eq. (11) is imaginary, so no virtual-particle velocity satisfies the stated positive-internal-energy condition; the derivation of Eq. (11) assumes mbary > mvirt, and for mvirt > mbary the 'remaining baryonic particle' in Eq. (10) has negative mass, making the energy bookkeeping unphysical. The statement that Eq. (12) applies only to large-angle scattering does not address this, because Eq. (11) is derived from the creation step, which is angle-independent. The test therefore relies on the rejection scheme and artificial heat conduction described in Section 2.2 and Section 3.3, which preferentially remove the highest-velocity scatterings. Consequently, the agreement in Fig. 11 does not demonstrate that the unequal-mass scattering kernel is accurately modeled, and the central claim of accurately accounting for the mass ratio is not supported. Please either run the test in a parameter regime satisfying Eq. (12) (or a stated small-angle analog), or quantify the rejected fraction and demonstrate that the result is insensitive to the rejection and heat-conduction corrections.","section":"Section 3.3, Eqs. (11)-(12)"},{"comment":"The analytic reference for the momentum-transfer test is an interpolation between two asymptotic regimes using a logistic weighting with hand-chosen parameters a=3.0 and b=0.1. The function f(X) = 1/(1 - exp(a(1-bX))) can take negative values for X < 1/b (for example, X=2 gives f about -0.1), which is unphysical because it reverses the sign of the drag contribution. Since the relative velocity in the simulation decreases toward about 1 km/s, corresponding to X about 2, the agreement in Fig. 7 is only as reliable as this interpolation. Please replace the reference with a full numerical solution or an exact solution, or justify the interpolation and demonstrate that the negative-f region has negligible effect on the comparison.","section":"Section 3.2, Eq. (25)"}],"minor_comments":[{"comment":"The claim that the scheme 'accurately accounts for' the mass ratio is stronger than what the current r=1000 test supports; please soften this claim or provide a valid unequal-mass test.","section":"Abstract and Section 6"},{"comment":"The word 'asymptomatically' should be 'asymptotically'.","section":"Section 3.3"},{"comment":"The typeset expression for xi' appears garbled in the text; please ensure the cube root is written unambiguously as xi' = (1/(2h)) * (Nidm/nmax)^(1/3).","section":"Section 2.4.1, Eq. (14)"},{"comment":"The legend labels 'Heat: DM Gas' and 'Heat: Gas DM' are ambiguous; please use arrows (e.g., 'Heat: DM -> Gas') to indicate the direction of heat flow.","section":"Figures 2 and 3"},{"comment":"The requirement that the thermalisation timescale of the baryons be small compared with the numerical time step is stated but never quantified; please provide an explicit criterion or discuss how to check it in practice.","section":"Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of Astronomy & Astrophysics and addresses an important gap in simulation methodology. The 'first scheme' claim should be checked against the literature for earlier or concurrent implementations of DM-baryon momentum exchange in cosmological codes; if the r=1000 test is withdrawn or substantially revised, that claim should be adjusted. The authors are honest about the scheme's limitations, which is a strength and should be preserved in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers something the field has needed: a practical scheme for including DM-baryon scattering in cosmological simulations. The virtual-particle construction is clever, and reusing the SIDM scattering routine is a sensible engineering choice. The test suite is thorough for a methods paper—heat conduction in both directions, momentum transfer, unequal-mass and isotropic scattering—and the SPH runs agree well with the analytic expectations. The implementation in OpenGadget3 is real, the limitations are discussed honestly, and the halo-collapse application shows the scheme is ready for physics use.\n\nThe soft spot is the r=1000 test, which is also the paper's headline claim. With N_bary=46656, N_DM=10^5, and equal total masses, m_bary/m_DM ≈ 2.14, so m_virt = 1000 m_DM ≈ 467 m_bary. That makes Eq. (11) have no real solution; the virtual particle cannot be created under the paper's own energy bookkeeping. The paper never acknowledges this, and the agreement in Fig. 11 may therefore come partly from the rejection scheme and artificial heat conduction rather than from the physical scattering kernel. The authors should either rerun with a configuration satisfying m_bary > m_virt, or explain why the creation constraint is suspended in the frequent small-angle branch.\n\nTwo smaller issues: velocity-dependent cross-sections are claimed in the abstract but never actually tested—all benchmarks are velocity-independent. And the momentum-transfer comparison uses an analytic interpolation with hand-chosen weights (a=3.0, b=0.1), which weakens what would otherwise be a nice test. MFM energy non-conservation at the few-percent level is a limitation, not a flaw, but it means the scheme's accuracy is solver-dependent.\n\nOverall: this is a worthwhile methods paper that deserves a serious referee. The central machinery is sound, and the limitations are mostly acknowledged. The r=1000 test needs to be fixed or honestly reframed, and a velocity-dependent test should be added. I would recommend conditional acceptance after major revision.","headline":"A genuinely new and useful method for simulating DM-baryon scattering in N-body/hydro codes, but the headline r=1000 test seems to violate the scheme's own mass-ratio bookkeeping, so that capability claim needs to be revisited.","tokens_in":34666,"tokens_out":8760,"would_cite":true,"duration_ms":86430,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":null,"created_at":"2026-08-16T12:32:20.591868+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":null,"supporting_citations":[],"review_version":1}