{"id":"8affff70-cb18-4726-9541-c827af960570","arxiv_id":"2411.12418","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"SO2 emission is thermalised near the midplane at about 55 K, while CH3OH and H2CO are sub-thermally excited in higher disk layers, so their rotational temperatures underestimate the kinetic temperature.","lead":"Astronomers mapped the temperatures of three gas molecules in the lopsided dust ring around the young star Oph-IRS 48 and found that two of them are not in thermal balance with the surrounding gas. The results suggest the molecules sit at different heights in the disk and help explain why some species spread further around the ring than others.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"RADEX thermalization grids (Fig. 5) assume N=10^12 cm^-2, but the paper later argues CH3OH and H2CO are optically thick with beam-diluted column densities up to ~500 times higher; line trapping at such columns could thermalize the lines and invalidate the sub-thermal excitation claim.","rationale":"The reader identified the DALI model density structure as the weakest assumption. I agree that this is one important dependency, but the most load-bearing and immediately concerning element is the internal choice of column density in the RADEX thermalization grids. The paper explicitly adopts N=10^12 cm^-2 to guarantee optically thin emission, yet the same paper argues that the emission is likely optically thick and that the observed column densities are underestimated by a factor of ~500 due to beam dilution (Section 5.3.1). Radiative trapping at these higher column densities can thermalize lines at substantially lower volume densities, meaning the Figure 5 contours may be systematically shifted to higher densities than physically appropriate. If the corrected contours place CH3OH and H2CO in the thermalized regime at the DALI-predicted densities of the elevated layers, then the central claim is not merely uncertain but internally contradicted. The concrete test I propose is straightforward and decisive: recompute the thermalization grids at realistic column densities. This is a computational check that can be done without new observations, and it directly tests the logical bridge from the RADEX models to the sub-thermal interpretation. I therefore retain the reader's CONDITIONAL verdict, but with a sharper, more fundamental concern than the external model dependency.","tokens_in":30474,"tokens_out":6038,"duration_ms":60810,"concrete_test":"Re-run the RADEX calculations of Section 5.1 for A- and E-type CH3OH and ortho/para H2CO with column densities of N=10^14, 10^15, and 10^16 cm^-2 (spanning the observed and beam-dilution-corrected values) while keeping the same Tkin and n_H2 grids. Recompute the Delta-T contours in Figure 5 and compare the density at which Delta-T=10% to the DALI densities at z/r~0.17-0.25 (n_H2~10^7-10^8 cm^-3). If the thermalization density drops below the DALI densities for any of these columns, the sub-thermal excitation conclusion is invalid and the vertical layering interpretation must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that CH3OH and H2CO are sub-thermally excited (Section 5.1) rests on the RADEX grids in Figure 5, which are computed for a fixed column density of N=10^12 cm^-2, explicitly chosen to ensure optically thin emission. This is inconsistent with the paper's own later argument (Section 5.3.1) that the emission is likely optically thick and that the observed column densities are beam-diluted by a factor of ~500. The measured peak column densities are already N(SO2)=4.8e13 cm^-2 and N(CH3OH)=1.2e14 cm^-2; after beam dilution, CH3OH would be ~6e16 cm^-2. At such high column densities, radiative trapping increases the radiation field within the lines and can thermalize the excitation at volume densities far below what the optically thin N=10^12 grid predicts. If the true thermalization boundary shifts to n_H2 < 10^7 cm^-3 for the emitting layers, then the conclusion that CH3OH and H2CO are sub-thermal at z/r~0.17-0.25 fails, and the derived vertical layering (SO2 deep, CH3OH/H2CO high) collapses. The scatter in the CH3OH rotational diagram would then need an alternative explanation (e.g., continuum oversubtraction, which the paper itself discusses in Section 5.3.2). This is a load-bearing internal inconsistency, independent of the DALI density structure, and it is not addressed in the paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analysis of ALMA Band 7 observations of SO2, CH3OH, and H2CO line emission in the asymmetric dust trap around Oph-IRS 48. Using 13, 22, and 7 transitions of the three species, the authors perform pixel-by-pixel rotational diagram analyses for SO2 and CH3OH and use line ratios with RADEX non-LTE calculations for H2CO. They derive rotational temperatures of ~55 K for SO2 and ~125 K for CH3OH, and temperature maps for H2CO in the range ~150-350 K. The well-behaved SO2 rotational diagram indicates thermalized emission near the midplane at densities above ~10^8 cm^-3, while the scattered CH3OH diagram and the H2CO line ratios are interpreted as sub-thermal excitation at elevated layers (z/r ~0.17-0.25), meaning the rotational temperatures underestimate kinetic temperatures. The authors also compute desorption, freeze-out, photodissociation, and turbulent mixing timescales using the DALI model of Leemker et al. (2023), and discuss implications for vertical layering, azimuthal extents, and a possible narrow radial 'sliver' for optically thick H2CO emission.","tokens_in":30814,"tokens_out":9136,"duration_ms":77816,"significance":"If the conclusions hold, the paper provides a valuable demonstration that multi-molecule, multi-transition ALMA analyses can constrain vertical thermal and chemical structure in low-mass planet-forming disks, and that non-LTE excitation is important for molecules such as CH3OH and H2CO in such disks. The study uses a large set of lines per molecule, applies a consistent pixel-by-pixel methodology, and presents maps of temperatures and column densities that can be compared with models. The paper is transparent about its limitations, including the reliance on a thermochemical model for gas densities and the ambiguity in optical depth. The proposed 'sliver' geometry and the predicted vertical ordering (SO2 deep, CH3OH/H2CO high) are falsifiable with higher angular resolution observations, which the authors explicitly identify. The analysis is based on public ALMA data and the methods are standard, making the work reproducible in principle.","major_comments":[{"comment":"The RADEX thermalization grids in Figure 5 are computed at a fixed column density of N=10^12 cm^-2 to ensure optically thin emission, yet Section 5.3.1 argues that the observed column densities are beam-diluted by a factor of ~500 and that the emission (especially H2CO, and likely CH3OH) is optically thick. At the implied column densities (e.g., N(CH3OH) ~6x10^16 cm^-2), line trapping can raise the radiation field within the lines and shift the thermalization boundary to volume densities well below those shown in Figure 5. The sub-thermal excitation conclusion for CH3OH and H2CO is therefore not robust against the column density uncertainty; the authors should re-run the RADEX grids over a range of column densities (or at the optically thick limits) and show that the key contours remain above the DALI densities, or else qualify the conclusion.","section":"5.1 and 5.3.1"},{"comment":"The H2CO temperature map is stated to be unconstrained at the peak-flux pixels: for the ortho and para line pairs (5,3-4,3,2 / 1,5-4,4 and 4,2-4,1 / 0,5-4,4), the observed ratios exceed the maximum RADEX grid values, and the peak position of the 5,0,5-4,0,4 transition falls in this unconstrained region (Appendix C). Consequently, the quoted T~150-350 K range and the associated claim that H2CO is sub-thermally excited at z/r~0.17-0.25 are not directly supported at the position of peak H2CO emission. The authors should either quantify the fraction of affected pixels or show that the sub-thermal conclusion still holds for the temperature range allowed by the unconstrained ratios (e.g., up to 500 K).","section":"4.3 and Appendix C"},{"comment":"The large scatter in the CH3OH rotational diagram is attributed to sub-thermal excitation, but the paper itself, in Section 5.3.2, presents evidence that the CH3OH line fluxes are affected by continuum oversubtraction or continuum blocking (see channel maps, Figure G.1). Such effects would also introduce scatter into the rotational diagram, independent of sub-thermal excitation. The paper does not quantify the expected magnitude of these systematic effects on the derived Trot. The sub-thermal interpretation would be strengthened by a demonstration that the residual scatter, after excluding or correcting pixels affected by the continuum feature, still leads to a similar conclusion.","section":"5.1 and 5.3.2"},{"comment":"All quantitative conclusions about sub-thermal excitation and the timescales rely on the gas density and temperature structure from the DALI model of Leemker et al. (2023), which is not independently validated for the IRS 48 disk. The thermalization boundary in Figure 5 is compared directly to the model's density profiles; if the actual densities in the CH3OH and H2CO emitting layers were higher by a factor of a few, both species could be thermalized. The paper would be considerably strengthened by a sensitivity analysis that perturbs the model densities (e.g., by 0.3-0.5 dex) and re-evaluates the crossing of the thermalization contours, and by a similar robustness check for the timescale comparisons in Section 6.","section":"5.1 and 6.1"}],"minor_comments":[{"comment":"The word 'azimtuhal' appears for 'azimuthal' in the abstract and in Section 6.2; please correct the spelling.","section":"Abstract and Section 6.2"},{"comment":"The equation for the dilution factor (DF^-1 = Omega2_source / (Omega2_beam + Omega2_source)) is not numbered; consider adding an equation number for clarity.","section":"Section 5.3.1"},{"comment":"The expression 's 2nsEb,i / pi2mi' for the vibrational frequency has ambiguous parentheses; please write it with explicit brackets, e.g., sqrt(2 n_s E_b,i / (pi^2 m_i)).","section":"Section 6.1, Eq. (1)"},{"comment":"The notation 'Delta V_line = sqrt(Delta V^2_thermal + Delta V^2_turbulence)' is clear, but consider defining Delta V_thermal explicitly as the thermal line width for each molecule in the text, rather than only giving a typical value.","section":"Section 3.2"},{"comment":"In the first paragraph, the phrase 'the temperature derived for these pixels is limited by the upper range of 500 K' could be clearer; for example, 'the ratios exceed the grid maximum, so only a lower limit of about 500 K can be inferred for these pixels.'","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The main technical concern is the internal inconsistency between the optically thin RADEX grid (N=10^12 cm^-2) used for the sub-thermal excitation conclusion and the paper's own later argument for optically thick, beam-diluted emission at column densities up to ~500 times higher. This issue is fixable with additional RADEX calculations and a sensitivity discussion, but it must be addressed before publication. The paper is otherwise well presented and the proposed vertical layering scenario is interesting and falsifiable. The reliance on a self cited thermochemical model is common in the field but should be acknowledged as a model-dependent step; a sensitivity test would improve confidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — quick take on arXiv:2411.12418 (Temmink et al., IRS 48 molecular line emission).\n\nThe new thing here is real: 13 SO2 and 22 CH3OH transitions from Cycle 8 ALMA data, analysed pixel-by-pixel, plus H2CO line ratios. The SO2 result is solid — well-behaved rotational diagram, T=55K, optically thin, likely from a deep layer. The CH3OH and H2CO sub-thermal excitation story is the headline, and it is the fragile part.\n\nThe authors do good housekeeping: they document blending, give upper limits, show uncertainty maps, and the timescale analysis is a thoughtful addition. The sliver scenario for optically thick H2CO is testable and appropriately flagged as tentative.\n\nThe soft spot is load-bearing. The RADEX grids in Figure 5 are computed at a fixed column N=10^12 cm^-2 to keep the lines optically thin. But later in the paper, from the H2CO/H2^13CO isotopologue ratio and their own optical depth grids, they conclude the emission is likely optically thick and beam diluted by a factor ~500. The actual column densities would then be ~10^16-17 cm^-2. At those columns, radiative trapping can thermalize the lines at much lower volume densities than the optically-thin grid predicts. The paper does not run RADEX at realistic columns, so the thermalisation boundary in Figure 5 may be misplaced. If the boundary falls below the emitting-layer densities, the sub-thermal conclusion collapses. The scatter in the CH3OH rotational diagram — the other pillar of the sub-thermal argument — can also be produced by optical depth or continuum oversubtraction, which the paper acknowledges elsewhere. So the two central claims rest on an internal inconsistency that is acknowledged piecewise but never resolved.\n\nSecond, the H2CO temperature map is unconstrained at the peak-flux pixels because the observed line ratios exceed the RADEX grid; the quoted 150-350K is honest, but the hottest region is where the action is.\n\nThe paper is still worth a serious referee: new data, careful analysis, and good questions. But the sub-thermal claim needs to be reworked with RADEX runs at actual columns and better H2CO temperature constraints. I would not accept as is; I'd send back for major revision. The SO2 part alone could be a short paper; the full story needs the excitation model fixed.","headline":"Careful new ALMA maps of IRS 48, but the sub-thermal excitation claim rests on RADEX grids that contradict the paper's own optically thick sliver scenario.","tokens_in":31424,"tokens_out":6219,"would_cite":true,"duration_ms":57409,"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 establishes that SO2 emission in the Oph-IRS 48 dust trap is thermalised at about 55 K near the midplane, while CH3OH and H2CO are sub-thermally excited in higher disk layers, so their rotational temperatures of roughly 125 K…","keywords":["astrochemistry","protoplanetary disks","dust traps","sub-thermal excitation","rotational diagram analysis","non-LTE radiative transfer","ALMA observations","IRS 48"],"falsifier":"Measure the kinetic temperature directly in the CH3OH and H2CO emitting layers, for example by observing a line with a much higher upper-level energy and critical density than those studied here, or by resolving the vertical emission height with channel maps at sub-0.1\" resolution and comparing with the rotational temperatures.","tokens_in":30278,"feed_emoji":"🔭","tokens_out":5363,"duration_ms":50824,"temperature":0.7,"pith_summary":"The paper aims to pin down where in the vertically layered Oph-IRS 48 disk the molecules SO2, CH3OH, and H2CO actually emit, and what that implies for the gas temperature and chemistry. Its central claim is that SO2 emission is thermalised in a deep layer near the midplane at about 55 K, while CH3OH and H2CO are sub-thermally excited in higher layers at z/r ≈ 0.17–0.25, so their LTE rotational temperatures of ~125 K and 150–350 K underestimate the true kinetic temperature. If this is right, low-mass disks like IRS 48 cannot be characterised with simple LTE rotational diagrams for these molecules, and the vertical ordering of molecular emission is CH3OH/H2CO highest, 13CO below, and SO2 deepest. The paper also argues that short photodissociation timescales explain the compact azimuthal extent of CH3OH and the wider extent of SO2, while H2CO's wider extent requires an additional gas-phase formation reservoir, and that the H2CO emission originates from a narrow radial sliver about 10 au wide at the inner edge of the dust trap.","feed_headline":"Two molecules in a dusty disk are colder than the gas around them","feed_subtitle":"In the IRS 48 dust trap, methanol and formaldehyde are sub-thermally excited; only SO2 traces the midplane.","key_machinery":"The load-bearing tools are the rotational diagram analysis (plotting upper-level column densities against upper-level energies and fitting a straight line to extract rotational temperature and total column density under LTE) and RADEX non-LTE grids that compute how much the rotational temperature deviates from the kinetic temperature across density and temperature space. The scatter in the CH3OH rotational diagram is the diagnostic sign of sub-thermal excitation, and the RADEX grids place the thermalisation boundary near n ≈ $10^{8}$ $cm^{-3}$ for SO2 and H2CO and above $10^{9}$ $cm^{-3}$ for CH3OH. The DALI thermochemical model of the IRS 48 dust trap supplies the gas density, gas temperature, and dust temperature structure that locates these boundaries in the disk and is also used to compute desorption, freeze-out, photodissociation, and turbulent-mixing timescales.","core_discovery":"Using 13, 22, and 7 transitions of SO2, CH3OH, and H2CO, the paper constructs pixel-by-pixel rotational diagrams and temperature maps. SO2 gives a well-behaved rotational diagram with T = 54.8 ± 1.4 K, which the paper reads as thermalised emission from a layer within about 5 au of the midplane where gas densities reach roughly $10^{8}$ $cm^{-3}$. CH3OH, by contrast, gives T = 125.5(+3.7/−3.5) K but with large scatter in the rotational diagram, and H2CO line ratios imply T ≈ 150–350 K; non-LTE RADEX grids show that at the model densities of the dust trap both CH3OH and H2CO are sub-thermally excited, meaning their derived rotational temperatures sit below the kinetic temperature. The paper further finds that the SO2 temperature map hints at a radial gradient but that no vortex-induced azimuthal temperature variation is present, and that photodissociation timescales, not turbulent mixing, best explain the observed azimuthal extents.","pith_inferences":["If the sub-thermal excitation picture is right, then existing and future LTE-based chemical surveys of low-mass disks may need re-interpretation: kinetic temperatures could be substantially higher than the rotational temperatures used to infer freeze-out, desorption, and volatile delivery to forming planets.","The sliver argument makes a sharp, testable prediction: high-resolution ALMA observations with beams smaller than about 8 au should resolve a narrow radial band of H2CO and CH3OH emission at the inner edge of the dust trap rather than a broad molecular layer.","The claim that H2CO is sustained by gas-phase formation implies that the azimuthal extent of H2CO should be spatially correlated with photodissociation products of CH3OH and H2O, such as CH3 and OH; mapping those radicals could confirm or refute the proposed reservoir.","The inferred ~2 K upper limit on vortex temperature imprints could be sharpened by observing a midplane tracer at higher spatial resolution, which would directly test whether the anticyclonic vortex motion is actually generating a detectable thermal signature."],"forward_implications":["For low-mass disks, LTE rotational-diagram temperatures for CH3OH and H2CO will systematically underestimate kinetic temperatures, so non-LTE analyses are needed to derive reliable gas temperatures and abundances.","SO2 traces the deepest molecular layer probed here, a region within about 5 au of the midplane (z/r < 0.1) at roughly 55 K, while CH3OH and H2CO trace elevated layers near z/r ≈ 0.17–0.25.","Photodissociation timescales of under a year in the elevated emitting layers explain why CH3OH is azimuthally compact and why SO2 can spread over roughly a quarter of an orbit, but they cannot explain H2CO's wider extent without invoking gas-phase formation from photodissociation products of CH3OH and H2O.","If H2CO is optically thick, as the low isotopologue ratio suggests, its emission must come from a narrow radial sliver about 10 au wide at the inner edge of the dust trap, and its measured column density is likely beam-diluted by a factor of order 500.","Any vortex-induced temperature variation in the midplane must be smaller than about 2 K, since the SO2 temperature map shows no azimuthal pattern despite its small pixel-to-pixel uncertainties."],"supporting_citations":[{"why":"Supplies the H2CO line observations and disk-integrated line ratios that the pixel-by-pixel temperature analysis extends.","marker":"van der Marel et al. (2021)"},{"why":"Reports the initial SO2 and CH3OH detections and the molecular inventory that this study draws its transition list from.","marker":"Booth et al. (2021)"},{"why":"Provides the Cycle 8 Band 7 data reduction, self-calibration, and continuum subtraction, and gives the H2CO/H2^13CO isotopologue ratio used in the optical-depth argument.","marker":"Booth et al. (2024)"},{"why":"Supplies the DALI thermochemical model of the dust-trap gas density and temperature structure used to locate thermalisation boundaries and compute all timescales.","marker":"Leemker et al. (2023)"},{"why":"Provides the RADEX non-LTE radiative transfer code used for the excitation grids and line-ratio calculations.","marker":"van der Tak et al. (2007)"},{"why":"Establishes that CH3OH lines with higher critical densities remain sub-thermally excited at densities near 10^7 cm^-3, which is the interpretive basis for the CH3OH scatter.","marker":"Johnstone et al. (2003)"},{"why":"Provides the standard rotational diagram method that underlies the LTE temperature and column density extractions.","marker":"Goldsmith & Langer (1999)"},{"why":"Supplies the 13CO J=6-5 data and the disk mass estimate used for the brightness temperature analysis and fractional abundance calculations.","marker":"van der Marel et al. (2016)"}],"fun_headline_variants":["Two molecules in a dusty disk are colder than the gas","Only SO2 traces the midplane in the IRS 48 dust trap","Sub-thermal methanol and formaldehyde reveal disk layering","Photodissociation, not turbulence, shapes molecular gas in IRS 48","Vortex leaves no mark on temperatures in Oph-IRS 48"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The DALI model of Leemker et al. (2023) correctly predicts the gas density and temperature structure of the IRS 48 dust trap, and it is used both to locate the thermalisation boundary for the RADEX comparison and to compute every desorption, freeze-out, photodissociation, and turbulent-mixing timescale.","fun_headline_variants_meta":{"raw":{"variants":["Two molecules in a dusty disk are colder than the gas","Only SO2 traces the midplane in the IRS 48 dust trap","Sub-thermal methanol and formaldehyde reveal disk layering","Photodissociation, not turbulence, shapes molecular gas in IRS 48","Vortex leaves no mark on temperatures in Oph-IRS 48"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000594,"raw_usage":{"total_tokens":2908,"prompt_tokens":1197,"completion_tokens":1711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":813,"completion_tokens_details":{"reasoning_tokens":1621}},"tokens_in":813,"tokens_out":1711,"duration_ms":12583,"temperature":1.0,"reasoning_tokens":1621,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:32:32.538179+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the kinetic temperature directly in the CH3OH and H2CO emitting layers, for example by observing a line with a much higher upper-level energy and critical density than those studied here, or by resolving the vertical emission height with channel maps at sub-0.1\" resolution and comparing with the rotational temperatures.","supporting_citations":[{"cited_title":"An ALMA molecular inventory of warm Herbig Ae disks: II. Abundant complex organics and volatile sulphur in the IRS 48 disk","cited_arxiv_id":"2402.04002","evidence_quote":"Provides the Cycle 8 Band 7 data reduction, self-calibration, and continuum subtraction, and gives the H2CO/H2^13CO isotopologue ratio used in the optical-depth argument."},{"cited_title":"S., van Dishoeck, E","cited_arxiv_id":null,"evidence_quote":"Supplies the DALI thermochemical model of the dust-trap gas density and temperature structure used to locate thermalisation boundaries and compute all timescales."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that CH3OH lines with higher critical densities remain sub-thermally excited at densities near 10^7 cm^-3, which is the interpretive basis for the CH3OH scatter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard rotational diagram method that underlies the LTE temperature and column density extractions."}],"review_version":1}