{"id":"3f7712ec-3d68-4727-af1b-f7f7ec3bd3f2","arxiv_id":"1908.05363","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The ratio of DLA absorption velocity width to host galaxy velocity dispersion decreases with impact parameter, and steep dark matter halo density profiles fit the data better than shallow ones, though no model can be excluded.","lead":"This paper uses the widths of quasar absorption lines and the emission lines of host galaxies to measure how fast gas moves at different distances from galaxies across cosmic time. The authors argue these motions trace dark matter halos, and that steeper dark matter density profiles fit the data better, though the scatter is large.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed steep-cusp preference is not yet established: it relies on treating Δv90 as σlos and on per-model vertical normalization, so the model ranking in Table 2 may not be a DM-profile measurement.","rationale":"The reader's weakest assumption—that Δv90 measures the same quantity as σlos—is indeed the load-bearing link, so I agree with that diagnosis. I would sharpen it: the manuscript calibrates the identification only at the median level (Section 3) and then compares radial shapes, for which no independent calibration exists. The paper's own Sections 6.2 and 6.3 limitations (outflows, groups, single-cloud contamination; \"we assume that the velocities represented by Δv90 is a measure of the projected velocity dispersion\") mark this as an assumption rather than a demonstrated result. The free vertical normalization of each model in Figures 1–2 compounds the problem, because a model with the wrong shape can still score well once scaled; all reduced chi-square values are far above unity, so the relative rankings are not a statistically strong preference. A simulation-based mock-spectrum test can settle the mapping question, and removing the free normalization would settle the shape-comparison question. Because the paper's own text (\"does not allow us to rule out any of the models\") is already more cautious than the abstract, the reader's CONDITIONAL verdict correctly captures the state of evidence; my stress-test does not move it. I give credit for the honest caveats and for the independent KS comparison with simulations in Section 5, although that comparison lacks impact parameters and redshift matching.","tokens_in":18709,"tokens_out":7000,"duration_ms":79589,"concrete_test":"Generate mock DLA spectra in a cosmological hydrodynamical simulation with known halo potentials (e.g., Bird et al. 2015 or TNG50), place sightlines at known impact parameters, measure Δv90 exactly as done for the real data, and compare to the true σlos computed from the same halos via Equations 1–3. If the median Δv90/σlos deviates from unity, or if that ratio depends on impact parameter, star-formation rate, or outflow velocity, the mapping used for Figures 1–2 is invalid; the model ranking in Table 2 would then not constrain the dark-matter density profile.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 sets the central identification, \"it is a good approximation to set Δv90 ≈ σlos\", but the calibration is only a median-level check (σlos ~ 100–130 km/s for log Mhalo = 11.7 at r = rs versus a median Δv90 = 141 km/s). The subsequent use in Figures 1–2 requires not just the median but a shape correspondence: the 90-percent velocity width of a handful of discrete absorbing clouds must track the isotropic Jeans σlos(R) of a smooth spherical halo at every impact parameter. The manuscript itself flags exactly the failure modes: outflows, group membership, and a single dominant cloud altering Δv90 are discussed in Section 6.2, and Section 6.3 states that the equivalence is an assumption. Additionally, every model curve in Figures 1 and 2 is normalized to minimize chi-square (see Section 3 and Table 2), so the comparison is of shapes only, not absolute predictions. With reduced chi-square values of 127–161 for all models in the scaled comparison, none of the models is a good fit, and the small ranking differences could be produced by the normalization or by a few high-ratio points at small b/rs. If Δv90 is contaminated by outflows or unrelated components, or if the normalization absorbs the model dependence, the conclusion that Dehnen/Jaffe/Einasto profiles fit \"better\" is not a robust measurement of inner dark-matter slopes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compiles a sample of 21 damped Lyman-alpha absorbers (DLAs) with spectroscopically confirmed host galaxies over 0.2 < z < 3.2, and studies how the ratio of the DLA velocity width Δv90 to the host galaxy emission-line velocity dispersion σem varies with projected impact parameter. The authors report that Δv90/σem decreases with impact parameter, and they compare the data with line-of-sight velocity dispersion profiles σlos(R) computed from a suite of dark matter halo density profiles (NFW, isothermal, Hernquist, Jaffe, Plummer, Dehnen, Einasto). The models are vertically normalized to minimize χ², and the paper concludes that steeper inner profiles (Dehnen, Jaffe, Einasto) fit the data better than shallower ones (NFW, isothermal), while cautioning that the scatter is large. The paper also compares Δv90 to halo virial velocities in numerical simulations, finding good agreement with a KS test, and argues that most DLAs are gravitationally bound to their host halos by comparing relative velocities to escape velocities.","tokens_in":19028,"tokens_out":6153,"duration_ms":59805,"significance":"If the central claim holds, the paper would provide a rare observational constraint on dark matter density profiles at intermediate redshifts using absorption-selected galaxies, complementing local kinematic studies. The authors deserve credit for assembling a unique sample of DLA host galaxies with both absorption and emission kinematics, for testing a broad range of standard halo profiles within a simple Jeans framework, and for honestly enumerating the assumptions and failure modes in Sections 6.2 and 6.3. The comparison with the Bird et al. (2015) simulations and the escape-velocity analysis are useful complementary checks. However, the headline conclusion about steeper dark matter profiles is currently not quantitatively established, because the model comparison is a shape-only test with heavily normalized and uniformly poor fits, and because the identification Δv90 ≈ σlos is validated only at the median level rather than as a point-to-point relation.","major_comments":[{"comment":"The model ranking in Table 2 is a shape-only comparison because every model curve is vertically normalized to minimize χ², as stated in Section 3 and the caption of Figure 1. The reduced χ² values for the scaled comparison in Table 2 range from 127 to 161, meaning no model is a good fit, and the differences between profiles are small relative to the overall badness of fit. The abstract's claim that models with steeper radial profiles 'provide a better fit' is therefore not established as a measurement of the dark matter inner slope; it is a statement about which normalized shape is least bad. The authors should fix the normalization using the absolute σlos calibration from Section 3, or apply a formal model comparison that accounts for the fitted normalization, and they should temper the abstract in line with Section 6.1's own admission that 'the currently known data sample does not allow us to rule out any of the models.'","section":"Section 3 / Table 2"},{"comment":"The central identification Δv90 ≈ σlos is validated only at the median level: σlos ~ 100–130 km/s for the median halo mass at r = rs versus a median Δv90 = 141 km/s. The shape comparison in Figures 1 and 2 requires Δv90(R) to track σlos(R) point by point, but Section 6.3 states that this is an assumption, and Section 6.2 lists mechanisms (outflows, group membership, a single dominant component) that can alter Δv90 without changing σlos. Because the dark-matter profile ranking depends entirely on this mapping, the authors should provide a direct test of the shape correspondence, for example by measuring Δv90 and σlos(R) in the Bird et al. (2015) simulations with known impact parameters, or by demonstrating that Δv90 scaled by host properties correlates with the expected σlos at the measured impact parameter.","section":"Section 3 / Section 6.3"},{"comment":"The paper asserts that Δv90/σem decreases with impact parameter, but no quantitative significance is reported. Please provide a rank correlation coefficient (e.g., Spearman) with a p-value, and ideally a fit of log(Δv90/σem) versus log(b) with uncertainties. The large scatter and the leverage of a few high-ratio points at small impact parameters (e.g., DLA0918+1636 at b = 2.0 kpc with Δv90/σem ≈ 16) make such a test necessary to establish the observational foundation for the subsequent model comparison.","section":"Section 2 / Figure 1"},{"comment":"The x-axis scaling of the data in Figure 2 depends on halo masses and scale radii derived from abundance matching and concentration–mass relations. Section 6.2 notes that individual halo masses can be underpredicted by up to a factor of 10, which shifts log(b/rs) by about –0.4 dex. This systematic uncertainty is not propagated into the χ² values in Table 2 and could affect the relative ranking of models. The authors should show how the ranking changes when rs is varied within the abundance-matching scatter, or at least discuss the sensitivity of the Table 2 values to these systematic shifts.","section":"Section 4.1 / Section 6.2"}],"minor_comments":[{"comment":"Several systems lack σem measurements (e.g., 0439–433, 0738+313, 1127–145). Please clarify explicitly which systems enter the Δv90/σem analysis and which enter the virial-velocity comparison, and state the sample size used for each figure.","section":"Table 1"},{"comment":"The data points are plotted without error bars. Please either add representative error bars on Δv90, σem, and impact parameter, or state in the captions the typical uncertainties and why they are omitted.","section":"Figures 1–3"},{"comment":"The phrase 'computed from Equation 3 to represent the tangential velocity' appears to refer to the wrong equation, since Equation (3) is the radial velocity dispersion from the Jeans equation; please correct the cross-reference.","section":"Section 5.2"},{"comment":"The text states that outflows increase both Δv90 and σem, which would partly cancel in the ratio Δv90/σem; the abstract's statement that outflows 'may cause an increased scatter' should be reconciled with this cancellation or made more precise.","section":"Section 6.2"},{"comment":"Equations (14) and (15) give different scaling exponents for rs with halo mass and redshift; please state explicitly which concentration–mass relation is used to compute the scale radii in Table 1 and which variant is used in each panel of Figure 2, to avoid ambiguity.","section":"Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript depends on the companion paper Møller & Christensen (2019, submitted) for the sample compilation and for the claimed scaling relation. Given that this paper is not yet published, the editors may wish to ask the authors to include the essential sample details here or to update the reference status. The paper's scope is appropriate for MNRAS, but the abstract overstates the strength of the dark-matter profile conclusion relative to the quantitative analysis presented; the revisions listed in the major comments should bring the claims in line with the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper is worth reading for the empirical Δv90/σem–impact parameter trend and for the honest presentation of its caveats, but the claimed preference for steep dark-matter profiles is not established. The model comparison in Fig. 2 and Table 2 normalizes every profile to minimize χ², so it only tests shape, and the reduced χ² values (127–161) mean no model actually fits. The ranking between Dehnen/Jaffe/Einasto and NFW/isothermal is a ranking of shapes with a fitted vertical offset, not a clean measurement of inner slope.\n\nWhat is genuinely new: the radial scaling relation itself, introduced in the companion paper, is now confronted with analytic σ_los(R) from a suite of halo profiles and with virial velocities from simulations. That is a sensible next step, and the KS test against Bird et al. (P=0.98) is encouraging, though the authors themselves note it could be coincidence given their metal-rich, massive-host selection. The escape-velocity argument that most DLA gas is bound to the host halo is straightforward and useful.\n\nThe soft spots are the ones the stress-test flags. Section 3 asserts Δv90 ≈ σlos based on a median comparison (σlos ~ 100–130 versus median Δv90 = 141 km/s). That gives a rough normalization check, but the figures use σlos(R) as a per-impact-parameter prediction. For that to work, the 90-percent velocity width of a handful of clouds must track the isotropic Jeans dispersion of a smooth spherical halo at every radius. The paper itself lists outflows, groups, and single dominant clouds as failure modes (Sec. 6.2), and admits the equivalence is an assumption (Sec. 6.3). With 21 points and large scatter, the steep-profile preference in the abstract is stronger than the text's own \"cannot rule out any model.\"\n\nAlso, reliance on Møller & Christensen (2019, submitted) for sample construction means key details are not independently checkable from this paper alone. That is not fatal, but it matters for a claim this subtle.\n\nWho it is for: people working on DLA-host connections and galaxy-halo scaling relations. The empirical trend is a useful constraint even if the DM-slope conclusion is provisional. A serious referee should engage because the idea is new and testable; the authors list the right next steps (per-object error bars, normalization-free predictions, simulations with impact parameters). I would accept it for review with a request to temper the abstract and either add error bars or drop the model ranking.","headline":"Interesting first attempt to turn DLA line widths into a dark-matter profile probe, but the steep-cusp preference rests on a loose Δv90≈σlos identification and per-model normalization, so the model ranking is not yet a measurement.","tokens_in":19596,"tokens_out":3272,"would_cite":true,"duration_ms":26886,"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":"Damped Lyman-α absorbers trace their host halos' dark matter potentials, with velocity-width ratios favoring steep inner density profiles over shallow ones.","keywords":["Damped Lyman-alpha absorbers","quasar absorption lines","dark matter halo profiles","galaxy velocity dispersion","impact parameter","circumgalactic medium","high-redshift galaxies","halo abundance matching"],"falsifier":"Measure $\\Delta v_{90}$ for a DLA with several quasar sightlines through the same halo, for example a gravitationally lensed quasar or a close quasar pair, and check whether the per-sightline widths follow the $\\sigma_{\\rm los}(R)$ curve of a single spherical halo at the corresponding impact parameters; if the widths are set by outflow or unbound gas, they will scatter independently of $R$ and the steep-profile preference will disappear. A more direct test is to compare a DLA host with a measured rotation curve at $z\\sim0.7$: the observed $\\Delta v_{90}$ at the impact parameter should equal the $\\sigma_{\\rm los}$ computed from the rotation curve, not exceed it by the outflow velocity scale.","tokens_in":18521,"feed_emoji":"🔭","tokens_out":11899,"duration_ms":95895,"temperature":0.7,"pith_summary":"This paper tries to establish that the widths of quasar absorption lines in damped Lyman-α systems (DLAs) trace the gravitational potential of the host galaxy's dark matter halo along a single line of sight. Normalizing the absorption velocity width $\\Delta v_{90}$ by the host galaxy's emission-line velocity dispersion $\\sigma_{\\rm em}$ removes the stellar-mass dependence, so the ratio can be compared with the projected line-of-sight velocity dispersions $\\sigma_{\\rm los}(R)$ predicted by different dark matter halo density profiles. The paper reports that $\\Delta v_{90}/\\sigma_{\\rm em}$ decreases with impact parameter and is better matched by models with steep inner density slopes (Dehnen, Jaffe, Einasto) than by shallower profiles such as NFW and isothermal, although the scatter is large. If correct, the result turns absorption-line surveys into a redshift-resolved probe of the inner structure of dark matter halos, and it implies that most DLA gas remains gravitationally bound to its host halo.","feed_headline":"Absorbers reveal steep dark matter halos","feed_subtitle":"Normalized DLA velocity widths favor cuspy halo profiles over NFW across redshifts 0.2 to 3.2.","key_machinery":"The central object is the mapping $\\Delta v_{90}\\approx\\sigma_{\\rm los}(R)$, where the absorption-line velocity width is treated as the line-of-sight velocity dispersion of a spherical dark matter halo at the projected impact parameter $R$. The machine starts with the Jeans equation and expresses $\\sigma_{\\rm los}(R)$ as an integral over the halo density profile, parameterized as a double power law $\\rho(r)=\\rho_0\\,(r/r_s)^{-\\gamma}[1+(r/r_s)^\\alpha]^{-(\\beta-\\gamma)/\\alpha}$, whose special cases are the Hernquist, Jaffe, Plummer, NFW, and isothermal profiles, plus the Dehnen and Einasto families. The comparison is placed on a common footing by normalizing each DLA's $\\Delta v_{90}$ by its host's $\\sigma_{\\rm em}$ and each impact parameter by the halo scale radius $r_s$ computed from abundance-matching halo masses, so that galaxies spanning nearly three orders of magnitude in stellar mass can be compared with one predicted profile at a time. The ratio $\\Delta v_{90}/V_{\\rm vir}$ with virial velocities computed from $V_{\\rm vir}^3=10GM_{\\rm halo}H(z)$ provides a second, simulation-based comparison.","core_discovery":"Using a sample of 21 DLAs at $0.2