{"id":"12e1f0fc-7080-439f-86ff-88069ad00cd9","arxiv_id":"2411.08956","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The offset between lensing-based and simulated concentration-mass relations can be largely reconciled by adopting more flexible dark matter profiles and more extreme stellar initial mass functions.","lead":"Strong-lensing studies of galaxy haloes find higher dark matter concentrations than simulations predict, but this work shows the gap mostly disappears when the assumed dark matter profile shape and stellar mass function are allowed to vary. The result points to model assumptions, not new physics, as the source of the tension, with the inner mass profile as the key unknown.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reconciliation rests on an unconstrained GNFW slope; Sect. 4.4 finds no preference for non-NFW profiles, so the Conclusion attributes the c-M offset to a parameter the data do not pin down.","rationale":"The reader's verdict is CONDITIONAL, and my stress test does not move it: the paper is a careful sensitivity analysis of how modeling choices (radial coverage, enclosed-mass distortion, GNFW slope, IMF) shift the inferred c-M relation. The common PixeLens treatment of observed and EAGLE systems is a real strength, and the resolution/IMF sensitivity experiments are informative. The non-parametric R90/R50 comparison and the residual patterns are suggestive, but the authors themselves note the KS test ignores errors and the GNFW fits do not prefer a non-NFW slope (Sect. 4.4). The strongest version of the conclusion - that the simulated/observed offset is 'largely attributed to differences in the inner mass profiles' - requires that the inner slopes of real lens haloes are actually different from NFW, or at least that the data identify a different gamma for the lens sample than for EAGLE. The paper does not show this; it shows that if one assumes a shallower or steeper inner slope, the derived c and M move by about 0.3 dex per 0.5 step in gammaGNFW. That is a model-dependence result, not a measurement of the inner profile. Since the reader already captured this as the weakest assumption, no verdict change is needed. The proposed concrete test - allowing gamma free and comparing posterior distributions - would settle whether the reconciliation is physical or purely a freedom-of-fit artifact.","tokens_in":18332,"tokens_out":5917,"duration_ms":57854,"concrete_test":"For every lens realization and EAGLE projection, fit the DM residual profiles with GNFW allowing gamma to vary freely over [0, 2] and compute DeltaBIC (or Delta ln Z) relative to NFW; then compare the inferred gamma distributions (posterior mean and credible interval) of the lens and EAGLE samples. Also recompute the c-M relation using each sample's maximum-a-posteriori gamma. If the lens gamma posteriors are consistent with 1 and overlap the EAGLE posteriors, while the c-M offset at the best-fit gamma remains, the reconciliation is an artifact of unconstrained GNFW freedom rather than evidence for different inner profiles; the Conclusion should be downgraded to a sensitivity analysis.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Conclusion) is that the c-M discrepancy between lenses and EAGLE 'can largely be attributed to differences in the inner mass profiles and the assumed dark matter density profile.' The load-bearing support would be evidence that the lens and simulated DM residual profiles actually have different inner logarithmic slopes. That evidence is missing. In Sect. 4.4 the authors fit GNFW models with gammaGNFW in {0, 1, 1.5, 2} and find that gammaGNFW = 0 and 1 have similar goodness-of-fit to NFW; they state 'we find no evidence that an NFW is a worse representation of the dark matter profiles.' Consequently the strong monotonic shift of log c with gammaGNFW shown in Fig. 5 is a demonstration of the profile-concentration degeneracy, not an inference that the lenses have shallower or steeper cores. The same applies to the IMF: bottom-heavy slopes are 'extreme choices' not preferred by the data, and several lenses rule them out entirely. Therefore the observed tension is not actually resolved; it is absorbed by a parameter (gammaGNFW, or IMF slope) that the lens data do not constrain. The paper is honest about this ('The large error bars, however, prevent a definitive conclusion'), but the Conclusion drops that hedging. The central claim as stated thus overstates the evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates possible sources of the offset between strong-lensing-based and EAGLE-simulation concentration–virial mass (c-M) relations, using 18 lens systems, PixeLens mass reconstructions, stellar mass maps with two IMF choices, and orientation-averaged EAGLE projections. It sequentially studies the effect of inner radial resolution (§4.1), distortion of the enclosed mass profile (§4.2), a non-parametric concentration R90/R50 (§4.3), generalized NFW profile slopes (§4.4), and IMF slopes (§4.5). The authors conclude that much of the observed tension can be attributed to differences in the inner mass profiles and to the assumed dark matter density profile, with a smaller methodological effect from fitting NFW to poorly constrained inner profiles.","tokens_in":18621,"tokens_out":4261,"duration_ms":40128,"significance":"If the central claim were established, the paper would be significant: the long-standing lensing-vs-simulation c-M tension would be reinterpreted as a modeling artifact rather than as evidence for new physics. The paper has real strengths: it applies the same pipeline to observed and simulated lenses, it includes a non-parametric concentration measure, it is unusually candid about its limitations, and the distortion experiment in §4.2 is a useful quantitative sensitivity analysis. However, as detailed below, the evidence does not support the causal attribution in the Conclusion; the paper is better read as a systematic exploration of how much the c-M relation can move under plausible modeling choices than as a demonstration that the offset is actually caused by inner-profile differences.","major_comments":[{"comment":"The central conclusion that variations in the c-M relation 'can largely be attributed to differences in the inner mass profiles and the assumed dark matter density profile' is not supported by the GNFW analysis. In §4.4 the authors state that 'we find no evidence that an NFW is a worse representation of the dark matter profiles that we produced with our method' and that γGNFW=0 and 1 give similar goodness-of-fit; the strong monotonic shift of log c with γGNFW in Fig. 5 therefore demonstrates a profile-concentration degeneracy (the scale radius rs is redefined by each GNFW model), not that the lens haloes actually have non-NFW inner slopes. A reconciliation achieved by letting an unconstrained slope float is a possibility, not an attribution; the Conclusion should be reframed accordingly, and ideally accompanied by constraints on γGNFW from the fit residuals.","section":"§4.4 and Conclusion"},{"comment":"The distortion experiment shows that a family of transformations f(x,L) that preserve M(<Rlens) can spread NFW-fitted c-M values along the observed elongation, but the logit parameters (a=1, k=10, x0=0.5) and the L range are chosen ad hoc, and the experiment does not test whether the actual PixeLens profiles of the 18 lenses have this type of distortion. As written, the experiment is a sensitivity demonstration ('we recreated the uncertainty of the ensemble of lens models'), not evidence about the real inner profiles. Without a quantitative comparison between the distorted EAGLE profiles and the reconstructed lens profiles over the same radial range, this experiment cannot support the causal attribution in the Conclusion.","section":"§4.2, Eq. (1)"},{"comment":"The two-sample KS test used to claim that lens and EAGLE samples are 'drawn from different populations' ignores the quoted measurement uncertainties; the authors admit that 'the error bars are not taken into account for the KS test.' Given the large asymmetric errors in Table 2 (e.g., c=43.97+12.1−0.65 for J1402 and R90/R50 upper errors reaching 7.65 for J1538), a KS test on point estimates is not a reliable basis for the 'tentative evidence' that the discrepancy traces to inner-profile differences. The paper should implement an error-aware comparison (e.g., Monte Carlo sampling of the posterior distributions or a generalized distance test) or explicitly downgrade this item to a hypothesis-generating observation.","section":"§4.3, Fig. 4"},{"comment":"The IMF-based reconciliation rests on 'extreme choices' of ΓBM and µ2PL, and the authors note that bottom-heavy 2PL slopes can be ruled out for J0037 (µ>1.8) and J0044 (µ>1.5), with only J0946 showing the trend 'with sufficient significance.' Since the lens data disfavor the very IMF slopes that produce the largest reduction in c, the IMF channel cannot be cited as evidence that the observed c-M relation is reconciled with simulations; at most it provides an upper bound on how much the IMF can contribute. The Conclusion should not present this channel as support for the central attribution claim.","section":"§4.5, Figs. 6–9"}],"minor_comments":[{"comment":"The function defined in Eq. (2) is a logistic (sigmoid) function, not a logit; rename it or adjust the notation to avoid a terminological error.","section":"Eq. (2)"},{"comment":"The horizontal axis labels in Figs. 6–9 read 'γNFW' while the text and Fig. 5 use 'γGNFW'; unify the notation.","section":"Figs. 6–9"},{"comment":"The projected mass expressions are written with 'M∼' and leave κs and the numerical coefficients unspecified; define all symbols or state explicitly that the normalization is absorbed.","section":"Eqs. (4)–(7)"},{"comment":"The asymmetric errors in column 3 are not accompanied by a confidence level; specify whether they are 68% or 90% intervals.","section":"Table 1"},{"comment":"The phrase 'changing the slope of a generic NFW profile' is ambiguous; use 'generalized NFW (GNFW) profile' consistently from the abstract onward.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its limitations, but the abstract and Conclusion present the reconciliation as a demonstration while the body repeatedly hedges it as a possibility. The paper would be a valuable systematic sensitivity study if reframed accordingly, with the causal claim replaced by a statement about how much the c-M relation can shift under flexible profile and IMF choices. The central issue is inference strength rather than technical error, so major revision is appropriate. No concerns about citation practice or novelty disclosure beyond the usual expectation that the authors cite the earlier papers on which the lens sample and method are based."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a careful, systematic study of why strong-lensing c-M estimates sit above simulation predictions. What is new is the combination: varying the GNFW inner slope, distorting the enclosed mass profile inside the Einstein radius, and pushing IMF parameters, all within the same pipeline applied to 18 real lenses and EAGLE projections. The non-parametric R90/R50 concentration is a genuinely fresh diagnostic, and it gives tentative support to the idea that the offset is tied to inner-profile differences rather than to NFW fitting alone. The authors are also admirably candid: they report success rates, rule out bottom-heavy IMFs for some lenses, and admit the KS test ignores error bars and that their own fits do not prefer the non-NFW slopes.\n\nThat last point is the soft spot, and it is load-bearing. The Conclusion says the c-M discrepancy 'can largely be attributed' to differences in inner mass profiles and the assumed DM density profile. But in Sect. 4.4 they state plainly that they find no evidence an NFW is a worse representation of the profiles. The shift in log c with gamma_GNFW is then a demonstration of the profile-concentration degeneracy, not evidence that real lenses have shallower or steeper cores. The R90/R50 KS test points the right way, but with error bars ignored and uncertainties large, it is tentative at best. The distortion experiment is a possibility proof: allowing more freedom in the inner slope spreads the c-M cloud, which is not the same as showing that freedom is realised in the lenses. The paper itself hedges in the Discussion ('the large error bars prevent a definitive conclusion'), but the Conclusion drops that hedge.\n\nThere is also mild circularity: since cvir is defined through fitted rs, changing the profile family redefines the quantity being measured. The non-parametric comparison is the right check against that, but under-powered here.\n\nSo: a solid, honest systematics study, with a headline that overstates. The reader who takes away 'the tension can be modelled away with flexible profiles' is getting the right message; the reader who thinks the tension is resolved is getting more than the evidence says. Worth citing for the method and the diagnostic; worth sending to a serious referee, with a clear request to soften the causal claim and do proper error treatment on the R90/R50 comparison.","headline":"Honest systematics exploration whose 'reconciliation' headline outruns the data.","tokens_in":19157,"tokens_out":3635,"would_cite":true,"duration_ms":30429,"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 offset between lensing-based and simulated dark-matter halo concentrations can be largely attributed to assumptions about the inner density profile and the stellar initial mass function, not to new physics.","keywords":["strong gravitational lensing","dark matter haloes","concentration-mass relation","NFW profile","generalized NFW profile","initial mass function","EAGLE simulation","non-parametric concentration"],"falsifier":"A decisive test would be to measure the inner dark matter slope of a few lens galaxies independently, for example with spatially resolved stellar kinematics inside the scale radius, and check whether the best-fit GNFW slope from lensing agrees. Alternatively, if a c-M offset persisted when the GNFW inner slope is free and the IMF is allowed to vary over independently constrained ranges, the claimed reconciliation would fail. A simpler check: determine whether the extreme bottom-heavy IMFs (e.g., $\\mu_{2PL} > 1.5$) required to shift the c-M relation are excluded by stellar-population or kinematic constraints on the same galaxies.","tokens_in":18089,"feed_emoji":"🔭","tokens_out":6246,"duration_ms":53605,"temperature":0.7,"pith_summary":"This paper attempts to establish that the long-standing tension between observed (strong-lensing) and simulated concentration-virial mass relations is largely a modeling artifact. By fitting a generalized NFW profile with a free inner slope to the dark matter component of 18 lens galaxies, and by varying the stellar IMF, the authors show that the two relations can be reconciled. If true, the high concentrations found in lensing studies do not signal new dark matter physics; they reflect the flexibility of the assumed mass profile and the stellar-to-dark matter decomposition. The paper also introduces a non-parametric concentration measure, the projected radius enclosing 90% versus 50% of the mass, which tentatively connects the offset to differences in the inner radial region of haloes, though large uncertainties prevent a definitive conclusion.","feed_headline":"Inner-profile shape reconciles lensing and simulated haloes","feed_subtitle":"Allowing the dark-matter inner slope to vary or adopting extreme IMFs removes the observed offset, the authors argue.","key_machinery":"The central object is the generalized NFW (GNFW) density profile $\\rho(r) \\propto (r/r_s)^{-\\gamma_{\\mathrm{GNFW}}} (1 + r/r_s)^{\\gamma_{\\mathrm{GNFW}}-3}$, whose inner slope $\\gamma_{\\mathrm{GNFW}}$ is free. Fitting this to the enclosed dark matter mass profile (total lensing mass minus stellar mass) changes the significance of the scale radius $r_s$; with cuspier profiles the fitted scale radius drifts outward and the concentration $r_{\\mathrm{vir}}/r_s$ falls. A secondary mechanism is the distortion function applied to the inner enclosed mass profile while keeping the mass within the lens radius fixed, which mimics reconstruction uncertainty and produces the elongated scatter seen in observed lenses. Finally, the non-parametric concentration $R_{90}/R_{50}$ is used to check that the parametric NFW fit itself is not the root cause.","core_discovery":"The central claim is that variations in the c-M relation between simulated and observed dark matter haloes can largely be attributed to differences in the inner mass profiles and the assumed dark matter density profile. The authors demonstrate that when the inner slope $\\gamma_{\\mathrm{GNFW}}$ of a generalized NFW profile is treated as a free parameter between 0 and 2, the fitted concentrations drop as the profile becomes cuspier (higher $\\gamma$), while virial masses rise; for $\\gamma=0$ and $\\gamma=1$ the goodness-of-fit is similar to the standard NFW. Along the same lines, extreme bottom-heavy stellar IMFs systematically lower concentrations. The paper argues that these choices, rather than new physics, can remove much of the offset between lensing and simulation.","pith_inferences":["The same GNFW/IMF framework could be applied to cluster-scale lenses, where the inner profile is better resolved, to test whether the reconciliation persists or the offset reappears.","A direct measurement of the inner slope $\\gamma_{\\mathrm{GNFW}}$ from combining lensing with stellar kinematics would break the degeneracy and either confirm or refute the modeling-artifact explanation.","The distortion-function experiment suggests that small systematic errors in inner-profile reconstruction can create the observed scatter; a full forward-modeling test, simulating lens images and reconstructing them with the same pipeline, would quantify this bias directly."],"forward_implications":["If the reconciliation holds, the higher concentrations found in strong-lensing studies are not evidence for exotic dark matter but follow from assuming a standard NFW profile and a canonical IMF.","The c-M relation measured from lenses should be reported as a function of the assumed GNFW inner slope; quoting only the NFW-based point may be misleading.","Cuspier inner profiles, whether from baryonic contraction or substructure, will shift galaxies to lower apparent concentrations at fixed virial mass.","The non-parametric $R_{90}/R_{50}$ measure provides a way to compare observed and simulated haloes without profile assumptions, though its current error bars limit its discriminating power."],"supporting_citations":[{"why":"Supplies the NFW profile, the baseline density model whose scale radius defines concentration in the c-M relation.","marker":"Navarro et al. 1997"},{"why":"Provides the observed lensing c-M relation and the simulation comparison baseline that this paper aims to reconcile.","marker":"Leier et al. 2022"},{"why":"Supplies the generalized NFW (GNFW) profile family used to vary the inner density slope.","marker":"Keeton 2001"},{"why":"Supplies the EAGLE simulation haloes used as the simulated comparison sample.","marker":"Schaye et al. 2015"},{"why":"Supplies the pixelised lens reconstruction method used to derive total mass profiles of the observed lenses.","marker":"Saha & Williams 2003"},{"why":"Supplies the lens sample, stellar mass maps, and the bimodal and two-power-law IMF choices used in the stellar-to-dark-matter decomposition.","marker":"Leier et al. 2016"},{"why":"Documents the earlier lensing-based c-M relation whose offset from simulations motivates this reconciliation study.","marker":"Leier et al. 2012"},{"why":"Provides an observational c-M baseline from X-ray clusters used for comparison in the c-M plane.","marker":"Buote et al. 2007"}],"fun_headline_variants":["Varying inner slope fixes halo concentration-mass mismatch","Lensing-simulation offset resolved by free inner dark matter slope","Extreme IMFs or adjustable inner slope reconcile halo scaling","Concentration-mass relation unified by profile shape choices","Dark matter inner slope plasticity bridges lensing and theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reconciliation depends on the assumption that the real dark matter haloes of lens galaxies are well represented by the GNFW family with inner slope between 0 and 2; the data themselves do not prefer these alternative profiles over a standard NFW, so if the true inner profiles differ in a way not captured by this family the reconciliation could be artificial.","fun_headline_variants_meta":{"raw":{"variants":["Varying inner slope fixes halo concentration-mass mismatch","Lensing-simulation offset resolved by free inner dark matter slope","Extreme IMFs or adjustable inner slope reconcile halo scaling","Concentration-mass relation unified by profile shape choices","Dark matter inner slope plasticity bridges lensing and theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000886,"raw_usage":{"total_tokens":3827,"prompt_tokens":946,"completion_tokens":2881,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":2814}},"tokens_in":562,"tokens_out":2881,"duration_ms":20420,"temperature":1.0,"reasoning_tokens":2814,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:13:11.287176+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to measure the inner dark matter slope of a few lens galaxies independently, for example with spatially resolved stellar kinematics inside the scale radius, and check whether the best-fit GNFW slope from lensing agrees. Alternatively, if a c-M offset persisted when the GNFW inner slope is free and the IMF is allowed to vary over independently constrained ranges, the claimed reconciliation would fail. A simpler check: determine whether the extreme bottom-heavy IMFs (e.g., $\\mu_{2PL} > 1.5$) required to shift the c-M relation are excluded by stellar-population or kinematic constraints on the same galaxies.","supporting_citations":[{"cited_title":"2022, MNRAS, 510, 24","cited_arxiv_id":null,"evidence_quote":"Provides the observed lensing c-M relation and the simulation comparison baseline that this paper aims to reconcile."},{"cited_title":"& Williams, L","cited_arxiv_id":null,"evidence_quote":"Supplies the pixelised lens reconstruction method used to derive total mass profiles of the observed lenses."},{"cited_title":"2016, MNRAS, 459, 3677–3692","cited_arxiv_id":null,"evidence_quote":"Supplies the lens sample, stellar mass maps, and the bimodal and two-power-law IMF choices used in the stellar-to-dark-matter decomposition."},{"cited_title":"2012, MNRAS, 424, 104–114","cited_arxiv_id":null,"evidence_quote":"Documents the earlier lensing-based c-M relation whose offset from simulations motivates this reconciliation study."},{"cited_title":"A., Gastaldello, F., Humphrey, P","cited_arxiv_id":null,"evidence_quote":"Provides an observational c-M baseline from X-ray clusters used for comparison in the c-M plane."}],"review_version":1}