{"id":"6286e53f-d327-49fd-8e26-aa45dcd241b9","arxiv_id":"2501.09077","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Dark matter constitutes 53 to 61 percent of the mass within the effective radii of two z~6 quasar hosts, implying very massive dark matter halos near 10^12.5 solar masses.","lead":"Using ALMA observations of glowing carbon gas, astronomers traced the rotation of two ancient quasar host galaxies out to nearly 8 kiloparsecs. They find dark matter makes up about 53 to 61 percent of the mass inside the central regions, implying these early galaxies sit in very heavy dark matter halos.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted fDM(R<Re) uses the gas half-mass radius, but the lower-redshift comparison uses the stellar effective radius; re-evaluating at the stellar Re could erase the 'DM dominated' claim.","rationale":"I read the paper in good faith and credit its genuine strengths: simultaneous fitting of high- and low-resolution data, the JvM correction for non-Gaussian beams, cross-checking with 3DBarolo, and Monte-Carlo variation of fixed parameters. The concentration concern identified by the Reader is real, but it mostly affects the extrapolated halo mass and not the measured rotation-curve decomposition. The more load-bearing issue is the definition of Re in fDM. The abstract and Figure 7 imply a comparison to lower-redshift fDM(R<Re) determinations, which use stellar effective radii, while the present analysis uses the gas half-mass radius. This inflates fDM relative to the comparison standard and directly bears on the central claim that these z~6 quasars are dark-matter dominated and significantly different from lower-redshift systems. The issue is concrete, testable with the existing best-fit models, and not a matter of outside consensus. Therefore the appropriate verdict remains CONDITIONAL: the kinematic analysis is a real step forward, but the headline dark-matter fractions need to be re-derived at the standard stellar effective radius before the low-redshift comparison can be accepted.","tokens_in":39326,"tokens_out":6059,"duration_ms":69492,"concrete_test":"Using the best-fit DysmalPy mass profiles from Table 3, compute the enclosed DM fraction at R = R_star for both objects: fDM(R_star) = M_DM(<R_star) / [M_DM(<R_star) + M_baryon(<R_star)]. If either value drops below 0.5, the claim that both quasars are dark-matter dominated within their effective radius is not supported. As a second check, recompute the Fig. 7 comparison using the same Re definition (stellar effective radius) as Nestor Shachar et al. (2023) and quantify the resulting shift. This requires no new ALMA data.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central comparison in Fig. 7 compares fDM(R<Re) = 0.61 and 0.53 with z~2 measurements where Re is the stellar effective radius. In this paper, Table 3 defines Re as the half-mass radius of the gas component. The fits give R_gas = 2.77 kpc (P009-10) and 2.57 kpc (J2318-3029), while the fitted stellar effective radii are 1.96 kpc and 1.30 kpc. Because the dark-matter fraction increases with radius, fDM at the stellar Re is systematically lower than the quoted values. For J2318-3029, with R_star about half of R_gas, fDM(R_star) could fall below 0.5, in which case the headline that both systems are dark-matter dominated would not hold, and the statement that these fractions are 'significantly larger than those extrapolated from lower redshift studies' would be weakened. The halo masses, which are derived from the DM mass enclosed at this same gas radius, inherit the same offset. This concern is distinct from the concentration issue in Sec. 5.2: fixing c=3.5 affects the halo mass extrapolation by up to ~2 dex but does not change fDM; the radius definition directly changes the headline fDM values themselves.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents ALMA [C II] observations of two z~6 quasar host galaxies (P009-10 and J2318-3029), combining high- and low-resolution data to recover extended emission. Using two independent modeling tools (DysmalPy and 3DBarolo), the authors derive rotation curves out to ~6-8 kpc and decompose the mass into baryonic and dark matter components. They report dark matter fractions fDM(R<Re)=0.61(+0.08,-0.08) and 0.53(+0.20,-0.23), infer halo masses log Mh~12.5-12.8 Msun, and argue that these quasars reside in the most massive halos at their epoch, with SMBH masses aligned with a local MBH-Mh relation rather than the MBH-Mstar relation.","tokens_in":39686,"tokens_out":3353,"duration_ms":34792,"significance":"If the central results hold, this would be one of the first direct dynamical measurements of dark matter content in quasar host galaxies at z~6, providing unique constraints on early structure formation and SMBH-host co-evolution. The paper has notable strengths: it makes use of both compact and extended ALMA configurations, applies a JvM beam correction, cross-checks the kinematics with two independent modeling codes, and includes Monte Carlo robustness tests showing that fDM(R<Re) is largely stable against variations in Sersic index, axis ratio, and halo concentration. However, the headline comparison with lower-redshift studies is undermined by an inconsistency in the definition of effective radius, and the halo mass inference depends strongly on a fixed NFW concentration. These issues affect the paper's main astrophysical claims, though they appear addressable with additional analysis.","major_comments":[{"comment":"The quoted fDM(R<Re) uses the gas half-mass radius, as stated in the note to Table 3 ('Re denotes the half mass radius of the gas component'), whereas the lower-redshift comparison in Fig. 7 (e.g., Nestor Shachar et al. 2023) is based on the stellar effective radius. The fitted values are Rgas=2.77 kpc and Rstar=1.96 kpc for P009-10, and Rgas=2.57 kpc and Rstar=1.30 kpc for J2318-3029. Since the dark matter fraction increases with radius, the values at the stellar effective radius will be lower than the quoted numbers; for J2318-3029, where Rstar is roughly half of Rgas, fDM(Rstar) could fall below 0.5. This would weaken or invalidate the claim that both systems are dark-matter-dominated and that the fractions are significantly larger than the lower-redshift extrapolation. Please recompute fDM at the stellar effective radius (or provide a radial profile of fDM) and restrict the comparison in Fig. 7 to a consistent definition of Re.","section":"Table 3 and Fig. 7"},{"comment":"The virial halo masses quoted in Table 3 and the abstract (log Mh~10^12.5-10^12.8 Msun) are not direct measurements but are extrapolated from the fitted fDM(R<Re) under the assumption of an NFW profile with a fixed concentration c=3.5. Fig. 9 shows that log Mh varies by roughly 2 dex when c is allowed to range from 1.5 to 5.5, and the text acknowledges this. The conclusion that these quasars reside in 'the most massive halos at these redshifts' is therefore contingent on the adopted concentration and is not a robust inference from the data alone. I recommend marginalizing over c with a physically motivated prior (e.g., from Dutton & Maccio 2014 or Diemer & Kravtsov 2015) or, at minimum, presenting the halo mass as a function of c and softening the abstract/conclusion statements accordingly.","section":"Sec. 5.2 and Fig. 9"},{"comment":"The filtering of 'non-circular' components in P009-10 is a potentially large intervention: residual Gaussian components are subtracted from the data cube before the kinematic fitting. The assumption that these components are non-gravitational and can be removed is not directly tested. If the removed emission actually traces part of the gravitational potential (e.g., a merger or an infalling clump), the derived fDM(R<Re) could be biased. The virial-theorem consistency check in Sec. 5.7 indicates that the total dynamical mass within ~3 kpc is consistent with the rotating-disk model, but it does not validate the dark matter fraction itself. Please quantify how the filtering changes the fitted fDM (for example, by running the DysmalPy fit on the unfiltered cube and comparing the results), or provide additional justification that the removed components do not affect the mass decomposition.","section":"Sec. 3.1"}],"minor_comments":[{"comment":"The abstract and Sec. 4.5 state fDM(R<Re) without specifying that Re is the gas half-mass radius; please define this explicitly in both places, and ideally add a note whenever referring to 'effective radius' to avoid confusion with the stellar effective radius used in the literature comparison.","section":"Abstract and Sec. 4.5"},{"comment":"The caption should state which effective radius is used for each dataset (stellar for lower-redshift points, gas for the present work) so that the reader can immediately see the potential inconsistency.","section":"Fig. 7 caption"},{"comment":"The dark matter fraction fDM(R<Re) is a fitted parameter with a flat prior [0,1], and the halo mass is derived from it; phrases such as 'we find' or 'the dynamic measurements indicate' (e.g., Sec. 4.5) could be phrased as 'we fit' to avoid implying an independent measurement.","section":"Sec. 3.2"},{"comment":"The paper acknowledges in Sec. 5.7 that 'we cannot entirely rule out the possibility that the observed kinematics are affected by these effects,' yet the Conclusions present the dark-matter-dominated result without this caveat; please carry the caveat through to the summary.","section":"Sec. 5.7 and Conclusions"},{"comment":"There is a typo: 'hugh fraction' should be 'huge fraction'.","section":"Sec. 5.2"}],"recommendation":"major_revision","confidential_remarks":"The radius-definition issue is the most serious technical point: it directly affects the headline fDM values and the comparison with the literature. It is fixable if the authors can recompute fDM at the stellar effective radius or present a radial profile, but it needs to be done before publication. The concentration dependence of the halo mass is acknowledged but should be moved more prominently into the abstract and conclusions, since the 'most massive halo' claim is a key selling point. The paper is otherwise solid in its data handling and kinematic modeling; I would be willing to see a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you the short version. This paper combines high- and low-resolution ALMA data on two z>6 quasar hosts, fits [CII] kinematics with two independent codes, and derives dark matter fractions within the effective radius. The data work is genuinely good: the JvM correction is handled carefully, the modeling with DysmalPy and 3DBarolo agree, and they show the fDM result is stable under several prior choices. This is the first time anyone has done a mass decomposition at z~6 from gas kinematics, so it's a real step forward.\n\nThe problem is the definition of Re. Table 3 defines Re as the half-mass radius of the gas component, but the comparison points from Nestor Shachar et al. and Rizzo et al. use the stellar effective radius. Your stress-test is right: for these two galaxies, the gas Re is roughly twice the stellar Re, so the quoted fDM values of 0.61 and 0.53 are measured at larger radii than the lower-z values. Because fDM rises with radius, the comparison in Figure 7 is not apples-to-apples. Re-evaluating at the stellar Re could easily push J2318-3029 below 0.5, which would kill the 'DM-dominated' headline for that object. The halo masses inherit the same offset.\n\nThe other soft spots are less severe but worth noting. The halo masses rely on a fixed NFW concentration c=3.5, and their own Figure 9 shows log Mh changes by ~2 dex across c=1.5 to 5.5. So the claim that these are the most massive halos at z~6 is only as good as that assumption. Also, the sample is two out of five; one of the two needed the non-circular component removed before fitting. That is a selection effect they acknowledge, but it means the result may not generalize.\n\nNone of this is fatal to the paper's core value. The rotation curves and the gas-phase mass decomposition are solid, and the qualitative conclusion that these systems have substantial dark matter at a few kpc seems robust. What needs fixing is the radius definition, a consistent comparison, and a more honest error budget on the halo masses. If the authors do that, this becomes a useful anchor for high-z galaxy evolution.\n\nI'd send it to a good referee. A serious referee will catch the Re issue, so the authors should address it in revision. It deserves review, not rejection, because the measurements are new and the analysis is careful.\n\nI'd probably cite the rotation curves even now, but I'd wait to cite the fDM values until they are recomputed at a consistent radius. For a reading group, it's a good case study in how subtle definitions can change a headline claim.","headline":"First z~6 fDM estimates from gas kinematics are worth taking seriously, but the radius definition mismatch with the comparison sample means the headline 'DM dominated' claim is not yet supported.","tokens_in":40268,"tokens_out":5057,"would_cite":true,"duration_ms":47223,"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":"Two quasars at redshift six appear to live in dark-matter-dominated host galaxies, with dark matter fractions near 0.6 and 0.5 within the effective radius.","keywords":["quasar host galaxies","dark matter fraction","rotation curves","[C II] kinematics","z~6 quasars","ALMA","supermassive black holes","galaxy evolution"],"falsifier":"A JWST measurement of the stellar light of P009-10 and J2318-3029 that puts their stellar masses near $10^{11}\\,M_\\odot$ rather than the fitted $\\sim10^{10.5}\\,M_\\odot$ would lower the inferred dark matter fractions below the quoted values and falsify the claim that these hosts are dark-matter dominated. An independent constraint on the halo concentration from the full rotation curve shape would also settle whether the virial masses are actually $\\gtrsim10^{12.5}\\,M_\\odot$.","tokens_in":39108,"feed_emoji":"🌌","tokens_out":6285,"duration_ms":58586,"temperature":0.7,"pith_summary":"This paper uses ALMA observations of the [C ii] 158 micron line to measure gas rotation in two quasar host galaxies at $z\\gtrsim6$, combining high- and low-resolution data to trace the velocity field from the inner galaxy out to roughly 8 kiloparsecs. The authors argue that both systems are rotating disks with $V_{\\rm rot}/\\sigma\\approx2$, and that decomposing the rotation curves into stars, gas, and a dark halo yields dark matter fractions within the effective radius of $f_{\\rm DM}(R<R_e)=0.61^{+0.08}_{-0.08}$ and $0.53^{+0.21}_{-0.23}$. These fractions are about twice what lower-redshift trends would predict, which the paper reads as evidence that these first quasars live in unusually massive dark matter halos, around $10^{12.5}$ to $10^{12.8}\\,M_\\odot$. If correct, the measurement would tie the early growth of supermassive black holes to their dark matter halos and provide a direct kinematic anchor for dark matter assembly in the first billion years.","feed_headline":"Two z~6 quasar hosts are dark-matter dominated","feed_subtitle":"ALMA rotation curves give dark matter fractions of 0.6 and 0.5, pointing to the most massive halos of that era.","key_machinery":"The central object is the rotation curve of the [C ii]-emitting gas, built from ALMA observations at two angular resolutions. The machinery is forward modeling with DysmalPy, which constructs a three-component mass model (stellar bulge, gaseous disk, and an NFW dark matter halo; NFW is the Navarro-Frenk-White density profile) and predicts the observed data cube, plus 3DBarolo, a tilted-ring model that recovers rotation velocities non-parametrically. The key step is that the low-resolution data reach radii of 6-8 kpc where the dark halo dominates the circular velocity, so the dark matter fraction inside the effective radius can be determined rather than treated as a nuisance parameter.","core_discovery":"On the paper's own terms, the discovery is that two $z\\sim6$ quasar hosts are dark-matter-dominated systems on kiloparsec scales. Using [C ii] emission as a dynamical tracer, the authors recover extended rotation curves that stay flat or keep rising out to 6-8 kpc, well beyond the likely stellar distribution, in contrast to the declining curves found for many massive star-forming galaxies at $z\\sim2$. Forward modeling of the mass distribution gives $f_{\\rm DM}(R<R_e)=0.61^{+0.08}_{-0.08}$ for P009-10 and $0.53^{+0.21}_{-0.23}$ for J2318-3029, with inferred halo masses of $\\sim10^{12.85}$ and $\\sim10^{12.50}\\,M_\\odot$ under an NFW profile with concentration 3.5. The authors show these values are stable under changes in S\\'ersic indices, disk thicknesses, and inclination, and that adding low-resolution data is what breaks the degeneracy between baryonic mass and dark matter fraction. They further find that the black hole masses, while roughly ten times above the local $M_{\\rm BH}$-stellar-mass relation, sit closer to the local $M_{\\rm BH}$-halo-mass relation, suggesting halo mass rather than stellar mass may set the scale for the first supermassive black holes.","pith_inferences":["If the fixed NFW concentration of 3.5 is relaxed, the paper's own tests show the virial halo mass shifts by roughly two orders of magnitude, so the 'most massive halo' conclusion is weaker than the dark matter fraction measurement itself.","The same high-plus-low-resolution [C ii] approach could be applied to non-quasar galaxies at $z>6$, testing whether these high dark matter fractions are intrinsic to massive halos or a selection effect of quasar environments.","Direct JWST imaging of the host starlight would replace the assumed stellar S\\'ersic index and thickness with measured values, turning the quoted uncertainties on $f_{\\rm DM}(R<R_e)$ into a sharper test.","A larger sample spanning a range of quasar luminosity would show whether the alignment with the local $M_{\\rm BH}$-$M_{\\rm halo}$ relation is a genuine evolutionary link or a consequence of selecting the most luminous systems."],"forward_implications":["The dark matter fraction at $z\\sim6$ does not follow the extrapolated decline from cosmic noon; massive quasar hosts can be dark-matter dominated already at this epoch.","Halo masses of $\\sim10^{12.5}$-$10^{12.8}\\,M_\\odot$ make these quasars tracers of the most massive halos at $z\\sim6$, testable through galaxy overdensity and clustering.","The offset from the local black hole-stellar mass relation combined with the proximity to the black hole-halo mass relation suggests the halo, not the stellar bulge, is the regulating reservoir for early black hole growth.","Deep low-resolution ALMA data are essential: with high-resolution data alone, the dark matter fraction of J2318-3029 is essentially unconstrained.","Flattened or rising rotation curves at large radius imply that pressure support must be removed via asymmetric drift correction before computing circular velocities."],"supporting_citations":[{"why":"Supplies the mass-concentration-redshift relation used to fix the halo concentration at c=3.5 and to interpret the inferred halo mass.","marker":"Dutton & Macciò 2014"},{"why":"Provides the DysmalPy forward-modeling machinery for fitting multi-component mass models to kinematic data.","marker":"Price et al. 2021"},{"why":"Provides 3DBarolo, the non-parametric tilted-ring code used to cross-check the rotation curves.","marker":"Di Teodoro & Fraternali 2015"},{"why":"Establishes the lower-redshift comparison for fDM(R<Re) and the asymmetric drift treatment for pressure-supported disks.","marker":"Genzel et al. 2020"},{"why":"Gives the z~2 dark matter fraction trend from which the paper's values deviate by a factor of two.","marker":"Nestor Shachar et al. 2023"},{"why":"Supplies the black hole masses and circular velocities of z~6 quasars used for the MBH-halo comparison.","marker":"Neeleman et al. 2021"},{"why":"Provides z~4 dusty star-forming galaxy dark matter fractions and rotation curve methods used as comparison.","marker":"Rizzo et al. 2021"},{"why":"Provides the semi-analytic Vrot/sigma redshift evolution and Toomre-stability prediction the paper compares against.","marker":"Wisnioski et al. 2015"},{"why":"Describes the JvM correction applied to the ALMA cubes, which recovers the diffuse extended emission used for the outer rotation curves.","marker":"Czekala et al. 2021"}],"fun_headline_variants":["Dark matter dominates two z~6 quasar hosts","z~6 quasar hosts are dark-matter dominated","Early quasar galaxies found dark-matter heavy","Two quasar hosts at z~6 have massive dark halos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the dark matter halo follows a standard NFW density profile with a concentration fixed to 3.5; if the real concentration is higher or lower, the inferred total halo mass changes by up to two orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter dominates two z~6 quasar hosts","z~6 quasar hosts are dark-matter dominated","Early quasar galaxies found dark-matter heavy","Two quasar hosts at z~6 have massive dark halos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000428,"raw_usage":{"total_tokens":2315,"prompt_tokens":1195,"completion_tokens":1120,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":811,"completion_tokens_details":{"reasoning_tokens":1054}},"tokens_in":811,"tokens_out":1120,"duration_ms":10454,"temperature":1.0,"reasoning_tokens":1054,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:10:34.378165+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A JWST measurement of the stellar light of P009-10 and J2318-3029 that puts their stellar masses near $10^{11}\\,M_\\odot$ rather than the fitted $\\sim10^{10.5}\\,M_\\odot$ would lower the inferred dark matter fractions below the quoted values and falsify the claim that these hosts are dark-matter dominated. An independent constraint on the halo concentration from the full rotation curve shape would also settle whether the virial masses are actually $\\gtrsim10^{12.5}\\,M_\\odot$.","supporting_citations":[],"review_version":1}