{"id":"3fbb2b4e-b364-4235-a483-5fdc11406f20","arxiv_id":"2411.08958","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"New ALMA observations of a z about 2.6 AGN-host galaxy reveal a cold, rotation-supported molecular disk with a flat rotation curve, several non-circular gas structures, and dynamical masses that disagree with photometric estimates.","lead":"Astronomers mapped the cold gas in PKS 0529-549, a distant galaxy hosting a radio-loud black hole, using ALMA. They found a surprisingly orderly rotating disk, complete with gas tails and an anomalous central structure, and they show that standard mass estimates for such galaxies may be off by factors of several.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Kinematic center fixed by hand and Vrad=0 make disk parameters and mass limits conditional on non-circular contamination being negligible, which the current fits do not test.","rationale":"After reading the paper, the main load-bearing assertion is that one can extract a clean rotation curve and velocity dispersion from a system that the authors themselves show contains a high-velocity anomalous component and two tails. The 3D modeling is competent and the PV residuals are encouraging, but the manual fixing of the kinematic center and the Vrad=0 prior are choices made precisely because the data are not clean. Because the inner two of five rings overlap the anomalous structure, the fitted Vrot and sigma are not independent of the non-circular features. This directly propagates into Vc and thus into the claimed hard upper limits on gas and stellar mass. No internal inconsistency is alleged; the point is that the decisive test (masking the non-circular components and freeing the center/Vrad) has not been reported. I therefore suggest a conditional acceptance: the results are credible and well presented, but the central numbers should be shown not to depend on the contested assumptions. This aligns with, though sharpens, the reader's weakest_assumption, so agreement_with_reader is agree.","tokens_in":25837,"tokens_out":19805,"duration_ms":203942,"concrete_test":"Re-run 3DBarolo on the observed cube after masking the anomalous central structure (LoS -501 to -346 km/s, R~0.1-0.3 arcsec) and the SW/E tails, with all seven parameters free (center, Vsys, PA, i, Vrot, sigma_v, and Vrad). Compare the resulting Vrot(R), sigma_v, and Vc(R) with Fig. 4 and Fig. 6. If the rotation curve and Vrot/sigma shift by more than the quoted uncertainties, the current disk parameters and the mass-model upper limits are biased by non-circular contamination; if they are unchanged, the central claim is robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central chain of inference is the tilted-ring model in Sections 4.1-4.2. Because free 3DBarolo fits did not converge to a center on the minor axis, the kinematic center was fixed by hand; Vrad was then set to zero and the vertical scale height fixed. The non-circular structures are treated as residuals, but the anomalous central component (R~0.1-0.3 arcsec, LoS velocities -501 to -346 km/s) overlaps the innermost two of the five fitted rings. If this component or the tails contribute real flux within the fitted mask, the fitted Vrot and sigma_v are biased; those values propagate through the asymmetric-drift correction (Sec. 4.3) into Vc, and the mass-model upper limits in Sec. 5 rest directly on Vc. The approaching/receding side agreement is a good check but not a decisive one, because both sides are fit with the same fixed geometry and with Vrad fixed to zero. Thus the quoted Vrot/sigma = 6±3 and the factor-of-three mass discrepancies are conditional on the non-circular emission being small perturbations, which is precisely what the analysis assumes rather than demonstrates.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new ALMA [C I] (2-1) observations of PKS 0529-549, a radio-loud AGN-host galaxy at z ≃ 2.6, at roughly 1.5 kpc spatial resolution. Using the 3D tilted-ring code 3DBarolo, the authors model the [C I] kinematics and conclude that the molecular gas is in a dynamically cold, rotation-supported disk with Vrot/σv = 6 ± 3 and a flat rotation curve out to ~3.3 kpc. They also identify several non-circular components: a kinematic anomaly near the galaxy center, a South-West gas tail, and a weaker Eastern tail, together comprising at least 12% of the total [C I] flux. The rotation curve is used to construct mass models with different combinations of gas, stars, and a Navarro-Frenk-White dark matter halo. The single-component models give hard upper limits on gas and stellar masses that are inconsistent with photometric estimates: the stars-only limit is about a factor of three below the SED-based stellar mass, and the gas-only limit is below gas masses inferred from [C I] and dust. The paper discusses possible reasons, including non-equilibrium dynamics, an unresolved two-galaxy line of sight, and systematic problems in photometric mass estimates.","tokens_in":26089,"tokens_out":4632,"duration_ms":48054,"significance":"If the kinematic modeling is trustworthy, this is a valuable addition to the sparse sample of high-resolution, multi-tracer kinematic studies of galaxies at cosmic noon. The paper shows that a vigorous starburst/AGN host can still contain a dynamically cold molecular disk, in line with recent ALMA results, and that non-circular structures can coexist with regular rotation. The mass-model comparison highlights a real tension between dynamical and photometric mass estimators, which is an important problem for high-redshift galaxy studies. The authors are appropriately cautious in interpreting the discrepancies and lay out concrete observational tests (HST/JWST imaging, multi-line data) to resolve them. The main limitation is that the central kinematic results rest on assumptions that are not fully stress-tested, as detailed below.","major_comments":[{"comment":"The kinematic center is fixed by hand because free 3DBarolo fits do not converge to the minor axis, and the radial velocity Vrad is set to zero. These two choices are consequential: a non-zero Vrad or an offset center can change Vrot and σv, and both are used as inputs to the rotation curve and mass models. The paper should demonstrate that the derived disk parameters are robust against reasonable changes in these assumptions. Concretely, re-run the fits with Vrad as a free parameter (or test whether a non-zero Vrad improves the fit and changes Vrot), and repeat the fits with the center shifted by ±0.1–0.2 arcsec along the minor axis to see how Vrot and σv respond. Without this, the quoted Vrot/σv = 6 ± 3 remains conditional on the very assumptions that non-circular motions are negligible and that radial motions are absent.","section":"§4.1"},{"comment":"The 3DBarolo fit is performed on a masked cube that includes the anomalous central component (R ≃ 0.1–0.3 arcsec at LoS velocities −501 to −346 km s⁻¹) and the two tails. Because the innermost two or three fitted rings overlap the anomalous structure, the fitted Vrot and σv may be biased by real flux that the axisymmetric model cannot represent. The paper should quantify this by re-fitting after explicitly masking the anomalous component (e.g., excluding voxels blueward of −300 km s⁻¹) and the tails, and by showing that the resulting rotation curve and velocity dispersion are consistent within errors. The current statement that the non-circular components are 'residuals' (Section 4.4) is not a substitute for such a test, since the mask used for the fit includes them.","section":"§4.2 and §4.4"},{"comment":"The MCMC likelihood uses only the statistical uncertainties on the circular velocity Vc. The systematic uncertainties introduced by the fixed kinematic center, the fixed Vrad, the fixed vertical scale height, and the choice of a constant σv are not propagated into the mass-model upper limits. Given that the factor-of-three discrepancy between the dynamical stellar-mass upper limit and the SED-based mass is a central claim, the paper should provide a systematic error budget for the rotation curve. As a minimum, estimate how much Vc would change under alternative (yet reasonable) kinematic assumptions—e.g., non-zero Vrad, center offsets, or a 2× thicker disk—and propagate that change into the inferred masses. If the factor-of-three discrepancy persists under these perturbations, the claim is much stronger.","section":"§5, Eq. (2)"}],"minor_comments":[{"comment":"The chemical symbol 'C I' is typeset inconsistently as '[Ci]' and '[C i]' (also in the abstract and Section 2). Please unify the notation, e.g., '[C I]'.","section":"Throughout"},{"comment":"The uncertainty formula uses the abbreviation NMAD, which is never defined; please spell out 'normalized median absolute deviation' and explain why this particular combination is used.","section":"Eq. (1)"},{"comment":"The combined uncertainty on Vrot/σv = 6 ± 3 is quoted without showing the propagation of the individual errors on Vrot and σv. Please state the assumed errors on Vrot and σv used in this calculation.","section":"§4.2"},{"comment":"The note says 'random errors of the fluxes are less than 1%', but the uncertainties shown are 10% of the flux. If this is a calibration error, please label it as such and use the term 'systematic' rather than 'random'.","section":"Table 3"},{"comment":"The term 'Renzograms' (a contraction of Renzini + diagrams) is used without definition or a citation; please add a brief explanation or a reference to the original definition.","section":"§4.4"},{"comment":"The statement that 1.1 × 10¹¹ M⊙ is a 'hard upper limit' on the stellar mass should be qualified, because it depends on the assumed spherical geometry; the flattened disk+bulge model in the appendix gives a lower value. Consider stating that this is an upper limit under the spherical assumption.","section":"§5.2.2"},{"comment":"The received/accepted dates in the manuscript header ('Received September 15, 1996; accepted March 16, 1997') are clearly template artifacts and should be corrected before final submission.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of A&A and presents a valuable dataset. The main technical concern is that the kinematic modeling assumptions (fixed center, Vrad = 0, mask including non-circular features) are not demonstrated to be innocuous for the derived rotation curve and mass limits. These issues are addressable with additional robustness tests, so major revision seems appropriate rather than rejection. I would also flag that the paper's mass discrepancies rest on a comparison with published photometric masses from a single SED-fitting code (De Breuck et al. 2010); the authors' own preliminary CIGALE fits (Section 6.2.2) show that the SED-based mass is highly model-dependent, so the factor-of-three claim should be framed as a systematic tension rather than a definitive contradiction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a strong single-object kinematics paper, and the headline claims are mostly believable. The new ALMA [CI](2-1) cube at 0.18\" resolves the disk in PKS 0529-549, shows Vrot/sigma = 6±3 with a flat rotation curve to 3.3 kpc, and reveals a SW tail, a probable E tail, and a central anomalous structure. That is new relative to Lelli et al. 2018, and it is presented with the usual supporting evidence: channel maps, a PV diagram, separate fits to approaching and receding sides, and a clean residual map.\n\nWhat the paper does well: the kinematic modeling is careful, the mass models are clearly specified, and the authors are honest about degeneracies. The dynamical masses are framed as upper limits, and the comparison with SED/CO/[CI] masses is an external benchmark, not a self-fulfilling fit. The discussion of possible explanations for the mass discrepancy—SED systematics, non-equilibrium kinematics, or two galaxies in projection—is genuinely balanced.\n\nThe soft spots are real but not fatal. The stress-test note is on target: the kinematic center is fixed by hand because free fits do not converge, Vrad is set to zero, and the anomalous central component overlaps the innermost fitted rings. If that component or the tails contribute flux inside the mask, Vrot and sigma_v are biased, and those feed directly into the circular velocities and mass limits. I would add that the residual map shows negative residuals along the minor axis, which is an artifact of azimuthal averaging, so the non-circular flux fraction of 12% is a lower limit. That said, the approaching/receding consistency is a decent check, and the authors explicitly discuss the possibility that the disk is not in equilibrium. The mass discrepancy is therefore conditional on assumptions the paper mostly names, even if it cannot test them.\n\nOne more thing: the SED comparison uses the old De Breuck et al. (2010) stellar mass, and the authors' own preliminary CIGALE fits suggest a lower mass is possible. So the factor-of-3 to 6 discrepancy is suggestive, not established.\n\nBottom line: this paper deserves a serious referee. It is a high-quality single-object case study that will be useful to anyone working on high-z gas kinematics, AGN hosts, or disk-halo degeneracy. I would send it to review, with the expectation that the referee asks for a robustness test of the kinematic model—for example fitting with Vrad free or masking the anomalous component—but I would expect it to survive.","headline":"A careful single-object kinematics paper whose new ALMA [CI] data genuinely resolve the disk and non-circular structures, with the mass discrepancy claims real but conditional on kinematic assumptions the authors mostly name themselves.","tokens_in":26732,"tokens_out":3062,"would_cite":true,"duration_ms":31357,"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 molecular gas in z≈2.6 AGN-host galaxy PKS 0529-549 forms a dynamically cold, rotation-supported disk with a flat rotation curve to 3.3 kpc, and its rotation-curve mass limits contradict photometric mass estimates.","keywords":["dark matter","galaxies: active","galaxies: high-redshift","galaxies: kinematics and dynamics","galaxies: formation","galaxies: evolution","molecular gas kinematics","rotating disks"],"falsifier":"Deep rest-frame optical or near-infrared imaging (JWST or HST) of PKS 0529-549 would settle the central tension: a single stellar component with a photometric mass near $3\\times10^{11}\\,M_\\odot$ inside 3.3 kpc cannot coexist with an equilibrium flat rotation curve of roughly $300$ km s$^{-1}$, whereas resolving two separate stellar components along the line of sight would confirm the superposition scenario and invalidate the single-disk dynamical model. A complementary test is to map the anomalous central structure at higher resolution: a front-back velocity reversal along the minor axis would reveal the radial motions that the current model sets to zero.","tokens_in":25596,"feed_emoji":"🌌","tokens_out":21516,"duration_ms":305297,"temperature":0.7,"pith_summary":"This paper asks whether the molecular gas in PKS 0529-549 — a radio-loud active galaxy at redshift 2.6, observed in the epoch of peak star formation and black-hole growth — is in ordered motion. Using new ALMA observations of the [C i] (2-1) fine-structure line of neutral carbon, a tracer of molecular gas, at about 1.5 kpc resolution, it argues that the gas forms a dynamically cold, rotation-supported disk with $V_{\\rm rot}/\\sigma_{\\rm v}=6\\pm3$ (rotation about six times the gas turbulence) and a flat rotation curve out to 3.3 kpc. That order is surprising in a galaxy with such extreme star formation and AGN activity, and it coexists with clear disorder: a kinematically anomalous central structure and two gas tails that together carry at least 12% of the line flux and may record a past merger. Turning the rotation curve into mass models gives hard upper limits on stellar and gas mass that clash with photometric estimates, leaving open whether the photometric masses, the equilibrium assumption, or a pair of galaxies along the line of sight is at fault.","feed_headline":"Gas disk spins six times faster than it jitters at z≈2.6","feed_subtitle":"New ALMA maps of PKS 0529-549 reveal merger-like gas tails and dynamical masses that don't match photometric estimates.","key_machinery":"The argument is carried by 3D kinematic modeling of the ALMA line cube rather than by 2D moment maps. The tilted-ring code 3DBarolo — which fits concentric rotating rings directly to the three-dimensional data cube, convolving the model with the telescope beam — corrects for beam smearing and for the intensity-weighting bias that flattens apparent rotation curves; each ring carries geometric parameters (center, position angle, inclination, systemic velocity) and physical ones (rotation velocity, velocity dispersion, surface density), with the kinematic center fixed by hand because free fits do not converge. An asymmetric-drift correction converts the fitted rotation velocity into the circular velocity $V_{\\rm c}(R)$ using the radial profiles of gas surface density and velocity dispersion. That circular-velocity curve is then fed into MCMC mass models whose model velocities combine a Sérsic-profile stellar component, a disk gas component, and an optional NFW dark-matter halo, with $\\Lambda$CDM scaling relations used as priors.","core_discovery":"The central claim is that PKS 0529-549 contains a regular, dynamically cold rotating disk of molecular gas, with $V_{\\rm rot}/\\sigma_{\\rm v}=6\\pm3$ and a flat rotation curve that yields a total dynamical mass of about $10^{11}\\,M_\\odot$ within 3.3 kpc — an orderly disk inside a galaxy whose star formation rate of order $10^3\\,M_\\odot\\,\\rm yr^{-1}$ and radio-loud AGN make such order unexpected. The same data show that the disk coexists with non-circular motions: a kinematically anomalous structure at about 2 kpc from the center, a South-West gas tail, and a weaker Eastern tail, together at least 12% of the [C i] (2-1) flux, plausibly the leftovers of a major merger. Mass models fitted to the asymmetric-drift-corrected circular velocity curve set hard dynamical upper limits — $M_{\\rm gas}\\simeq 8.6\\times10^{10}\\,M_\\odot$ and $M_\\star\\simeq1.1\\times10^{11}\\,M_\\odot$ — but the gas-only model cannot reproduce the inner rotation curve, and the stellar limit sits about a factor of three below the SED-based stellar mass of $3\\times10^{11}\\,M_\\odot$. Models with and without a dark-matter halo fit equally well, leaving the dark-matter content unconstrained, and the paper concludes that the discrepancy could stem from the photometric masses, from a disk not yet in dynamical equilibrium, or from a second galaxy hiding along the line of sight.","pith_inferences":["A plausible generalization the paper does not make: if cold molecular disks persist inside such violently active hosts, AGN feedback at cosmic noon may act mainly on the ionized gas phase while leaving the cold gas disk intact, which would sharpen how feedback models are judged.","The superposition scenario, if confirmed in this system, would imply that a fraction of high-redshift 'single' galaxies with AGN are chance alignments, so any kinematic mass measurement on unresolved targets would be systematically suspect.","A testable extension the five-ring model cannot address: higher-resolution mapping of the anomalous central structure should show a minor-axis velocity reversal if it is an inflow related to the South-West tail, as the paper suggests."],"forward_implications":["If the measured ratio is right, extreme star formation and a radio-loud AGN do not necessarily destroy ordered rotation in the molecular gas of a cosmic-noon galaxy.","The dynamical upper limits imply that SED-based stellar masses and standard [C i]- and dust-based gas masses can exceed what the gravitational potential allows, unless the disk is out of equilibrium or two galaxies are aligned along the line of sight.","Because baryons-only and baryons-plus-DM models fit the rotation curve equally well, dark-matter fractions at high redshift cannot be measured securely until baryonic masses are known to about 25 percent precision.","Non-circular structures carrying at least 12% of the line flux would be smeared away or misread at lower resolution, so full 3D line-cube analysis is necessary to separate rotation, merger remnants, and possible outflows in high-z galaxies."],"supporting_citations":[{"why":"Supplies the earlier low-resolution [C i] (2-1) detection of the rotating disk and the initial geometric values (position angle 75 degrees, inclination 50 degrees) that the new fits confirm.","marker":"Lelli et al. 2018"},{"why":"Provides the 3DBarolo tilted-ring code used to fit the rotating disk directly to the three-dimensional [C i] cube with beam smearing accounted for.","marker":"Di Teodoro & Fraternali 2015"},{"why":"Supplies the disk-instability model prediction of the rotation-to-dispersion ratio that the measured value of 6 ± 3 is compared against.","marker":"Wisnioski et al. 2015"},{"why":"Supplies the very high star formation rate of about one thousand solar masses per year, which makes the cold orderly disk unexpected.","marker":"Falkendal et al. 2019"},{"why":"Supplies the SED-based stellar mass (3×10^11 M_sun) that the dynamical upper limit falls below by a factor of about three.","marker":"De Breuck et al. 2010"},{"why":"Supplies the CO (4-3) and [C i] (1-0) luminosities and dust-based gas masses that the dynamical gas limits are compared against.","marker":"Huang et al. 2024"},{"why":"Supplies the [O iii] velocity field and its redshift offset of about 350 km/s from [C i], the key evidence for the two-galaxies scenario.","marker":"Nesvadba et al. 2017"},{"why":"Provides the abundance-matching stellar-to-halo mass relation used as a prior on the dark-matter halo mass in the baryons-plus-DM model.","marker":"Legrand et al. 2019"},{"why":"Provides the mass-concentration relation for NFW halos used as a prior in the baryons-plus-DM model.","marker":"Dutton & Macciò 2014"}],"fun_headline_variants":["Cosmic noon galaxy spins fast but hides mass mismatch","AGN host at z≈2.6 has orderly gas disk, messy edges","Rotating gas disk in distant AGN galaxy challenges mass models","PKS 0529-549: fast disk, slow tails, and mass mystery"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The kinematic modeling assumes a single, axisymmetric, flat rotating disk with zero radial motions and a hand-fixed center (free fits of the center do not converge), so if the anomalous central structure or the gas tails belong to the disk, or if radial flows are significant, the fitted rotation speeds and the mass limits derived from them would be biased.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic noon galaxy spins fast but hides mass mismatch","AGN host at z≈2.6 has orderly gas disk, messy edges","Rotating gas disk in distant AGN galaxy challenges mass models","PKS 0529-549: fast disk, slow tails, and mass mystery"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000731,"raw_usage":{"total_tokens":3410,"prompt_tokens":1225,"completion_tokens":2185,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":841,"completion_tokens_details":{"reasoning_tokens":2119}},"tokens_in":841,"tokens_out":2185,"duration_ms":15568,"temperature":1.0,"reasoning_tokens":2119,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:12:49.862100+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Deep rest-frame optical or near-infrared imaging (JWST or HST) of PKS 0529-549 would settle the central tension: a single stellar component with a photometric mass near $3\\times10^{11}\\,M_\\odot$ inside 3.3 kpc cannot coexist with an equilibrium flat rotation curve of roughly $300$ km s$^{-1}$, whereas resolving two separate stellar components along the line of sight would confirm the superposition scenario and invalidate the single-disk dynamical model. A complementary test is to map the anomalous central structure at higher resolution: a front-back velocity reversal along the minor axis would reveal the radial motions that the current model sets to zero.","supporting_citations":[{"cited_title":"2018, MNRAS, 479, 5440","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier low-resolution [C i] (2-1) detection of the rotating disk and the initial geometric values (position angle 75 degrees, inclination 50 degrees) that the new fits confirm."},{"cited_title":"corner plots","cited_arxiv_id":null,"evidence_quote":"Supplies the disk-instability model prediction of the rotation-to-dispersion ratio that the measured value of 6 ± 3 is compared against."},{"cited_title":"D., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the very high star formation rate of about one thousand solar masses per year, which makes the cold orderly disk unexpected."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the [O iii] velocity field and its redshift offset of about 350 km/s from [C i], the key evidence for the two-galaxies scenario."},{"cited_title":"J., Davidzon, I., et al","cited_arxiv_id":null,"evidence_quote":"Provides the abundance-matching stellar-to-halo mass relation used as a prior on the dark-matter halo mass in the baryons-plus-DM model."}],"review_version":1}