{"id":"9e7eb6d9-7232-4322-947e-944355b4b362","arxiv_id":"2501.12636","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A new semianalytic model in Galacticus predicts that decaying dark matter flattens halo density profiles and suppresses subhalo mass functions, matching N-body simulations.","lead":"This paper presents a fast, semianalytic model for how decaying dark matter changes the inner structure of halos and the number of surviving subhalos. It matches expensive N-body simulations well enough to scan dark matter decay parameters cheaply.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uncalibrated mass-loss heating parameter gamma is degenerate with tau and vk, so the model's parameter-space exploration claim rests on an untested assumption until gamma is calibrated or marginalized.","rationale":"The paper is a careful and useful semianalytic extension of Galacticus, with public code and a clear exposition of the DDM physics. The SHMF comparison to Mau et al. (2022) is encouraging, and the model's modularity is a genuine strength. However, the claim that the model 'enables efficient and accurate exploration of DDM parameter space' is weakened by the free parameter gamma in Eq. (15). Unlike the shell-crossing prescription (Eq. 7), which is an acknowledged numerical artifact with an explicit (though not shown as a figure) comparison to N-body profiles, gamma directly controls the dominant heating channel for exactly the halos most affected by DDM (Vmax<vk; Sec. III D). The paper itself demonstrates that gamma is partially degenerate with the physical DDM parameters tau and vk (Appendix A), and it defers calibration. Without a quantitative posterior on gamma, the model's predictions for arbitrary (tau, vk) are not predictive in the sense the abstract claims; they are conditional on an arbitrary choice. The single-host N-body SHMF comparison is too coarse to break the degeneracy. This does not invalidate the model's qualitative behavior, but it means the central claim of accurate parameter-space exploration is not yet supported. The shell-crossing concern identified by the reader is real and should be addressed, but it is secondary: the paper already acknowledges the artifact and the SHMF comparison suggests it does not dominate the results at the tested mass scale. The gamma issue is more load-bearing because it threatens the model's primary purpose—deriving constraints on DDM parameters—unless gamma is calibrated or marginalized. Therefore, the reader's conditional verdict should stand, with the condition extended to include gamma calibration.","tokens_in":15970,"tokens_out":7538,"duration_ms":76099,"concrete_test":"Perform a joint Bayesian calibration of gamma, tau, and vk against the N-body SHMF data from Ref. [42] (and, if available, isolated-density profiles from Ref. [34]) using the Galacticus DDM model, with a nested-sampling or MCMC sampler. Evaluate the marginalized posterior on gamma and its correlation with tau and vk. If the 68% credible interval on gamma is wider than about 0.2, or if gamma and tau/vk are strongly degenerate, the fiducial gamma=0.5 is not validated, and the model's parameter-space exploration claim requires either calibration of gamma or explicit marginalization over it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The model's central claim is that it enables efficient and accurate exploration of DDM parameter space (Abstract, Sec. V). This requires the mass-loss heating efficiency parameter introduced in Eq. (15), gamma, to be either known or marginalized over. The paper sets gamma=0.5 using a virial-energy heuristic and defers calibration to future work (Appendix A). For halos with Vmax<vk, mass-loss heating dominates over velocity-kick heating (Sec. III D, Fig. 4), and gamma controls the degree of inner density flattening and thus subhalo tidal disruption. Appendix A shows gamma is partially degenerate with tau and vk (Fig. 7), meaning different (gamma, tau, vk) combinations can produce similar density profiles. The SHMF comparison in Fig. 5, which the paper cites as evidence that gamma=0.5 is 'reasonable,' is based on a single N-body host with large Poisson uncertainties and does not tightly constrain gamma. Without a posterior on gamma, predictions for unsimulated DDM models carry an unquantified systematic error, so the model cannot yet be claimed to enable accurate DDM parameter constraints.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a semianalytic model for two-body decaying dark matter (DDM) implemented in the Galacticus framework. The model combines adiabatic heating from daughter-particle velocity kicks with mass-loss heating and direct mass reduction, parametrized by decay lifetime τ and kick velocity vk plus one efficiency parameter γ for mass-loss heating. The authors apply the model to isolated halos and subhalo populations, predicting suppressed and flattened inner density profiles, a bend in the Rmax–Vmax relation for halos with Vmax ≲ vk, and a mass-dependent suppression of the subhalo mass function. They compare the subhalo predictions to cosmological zoom-in N-body results from Mau et al. (2022) and claim general consistency with isolated and cosmological DDM simulations, arguing that the model enables efficient exploration of DDM parameter space.","tokens_in":16214,"tokens_out":7825,"duration_ms":78927,"significance":"If the model's accuracy claims hold, this is a valuable tool for DDM phenomenology: it is open-source, modular, computationally cheap relative to N-body simulations, and it makes falsifiable predictions (e.g., the Rmax–Vmax bend and the subhalo mass function suppression). The manuscript is commendable for separating velocity-kick heating, mass-loss heating, and direct mass loss, and for testing against external N-body benchmarks without fitting the benchmark data. However, the validation evidence is currently thinner than the central claims: the mass-loss efficiency γ is uncalibrated, the density-profile consistency with Ref. [34] is asserted but not shown, and the shell-crossing prescription produces a feature absent from N-body profiles. These gaps affect the advertised ability to accurately explore parameter space, but they are addressable in revision.","major_comments":[{"comment":"The mass-loss heating efficiency γ is effectively uncalibrated. For halos with Vmax <~ vk, mass-loss heating dominates the profile evolution (Sec. III D and Fig. 4), so γ controls the inner density flattening and hence the subhalo disruption that drives the SHMF suppression in Sec. IV. Appendix A shows that γ is partially degenerate with τ and vk (Fig. 7), and the paper states that calibrating γ is beyond its scope. The SHMF agreement in Fig. 5 is based on a single N-body host with large Poisson uncertainties and does not tightly constrain γ. Consequently, the Abstract and Sec. V claim that the model enables 'efficient and accurate exploration' of DDM parameter space is not yet supported; the authors should either calibrate γ against simulations, marginalize over it in predictions, or soften the accuracy claim.","section":"Appendix A; Eq. (15)"},{"comment":"The shell-crossing prescription is load-bearing. Equation (7) determines where the heating energy ratio is frozen, and the paper itself notes that the resulting sharp feature in predicted density profiles is not visible in N-body simulations (Sec. V). This feature also creates an apparent bifurcation in the Rmax–Vmax relation (Sec. III C). Because the frozen heating ratio below rc sets the inner profile response, the absence of this feature in N-body profiles raises a correctness risk for the predicted inner flattening and the resulting subhalo disruption rates. A direct comparison with Ref. [34] profiles, or a revised shell-crossing treatment that does not produce the artifact, is needed before this part of the model can be regarded as validated.","section":"Eq. (7); Secs. III B and V"},{"comment":"The consistency of the predicted density and velocity dispersion profiles with Ref. [34] is asserted but not demonstrated: the text says 'we have checked' but provides no comparison plot or quantitative metric. Since this comparison is part of the central claim that the model matches isolated N-body simulations, the evidence should be shown explicitly, or the claim should be qualified.","section":"Sec. III B"},{"comment":"The SHMF validation rests on a single N-body host, and the quoted Poisson uncertainties are large. The agreement within these uncertainties is encouraging, but it does not by itself validate the model across the parameter space advertised in the Abstract. The authors should either present additional simulation comparisons where available or explicitly restrict the validation claim to the tested parameter range.","section":"Sec. IV B; Fig. 5"}],"minor_comments":[{"comment":"The velocity distribution is written as p(v, θ|r, s) although the text states it is independent of θ after assuming isotropy; consider using p(v|r,s) and defining θ only in the kick calculation.","section":"Eq. (8)"},{"comment":"The caption of the lower-right circular-velocity panel appears to repeat 'M = 10^9 M⊙' instead of giving the 10^10 M⊙ case; please check and correct.","section":"Fig. 2"},{"comment":"The statement that setting the second-order energy perturbation coefficient to zero is justified because this term is 'most relevant for cuspy halos' would benefit from a brief quantitative justification or a reference to a comparison.","section":"Sec. IV A"},{"comment":"The explanation that the shell-crossing feature is absent from N-body profiles 'may be due to anisotropic accretion in a cosmological environment' is speculative; consider framing it as an open question or supporting it with a test.","section":"Sec. V"},{"comment":"The label 'Mau + 2022' in the figure should be 'Mau et al. (2022)' to match the text and reference list.","section":"Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the underlying model is promising. The main gap is between the advertised accuracy for parameter-space exploration and the current state of validation, particularly the uncalibrated γ parameter and the shell-crossing artifact. I would like the revision to address these points explicitly rather than only deferring them to future work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look. This is a genuine step forward for semianalytic DDM modeling: it cleanly separates velocity-kick heating from mass-loss heating, derives the injection energy from decay kinematics, predicts a bend in the Rmax–Vmax relation, and generalizes subhalo suppression beyond the two kick velocities that N-body simulations have actually run. The implementation is in the public Galacticus code, and the validation against the Mau et al. subhalo mass functions is the right kind of check. The paper is also unusually candid about its own limitations — the shell-crossing feature, the uncalibrated gamma, the single-host SHMF comparison. That honesty earns credit.\n\nThe soft spots are real, though. The mass-loss heating efficiency gamma is set to 0.5 by a virial heuristic and never calibrated. Appendix A explicitly shows that gamma is partially degenerate with tau and vk, and the Fig. 5 comparison is one N-body host with Poisson errors — not enough to pin gamma down. So the abstract's claim that the model enables \"accurate\" exploration of DDM parameter space is ahead of the evidence. The model can efficiently generate predictions, but those predictions carry an unquantified systematic error until gamma is constrained or marginalized. That is a load-bearing caveat for anyone planning to use this to derive competitive limits.\n\nSecond, the shell-crossing prescription in Eq. (7) produces a sharp density feature that the authors themselves note is absent from N-body profiles. They hand-wave it to anisotropic accretion, but it sits exactly where the heating model stops being conservative. It does not obviously wreck the SHMF comparison, but it weakens the case that the inner-profile flattening is quantitatively right.\n\nThird, the claim of consistency with the Peter et al. profiles is asserted in the text but no comparison figure is shown. For a paper whose validation section is otherwise the main selling point, that omission is conspicuous.\n\nWho is this for? Modelers who need fast DDM subhalo predictions for satellite, Lyman-alpha, or lensing work, and who want a public code base to modify. It deserves a serious referee; I would accept it with the expectation of major revision — specifically, calibrate or marginalize gamma, show the profile comparison to Peter et al., and either fix or isolate the shell-crossing artifact so it cannot be mistaken for a physical core.","headline":"A useful, honest semianalytic DDM model with real code and real validation, but the headline claim of accurate parameter-space exploration rests on an uncalibrated gamma and an acknowledged shell-crossing artifact.","tokens_in":16708,"tokens_out":1469,"would_cite":true,"duration_ms":16947,"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":"Semianalytic model reproduces decaying-dark-matter halo simulations","keywords":["decaying dark matter","subhalo mass function","semianalytic model","adiabatic heating","halo density profiles","Rmax-Vmax relation","Galacticus","N-body simulations"],"falsifier":"A high-resolution cosmological DDM zoom-in simulation with $\\tau = 10$ Gyr and $v_k = 40$ km/s for a $10^{10}\\,M_\\odot$ halo would settle the matter: if its density profile interior to $r/R_{\\rm vir} \\approx 0.2$ does not flatten, or if its subhalo population shows a different mass-dependent survival ratio than predicted, the shell-crossing and heating model is falsified. A cheaper check is to measure whether the sharp feature at the shell-crossing radius $r_c$ appears in N-body profiles at all; the paper itself notes that it does not.","tokens_in":15767,"feed_emoji":"🌌","tokens_out":6307,"duration_ms":57791,"temperature":0.7,"pith_summary":"The paper claims that the effects of two-body decaying dark matter on halo structure can be captured by a semianalytic model that combines two physical ingredients: velocity kicks, which heat the dark matter, and mass loss, which reduces the mass of decay products. The model is implemented in the authors' open-source structure-formation code and predicts that decays flatten and lower the inner density profiles of halos, and that the resulting fragile subhalos are more easily tidally disrupted, suppressing the subhalo mass function relative to cold dark matter. The authors show that these predictions match results from both isolated and cosmological DDM N-body simulations without tuning the new heating-efficiency parameter. Because the calculation runs at a tiny fraction of the cost of simulations, it opens up parameter-space exploration for DDM, which matters for upcoming small-scale structure observations such as Milky Way satellite counts.","feed_headline":"Fast model matches dark-matter decay halo simulations","feed_subtitle":"Predicts flattened inner profiles and suppressed subhalo counts, opening cheap scans of decay lifetime and kick velocity.","key_machinery":"The model generalizes an adiabatic heating scheme originally developed for tidal heating: a halo is decomposed into spherical shells, and conservation of energy relates each shell's final radius to the energy injected by decays. The injected specific energy has two parts: a velocity-kick term computed by averaging over a truncated Maxwell-Boltzmann velocity distribution, keeping only daughter particles that remain bound, and a mass-loss term proportional to $GM(<r)/r$ with an efficiency parameter $\\gamma$. Shell crossing is handled by finding the radius $r_c$ where $dr_f/dr_i = 0$ and freezing the heating energy ratio at $r_c$ for all smaller radii; this is what produces the sharp feature in the predicted density profiles. This machinery is what converts the two microphysical DDM parameters (lifetime and kick velocity) into density profiles and, through merger trees and a tidal evolution model, into subhalo populations.","core_discovery":"DDM halos are not just lighter copies of CDM halos. The paper establishes that, in a two-body DDM model with a lifetime comparable to the age of the Universe, the velocity kick from each decay heats the halo and causes daughter particles to escape; combining this heating with the mass lost to the daughter particles flattens and suppresses the inner density profile. The strength of the effect is set by the ratio of the kick velocity to the halo's internal velocity: halos with $V_{\\rm max}$ above roughly $v_{\\rm kick}$ are affected mainly by kick heating, while lower-mass halos with $V_{\\rm max}$ below $v_{\\rm kick}$ are dominated by mass-loss-induced heating and can even be completely unbound. Because the heated, lower-density DDM subhalos are more easily stripped and tidally disrupted than cuspy CDM subhalos, the subhalo mass function is suppressed in a mass-dependent way. The central quantitative claim is that this semianalytic treatment reproduces the density, circular velocity, and velocity dispersion profiles of DDM N-body simulations, as well as the subhalo mass functions and radial distributions from cosmological zoom-in simulations, with one efficiency parameter (set to 0.5) left effectively uncalibrated.","pith_inferences":["If the model is right, existing constraints from Milky Way satellite counts at $v_k = 20$ and $40$ km/s can be extended into a continuous exclusion curve in the lifetime–kick plane, and the same pipeline can forecast sensitivity for upcoming surveys.","The shell-crossing feature the authors identify—present in their profiles but absent from N-body simulations—suggests that the fixed-heating-energy prescription inside $r_c$ may overstate the sharpness of the transition; if so, the precise disruption rates of subhalos near that radius could shift, though the overall SHMF suppression is validated by comparison.","The $\\gamma$ parameter measuring mass-loss heating efficiency is degenerate with lifetime and kick velocity for low-mass density profiles; subhalo statistics are the stated way to break this degeneracy, so a calibration run against higher-resolution DDM simulations would sharpen all parameter constraints.","The same energy-injection machinery could be adapted to other beyond-CDM scenarios, such as self-interacting or annihilating dark matter, as long as the deposited energy as a function of radius and time can be written down."],"forward_implications":["Decaying dark matter suppresses the subhalo mass function relative to CDM in a mass-dependent way, with the strongest suppression at low masses; for a $10^{12}\\,M_\\odot$ host the surviving fraction ranges from roughly a third for $\\tau=10$ Gyr, $v_k=20$ km/s to most subhalos surviving for $\\tau=80$ Gyr, $v_k=40$ km/s.","The $R_{\\rm max}$–$V_{\\rm max}$ relation bends away from CDM for low-mass DDM halos, providing a signature that can be distinguished from self-interacting dark matter, whose core-collapsing branch shifts in the opposite direction.","Halo evolution splits into two regimes: halos with $V_{\\rm max} \\gtrsim v_k$ are heated mainly by velocity kicks, while halos with $V_{\\rm max} \\lesssim v_k$ are dominated by mass-loss heating and can unbind, making DDM effects on dwarf-scale structure strongly mass-dependent.","The semianalytic model can generate DDM predictions at any lifetime and kick combination, not just the discrete values covered by existing N-body runs, enabling constraints from Milky Way satellite observations to be evaluated continuously in parameter space.","The model's computational cost is a small fraction of N-body simulation cost, so it can be used to derive constraints from upcoming small-scale structure surveys over a wide range of DDM parameters."],"supporting_citations":[{"why":"The open-source structure-formation code in which the DDM heating and mass-loss models are implemented, providing merger trees and halo evolution infrastructure.","marker":"[45]"},{"why":"Supplies the adiabatic heating framework (shell energy conservation) that the DDM energy injection extends.","marker":"[46]"},{"why":"Provides the shell-crossing treatment ($dr_f/dr_i = 0$ and frozen heating ratio) used to handle non-monotonic shell behavior.","marker":"[59]"},{"why":"Isolated DDM N-body simulations whose density and velocity dispersion profiles the model is checked against.","marker":"[34]"},{"why":"Earlier DDM density-profile fitting and Milky Way satellite constraints whose simulation limitations motivate the semianalytic approach.","marker":"[41]"},{"why":"Cosmological DDM zoom-in simulations and subhalo mass function measurements that provide the primary validation data for subhalo population predictions.","marker":"[42]"},{"why":"Supplies the finite-resolution artificial disruption model and tidal disruption framework needed to compare subhalo populations with simulations.","marker":"[48]"},{"why":"Provides the tidal heating model (with second-order coefficient) used for CDM and DDM subhalo evolution.","marker":"[49]"}],"fun_headline_variants":["Fast model for decaying dark matter halos matches simulations","Semianalytic DDM model reproduces halo flattening and subhalo loss","Decaying dark matter halos get a fast semianalytic treatment","Efficient semianalytic DDM model fits N-body halo simulations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The shell-crossing prescription—finding the radius where $dr_f/dr_i = 0$ and then holding the heating energy fixed at that radius for all smaller radii—is what produces the predicted inner flattening and the subhalo disruption rates; if that treatment is wrong, the model's central profile and subhalo mass function predictions lose their support.","fun_headline_variants_meta":{"raw":{"variants":["Fast model for decaying dark matter halos matches simulations","Semianalytic DDM model reproduces halo flattening and subhalo loss","Decaying dark matter halos get a fast semianalytic treatment","Efficient semianalytic DDM model fits N-body halo simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000377,"raw_usage":{"total_tokens":2040,"prompt_tokens":1011,"completion_tokens":1029,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":951}},"tokens_in":627,"tokens_out":1029,"duration_ms":10107,"temperature":1.0,"reasoning_tokens":951,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:58:31.722263+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution cosmological DDM zoom-in simulation with $\\tau = 10$ Gyr and $v_k = 40$ km/s for a $10^{10}\\,M_\\odot$ halo would settle the matter: if its density profile interior to $r/R_{\\rm vir} \\approx 0.2$ does not flatten, or if its subhalo population shows a different mass-dependent survival ratio than predicted, the shell-crossing and heating model is falsified. A cheaper check is to measure whether the sharp feature at the shell-crossing radius $r_c$ appears in N-body profiles at all; the paper itself notes that it does not.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The open-source structure-formation code in which the DDM heating and mass-loss models are implemented, providing merger trees and halo evolution infrastructure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the adiabatic heating framework (shell energy conservation) that the DDM energy injection extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the shell-crossing treatment ($dr_f/dr_i = 0$ and frozen heating ratio) used to handle non-monotonic shell behavior."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Isolated DDM N-body simulations whose density and velocity dispersion profiles the model is checked against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier DDM density-profile fitting and Milky Way satellite constraints whose simulation limitations motivate the semianalytic approach."},{"cited_title":"Mau et al","cited_arxiv_id":null,"evidence_quote":"Cosmological DDM zoom-in simulations and subhalo mass function measurements that provide the primary validation data for subhalo population predictions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-resolution artificial disruption model and tidal disruption framework needed to compare subhalo populations with simulations."},{"cited_title":"Du et al., Phys","cited_arxiv_id":null,"evidence_quote":"Provides the tidal heating model (with second-order coefficient) used for CDM and DDM subhalo evolution."}],"review_version":1}