{"id":"aeae44cc-3aba-4f87-a9a5-2d1049b257ce","arxiv_id":"2507.01336","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Self-consistent simulations show that planet-carved gaps heat the disk midplane by tens of Kelvin, while dust rings cool it, shifting and multiplying the icelines of water, CO2, and CO compared to smooth-disk models.","lead":"This paper uses computer simulations of planet-forming disks to show that a forming giant planet heats the gas in the gap it carves out, shifting where volatile ices like CO can exist. Adding realistic dust grains to the models changes those temperatures by about 10 Kelvin and creates cooler dust rings, which matters for interpreting ALMA observations and for where planets can form.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §2.3 assumption T_gas = surface-area-averaged T_dust is load-bearing: the Model D 'shallower gap, +10 K' result in §3.1.2 has the sign that reverses if gas and dust decouple in the gap, as the authors acknowledge in §4.4.","rationale":"The reader and I converge on the same premise. The central claim has two layers: (i) planet-carved gaps are locally heated by tens of K relative to a smooth disk; (ii) adding multiple dust species makes gaps shallower and ~10 K hotter than gas-only feedback. Layer (i) is supported by a clear physical mechanism — reduced optical depth in the gap lets stellar radiation penetrate deeper — and is internally consistent with the MCRT maps, so I do not see a reason to attack it. Layer (ii) is the genuinely new quantitative result, and it rests entirely on converting dust temperatures to the gas temperature that controls the hydrostatic scale height. The paper's own §4.4 acknowledges the Facchini et al. (2018) counterexample. In a deep gap, the dust surface area is depleted, the dust-gas collisional coupling time lengthens, and gas can be colder than dust; if realized in these simulations, the sign of the feedback flips: lower T_gas → smaller h/r → deeper gap, which would erase or reverse the headline Model D vs Model G difference. The 10 K/25% numbers in §3.1.2 are therefore the least secure quantitative claims. I also considered the absence of viscous heating and possible non-convergence at low α; both are real caveats, but they affect the absolute temperature normalization and outer-disk details more than the sign of the multi-dust feedback, and the paper discusses them explicitly. The T_gas = T_dust equality is not tested at all. Because the concern is acknowledged but unquantified, the appropriate disposition is the same conditional acceptance the reader gave: the central scenario is plausible, but the headline quantitative comparison should be accompanied by either a two-temperature check or an explicit statement of its dependence on this equality. My recommended verdict is unchanged. I credit the paper for including a dust-species-number convergence test and a puffed-up-rim control, which strengthen the dust-ring cooling result, but neither test probes the gas-dust temperature coupling.","tokens_in":16325,"tokens_out":7110,"duration_ms":81112,"concrete_test":"Post-process the final snapshot of the 100 M_E at 10 au, α=1e-3 Model D run in a two-temperature calculation: solve the local gas energy balance from dust-gas collisional heating/cooling (e.g., the Facchini et al. 2018 framework or RADMC-3D gas-temperature mode) using the same gas density, dust densities, and opacities, instead of imposing Eq. (2). Compare the resulting T_gas(r,z) with the area-averaged T_dust in the gap midplane. If T_gas/T_dust < ~0.9 in the gap, recompute h/r and the iterated gap depth; if the Model D gap is no longer shallower than Model G, the §3.1.2 conclusion is not robust and the paper should report a bracketed range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's new quantitative claims — Model D gaps ~1.5× shallower and ~10 K hotter than Model G (§3.1.2, Fig. 3) — are produced by feeding the MCRT dust temperature back into the HD as the gas temperature: Eq. (2) averages over dust grain surface area, Eq. (3) vertically averages, and §2.3 states T_dust = T_gas. The authors themselves flag in §4.4 that this equality can break: Facchini et al. (2018) find T_gas/T_dust < 1 in deep-gap midplanes because the reduced dust surface area weakens collisional coupling. Since the aspect ratio h/r is set by T_gas, a cooler-than-dust gas in the gap lowers h/r and deepens the gap, which could counteract or reverse the headline comparison. The 10 K gap-temperature enhancement is the stated cause of the shallower gap; if T_gas < T_dust in the gap, both the magnitude and sign of that feedback are uncertain. The ice-line predictions also use T_mid (dust) as the gas-phase tracer, so volatile sublimation/freeze-out locations inherit the same decoupling error in low-density gap regions. This limitation is acknowledged but not quantified or propagated, making it the most load-bearing untested step in the central argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an iterative hydrodynamics + Monte Carlo radiative transfer framework (extending C23) in which gas and multiple dust species evolve, and the resulting dust temperature is fed back as the gas temperature. The authors show that planet-carved substructures alter the disk temperature relative to smooth disks: gap midplanes are warmer by tens of kelvin, dust-trap rings are cooler by several kelvin, and volatile icelines shift or multiply. They further compare the multi-dust model (Model D) with a gas-only/0.1-micron-dust model (Model G), reporting that Model D produces ~1.5x shallower gaps and ~10 K (~25%) higher gap temperatures, and they study the effect of alpha viscosity on temperature and icelines. An ALMA observational strategy is proposed.","tokens_in":16588,"tokens_out":5079,"duration_ms":51307,"significance":"The iterative coupling of hydrodynamics and radiative transfer is a genuine strength: the temperature changes are emergent outputs, not imposed by fitted parameters, and the paper extends the prior gas-only approach by including multiple dust sizes with dynamics in both HD and MCRT. If the quantitative results hold, they would imply that smooth-disk temperature assumptions are inadequate for interpreting gap/ring observations and volatile distributions. The comparison with C23 as a baseline is appropriate and non-circular. However, the central quantitative claims rely on the assumed equality between dust and gas temperatures in gap regions, which the authors themselves identify as fragile, and on the neglect of viscous heating; these issues must be addressed before the specific magnitudes (10 K, 1.5x, tens of kelvin) can be accepted.","major_comments":[{"comment":"The assumption T_gas = surface-area-averaged T_dust is load-bearing for the central Model D versus Model G result in Section 3.1.2. The authors note in Section 4.4 that Facchini et al. (2018) find T_gas/T_dust < 1 in deep-gap midplanes because the reduced dust surface area weakens collisional coupling. Since the aspect ratio h/r is set by T_gas, a cooler gas temperature in the gap would reduce h/r and deepen the gap, counteracting or even reversing the claimed shallower-gap trend. The magnitude and sign of the feedback are therefore uncertain. Please quantify this uncertainty, for example by repeating the iteration with a parametric T_gas/T_dust ratio or by implementing a two-temperature prescription in the gap region.","section":"2.3, Eq. (2); 4.4"},{"comment":"The omission of viscous heating may change the sign of the gap temperature perturbation in the inner disk. The authors cite Broome et al. (2023) showing that with viscous heating, T_mid in a Jovian gap at ~3 au can decrease by 20-30% relative to a gap-free model, which is opposite in sign to the stellar-heating-only result presented here. Since the abstract's 'midplane temperatures in gaps can increase by tens of kelvin' is a headline claim, the radial range over which this claim holds should be stated explicitly, and the qualitative conclusion should be qualified accordingly, particularly for the r_p = 4 au cases.","section":"3.1.2, Fig. 3; 4.4"},{"comment":"The claim that Model D gaps are about 1.5 times shallower than Model G gaps is presented as a general statement but is based on one representative case (100 M_Earth at 10 au). The authors note exceptions for eccentric-gap cases with 3 M_J, where Model G gaps are actually shallower than Model D. This exception is not quantified and weakens the generality of the headline comparison. Please provide gap-depth statistics across the full parameter grid (Table 1) or restrict the claim to the cases for which it holds.","section":"3.1.2"},{"comment":"The dust-ring cooling of 'several kelvin' is close to the expected uncertainty from dust opacity and grain-size choices: the 4-species versus 8-species test in Section 3.1.2 shows temperature differences up to 15% in most regions. The claim that dust rings 'create volatile freeze-out regions' should be softened or supported by a sensitivity test on the grain-size distribution, since the magnitude of the cooling is near the model's internal scatter.","section":"3.1.3, Fig. 5"}],"minor_comments":[{"comment":"The heading 'Prossessing between radiative transfer and hydro' contains a typo; it should be 'Processing'.","section":"2.3 heading"},{"comment":"The caption refers to the 'nineth row' of the table; this should be 'ninth row'.","section":"Table 1 caption"},{"comment":"The word 'locaitons' in the parameter-list sentence should be 'locations'.","section":"Section 2"},{"comment":"The heading 'Limits of our model' is spelled 'limtis' in the text; please correct the typo.","section":"Section 4"},{"comment":"The captions contain typos: 'differemt' and 'Miplane' in Fig. 6, and 'Comparions' in Fig. 7.","section":"Figures 6 and 7 captions"},{"comment":"The statement that the radiative transfer simulations do not include viscous heating appears in the results section; it would be better placed in Section 4.4 where the limitation is discussed in detail.","section":"3.2.1"},{"comment":"Please clarify that Eq. (3) produces a vertically averaged, density-weighted temperature that is not necessarily equal to the midplane temperature used later for iceline analysis; the difference between T_iterate and T_mid should be stated explicitly.","section":"2.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the iterative framework is a useful contribution. The main concern is that the headline quantitative claims (shallower gaps, 10 K differences, tens-of-K gap heating) rest on assumptions that are acknowledged in Section 4.4 but not quantitatively tested. I recommend major revision rather than rejection because the framework is promising and the limitations are disclosed; however, the authors should either provide additional simulations that relax the T_gas = T_dust assumption and include viscous heating, or substantially soften the quantitative claims in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. Verdict: this is a solid, incremental extension of Chen et al. 2023, worth engaging despite one load-bearing simplification.\n\nWhat is actually new: C23 ran gas-only HD and assumed a single 0.1 micron dust well mixed in MCRT. Here they put four dust sizes (0.1 micron to 1 mm) into both FARGO3D and RADMC-3D, let the dust drift and settle, and feed a surface-area-averaged dust temperature back into the gas hydrodynamics. The result is that multi-dust models give gas gaps about 1.5 times shallower and gap temperatures about 10 K higher than gas-only Model G. They also show mm dust rings cool the midplane by a few Kelvin and act as freeze-out zones, and that the ice-line response to viscosity is non-monotonic. The 8-species convergence run is a nice check; the paper is clear and the limitations section is honest, including the T_gas = T_dust caveat and the neglect of viscous heating.\n\nThe soft spot is exactly the one you flagged. Equation (2) sets T_gas equal to the dust surface-area average; that equality is what produces the +10 K and the shallower gap. In a deep gap, Facchini et al. (2018) find T_gas/T_dust < 1 in the midplane because the dust surface area drops. If the gas is cooler than the dust, the aspect ratio drops and the gap deepens, which could reverse the Model G vs Model D comparison. The authors acknowledge this in Section 4.4 but do not quantify it, so the central quantitative claim is conditional. The main qualitative point -- that a giant planet's gap can heat the midplane by tens of Kelvin and shift or multiply icelines -- does not depend on that equality and is consistent with earlier work, so it likely stands.\n\nThe other caveats are minor: no viscous heating (they discuss it, and it would mainly affect the inner few au), no photon-noise error bars, azimuthal averaging, and data only on request. The viscosity study is a useful parameter sweep, and the ALMA observability section is practical.\n\nWho is it for: people working on planet-disk interaction, disk thermal structure, and volatile delivery. It is a serious simulation paper with clear logic; send it to a good referee. In revision I'd ask for a sensitivity test with a decoupled gas temperature (or at least a simple correction) and a public data release of the temperature maps.\n\nNet: cite it, engage with it, but treat the 10 K Model D vs G numbers as provisional.","headline":"A careful multi-dust extension of C23 that makes the disk thermal-structure feedback richer; the headline 10 K gap-temperature shift is plausible but hinges on an acknowledged, unquantified T_gas = T_dust assumption.","tokens_in":17155,"tokens_out":2470,"would_cite":true,"duration_ms":28503,"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":"Planet-carved gaps heat protoplanetary disk midplanes by tens of Kelvin, shifting and multiplying volatile icelines while the heating feeds back to make the gaps shallower.","keywords":["protoplanetary discs","planet-disc interactions","radiative transfer","hydrodynamics","dust trapping","icelines","volatile distribution","disk temperature"],"falsifier":"A resolved measurement of gas temperature inside a planet-carved gap — for example from CO rotational line ratios or HD emission — that shows $T_{\\rm gas}$ below $T_{\\rm dust}$ in the gap midplane, or a hydrodynamics run that couples gas and dust temperatures separately and produces gaps deeper rather than shallower than a gas-only model, would refute the central feedback claim. The paper itself notes that Facchini et al. (2018) already found $T_{\\rm gas}/T_{\\rm dust} < 1$ in gap midplanes, so this test is the natural next step.","tokens_in":16087,"feed_emoji":"🪐","tokens_out":8585,"duration_ms":81639,"temperature":0.7,"pith_summary":"This paper argues that a protoplanetary disk carved by a giant planet is not thermally smooth: the planet's gap heats the midplane by tens of Kelvin, dust rings cool by several Kelvin, and the radii at which water, carbon dioxide, and carbon monoxide freeze out shift and multiply. To show this, the authors iterate between hydrodynamical simulations with four dust grain sizes and Monte Carlo radiative transfer, letting the temperature change the disk scale height and feed back on the gap depth. They find that including multiple dust species makes gas gaps about 1.5 times shallower and gap temperatures about 10 K (~25%) higher than in gas-only models, because millimeter grains trap cooling radiation. If correct, the disk's thermal structure is coupled to its substructure in a way that smooth-disk models miss, with consequences for where volatiles freeze out and for how molecular emission from gaps is interpreted.","feed_headline":"Planet-carved gaps heat disks and shift their icelines","feed_subtitle":"Heated gaps hollow out less deeply, and volatile ice lines shift and multiply, changing where ices freeze and planets assemble.","key_machinery":"The mechanism that carries the argument is an iteration loop between the FARGO3D hydrodynamics code and the RADMC-3D Monte Carlo radiative transfer code, run to 2000 planetary orbits with temperature feedback every 100 orbits. The central identity is Eq. (2): the gas temperature fed into the next hydrodynamics step is the dust surface-area-averaged temperature over the four grain sizes, computed from the dust temperatures and number densities in each cell. This identity is what couples the thermal structure to the grain-size-dependent opacity and dust distribution: in the gap, sparse dust lets stellar photons penetrate and heat the midplane, raising the aspect ratio and shallowing the gap; at the dust ring, millimeter grains raise the optical depth and cool the midplane, lowering the sublimation temperature of volatiles.","core_discovery":"The central claim is that planet-induced substructures and disk temperature must be computed self-consistently: a giant planet's gap heats the surrounding gas and dust significantly, and that heating raises the disk aspect ratio, which in turn limits how deep the gap can grow. In the authors' simulations, a 3 Jupiter-mass planet at 30 au raises the midplane temperature in the gap from about 30 K to about 60 K, while a 100 Earth-mass planet at 10 au raises it by about 10 K; dust rings at the outer gap edge, formed by pressure-bump trapping of millimeter grains, cool the midplane by several Kelvin and act as volatile freeze-out regions. Including multiple dust species (0.1 μm to 1 mm) rather than a single well-mixed grain size yields gas gaps roughly 1.5 times shallower and gap temperatures about 10 K (25%) higher than the gas-only model, because the larger grains make the disk marginally optically thick at long wavelengths and slow the escape of cooling radiation. The overall midplane ice distribution of H2O, CO2, and CO is similar in both models, but the hot-gap/cold-ring configuration creates a volatile sublimation and re-freeze-out cycle.","pith_inferences":["A consequence not explored in the paper: because the feedback strength depends on the grain size distribution, disks with ongoing dust growth or fragmentation may oscillate between shallow and deep gap states as grain sizes evolve, testable with multi-epoch observations of the same disk.","The shallow-gap result in the multi-dust model hinges on the gas-dust temperature equality; if the two thermally decouple in deep gaps, the gap-depth comparison between Model G and Model D could reverse, so the quantitative claim should be re-checked with a two-temperature coupling.","The proposed ALMA test — comparing CO iceline extent in a gap-hosting disk against a smooth-disk counterpart — could be applied to existing disk surveys, since the predicted outward shift of the CO iceline is a large, resolvable effect."],"forward_implications":["Volatile icelines in a planet-structured disk are not single, smoothly located radii: the CO iceline can multiply and shift outward, so comparing observed ring positions to iceline predictions from smooth temperature profiles is unreliable.","A deep H$_2$ gap located beyond the smooth-disk CO iceline can appear as a bright CO molecular ring rather than an emission gap in ALMA observations, complicating the interpretation of molecular emission gaps.","Dust rings at pressure bumps cool the midplane by a few Kelvin and act as volatile freeze-out zones, while the adjacent hot gap sublimates the same volatiles, creating a cycle that can concentrate solids for planetesimal growth.","Disk viscosity affects midplane temperature through turbulent dust mixing, but gap opening counteracts this, so iceline locations and the number of ice regions do not vary monotonically with $\\alpha$.","The predicted temperature changes and outward CO iceline shift in a gap-hosting disk are testable with ALMA CO intensity maps or spectral lines at roughly 30 au resolution."],"supporting_citations":[{"why":"Provides the previous gas-only iteration model (Model G) that this paper extends by adding multiple dust species; supplies the baseline gap-depth and temperature results being compared.","marker":"Chen et al. (2023)"},{"why":"Supplies the multi-species FARGO3D hydrodynamics framework with multiple dust grain sizes, including dust feedback and diffusion, used for all Model D runs.","marker":"Weber et al. (2019)"},{"why":"Provides the steady-state dust vertical distribution (their Eq. 19) used to extend 2D dust surface densities into 3D for the radiative transfer.","marker":"Fromang & Nelson (2009)"},{"why":"Cited as the key limitation: shows T_gas/T_dust < 1 in gap midplanes, questioning the paper's assumption that gas temperature equals surface-area-averaged dust temperature.","marker":"Facchini et al. (2018)"},{"why":"Establishes the dust-trapping mechanism at pressure maxima that produces the mm dust rings whose cooling is a central result.","marker":"Pinilla et al. (2012a,b)"},{"why":"Previous finding that dust rings can act as volatile freeze-out regions, which this paper's hot-gap/cold-ring configuration invokes.","marker":"Alarcón et al. (2020)"},{"why":"Independent simulation work that found a temperature drop at dust rings with an assumed Gaussian density profile; this paper confirms the drop with self-consistent hydrodynamics.","marker":"Zhang et al. (2021)"},{"why":"Provides comparison for radiative transfer with prescribed gap density profiles and viscous heating, used to frame the temperature deviations in gaps.","marker":"Broome et al. (2023)"}],"fun_headline_variants":["Planet gaps heat midplane, shifting and multiplying icelines","Hot gaps and cold dust rings reshape disk volatile ice zones","Self-consistent simulations show planet-carved gaps warm disks","Giant planets cause tens of Kelvin gaps, moving ice lines"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire feedback loop rests on treating the surface-area-averaged dust temperature as equal to the gas temperature in every cell; if the gas in a deep gap is actually cooler than the dust, the gap would deepen instead of shallowing, reversing the paper's headline result about shallower gaps in the multi-dust model.","fun_headline_variants_meta":{"raw":{"variants":["Planet gaps heat midplane, shifting and multiplying icelines","Hot gaps and cold dust rings reshape disk volatile ice zones","Self-consistent simulations show planet-carved gaps warm disks","Giant planets cause tens of Kelvin gaps, moving ice lines"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00058,"raw_usage":{"total_tokens":2811,"prompt_tokens":1100,"completion_tokens":1711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":1654}},"tokens_in":716,"tokens_out":1711,"duration_ms":16168,"temperature":1.0,"reasoning_tokens":1654,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:54:16.699770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A resolved measurement of gas temperature inside a planet-carved gap — for example from CO rotational line ratios or HD emission — that shows $T_{\\rm gas}$ below $T_{\\rm dust}$ in the gap midplane, or a hydrodynamics run that couples gas and dust temperatures separately and produces gaps deeper rather than shallower than a gas-only model, would refute the central feedback claim. The paper itself notes that Facchini et al. (2018) already found $T_{\\rm gas}/T_{\\rm dust} < 1$ in gap midplanes, so this test is the natural next step.","supporting_citations":[{"cited_title":"Self-consistent ring model in protoplanetary disks: temperature dips and substructure formation","cited_arxiv_id":"2110.00858","evidence_quote":"Independent simulation work that found a temperature drop at dust rings with an assumed Gaussian density profile; this paper confirms the drop with self-consistent hydrodynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides comparison for radiative transfer with prescribed gap density profiles and viscous heating, used to frame the temperature deviations in gaps."}],"review_version":1}