{"id":"6215cfd3-9266-4418-8856-5b15c621a833","arxiv_id":"1908.02508","paper_version":3,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Including scatter in the ULDM core-halo mass relation can make ULDM halos less cuspy than NFW up to 10^12 M_sun, but SPARC rotation curves do not discriminate between ULDM and CDM.","lead":"This astrophysics paper tests whether ultralight dark matter (ULDM) really solves the 'core-cusp' problem, the mismatch between simulated cuspy dark matter halos and observed cores in dwarf galaxies. It finds that allowing realistic scatter in the ULDM core-halo relation creates profiles that can be less dense than Cold Dark Matter halos, but current galaxy rotation data cannot clearly favor either model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that ULDM profiles can be less cuspy than NFW up to 10^12 Msun, and hence that the comparison is inconclusive, depends on an uncalibrated ±50% scatter in the core-halo mass relation.","rationale":"The reader's weakest_assumption identifies the same point: the ±50% scatter in the ULDM core-halo mass relation is not empirically calibrated, and the paper itself disclaims universal application of that relation. My concern is that this assumption is load-bearing for the two most concrete claims: (1) that feasible ULDM profiles exist up to 10^12 Msun which are less cuspy than NFW, and (2) that the two models are observationally indistinguishable in the SPARC comparison. If the true scatter is smaller, the blue band in Figure 2 narrows, and at the high-mass/high-particle-mass corner (Mvir = 10^12 Msun, m22 = 2.5) the ULDM profiles may be systematically more peaked than NFW, reproducing the earlier conclusion of Ref. [37] that ULDM can exacerbate the core-cusp problem. That would not destroy the paper's meta-point about baryonic feedback and limited data, but it would weaken the stated finding that ULDM remains viable up to 10^12 Msun and that current semi-analytic comparisons are inherently inconclusive. The paper is honest about this limitation and explicitly says more work is needed, so a rejection is not warranted. A conditional acceptance is appropriate: the authors should either calibrate the scatter from the underlying simulations or present the narrow-scatter/zero-scatter case as a robustness check, and soften the corresponding claims if the crossing disappears.","tokens_in":13465,"tokens_out":5091,"duration_ms":59450,"concrete_test":"Recompute Figure 2 using the actual halo-to-halo scatter in the core-halo mass relation measured from the Schive et al. simulations (or from a cosmological ULDM simulation such as Veltmaat et al. 2018), with the same percentile coverage as the NFW concentration band (e.g., 2σ). Specifically, test whether the 2σ ULDM central-density band at Mvir = 10^12 Msun and m22 = 2.5 dips below the NFW 2σ band at radii 10^-3 to 10^-2 rvir. Also run the zero-scatter case (Mc = Mcp exactly). If the blue band no longer crosses the red band, the 'less cuspy' and 'inconclusive' conclusions are artifacts of the arbitrary ±50% scatter.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central supporting claim—that for Mvir up to 10^12 Msun there exist 'feasible' ULDM profiles with central densities below NFW, making the model comparison inconclusive—rests on taking Mc = Mcp ±50% as the plausible range (Section 2, Eqs. 11-12). This range is not derived from a measured distribution of core masses. The paper itself states that 'universal application of the core-halo mass relation cannot be fully justified' and that the small sample in Ref. [38] precludes detailed statistical analysis. Moreover, the ±50% band is not matched to the ±2σ band used for NFW concentrations, so the blue and red shaded regions in Figure 2 are not commensurate containment intervals. If the true scatter is smaller, or if the core-halo relation extrapolates differently to Mvir = 10^12 Msun (outside the 10^8-10^11 Msun range probed by [38]), the blue region shrinks. At m22 = 2.5 and Mvir = 10^12 Msun, the ULDM band could lie entirely above the NFW band, reversing the 'less cuspy' claim and making the simplified models potentially discriminating rather than inconclusive. The headline conclusion that simplified semi-analytic DM-only models 'cannot meaningfully test' these scenarios is therefore partly an inference from an uncalibrated scatter assumption, not from a robust statistical analysis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the core-cusp problem as a discriminator between ultralight dark matter (ULDM) and standard cold dark matter. Using the piecewise semi-analytic ULDM halo model of Robles, Bullock, and Boylan-Kolchin, the authors allow for scatter in the Schive et al. core-halo mass relation, and compare the resulting density profiles and rotation-curve predictions with NFW halos for dwarf-to-LMC-mass halos (10^11 to 10^12 M_sun). They find that, for some choices of halo mass and ULDM particle mass, the allowed scatter yields ULDM profiles that are less cuspy than their NFW counterparts at observationally accessible radii, in contrast to the no-scatter case. Comparisons with SPARC rotation curves are presented as qualitative and inconclusive. The paper concludes that simplified semi-analytic DM-only models cannot meaningfully test the core-cusp problem and that robust conclusions will require better small-radius data, tighter particle-mass constraints, and simulations with baryonic feedback.","tokens_in":13772,"tokens_out":5128,"duration_ms":59260,"significance":"If taken as a cautionary statement, the paper has a worthwhile message: the ULDM core-halo relation, including its scatter, is not yet well enough characterized to make strong model-selection claims from DM-only semi-analytic profiles. The paper is transparent about its limitations, explicitly warns against over-interpreting the core-halo relation, and refrains from over-claiming a detection or exclusion. Its parameterization and the appendices on the core-halo scaling and SPARC errors are useful. The main caveat is that the central 'feasible profiles' result depends on an assumed, uncalibrated scatter range, and the comparison band for NFW is a different statistical containment interval; these issues need to be addressed before the main conclusion is fully persuasive.","major_comments":[{"comment":"The central existence claim—that ULDM profiles less cuspy than NFW remain feasible up to Mvir = 10^12 M_sun—is built directly on the assumed ±50% scatter in the core mass Mc. The paper itself notes that the small sample in Schive et al. precludes a detailed statistical analysis and that universal application of the core-halo relation cannot be fully justified. Because the lower envelope of the blue band in Figure 2 is set by this assumed range, the claim is conditional on an uncalibrated quantile of the core-mass distribution. To make the conclusion robust, the authors should either calibrate the scatter from available simulations/observations or, at minimum, show a sensitivity study with smaller scatter values (e.g., ±25%, ±10%) and state explicitly at what scatter the 'less cuspy' crossover disappears for m22 = 2.5 and Mvir = 10^12 M_sun.","section":"Section 2, Eqs. (11)-(12) and Figure 2"},{"comment":"The blue ULDM band (±50% 'possible' scatter in Mc) and the red NFW band (±2σ concentration scatter) are not commensurate containment intervals. A ±50% range could correspond to a very different probability content than the nominal 2σ range, depending on the (unknown) distribution of core masses. As a result, the overlap or separation of the bands in Figure 2 is not a statistical statement about the relative plausibility of the two models. The figure and the surrounding text should either use the same containment level for both models, with an explicitly stated error distribution, or be labeled clearly as illustrative envelopes with different quantile meanings.","section":"Figure 2 and Section 2"},{"comment":"The observational comparison is qualitative and uses inconsistent parameter choices across the two figures. In particular, Figure 4 adopts m22 = 0.1, which the authors themselves note is in tension with Lyman-alpha forest and other constraints, and the best-fitting curves in Figures 3 and 4 use different Mvir values. The conclusion that 'neither theoretical ULDM nor CDM model can reliably reproduce the data across a broad range' is therefore not demonstrated by the presented fits: the m22 = 0.1 fit lies outside the preferred 0.8-2.5 m22 window, and no quantitative model-selection statistic is applied, although BIC is discussed in Section 4. At minimum, the paper should state more carefully that this is an illustrative demonstration of non-uniqueness, not an empirical test, and should consider whether the qualitative conclusion survives if m22 is restricted to the preferred range.","section":"Section 3, Figures 3-4, and Section 4"}],"minor_comments":[{"comment":"The colors identifying ULDM and NFW are swapped relative to Figure 2: the captions for Figures 3 and 4 say the shaded blue region is NFW and the shaded red region is ULDM, whereas in Figure 2 blue is ULDM and red is NFW. Please make the color scheme consistent across all figures.","section":"Figures 3 and 4"},{"comment":"The text contains a typographical error, 'UDLM' should read 'ULDM' in the sentence beginning 'the relative merits of the UDLM and NFW profiles'.","section":"Section 4, first paragraph"},{"comment":"The sentence 'This spread of data may suggest that grouping galaxies based on asymptotic velocities alone is an insufficient method of characterisation' is repeated nearly verbatim in the same paragraph. Please remove the duplicate.","section":"Appendix C"},{"comment":"The phrase 'line-of-site velocities' should be 'line-of-sight velocities'.","section":"Appendix C"},{"comment":"Several references contain corrupted accented characters (e.g., [2] 'Pseudo-Goldsone-boson', [3] 'C. YÃĺche', [26] 'MacciÚ'). Please clean up the bibliography encoding.","section":"References"},{"comment":"Error bars are omitted from the main rotation-curve figures and deferred to Appendix C; a sentence near the figures indicating this deferral, with a brief statement of the typical error size, would help the reader judge the fits without flipping to the appendix.","section":"Section 3, Figures 3-4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a fair and honest cautionary study, and the topic is appropriate for a general astrophysics journal. The main concern is that the headline conclusion is an inference from an assumed scatter range rather than from a calibrated distribution; this is fixable with a sensitivity analysis and a more careful framing of the comparison bands. I would not recommend rejection, but the revision should directly address the scatter dependence and the incommensurate statistical intervals before the paper is accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is explicit treatment of scatter in the ULDM core-halo mass relation, something Robles et al. left out. Kendall and Easther show that with ±50% scatter in core mass, you can get ULDM profiles less cuspy than NFW for halos up to 10^12 Msun, which reverses the earlier claim that ULDM makes the core-cusp problem worse at the high-mass end. The SPARC comparison is deliberately qualitative, and they are upfront that neither model fits well and that baryonic physics and data limitations preclude strong statements. That honesty is the paper's best feature.\n\nThe soft spot is exactly what your stress-test note flags: the ±50% scatter is not derived from a measured distribution. The paper itself says universal application of the core-halo relation cannot be fully justified and that the Schive et al. sample is too small for statistics. So the size of the blue region in Figure 2 is essentially an assumption, not a result. And it is not commensurate with the red 2σ NFW band, since ±50% is not a stated confidence interval. If the true scatter is smaller, or the relation extrapolates differently above 10^11 Msun, the 'less cuspy' profiles for m22 = 2.5 and Mvir = 10^12 can disappear, and the comparison reverts to Robles et al. That is a real weakness in the paper's most concrete claim.\n\nBut the broader conclusion—that simplified semi-analytic DM-only models, combined with current SPARC data, cannot meaningfully test the core-cusp problem—does not depend on the exact scatter. The systematic uncertainties (baryonic feedback, stellar mass-to-light ratio, particle mass, limited radial coverage) are real and acknowledged. So the stress-test concern narrows the headline but does not sink the paper. What the authors should do in revision is either calibrate the scatter from some actual distribution, or explicitly label the ±50% band as illustrative and compare like with like (same containment probability for both models).\n\nThis paper deserves a serious referee. It is a useful caution against over-interpreting semi-analytic DM-only fits, and it highlights a genuine gap in the ULDM literature. With the scatter band properly framed, it would be a solid, if modest, contribution. I would send it out, with instructions to push on the calibration issue.","headline":"A modest and honest cautionary paper whose specific 'ULDM can be less cuspy up to 10^12 Msun' claim leans on an uncalibrated scatter assumption, but whose broader 'these simplified models cannot meaningfully discriminate' conclusion survives that weak spot.","tokens_in":14318,"tokens_out":1530,"would_cite":true,"duration_ms":19554,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d"],"model":"deepseek-v4-flash","headline":"With realistic scatter, ultralight dark matter can match cold dark matter's cusps, so current dwarf data cannot decide between the two models.","keywords":["core-cusp problem","ultralight dark matter","fuzzy dark matter","solitonic cores","NFW profile","rotation curves","core-halo mass relation","dark matter halos"],"falsifier":"A determination of the true scatter in the ULDM core-halo mass relation from higher-resolution simulations—showing, say, that core masses for a given virial mass vary by less than $\\pm 20\\%$—would shrink the band of feasible ULDM profiles, potentially restoring the earlier conclusion that ULDM worsens the core-cusp problem at high dwarf masses. Alternatively, a single well-measured dwarf rotation curve at $M_{\\mathrm{vir}} \\sim 10^{12}\\,M_\\odot$ that reaches inside the radii where ULDM and NFW predictions separate and cleanly excludes one model band would settle the direction of the discrepancy.","tokens_in":13258,"feed_emoji":"🌌","tokens_out":9244,"duration_ms":83858,"temperature":0.7,"pith_summary":"The paper argues that simplified, dark-matter-only analytic models cannot yet decide whether ultralight dark matter (ULDM) or cold dark matter (CDM) better matches the inner regions of dwarf galaxies, because both models can be adjusted to fit subsets of the data. It revisits the core-cusp problem, the mismatch between the steep 'cusp' of simulated CDM halos and the flatter 'cores' inferred from observations. The authors show that once a ±50% scatter is allowed in the ULDM core-halo mass relation, there exist feasible ULDM halos up to $10^{12}\\,M_\\odot$ that are less cuspy than comparable NFW halos, contradicting earlier claims that ULDM makes the problem worse at these masses. However, fits of these models to the observed rotation curves do not strongly favor either model, so the paper concludes that robust tests require more comprehensive observations and simulations that include baryonic feedback.","feed_headline":"Scatter erases ULDM's cusp advantage over CDM","feed_subtitle":"With ±50% scatter in ULDM's core-halo relation, both models fit dwarf rotation curves, so data can't decide.","key_machinery":"The key machinery is the piecewise semi-analytic ULDM halo profile: an inner solitonic core, the ground state of the Schrödinger–Poisson system, matched at a transition radius $r_\\alpha = \\alpha r_c$ to an outer NFW-like tail. The core's central density $\\rho_c$ and radius $r_c$ are tied to the halo's virial mass by the core-halo mass relation, and the paper treats that relation as a statistical band by allowing a $\\pm 50\\%$ scatter in core mass $M_c$ (equivalently a range in the soliton scaling parameter $\\gamma$ from $\\gamma_p/4$ to $9\\gamma_p/4$). This scatter, combined with the standard scatter in NFW concentration, generates the shaded profile bands into which observations are compared via converted circular-velocity curves.","core_discovery":"The central discovery is that the earlier conclusion that ULDM worsens the core-cusp problem for dwarf galaxies with virial masses above roughly $10^{11}\\,M_\\odot$ is not robust. Within the piecewise soliton-plus-NFW parameterization, the ULDM core-halo mass relation is not a single line but a band: allowing a $\\pm 50\\%$ scatter in the soliton core mass for a given virial mass produces ULDM profiles whose central densities can be lower than their NFW counterparts for halos up to $10^{12}\\,M_\\odot$. When these profiles are converted to rotation curves and compared to the observational database, neither the ULDM band nor the NFW concentration-scatter band provides a convincing fit across the full radial range. The paper therefore concludes that semi-analytic dark-matter-only comparisons, with limited data and without baryonic feedback, cannot meaningfully discriminate between ULDM and CDM, and that stronger claims about the core-cusp problem in either direction are premature.","pith_inferences":["The paper's conclusion likely generalizes beyond ULDM: any dark-matter model with a cored central profile and a free core-size parameter will face the same degeneracy when fitted to a handful of inner rotation-curve points.","A physically motivated scatter in the core-halo mass relation could be extracted from future ULDM cosmological simulations; the parameterization introduced here offers a direct template for that measurement.","Because baryonic feedback can also produce flatter cores, the cleanest discriminator between ULDM and CDM may come from scaling relations, such as core mass versus halo mass, rather than from individual profile fits.","The paper's caution implies that single-galaxy claims of detected ULDM soliton cores should be treated as provisional until the scatter distribution and baryonic influence are quantified."],"forward_implications":["Earlier claims that ULDM worsens the core-cusp problem for dwarf halos above about $10^{11}\\,M_\\odot$ are not robust once scatter in the core-halo mass relation is included.","The rotation-curve data used here do not favor either ULDM or CDM; both models can fit subsets of galaxies with different parameter choices, so no model-selection conclusion is warranted from these data.","Discriminating between ULDM and CDM on small scales will require rotation curves covering a much wider radial range, better constraints on the ULDM particle mass, and simulations that include baryonic feedback.","The ULDM particle mass cannot be pinned down by current rotation-curve fits because the fits are degenerate with core-mass scatter and transition-radius choices."],"supporting_citations":[{"why":"Supplies the core-halo mass relation and the ±50% core-mass scatter used to build the ULDM profile band.","marker":"[38]"},{"why":"The earlier semi-analytic ULDM-versus-NFW comparison whose conclusion about worsened cusps above about $10^{11}\\,M_\\odot$ this paper revisits.","marker":"[37]"},{"why":"Provides the observed rotation curves used to compare the theoretical profile bands.","marker":"[21]"},{"why":"Gives the NFW concentration parameter distribution and its scatter used for the CDM band.","marker":"[26]"},{"why":"Shows temporal fluctuations in soliton core densities of up to a factor of two, supporting the treatment of the core-halo relation as a statistical band.","marker":"[45]"},{"why":"Demonstrates that baryonic feedback can create cores in CDM halos, motivating the paper's caution that dark-matter-only models cannot test the small-scale predictions.","marker":"[5]"}],"fun_headline_variants":["ULDM vs CDM: scatter blurs core-cusp verdict","Core-cusp test: ULDM and CDM tie with scatter","Scatter in ULDM halos undermines cusp solution","Both dark matter models fit dwarfs when scatter allowed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument that ULDM can avoid cuspy centers for halos up to $10^{12}\\,M_\\odot$ rests on the assumed $\\pm 50\\%$ scatter in the core-halo mass relation, a spread that the paper acknowledges cannot be fully justified from the limited simulations currently available.","fun_headline_variants_meta":{"raw":{"variants":["ULDM vs CDM: scatter blurs core-cusp verdict","Core-cusp test: ULDM and CDM tie with scatter","Scatter in ULDM halos undermines cusp solution","Both dark matter models fit dwarfs when scatter allowed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000743,"raw_usage":{"total_tokens":3369,"prompt_tokens":1053,"completion_tokens":2316,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":2245}},"tokens_in":669,"tokens_out":2316,"duration_ms":16108,"temperature":1.0,"reasoning_tokens":2245,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:41:31.938719+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A determination of the true scatter in the ULDM core-halo mass relation from higher-resolution simulations—showing, say, that core masses for a given virial mass vary by less than $\\pm 20\\%$—would shrink the band of feasible ULDM profiles, potentially restoring the earlier conclusion that ULDM worsens the core-cusp problem at high dwarf masses. Alternatively, a single well-measured dwarf rotation curve at $M_{\\mathrm{vir}} \\sim 10^{12}\\,M_\\odot$ that reaches inside the radii where ULDM and NFW predictions separate and cleanly excludes one model band would settle the direction of the discrepancy.","supporting_citations":[{"cited_title":"Schive, M.-H","cited_arxiv_id":null,"evidence_quote":"Supplies the core-halo mass relation and the ±50% core-mass scatter used to build the ULDM profile band."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The earlier semi-analytic ULDM-versus-NFW comparison whose conclusion about worsened cusps above about $10^{11}\\,M_\\odot$ this paper revisits."},{"cited_title":"Lelli, S","cited_arxiv_id":null,"evidence_quote":"Provides the observed rotation curves used to compare the theoretical profile bands."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the NFW concentration parameter distribution and its scatter used for the CDM band."},{"cited_title":"Veltmaat, J","cited_arxiv_id":null,"evidence_quote":"Shows temporal fluctuations in soliton core densities of up to a factor of two, supporting the treatment of the core-halo relation as a statistical band."}],"review_version":1}