{"id":"52fdc28a-a4b5-477a-9e42-ef2026876a28","arxiv_id":"2607.25304","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"In IRAS 06530-0213, HC3N and C4H emission is concentrated in an inner equatorial zone, and radiative-transfer modeling indicates HC3N's abundance, not just its excitation, peaks there.","lead":"New radio maps of the dying star IRAS 06530-0213 show that carbon-chain molecules HC3N and C4H are packed into a compact equatorial ring, while carbon monoxide fills a much larger barrel-shaped shell. This suggests the dense inner belt of this protoplanetary nebula is a chemical factory for carbon chains, possibly related to the unexplained 21-micron infrared feature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"1D spherical RATRAN model cannot uniquely separate abundance gradients from barrel geometry; the 'requires' claim is underdetermined.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the 1D spherical RATRAN model may not adequately represent the equatorial sector of a barrel-like, likely non-spherical source. This concern is well founded. The paper's own language in Sect. 4.1 acknowledges the model is not a unique reconstruction. The constant-abundance HC3N model failing to match the data is real evidence, but it is evidence only within the adopted spherical geometry; it does not by itself rule out geometric or density-structure explanations for the compact HC3N distribution. The proposed 3D radiative transfer test is the direct, decisive check: if a barrel geometry with constant HC3N abundance reproduces the observations, the paper's central 'requires' claim would need to be softened; if it still fails, the abundance-gradient interpretation is strongly supported. Since the reader already assigned CONDITIONAL verdict on this basis, my stress-test does not change that verdict. I see no reason to upgrade to ACCEPT (the test is not yet done) nor to downgrade to REJECT (the observational results are new and the constant-abundance spherical model does fail, providing partial support). The paper would be strengthened by running the 3D test or by explicitly reducing the 'requires' language to 'is consistent with' or 'suggests' in the abstract and conclusion.","tokens_in":20548,"tokens_out":3466,"duration_ms":40679,"concrete_test":"Run a 3D radiative transfer model (e.g., RADMC-3D or LIME) with a barrel/torus density structure that reproduces the CO and 13CO moment-zero maps and PV diagrams (including the equatorial-sector profiles). Keep the HC3N abundance constant with radius and compute the predicted HC3N J=24–23 equatorial radial profile, convolved with the same beam. If this constant-abundance model matches the observed compact HC3N profile within the noise, the abundance-gradient claim is not required; if it still overextends beyond ~0.8′′, the abundance gradient is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that reproducing the HC3N morphology 'requires' an inwardly concentrated abundance distribution (Sect. 5, abstract)—rests on the 1D spherical RATRAN model in Sect. 4.1. The model adopts a single power-law density n_H2(r) ∝ r^−2, a homologous expansion, and a CO-constrained shell (R_in = 0.45′′, R_out = 0.87′′), then compares the predicted radial profile with observed equatorial-sector profiles. The load-bearing weakness is that this 1D spherical geometry may not represent the equatorial sector of a clearly barrel-like, non-spherical PPN. In the real source, the observed emission at large projected radii passes through the thin barrel walls and the surrounding low-density halo; in the spherical model, the line of sight at those projected radii traverses full spherical shells, which can overproduce emission for any molecule with constant abundance. The CO and 13CO fits do not remove this ambiguity because CO is optically thick and the shell parameters were tuned to CO. Thus the failure of the constant-abundance HC3N model could reflect the assumed spherical filling of the outer regions rather than a true depletion of HC3N. The authors explicitly concede (Sect. 4.1): 'This one-dimensional model is not designed to yield a unique reconstruction of the source's intrinsic non-spherical structure.' Without a test that uses a geometry consistent with the observed barrel/EDE, the conclusion that abundance—rather than density structure or viewing geometry—is concentrated inward is not uniquely established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents new NOEMA interferometric observations of the 21 micron protoplanetary nebula IRAS 06530-0213, reporting the first spatially resolved detections of HC3N (J=24-23 and 26-25), C4H (J=47/2-45/2 and 45/2-43/2), and SiC2. The authors compare these maps with previously published CO and 13CO data and find that HC3N and C4H are concentrated in a compact inner, low-latitude torus-like structure, whereas CO and 13CO trace a larger barrel-like nebula. Using LTE column-density estimates, they derive fractional abundances and compare them with other carbon-rich evolved stars. To interpret the different spatial distributions, they perform 1D radiative transfer modeling with RATRAN. A constant-abundance model reproduces CO and 13CO but overproduces HC3N at large radii; a variable HC3N abundance peaked near 0.64 arcsec reproduces the observed radial profile. The paper concludes that HC3N must possess a more inwardly concentrated abundance distribution than CO/13CO, indicating that the inner equatorial density enhancement is chemically favorable for carbon-chain formation. The paper also compares Spitzer/IRS spectra of IRAS 06530-0213 and IRAS 23304+6147 and discusses possible links between the inner EDE and the 15.8 and 21 micron features.","tokens_in":20959,"tokens_out":6386,"duration_ms":70544,"significance":"The observational core of this paper is valuable. It adds a spatially resolved molecular inventory for a rare 21 micron PPN, and the compact HC3N/C4H morphology relative to CO/13CO is a robust, interesting observational result that will be useful for future studies of chemical differentiation in post-AGB envelopes. The abundance and column-density measurements, while model-dependent, expand the small sample of 21 micron sources with interferometric molecular maps. The central interpretive claim, however, rests on a simplified 1D spherical RATRAN model with several hand-set parameters and a fitted abundance profile. As the authors themselves note, this model is not designed to reconstruct the non-spherical structure. The claim that the observed morphology 'requires' an inwardly concentrated abundance distribution therefore goes beyond what the modeling actually demonstrates, and the paper would need to break the geometry-abundance degeneracy before the strong version of the conclusion can be accepted.","major_comments":[{"comment":"The central claim that reproducing the observed HC3N morphology 'requires' an inwardly concentrated abundance distribution is not established by the 1D spherical RATRAN model. In the true source, the outer equatorial sectors contain a thin barrel wall and tenuous lobes/halo, whereas in the spherical model lines of sight at those projected radii pass through full spherical shells, which will overproduce emission for any molecule with constant abundance. The failure of the constant-abundance HC3N model could therefore reflect the assumed spherical filling of the outer regions rather than a true chemical depletion at large radii. Moreover, the variable abundance profile is fitted to the very radial profile that defines the discrepancy, so the conclusion partly reduces to model input. The authors explicitly concede in §4.1 that the 1D model is not designed to yield a unique reconstruction. I","section":"§4.1, Fig. 10, Abstract, §5"},{"comment":"The physical structure adopted for all species is constrained only by fitting the CO radial profile with a constant CO abundance and hand-set parameters (n_H2,out = 2.0e5 cm^-3, T_out = 60 K, CO shell 0.45-0.87 arcsec). CO is optically thick, so its normalized radial profile is not a direct tracer of the volume density distribution, and many combinations of density normalization, temperature, and shell boundaries could match the same profile. The paper does not report an exploration of the acceptable parameter range, and it does not use 13CO (which is optically thin and is reproduced with constant abundance) to independently constrain the inner density structure in the region where HC3N peaks. Without such sensitivity tests, the statement that 'excitation effects or optical depth alone cannot explain the central concentration' is not fully supported: a spherical model with a steeper inne","section":"§4.1, Fig. 10"},{"comment":"The derived fractional abundances rely on several stacked assumptions: LTE and optically thin emission, an excitation temperature T_ex = 24 K from only two HC3N lines observed at different array configurations, N(H2) estimated from 13CO with X(CO) = 8e-4 and 12C/13C = 30, and the same T_ex applied to C4H and SiC2. The resulting HC3N abundance in this source is roughly two orders of magnitude higher than that reported for other 21 micron sources (IRAS 23304+6147 and IRAS 22272+5435). This difference may partly reflect different observational beams, assumed excitation temperatures, and analysis methods rather than a genuine chemical difference. The comparison in §4.2 ('21 micron sources tend to show a higher HC3N/SiC2 ratio') should therefore be treated with caution, or reanalyzed with consistent assumptions across sources. This does not affect the morphological result, but it does affect","section":"§3.4, Table 2"}],"minor_comments":[{"comment":"The variable-abundance HC3N model is described as 'confined between R_in = 0.30 arcsec and R_out = 0.65 arcsec with a peak at r = 0.64 arcsec'; since the peak is essentially at the outer boundary, please clarify the adopted functional form (e.g., a power-law with a central hole, or a Gaussian ring) and whether the peak position is a fitted parameter or a fixed boundary.","section":"§4.1"},{"comment":"The lower panels show abundance profiles X(r) but the axis labels and the precise curves are not described in the text. Please add a short description of the CO, 13CO, and HC3N abundance profiles (constant versus variable) and specify whether the plotted X(r) for CO and 13CO are constants or also radial functions.","section":"Fig. 10"},{"comment":"The orientation language is sometimes ambiguous: e.g., HC3N 'extends along the southeast-northwest axis' while the barrel direction is described elsewhere as PA=35 deg. Please state explicitly how the HC3N/C4H torus orientation relates to the CO barrel axes defined in Sun et al. (2025b).","section":"§3.1-3.2"},{"comment":"The comparison of the 15.8 and 21 micron features between IRAS 06530-0213 and IRAS 23304+6147 is qualitative and does not account for slit position angle, possible extended emission, or continuum fitting details beyond the statement in the text. The discussion is appropriately cautious, but a brief caveat about these systematic effects should be added.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a solid observational dataset and a clearly presented new detection of compact HC3N/C4H emission in a rare 21 micron PPN. My concern is not with the observations but with the central interpretation: the abstract and conclusions state that the modeling 'requires' an abundance gradient, while the model itself is 1D spherical and the authors concede it is not unique. The abundance profile is effectively a fitted free parameter, so the strong claim is not yet warranted. I would be willing to accept the paper after the authors (1) add a concrete test or discussion that addresses the geometry-abundance degeneracy (e.g., an axisymmetric model or a sensitivity study using 13CO), and (2) soften the wording of the central claim to match the model's actual scope. The abundance ratio discussion in §3.4/§4.2 also needs a caveat about heterogeneous assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: solid observational paper, overreaching interpretive claim. The new NOEMA detections of HC3N, C4H, and SiC2 toward IRAS 06530-0213 are genuine, and the maps convincingly show that the carbon chains sit in a compact inner equatorial structure while CO/13CO trace a larger barrel. That morphological difference is the real result, and it's useful. Column densities and upper limits are reported cleanly.\n\nThe soft spot is the RATRAN step. A 1D spherical model with hand-set n_H2,out, T_out, shell boundaries, and a CO-normalized structure is used to argue that reproducing the HC3N radial profile 'requires' an abundance concentration. The fixed-abundance model fails, and the variable-abundance model fits, but that only shows a variable abundance within that assumed geometry. In a barrel-like source, line-of-sight projection through the thin outer walls and the central torus can produce the same compactness for a constant abundance. The authors acknowledge the model isn't unique (Sect 4.1), so the word 'requires' in the abstract and conclusions is stronger than the evidence supports. This is fixable by softening the wording or adding a crude axisymmetric/3D test. Either way, the morphological result doesn't depend on the model.\n\nI also think the discussion of the 15.8/21 micron profiles is speculative but clearly flagged, so no problem there. The LTE abundances are approximate but standard for this context.\n\nBottom line: worth publishing after moderate revision. Send it to a good referee; the observations deserve to be in the literature. I'd be happy to cite it for the maps.","headline":"New NOEMA carbon-chain maps are solid and useful; the claim that an HC3N abundance gradient is 'required' overreaches the 1D spherical model.","tokens_in":21493,"tokens_out":1946,"would_cite":true,"duration_ms":21863,"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 authors claim that in IRAS 06530-0213, HC3N and C4H are chemically concentrated in an inner equatorial torus, not merely excited there, and that this zone is favorable for carbon-chain chemistry.","keywords":["protoplanetary nebulae","circumstellar matter","carbon-chain molecules","HC3N","C4H","radiative transfer","21 micron feature","NOEMA"],"falsifier":"Map HC3N in at least three rotational transitions at matched high resolution to derive a resolved excitation-temperature map. If T_ex varies strongly with radius and a 3D non-LTE model with the actual barrel geometry and inclination reproduces the observed central concentration while keeping HC3N abundance constant, the abundance-gradient claim would collapse. Alternatively, resolving C4H in two or more transitions and finding that its radial profile also requires a constant abundance would weaken the chemical-differentiation interpretation.","tokens_in":20388,"feed_emoji":"🔭","tokens_out":5129,"duration_ms":47668,"temperature":0.7,"pith_summary":"The paper uses NOEMA interferometric maps of the protoplanetary nebula IRAS 06530-0213 to show that HC3N and C4H emission is concentrated in a compact equatorial torus, while CO and 13CO trace a larger barrel-like outflow. A one-dimensional RATRAN radiative transfer model with a constant HC3N abundance over-predicts the outward extent of the emission; only a model with an HC3N abundance peaking near 0.64 arcsec matches the data. The authors conclude that excitation and optical-depth effects alone cannot explain the central concentration, so the inner equatorial density enhancement must be a chemically distinct zone with elevated carbon-chain abundance. This matters because the 21-micron emission feature, whose carrier is unidentified, appears only in carbon-rich protoplanetary nebulae, and the inner zone may be where its carrier forms or survives.","feed_headline":"Carbon-chain gas hugs the inner torus of IRAS 06530-0213","feed_subtitle":"HC3N peaks near 0.64 arcsec, marking the inner equatorial zone as chemically active and a candidate home for the 21-micron carrier.","key_machinery":"The central tool is a one-dimensional RATRAN Monte Carlo radiative transfer model, run for CO, 13CO, and HC3N with a common power-law density profile, homologous expansion, and a temperature power law. CO fixes the physical structure; then 13CO is reproduced with a constant abundance while HC3N requires an ad hoc radial abundance profile confined between 0.30 and 0.65 arcsec with a peak at 0.64 arcsec. The model's role is to test whether the observed central concentration can be explained by excitation or optical depth; its failure for constant HC3N abundance is the evidence for a true abundance gradient.","core_discovery":"On the nebula's own terms: HC3N and C4H trace an inner, low-latitude torus with a deprojected expansion velocity of about 10.6 km/s, while CO and 13CO outline the full barrel. The constant-abundance model fails to reproduce the HC3N radial profile, and the variable-abundance model with a peak at r = 0.64 arcsec succeeds. Hence the observed morphology requires HC3N to have a more inwardly concentrated abundance distribution compared with CO and 13CO, implying that the inner equatorial regions provide favorable conditions for carbon-chain chemistry. The paper also finds that IRAS 06530-0213 shows stronger 15.8 and 21-micron features and a broader red wing on the 15.8-micron band than the compa","pith_inferences":["If the HC3N gradient is real chemical differentiation, higher-order cyanopolyynes such as HC5N and HC7N should be even more centrally concentrated; that is a testable prediction for future NOEMA or ALMA observations.","The inferred peak radius of 0.64 arcsec is tied to the assumed spherical geometry; a 3D model that includes the barrel walls and inclination could shift or sharpen this peak, so the quantitative value should be treated with caution.","The comparison between the two sources suggests a test: if the 21-micron carrier is produced in the inner EDE, then among 21-micron sources the feature strength should correlate with EDE density or compactness, not just with central star temperature."],"forward_implications":["The inner equatorial density enhancement in IRAS 06530-0213 is a chemically active zone where carbon-chain molecules form or survive, not just a region where they are easier to excite.","Because all five spatially resolved 21-micron sources show equatorial density enhancements, these structures may be a common prerequisite for the 21-micron feature, pending a larger sample.","The enhanced 15.8 and 21-micron features in IRAS 06530-0213 relative to IRAS 23304+6147 indicate that hotter, denser inner environments strengthen the bands associated with large carbonaceous material.","The non-Lorentzian red wing of the 15.8-micron band in IRAS 06530-0213 suggests anharmonic hot-band emission from vibrationally excited PAHs, implying a higher-temperature environment that could also boost carbon-chain chemistry."],"fun_headline_variants":["HC3N and C4H trace inner equatorial torus in IRAS 06530-0213","Carbon-chain chemistry peaks near 0.64 arcsec in IRAS 06530-0213","Inner equatorial zones of IRAS 06530-0213 favor carbon-chain chemistry","Carbon-chain gas concentrated in inner torus of IRAS 06530-0213"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The argument leans on a one-dimensional spherical RATRAN model with a single power-law density profile and a CO-constrained shell standing in for the equatorial sector of a barrel-shaped, likely non-spherical nebula; the paper itself states this model cannot uniquely reconstruct the intrinsic non-spherical structure, so if the true 3D geometry, inclination, or optical-depth gradients differ, the inferred HC3N abundance concentration could be an artifact of the assumed physics","fun_headline_variants_meta":{"raw":{"variants":["HC3N and C4H trace inner equatorial torus in IRAS 06530-0213","Carbon-chain chemistry peaks near 0.64 arcsec in IRAS 06530-0213","Inner equatorial zones of IRAS 06530-0213 favor carbon-chain chemistry","Carbon-chain gas concentrated in inner torus of IRAS 06530-0213"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001263,"raw_usage":{"total_tokens":5028,"prompt_tokens":786,"completion_tokens":4242,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":4145}},"tokens_in":530,"tokens_out":4242,"duration_ms":26491,"temperature":1.0,"reasoning_tokens":4145,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T02:48:44.550931+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Map HC3N in at least three rotational transitions at matched high resolution to derive a resolved excitation-temperature map. If T_ex varies strongly with radius and a 3D non-LTE model with the actual barrel geometry and inclination reproduces the observed central concentration while keeping HC3N abundance constant, the abundance-gradient claim would collapse. Alternatively, resolving C4H in two or more transitions and finding that its radial profile also requires a constant abundance would weaken the chemical-differentiation interpretation.","supporting_citations":[],"review_version":1}