{"id":"4ec2f0b3-6e87-44a9-b43a-9b84f4dc2f83","arxiv_id":"2412.09363","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A radio-selected, optically dark galaxy at z=4.82 is a massive, compact main-sequence galaxy with cold, optically thick dust and tentative evidence for an AGN.","lead":"Using ALMA, JWST, and X-ray data, astronomers confirmed a dusty, star-forming galaxy at redshift 4.82 that is invisible in most optical and near-infrared filters. The galaxy appears massive, compact, and shrouded in unusually cold, optically thick dust, with a tentative X-ray detection hinting at an active black hole.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The FIR SED is fit with a single-temperature MBB using only two robust submm points plus 3 mm; the same data yield a factor-of-2 lower L_IR in the MBB fit than in STARDUST, so the cold-T and SFR headline claims are not yet anchored.","rationale":"The central claims that set this paper apart are the surprisingly cold T=33 K, the optically thick FIR diagnosis, and the quantitative SFR/main-sequence characterization. The spectroscopic redshift and stellar mass are well supported by ALMA line detections and three independent SED codes, so those are not the weak point. The weak point is the FIR peak: the MBB fit uses only tail points, and the internal factor-of-two L_IR mismatch between Fig. B.1 and Fig. 4 is the clearest symptom. Two models can both pass through the few observed submm/mm points and differ in L_IR by 2x only because the peak is not sampled. This means the single-temperature T is a model-dependent extrapolation, not a direct measurement. The X-ray hint of an AGN raises the realistic possibility of a warm dust component, which would make the single-MBB T the temperature of the cold reservoir rather than the galaxy's effective dust temperature, weakening the 'coldest DSFG' comparison and making the SFR uncertain if AGN-heated dust contributes to L_IR. The paper's own diagnostics (Jin et al. 2022) are useful but use the same sparse photometry and the same area prior; they do not independently resolve the peak-luminosity ambiguity. The case study is valuable and likely correct in its essentials: the source is massive, dusty, and probably compact. But the quantitative headline claims need the FIR peak anchored, which is exactly the conditional concern the reader raised. I therefore recommend no change to the conditional verdict.","tokens_in":16424,"tokens_out":16928,"duration_ms":184321,"concrete_test":"Refit the same photometry (SCUBA-2 450/850, ALMA 2 mm and 3 mm, Herschel limits) with a two-temperature MBB and with a library SED code including an AGN component. Check whether the cold-component T shifts by more than ~5 K and whether total L_IR moves by more than ~30% when a warm component is added at a level consistent with the X-ray/AGN constraints. Also recompute tau_100 and lambda0 using twice the assumed emitting area; if T stays within 33+/-5 K, L_IR within 40% of the STARDUST value, and tau_100 remains >1, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline cold-T and SFR claims rest on an unconstrained FIR peak. The mercurius MBB fit (Fig. B.1, Sect. 3.3) uses SCUBA-2 450/850 um, tentative ALMA 2 mm, and 3 mm continuum, all on the long-wavelength tail (rest 77-515 um), with no data near the rest ~20-60 um peak of a T=33 K MBB. Consequently L_IR is strongly model-dependent: the MBB fit returns L_IR ~ 2.6e12 Lsun (SFR_IR ~ 390 Msun/yr) while STARDUST (Fig. 4) returns L_IR ~ 5.4e12 Lsun (SFR = 540), a factor ~2.1 discrepancy never reconciled in the paper. The abstract's 'main-sequence' and 'coldest DSFG' claims depend on these numbers. If a warmer AGN-heated component is present (as the tentative X-ray detection hints), the single-T MBB measures only the cold reservoir, the quoted T is not a luminosity-weighted dust temperature, and the optically-thick/lambda0 diagnosis is built on a model that underpredicts L_IR by 2x; the T-z comparison then is not apples-to-apples. The quoted +/-2 K and +/-0.2 errors do not include this systematic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the spectroscopic confirmation and panchromatic analysis of XS55, a radio-selected optically dark galaxy at z=4.8214, combining ALMA CO(5-4) and [CI](1-0) detections, JWST/NIRCam imaging, SCUBA-2 photometry, and X-ray data. The authors derive a stellar mass of M*=5e10 M_sun, a star formation rate of ~540 M_sun/yr, a dust mass of ~1.7e9 M_sun, and a dust temperature of ~33 K assuming optically thick dust, and they argue that the FIR dust is optically thick and surprisingly cold. They also report tentative X-ray emission and a compact central component, suggesting the presence of an obscured AGN. The abstract presents XS55 as a massive main-sequence galaxy that is one of the coldest DSFGs at z>4.","tokens_in":16828,"tokens_out":8651,"duration_ms":80905,"significance":"The spectroscopic redshift and the multi-wavelength census are valuable, given the small sample of spectroscopically confirmed optically dark DSFGs at z>4; the tentative X-ray detection adds to evidence for obscured AGN in this population. The morphological decomposition with JWST is a strength, and the paper makes productive use of standard fitting tools. However, the headline cold-dust and optically-thick claims rest on a sparse FIR SED and a model-dependent infrared luminosity, so the quantitative conclusions (T_dust, SFR, main-sequence status) are not yet at the level of robustness implied by the abstract. The paper would be more convincing if the FIR SED modeling and the consistency between the two SED codes were addressed.","major_comments":[{"comment":"The single-temperature modified blackbody fit is underconstrained: the available photometry (SCUBA-2 450 and 850 um, ALMA 3 mm, and a tentative 2 mm point) covers rest wavelengths of only 77-516 um, with no data near the expected rest-frame peak of a T=33 K modified blackbody (~110 um). The quoted T_dust uncertainties of +/-2 K are statistical only. The optically thick case sets the emitting area to the F444W half-light area (9.08 kpc^2) as an upper limit; this choice yields a lower limit on T_dust, so the 'cold' finding is conditional on this assumption and would weaken if the actual emitting area is smaller (e.g., the ALMA continuum upper limit Re < 1.76 kpc). The authors should demonstrate the robustness of T_dust and the optically-thick/thin diagnosis to (a) alternative SED models such as two-temperature or template-based fits, (b) the unknown FIR peak wavelength, and (c) the choice of emitting area.","section":"Section 3.3, Fig. B.1"},{"comment":"There is an unresolved factor-of-two discrepancy between the infrared luminosities from the two fitting approaches: the mercurius MBB fit gives L_IR ~ 2.6e12 L_sun and SFR_IR ~ 390 M_sun/yr, while the STARDUST panchromatic fit gives L_IR ~ 5.4e12 L_sun and SFR = 540 M_sun/yr. The abstract quotes the STARDUST SFR, while the cold-T and optically-thick claims come from the MBB fit. The text notes that 'the modified blackbody is not accounted for in the total SED fit' (Fig. 4 caption) but does not explain why the two models disagree or which L_IR should be adopted for the main-sequence and dust-temperature comparisons. This inconsistency is load-bearing: using the MBB SFR would reduce the main-sequence offset, while using the STARDUST L_IR implies that the single-temperature MBB is missing a substantial luminosity component, calling into question the cold-dust and optically-thick conclusions.","section":"Section 3.3, Fig. 4, Fig. B.1"},{"comment":"The gas-mass comparison that supports optically thick dust depends critically on the adopted alpha_CI and delta_gdr. The paper claims agreement between M_mol,thick ~ 1.7e11 M_sun and M_mol,[CI] ~ 1.8e11/alpha_CI with alpha_CI = 17, but the text also cites the high-z SMG calibration alpha_CI = 4.1 +/- 1.4 (Frias Castillo et al. 2024), which would give M_mol,[CI] ~ 4.3e10 M_sun, a factor of four lower. The choice of alpha_CI is therefore decisive for the first of the three optically-thick diagnostics. The paper should present the gas-mass comparison under both calibrations and state how robust the conclusion is to this uncertainty and to the assumed delta_gdr = 100.","section":"Section 4.1"},{"comment":"The comparison of T_dust with the Schreiber et al. (2018) main-sequence T-z relation may not be apples-to-apples, as that relation is presumably based on optically thin modified blackbody fits while XS55's T is from the optically thick fit. The paper should clarify the fitting conventions used for the comparison sample and, if appropriate, recompute the offset using a thin-fit temperature for XS55 or a consistently calibrated relation. Without this, the claim that XS55 lies 0.13 dex (4 sigma) below the relation is not well defined.","section":"Section 4.1, Fig. 5 (left)"}],"minor_comments":[{"comment":"The literature labels in Fig. 5 are inconsistent with the reference list: 'Marrone et al. 2017' should be 'Marrone et al. 2018' and 'Hodge et al. 2018' should be 'Hodge et al. 2019'.","section":"Fig. 5 captions"},{"comment":"The phrase 'As shown in see Fig. 2' should read 'As shown in Fig. 2'.","section":"Section 2.2"},{"comment":"The statement 'The diffuse component cannot be fit' is ambiguous; it likely means the diffuse component is not significantly detected in F277W. Please clarify.","section":"Section 3.2"},{"comment":"The description of the XMM-Newton detection as 'detected in the stacked soft, medium and hard X-ray XMM-Newton images with 3.1 sigma' is unclear; please specify whether this is a single combined image or three individual detections.","section":"Section 2.4"},{"comment":"The term 'optically dark' is used to describe XS55, but the source is detected in F277W and F444W; the definition (e.g., F150W dropout) should be stated precisely in the abstract or selection section to avoid confusion.","section":"Abstract / Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is well-written and presents a rich dataset; the spectroscopic redshift and multi-wavelength census are solid contributions. The main issue is the inconsistency between the mercurius MBB fit and the STARDUST panchromatic fit, which directly affects the headline SFR, T_dust, and main-sequence claims. I would not reject, but the authors need to reconcile the two models or clearly justify their separate use. The paper would also be strengthened by a systematic uncertainty analysis of the FIR SED fit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a solid case study of a single radio-selected, JWST-dark DSFG at z=4.82, with a secure ALMA redshift and a rich multi-wavelength dataset. The headline claims about cold dust and optically thick FIR emission are plausible but rest on sparse photometry and a factor-of-two disagreement between the two SED codes that the paper never reconciles. The redshift and basic stellar properties are safe.\n\nWhat is actually new: this is the first spectroscopic confirmation and panchromatic characterization of XS55. The ALMA detections of CO(5-4) and [CI](1-0) are convincing (9.8 and 3.9 sigma, respectively), and the authors do the right sanity check against the z=3.66 alternative. The JWST morphology fits are standard but well-executed, and the X-ray detection is honestly flagged as tentative. The gas-mass comparison between [CI] and dust mass is a nice cross-check, and the authors consider both thin and thick dust cases. The paper is careful about its limitations.\n\nThe soft spot is the FIR SED. The mercurius MBB fit uses only two robust submm points (SCUBA-2 450 and 850 um), a 3 mm continuum detection, and a tentative 2 mm point. For a T~33 K modified blackbody, the rest-frame peak is around 100-110 um, so the SCUBA-2 points do bracket the peak; the stress-test claim that the peak is entirely unobserved is not quite right. The real problem is that the MBB fit returns L_IR = 2.6e12 Lsun while STARDUST returns 5.4e12 Lsun, a factor of two that is never explained. That discrepancy propagates into SFR, dust mass, and the surface-density argument, so the cold-T and optically-thick claims carry a systematic error well larger than the quoted ±2 K. The abstract's \"coldest DSFG at z>4\" is therefore a bit strong for a single-component model on four photometric points. The gas-mass comparison also assumes a gas-to-dust ratio and [CI] calibration that favor the thick case; the authors mention the alternative alpha_CI, but the choice is not neutral.\n\nWho should read this: anyone working on optically dark / JWST-dark galaxies and the obscured AGN population. It is a useful addition to the small sample of spectroscopically confirmed DSFGs at z>4. It deserves serious refereeing, but the referee should push for a discussion of the L_IR discrepancy and a clearer caveat on the model dependence of T_dust. As it stands, the core redshift and mass are solid; the cold-dust claim is plausible but not yet anchored.","headline":"A solid case study of one radio-selected, JWST-dark DSFG at z=4.82 with a secure ALMA redshift; the cold-dust and optically thick claims are plausible but rest on sparse FIR photometry and a factor-of-two L_IR discrepancy between fitting codes that the paper never reconciles.","tokens_in":17511,"tokens_out":4546,"would_cite":true,"duration_ms":45123,"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 paper claims that XS55, an optically invisible radio-selected galaxy at z=4.82, is a massive main-sequence galaxy with surprisingly cold, optically thick dust and a likely buried active galactic nucleus.","keywords":["dusty star-forming galaxies","optically dark galaxies","high-redshift galaxy evolution","submillimeter galaxies","active galactic nuclei","dust temperature","JWST NIRCam","ALMA spectroscopy"],"falsifier":"Measure the dust continuum size with high-resolution ALMA at 3 mm or at 450 µm and test a two-temperature dust model: if the far-infrared SED requires a warm component or the emitting area exceeds the F444W half-light radius of 9.08 kpc$^2$, the cold, optically thick single-temperature interpretation fails.","tokens_in":16236,"feed_emoji":"🌌","tokens_out":11348,"duration_ms":104801,"temperature":0.7,"pith_summary":"XS55 is a galaxy at redshift 4.82 that is invisible in optical and near-ultraviolet images but bright at radio and submillimetre wavelengths — an 'optically dark' dusty star-forming galaxy, the kind thought to carry a large share of early massive-galaxy formation. The paper establishes a secure redshift from ALMA detections of the CO(5-4) and [CI](1-0) lines, and then assembles a panchromatic spectral energy distribution from JWST, infrared, submillimetre, and radio photometry. It finds that XS55 sits on the galaxy main sequence — a normal star-forming galaxy rather than a starburst — with stellar mass $(5\\pm1)\\times10^{10}\\,M_\\odot$, star-formation rate $540\\pm177\\,M_\\odot\\,\\mathrm{yr}^{-1}$, and a compact morphology: a small central component, an extended disk, and a likely companion. The paper's central physical claim is that the dust is optically thick in the far infrared and still cold, $T_\\mathrm{dust}=33\\pm2$ K, one of the coldest dusty star-forming galaxies found at $z>4$, and that tentative X-ray emission plus the compact central component indicate an active black hole buried in the dust. If right, XS55 is an example of an ordinary massive galaxy caught in a heavily obscured phase, visible only through the few photons that leak around the dust.","feed_headline":"Cold, thick dust hides a massive normal galaxy at z=4.82","feed_subtitle":"An optically invisible galaxy at z=4.82 is a massive main-sequence system hiding a likely black hole.","key_machinery":"The central object is a panchromatic spectral energy distribution built from sparse photometry: near-infrared detections in F277W and F444W with an F150W dropout, SCUBA-2 450 and 850 µm fluxes, ALMA 2 and 3 mm data, and 1.3 and 3 GHz radio points. The argument is carried by a modified blackbody model with emissivity index $\\beta_{\\rm IR}=2.0\\pm0.2$, fitted twice — once with dust assumed optically thin and once with a 'self-consistent' optically thick prescription in which the emitting area cannot exceed the 9.08 kpc$^2$ half-light area measured in F444W. Three independent diagnostics (gas mass from [CI](1-0) compared with gas mass from dust, opacity at 100 µm, and position on the infrared-luminosity-surface-density versus dust-temperature diagram) are what turn the fit into the claim that the dust is optically thick. The molecular-gas mass and the X-ray-based AGN interpretation complete the picture of all baryonic components.","core_discovery":"On the paper's own account, XS55 is not an extreme starburst but a massive main-sequence galaxy caught behind a dust veil: $M_\\ast=(5\\pm1)\\times10^{10}\\,M_\\odot$, $\\mathrm{SFR}=540\\pm177\\,M_\\odot\\,\\mathrm{yr}^{-1}$, within $2\\sigma$ of the main sequence at its redshift. The far-infrared SED is best described as optically thick with cold dust $T_\\mathrm{dust}=33\\pm2$ K, which resolves an apparent violation of the Stefan-Boltzmann limit on infrared surface brightness: the optically thin fit would require an implausibly large dust mass and would disagree with the gas mass measured from [CI] emission, whereas the optically thick fit with an emitting area capped by the F444W half-light radius (9.08 kpc$^2$) is consistent with both. The compact central component, tentatively detected in X-ray observations, suggests an active galactic nucleus, making XS55 a candidate for a buried black hole in a compact, cold, dusty massive galaxy at $z=4.82$.","pith_inferences":["The cold-dust, optically thick interpretation implies a high gas-phase metallicity (the local radiation field scales as star-formation efficiency divided by metallicity), so an independent metallicity measurement from far-infrared fine-structure or millimetre recombination lines is a testable consequence the paper does not carry out.","If XS55 is representative, X-ray and radio surveys may be systematically missing obscured AGN in the optically dark population; stacking the X-ray emission of a larger sample of similar radio-selected F150W dropouts would turn the single-source 3.1σ detection into a population constraint.","Because only two robust submillimetre bands anchor the far-infrared SED, the method could be extended with ALMA Band 8 or 9 photometry at 400–700 µm; detecting a second, warmer dust component would separate 'intrinsically cold' from 'optically thick' more cleanly."],"forward_implications":["A spectroscopically confirmed optically dark dusty star-forming galaxy at z≈4.8 can be a normal massive main-sequence galaxy rather than an extreme starburst, broadening the interpretation of optically dark selections.","The gas mass inferred from [CI](1-0) agrees with the optically thick dust model ($M_{\\rm dust}\\approx1.7\\times10^9\\,M_\\odot$) and favors it over the thin model, so self-consistent thick-dust SED fitting is needed to avoid overestimating dust and gas masses.","XS55's tentative X-ray detection and lack of radio excess imply that some AGN in optically dark dusty galaxies are visible only in X-rays, not in radio, so AGN fractions from radio-selected samples may be lower limits.","The compact stellar size, 2.6 times below the main-sequence mass-size relation, and the large cold dust mass together explain why XS55 is an F150W dropout, making compactness a predictor of optical darkness."],"supporting_citations":[{"why":"Provides the ALMA 3 mm line-scan method, the line-search algorithm, and the chance-probability check used to identify CO(5-4) and [CI](1-0) at z=4.8214.","marker":"Jin et al. 2019"},{"why":"Supplies the ultra-deep 3 GHz radio catalogue from which XS55 was selected as an optically dark radio source.","marker":"Algera et al. 2020"},{"why":"Gives the modified blackbody dust model and the gas-to-dust ratio assumed to convert dust mass to molecular gas mass.","marker":"Magdis et al. 2012"},{"why":"Provides the fitting code used for the optically thin and optically thick far-infrared SED fits.","marker":"Witstok et al. 2022"},{"why":"Defines the three optically thick dust diagnostics (gas-mass comparison, tau at 100 µm, and infrared-surface-density versus temperature) applied to XS55.","marker":"Jin et al. 2022"},{"why":"Supplies the panchromatic SED fitting code used to derive stellar mass and star-formation rate from the NIR-to-radio data.","marker":"Kokorev et al. 2021"},{"why":"Supplies the alpha_CI conversion used to derive molecular gas mass from [CI](1-0) for comparison with the dust-based gas masses.","marker":"Dunne et al. 2022"},{"why":"Defines the main-sequence relation used to classify XS55 as a massive main-sequence galaxy.","marker":"Schreiber et al. 2015"},{"why":"Provides the X-ray images and flux-conversion factors behind the tentative AGN detection.","marker":"Civano et al. 2016"},{"why":"Supplies the JWST/NIRCam survey data in which XS55 is detected in F277W and F444W but drops out in F150W.","marker":"Casey et al. 2023"}],"fun_headline_variants":["Dust-veiled galaxy at z=4.82 is massive and ordinary","Cold, thick dust hides a normal massive galaxy in early universe","Optically dark galaxy at z=4.82 hosts a buried black hole","A cold dark galaxy at z=4.82 hides an active nucleus","Massive normal galaxy discovered behind thick dust at z=4.82"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire cold-and-thick dust picture rests on a single-temperature model fitted to just a few submillimetre points, with the dust's emissivity and the size of the emitting region fixed by assumption; if the dust is a mixture of temperatures or the emitting area is larger than the adopted F444W half-light area, the quoted temperature and the optically thick conclusion would change.","fun_headline_variants_meta":{"raw":{"variants":["Dust-veiled galaxy at z=4.82 is massive and ordinary","Cold, thick dust hides a normal massive galaxy in early universe","Optically dark galaxy at z=4.82 hosts a buried black hole","A cold dark galaxy at z=4.82 hides an active nucleus","Massive normal galaxy discovered behind thick dust at z=4.82"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000504,"raw_usage":{"total_tokens":2563,"prompt_tokens":1150,"completion_tokens":1413,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":1330}},"tokens_in":766,"tokens_out":1413,"duration_ms":13630,"temperature":1.0,"reasoning_tokens":1330,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:06:41.866264+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dust continuum size with high-resolution ALMA at 3 mm or at 450 µm and test a two-temperature dust model: if the far-infrared SED requires a warm component or the emitting area exceeds the F444W half-light radius of 9.08 kpc$^2$, the cold, optically thick single-temperature interpretation fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the ultra-deep 3 GHz radio catalogue from which XS55 was selected as an optically dark radio source."}],"review_version":1}