{"id":"9ec9cd5a-377b-4d81-a134-9fc1c66b329f","arxiv_id":"1909.02573","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"First measurement of the dark matter halo mass-concentration relation below 10^9 solar masses, using flux ratios of 11 strongly lensed quasars.","lead":"Using 11 quadruple-image gravitational lenses, the authors measure the mass-concentration relation of dark matter halos on sub-galactic scales for the first time. Their constraint on the concentration of 10^8 solar mass halos is consistent with cold dark matter predictions, establishing strong lensing as a probe of dark matter structure.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Subhalo tidal-stripping approximation in §2.2 may bias the inferred c0; no test quantifies the shift.","rationale":"The reader's strongest claim is the constraint on c0, the normalization of the mass-concentration relation. The condition for this claim to hold is that the forward model correctly maps c0 to the lensing observables. The weakest place in that chain is the treatment of subhalo concentrations at infall without a full tidal-stripping model, exactly the assumption flagged by the reader and explicitly acknowledged in Section 2.2. I agree that this is the most load-bearing concern because it directly affects the density profiles of the objects that generate the flux-ratio perturbations and hence the inferred c0; a bias here would change the central result rather than merely widen its error bars. The paper does not provide any test of how sensitive c0 is to this approximation, and the deferred details in Gilman et al. (2019b) are not part of this submission. I also note a secondary concern that the posterior for beta is described as \"peaked towards higher values\" and may be truncated by the prior upper bound at 1.3, which could affect the marginal interval on c0; this is worth a robustness check but is less central than the tidal-stripping systematic. The paper otherwise makes a credible case: the Bayesian forward-modeling framework has been developed in prior work, the sample of 11 quads is real, the analysis marginalizes over nuisance parameters, and the quoted uncertainties are large enough that the result is consistent with CDM. Nothing in my reading requires changing the reader's conditional verdict; the result is plausible but conditional on the validity of the infall-concentration approximation and on the details in the companion paper. A concrete mock test comparing recovery when the fitting model is the same infall approximation but the simulated data include realistic tidal stripping would settle whether the concern lands.","tokens_in":8109,"tokens_out":11092,"duration_ms":116403,"concrete_test":"Generate mock lens populations with a known input c0 and with subhalo density profiles taken from a model that includes tidal stripping, for example subhalos in an N-body simulation or the Green & van den Bosch (2019) semi-analytic model. Run the same forward-model inference used here, which treats subhalos as NFW with concentration c(M, z_infall) and a truncation radius only, and compare the recovered c0 to the known input. If the input c0 lies outside the 68% credible interval of the recovered posterior, then the infall approximation is a load-bearing systematic and the quoted constraint c0=12^{+6}_{-5} is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.2 assigns subhalo concentrations at infall: \"we simply evaluate subhalo concentrations at the time of infall when the mass-concentration relation in Equation 4 is valid.\" The subhalo mass function in Equation 6 is used to render subhalos at the host redshift, so the mass m is naturally the bound mass at that epoch. Tidal stripping reduces that bound mass below the mass at infall while largely preserving the inner density profile, so the concentration of a stripped subhalo of present-day mass M is not equal to c(M, z_infall). Using the present-day mass with infall redshift therefore mis-models the density profile that produces the flux-ratio perturbations in Figure 2. Since c0 is inferred from these perturbations, a systematic mis-assignment of subhalo concentrations can bias the central result. The paper does not quantify the sensitivity of the c0 constraint to the infall assumption, nor does it resolve the ambiguity between infall and present-day mass when evaluating Equation 4. The closing paragraph of the paper admits that dedicated studies are \"necessary to refine theoretical predictions for the complicated processes that can shape the density profiles of dark matter subhalos in tidal fields,\" which underscores that this effect is currently unmodeled and therefore a load-bearing systematic in the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This letter uses 11 quadruply-imaged quasars to constrain, for the first time, the normalization and logarithmic slope of the cold dark matter mass-concentration relation on sub-galactic scales (10^6 to 10^10 Msun) at cosmological distances. The analysis assumes CDM, models both dark subhalos and line-of-sight halos, uses the Sheth-Tormen mass function and a power-law c(M,z) in peak height (Eq. 4), and performs Bayesian forward modeling with a summary statistic. The main result is c0 = 12^{+6}_{-5} at 68% confidence for 10^8 Msun halos at z=0, together with translated constraints c(10^7)=15^{+9}_{-8} and c(10^9)=10^{+7}_{-4}; these are reported as consistent with several literature mass-concentration relations.","tokens_in":8365,"tokens_out":5559,"duration_ms":64931,"significance":"If the modeling assumptions hold, this is the first observational measurement of halo concentrations below 10^9 Msun at cosmological distance and demonstrates that flux-ratio statistics carry information on subhalo inner density structure, not only subhalo abundance. The comparison with theoretical c-M relations is a genuine comparison rather than a circular fit: c0 is inferred from lensing data through the population model, and the literature relations are not used to set c0. The paper is honest about the large uncertainties and uses public numerical tools, a forward model that includes finite source size, and a simultaneous treatment of subhalos and line-of-sight halos. Its main vulnerabilities are the deferred inference machinery and the uncontrolled treatment of tidal stripping of subhalos.","major_comments":[{"comment":"The subhalo rendering draws present-day subhalo masses from the mass function in Eq. (6), but Eq. (4) is then evaluated at the infall redshift. For a subhalo that has lost mass to tidal stripping since infall while approximately preserving its inner density profile, the assigned concentration does not correspond to the drawn present-day mass. Because flux-ratio perturbations are sensitive to concentration (Figure 2), this systematic misassignment can bias the inferred c0. The paper does not quantify the shift. I request a sensitivity test that replaces the infall-concentration prescription with an explicit tidal-stripping correction, or that samples infall masses before stripping, and a statement of the resulting systematic error on the headline c0 value.","section":"Section 2.2, Eq. (4) and Eq. (6)"},{"comment":"The central posterior for c0 is obtained from an approximate likelihood and a forward model whose details, including the exact definition of the summary statistic, the likelihood estimator, and validation on simulated data, are deferred to a companion paper. Since the claim that the inferred c0 is unbiased depends entirely on this machinery, the letter should either include the summary-statistic definition and a concise posterior-recovery test, or explicitly cite the validation results from the companion paper and state any known biases of the approximate likelihood.","section":"Section 2.4 and Section 2 of Gilman et al. (2019b)"},{"comment":"The prior on zeta is quoted as a Gaussian with mean -0.25 and variance 0.05, but the paper also states that zeta is unconstrained. If the data do not constrain zeta, the prior is effectively the posterior and, through the (1+z)^zeta factor, can propagate into the inferred z=0 normalization c0 for lenses at nonzero redshift. Please specify whether 0.05 is the variance or the standard deviation, and show the robustness of c0 to a wider prior on zeta.","section":"Section 2.1, Eq. (4)"}],"minor_comments":[{"comment":"The phrase 'conclusively establish strong gravitational lensing by galaxies as perhaps the only probe' combines an absolute claim with a hedged qualifier; I suggest rewording to avoid the tension.","section":"Abstract and Section 1"},{"comment":"Please state explicitly whether the scatter of 0.1 dex is applied in log10 c and whether the same scatter is assumed for subhalos, since subhalo concentrations may have a different population scatter than field halos.","section":"Section 2.1, Eq. (4)"},{"comment":"The sentence 'We assume the population mean halo mass on log10(Mhalo) of 13.3' should be rephrased, and the mass units should be stated explicitly.","section":"Section 3"},{"comment":"The caption refers to a '30 pc background source'; please specify whether this is the radius or the full width of the source and state the source redshift in the caption text.","section":"Figure 2 caption"},{"comment":"The text says the parameter zeta is unconstrained; it would be clearer to state that the posterior on zeta equals its prior, given that it is not informed by the data.","section":"Section 2.4"}],"recommendation":"major_revision","confidential_remarks":"The paper's main result is interesting and the comparison to theory is not circular, but the tidal-stripping issue in Section 2.2 is a load-bearing systematic that is acknowledged but not quantified. The editor should also confirm that the companion paper (Gilman et al. 2019b) is available to the referee, since the likelihood and validation details are essential to evaluating the c0 constraint. I would support publication after the requested sensitivity tests and a clearer statement of the systematic uncertainty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline here is that this paper does what it says: it uses 11 quadruple lenses to put the first observational constraint on the mass-concentration relation below 10^9 solar masses. The central number, c0 = 12^{+6}_{-5} for a 10^8 M_sun halo at z=0, is a direct fit to flux-ratio perturbations, and the uncertainties are reported honestly. It lands where the CDM literature predicts. That is a real result.\n\nThe method is an extension of the authors' own forward-modeling framework, so the conceptual novelty is modest, but the measurement is new and the consistency with simulations is a meaningful step, not a circular fit. They include both subhalos and line-of-sight halos, marginalize over macromodel nuisance parameters, and account for finite source size. Good.\n\nThe soft spots are real but not fatal. First, the actual likelihood and prior details live in Gilman et al. 2019b, which is not part of this submission; that makes a full assessment hard. Second, and more concerning, subhalo concentrations are evaluated at infall using the field mass-concentration relation, with no tidal-evolution prescription. The paper says so plainly, but it does not test how the inferred c0 shifts if that assumption is wrong. The stress-test note is right that this could bias the central value; it just does not demonstrate the bias is large. Given the error bars are roughly a factor of 1.7, an unquantified systematic of that order would matter. I would want a robustness test before trusting c0 as a precision number.\n\nThe Discussion's phrase 'conclusively establish' overreaches. This is a first constraint with large uncertainties, not a conclusive establishment. Minor.\n\nNo code is shipped, but the public packages used are named, so reproduction is not impossible.\n\nWho should read this: anyone working on small-scale CDM, strong lensing, or halo concentration models. It opens a new observable window, and the constraints are worth taking seriously. I would cite it. It deserves peer review; the referee should ask for a sensitivity analysis on the infall-concentration assumption, or at least an explicit statement of the expected bias direction and magnitude.\n\nSend it to review.","headline":"First observational constraint on sub-10^9 M_sun halo concentrations, with honest uncertainties, but the infall-time subhalo concentration assumption is an unquantified systematic that deserves a robustness test.","tokens_in":8886,"tokens_out":2518,"would_cite":true,"duration_ms":28243,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The first sub-galactic measurement of the dark matter mass-concentration relation finds c=12 for 10^8 solar-mass halos.","keywords":["strong gravitational lensing","mass-concentration relation","cold dark matter","dark matter subhalos","flux ratios","quadruple-image lenses","NFW profile","sub-galactic scales"],"falsifier":"Run a cosmological simulation that follows full tidal stripping of subhalos and generate mock realizations of the eleven lenses; if the inferred $c_0$ from these mocks differs from $12^{+6}_{-5}$ by more than the quoted uncertainties, the infall-only concentration treatment underlying the result is falsified.","tokens_in":7863,"feed_emoji":"🔭","tokens_out":7055,"duration_ms":69539,"temperature":0.7,"pith_summary":"Using image positions and flux ratios from eleven quadruple-image gravitational lenses, the paper reports the first measurement of the dark matter mass-concentration relation on sub-galactic scales, for halo masses $10^6$--$10^{10}\\,M_\\odot$. Within the assumption of cold dark matter, the analysis jointly constrains the normalization and logarithmic slope of the relation together with the subhalo mass function, and finds that a $10^8\\,M_\\odot$ halo at $z=0$ has concentration $c = 12^{+6}_{-5}$ at 68% confidence. The inferred relation is consistent with theoretical predictions from CDM simulations, and establishes galaxy-scale strong lensing as a direct probe of dark halo structure across cosmological distance.","feed_headline":"Halo concentration at 10^8 solar masses measured: c=12","feed_subtitle":"Eleven lensed quasars give the first direct view of dark halo structure on sub-galactic scales.","key_machinery":"The central object is the parameterized mass-concentration relation $c(M,z) = c_0(1+z)^\\zeta [\\nu(M,z)/\\nu(10^8,0)]^{-\\beta}$, written in terms of the peak height $\\nu = \\delta_c/\\sigma(R,z)$, where $\\sigma^2$ is the variance of the linear matter density field and $\\delta_c = 1.686$ is the spherical-collapse threshold. The argument is carried by forward modeling of lens flux-ratio perturbations: a more concentrated halo has a larger magnification cross section (roughly 30% perturbations for $c = 22$ versus 10% for $c = 8$), so the population of subhalos and line-of-sight halos imprints the mass-concentration relation on observed flux ratios. A Bayesian forward-modeling procedure samples the dark-matter hyper-parameters jointly with lens nuisance parameters and compares simulated lens ensembles to the eleven lenses through a summary statistic, avoiding direct evaluation of an intractable likelihood.","core_discovery":"The paper's central result is the constraint on $c_0$, the normalization of the mass-concentration relation defined as the median concentration of a $10^8\\,M_\\odot$ halo at $z=0$. At 68% (95%) confidence the authors find $c = 12^{+6}_{-5}$ ($c = 12^{+9}_{-15}$), marginalized over the subhalo mass function normalization, the line-of-sight halo abundance, the slope of the subhalo mass function, and lens nuisance parameters. Converting the hyper-parameter posterior to physical concentrations, a $10^7\\,M_\\odot$ halo has $c = 15^{+9}_{-8}$ ($15^{+11}_{-18}$ at 95%) and a $10^9\\,M_\\odot$ halo has $c = 10^{+7}_{-4}$ ($10^{+14}_{-7}$ at 95%). These values are consistent with CDM predictions from the literature and provide the first observed anchor for halo concentrations on scales where halos are expected to be mostly or entirely dark.","pith_inferences":["If this measurement is right, cold dark matter remains viable at masses where warm or self-interacting dark matter models predict modified small-scale structure; a direct comparison of this measured relation with such model predictions would sharpen the test.","The authors evaluate subhalo concentrations at infall and do not model tidal stripping; a testable extension is to implement tidal-evolution prescriptions and check whether the inferred $c_0$ shifts outside the quoted interval.","Because $c_0$ and $\\Sigma_{\\rm sub}$ are covariant, independent substructure abundance measurements from the same lenses with deeper data could break the degeneracy and tighten the concentration constraint."],"forward_implications":["If the central measurement holds, CDM's predicted halo concentrations are validated down to $10^6$--$10^{10}\\,M_\\odot$, where halos are expected to be mostly dark.","The constraint at $10^8\\,M_\\odot$ anchors the density profiles of low-mass halos used in predictions for dark matter annihilation signals and for the abundance of faint satellite galaxies.","The result establishes flux-ratio strong lensing as a cosmological probe of dark matter structure on sub-galactic scales, complementary to dynamical and stellar probes.","The reported degeneracies imply that substantially larger lens samples are needed to pin down the logarithmic slope $\\beta$ and the subhalo mass function normalization $\\Sigma_{\\rm sub}$."],"supporting_citations":[{"why":"Defines the NFW density profile and the concentration parameter $c \\equiv r_{\\rm vir}/r_s$ used throughout the analysis.","marker":"Navarro et al. (1997)"},{"why":"Supplies the forward-modeling Bayesian inference method and the subhalo and line-of-sight mass function parameterizations adopted here.","marker":"Gilman et al. (2019b)"},{"why":"Provide narrow-line flux-ratio measurements for eight of the eleven lens systems, giving sensitivity to $10^6$--$10^9\\,M_\\odot$ halos.","marker":"Nierenberg et al. (2014, 2017, 2019)"},{"why":"Provides the radio flux-ratio data for the lens B0128+437.","marker":"Koopmans et al. (2003)"},{"why":"Provides the radio flux-ratio data for the lens MG0414+0543.","marker":"Katz et al. (1997)"},{"why":"Provides the radio flux-ratio data for the lens PG 1115+080.","marker":"Chiba et al. (2005)"},{"why":"The galacticus semi-analytic model supplies the infall-redshift distribution and host-halo scaling used to render subhalo concentrations and abundances.","marker":"Benson (2012)"},{"why":"The Sheth-Tormen mass function is the baseline for the line-of-sight halo abundance in Equation 3.","marker":"Sheth et al. (2001)"},{"why":"Provides the 0.1 dex scatter in the mass-concentration relation adopted for halo populations.","marker":"Dutton & Macciò (2014)"}],"fun_headline_variants":["First mass-concentration link from 11 lensed quasars","Sub-galactic halo concentrations measured for first time","Dark halo structure probed with 11 gravitational lenses","Lensed quasars reveal halo concentration at 10^8 solar masses","First constraints on sub-galactic halo concentration from lensing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on assuming that each subhalo's density profile is set by the field mass-concentration relation at infall and is not subsequently reshaped by tidal stripping; if tides substantially alter subhalo inner profiles, the inferred concentration normalization would be biased.","fun_headline_variants_meta":{"raw":{"variants":["First mass-concentration link from 11 lensed quasars","Sub-galactic halo concentrations measured for first time","Dark halo structure probed with 11 gravitational lenses","Lensed quasars reveal halo concentration at 10^8 solar masses","First constraints on sub-galactic halo concentration from lensing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000646,"raw_usage":{"total_tokens":3041,"prompt_tokens":1088,"completion_tokens":1953,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":1871}},"tokens_in":704,"tokens_out":1953,"duration_ms":13831,"temperature":1.0,"reasoning_tokens":1871,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:46:29.805635+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a cosmological simulation that follows full tidal stripping of subhalos and generate mock realizations of the eleven lenses; if the inferred $c_0$ from these mocks differs from $12^{+6}_{-5}$ by more than the quoted uncertainties, the infall-only concentration treatment underlying the result is falsified.","supporting_citations":[{"cited_title":"A., Moore C","cited_arxiv_id":null,"evidence_quote":"Provides the radio flux-ratio data for the lens MG0414+0543."}],"review_version":1}