{"id":"feec56ec-20d8-4a31-a4da-0c31b23d5ac6","arxiv_id":"2501.14858","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Monte Carlo post-processing of radiation hydrodynamical protoplanetary disk models shows typical 10 percent temperature agreement, with systematic FLD/MCRT discrepancies that significantly alter predicted NIR and submm fluxes.","lead":"The paper compares gas temperatures from radiation hydrodynamics simulations with flux-limited diffusion against Monte Carlo radiative transfer post-processing in protoplanetary disks with accreting planets. It finds typical 10 percent agreement but systematic discrepancies up to 40 percent near the photosphere, and shows these shift predicted near-infrared and submillimeter fluxes by up to 100 percent in some regions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Non-iterative MCRT opacity fixing (Eq. 12) is validated only for the reference model and can bias TMC near sublimation fronts, so the ~10% agreement claim is not yet securely established.","rationale":"The reader identified the non-iterative MCRT reference as the weakest assumption, and my reading agrees: the central claim's credibility depends on TMC being a faithful radiative-equilibrium solution, and the fixed-opacity approximation is the least secure part of that chain. The concern is concrete and testable. It is not fatal: the zone-agreement check provides some support, and the limitation is openly acknowledged, so conditional acceptance remains appropriate. The reader's other concerns (flux-error independence, missing code/data) are secondary and already justify the CONDITIONAL verdict; my concern strengthens the condition by requiring a quantitative sensitivity check of the reference solution. Honest non-finding was considered: the paper is careful, and the temperature-comparison methodology is otherwise sound. However, the absence of any iterative MCRT verification for the high-resolution, high-mass cases leaves a real, identifiable risk that the reported agreement level is partly an artifact of frozen opacities. Thus I recommend no change to the reader's CONDITIONAL verdict, with the added condition that the opacity-fixing bias be explicitly bounded.","tokens_in":32908,"tokens_out":8244,"duration_ms":74561,"concrete_test":"Run the MCRT temperature calculation with opacity iterations for the two most demanding models (N3 300 M⊕ A32, and the reference N1 300 M⊕ A32): after each iteration, update the cell-wise dust-to-gas ratio using κR(TMC,prev) instead of κR(TRHD) in Eq. (12), and recompute TMC; repeat until the cell classification and temperature converge. Compare the iterated TMC with the reported fixed-opacity TMC in the high-resolution region. If the mean relative change exceeds a few percent, or if the dust-dominated/sublimation classification shifts in more than ~2% of cells, the non-iterative reference is not reliable and the 10% agreement claim must be re-evaluated; if the change is negligible, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that RHD-FLD and MCRT temperatures agree to ~10% in the high-resolution region. This comparison treats TMC as the reference solution, but TMC is computed with opacities frozen at κR(TRHD) via the cell-wise dust-to-gas adjustment in Eq. (12). The approach is valid only if the radiative-equilibrium temperature is close enough to TRHD wherever opacity is temperature-sensitive. The paper's a posteriori check reports 98–100% zone-agreement, but only for the reference model and only as a coarse dust-dominated/sublimation classification. It does not verify that within the dust-dominated zone the Rosseland (or Planck/frequency-dependent) opacities are accurately set, nor does it bound the resulting temperature error. The RHD Rosseland opacity is steep near the silicate sublimation temperature (Tevp, Eq. 11, with ΔTevp=200 K), and the N3 300 M⊕ model reaches or exceeds Tevp in the planetary cell. For such cells, a 10% temperature offset in the reference solution could shift the opacity by a large factor, and the non-iterative MCRT would then produce a TMC that is systematically biased. Because the headline statement is a quantitative agreement level, an uncontrolled systematic in the reference solution is load-bearing: the 10% figure could either overstate or understate the true FLD-MCRT agreement. The authors themselves list consistent Planck opacities for evaporation regions as a future improvement (Sect. 4.2), confirming the issue is live.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper quantifies how closely gas temperatures from radiation-hydrodynamical (RHD) flux-limited-diffusion simulations of protoplanetary disks with accreting planets agree with temperatures recomputed by Monte Carlo radiative transfer (MCRT) post-processing. The authors compare 18 three-dimensional models spanning two planet masses, three resolutions, and three accretion timescales, embedding each local high-resolution RHD model in an axisymmetric global disk. They report typical relative temperature differences of about 10% in the high-resolution region, with larger differences at the photosphere and in far-outer regions, and they attribute these differences to known limitations of FLD (angle and frequency averaging, missing scattering, optical-depth limiter) and of MCRT (missing PdV work and internal-energy advection). Controlled axisymmetric experiments isolate specific mechanisms, and synthetic flux maps at VIS, NIR, and submillimeter wavelengths translate the temperature differences into observable flux errors.","tokens_in":33136,"tokens_out":10417,"duration_ms":94651,"significance":"The paper provides a valuable, quantitative benchmark for a widely used modeling pipeline: RHD simulations of planet-disk interaction followed by MCRT post-processing. Its strengths are the broad parameter scan (18 models, three resolutions, two planet masses, three accretion parameters), the deliberate matching of opacities and density fields between the two schemes, the controlled axisymmetric experiments that separate irradiation, viscosity, frequency dependence, and the optical-depth limiter, and the flux-error analysis that connects temperature discrepancies to observable quantities. If the central agreement figures survive scrutiny, the paper gives practical error estimates for synthetic observations of accreting planets and a clear roadmap for improving RHD temperature solvers.","major_comments":[{"comment":"The bullet 'For all cells in the high-resolution region, the temperature estimates agree across all tested resolutions and accretion parameters by |δT| ≲ 0.1−0.13' is contradicted by the paper's own reference-model results. In Fig. 3 (upper right panel), the midplane relative difference reaches about 1/3 near r ≈ 4–5 au, which is inside the high-resolution region, and the radial profile in Fig. 4 shows negative relative differences of order −0.2 to −0.3 in the green-shaded high-resolution region outside the planet. The abstract itself acknowledges values exceeding 40% at the photosphere. The 'all cells' claim should be replaced with a quantile-based statement (e.g., the 90th or 99th percentile), or the photosphere and outer regions should be explicitly excluded. As written, the contradiction between the headline summary and the displayed profiles weakens the central quantitative message.","section":"Section 4, first bullet list; Figs. 3 and 4"},{"comment":"The non-iterative opacity-fixing procedure is load-bearing: the MCRT temperature TMC is computed with opacities fixed at κR(TRHD) via the cell-wise dust-to-gas adjustment of Eq. (12), and the whole comparison treats TMC as the reference. The a posteriori validation is reported only for the reference model and only as a three-zone classification (dust-dominated, transition, gas-dominated). Because the Rosseland opacity changes by orders of magnitude across the 200-K sublimation transition (Eq. (11) and Fig. 1), a cell can be classified correctly while its opacity at TMC differs substantially from κR(TRHD); the N3 300M⊕ model reaches or exceeds the sublimation temperature in the planetary cell. The authors should provide a quantitative sensitivity test, for example an iterative MCRT run for one model or a subset of cells, or a map of |κR(TMC) − κR(TRHD)|/κR(TRHD) throughout the high-resolution region, and they should report the zone-agreement statistics for all 18 models rather than only the reference. Without this, the reported ~10% agreement could be biased by an uncontrolled systematic in the reference solution itself. The authors' own list of future work (Sect. 4.2: 'consistent Planck opacities for the dust-evaporation regions') confirms that the issue is live.","section":"Section 2.2, Eq. (12); Section 3.1.1; Section 4.2"},{"comment":"The claim that the observed temperature differences are 'systematic' rather than Monte Carlo noise is not backed by a quantitative noise estimate. The MCRT temperature calculation uses 10^8 photon packages, but no repeat-run or cell-wise variance estimate is provided. In the NIR flux maps (Sect. 3.3.2), the authors themselves state that differences outside the inner region are 'completely dominated by MC noise', yet the temperature-difference histograms in Figs. 5 and 8 are presented without error bars. A simple estimate of the MCRT temperature noise floor would substantiate the systematic-origin argument and make the 10% figure more robust.","section":"Section 3.1.1; Section 3.3.2"}],"minor_comments":[{"comment":"The stellar radius is given as 'R∗ = 2.5 L⊙'; this should be 'R∗ = 2.5 R⊙' to agree with Sect. 2.1.1.","section":"Section 2.2, first paragraph"},{"comment":"The word 'axisymetric' should be 'axisymmetric'.","section":"Section 2.1.1"},{"comment":"The quantity δT is an antisymmetric relative difference; the paper should state this explicitly when translating δT into 'percent' values, since δT ≈ ΔT/(2T) for small differences.","section":"Section 3.1.1, Eq. (13)"},{"comment":"The statement that NIR relative differences outside the inner region are 'completely dominated by MC noise' should be accompanied by a quantitative noise estimate for the flux maps, otherwise the reader cannot separate noise from genuine temperature-driven flux differences in the NIR.","section":"Section 3.3.2"},{"comment":"The phrase 'typical discrepancy' is used without a formal definition; specifying whether 'typical' means the mean, median, or standard deviation of the histogram would make the headline number reproducible.","section":"Section 4, first bullet list"}],"recommendation":"major_revision","confidential_remarks":"The paper is a careful and useful comparison, and the central methodology is sound in principle. The two main issues are fixable: (i) the 'all cells |δT|≲0.1–0.13' bullet is overstated relative to the paper's own figures and should be replaced by a quantile or region-specific statement; (ii) the non-iterative MCRT opacity freezing needs a quantitative sensitivity check beyond the reference-model zone-agreement statistic. The authors already identify related opacity-consistency improvements in Sect. 4.2, so the concern is acknowledged but not yet resolved. I would be comfortable with acceptance after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a genuinely useful error census, and the central numbers are probably about right in the bulk of the disk, but the '10 percent' headline is easy to misread. They define δT = (TRHD - TMC)/(TRHD + TMC), so δT ≈ 0.1 means about a 20–22% temperature difference in the usual sense. Keep that in mind when you see 'about 10 percent.'\n\nWhat's new: it is the first systematic parameter sweep (two planet masses, three resolutions, three accretion timescales, 18 models) that quantifies where FLD-based RHD and MCRT temperatures diverge for embedded accreting planets, and it connects those differences to flux errors at VIS/NIR/submm. The FLD flux-crossing diagnosis for midplane overheating is plausible and well supported by the controlled axisymmetric runs that switch irradiation and viscosity on and off. The authors are also honest that MCRT is not a 'true' reference: it lacks PdV work and internal-energy advection, which matters in spiral wakes.\n\nWhere it's soft: the non-iterative MCRT opacity fixing (Eq. 12) uses κR(TRHD) to set the dust-to-gas ratio per cell, and the a posteriori validation is only zone classification (dust vs. sublimation) for the reference model. The N3 300 M⊕ model reaches Tevp inside the planetary cell, where the Rosseland opacity is steep; a 10% temperature error there could shift the opacity a lot. The paper argues that cells crossing sublimation become optically thin and insensitive to the exact opacity, which mitigates the worry, but they don't demonstrate that for all 18 models. So the ~10% agreement claim is not bulletproof near sublimation fronts, though it's probably robust in the dust-dominated bulk.\n\nAlso minor: the abstract says flux errors are 'independently of the amount of gas piling up in the Hill sphere and the used model resolution,' but flux maps were only made for the A32 accretion timescale and N1/N3 resolutions, not the full grid. And no code or data are released, which hampers reproducibility.\n\nOverall: worth a serious referee. This delivers a quantitative error budget the community will use; it needs a clearer statement of the δT metric and a broader validity check for the non-iterative MCRT step. I'd bring it to a reading group.","headline":"A useful parameter-swept error census for RHD-vs-MCRT temperatures and fluxes in accreting-planet disks, but the '10 percent agreement' headline is a symmetric δT (closer to 20% in ordinary relative terms) and the non-iterative MCRT opacity check is thinner than the central claim needs.","tokens_in":33763,"tokens_out":5453,"would_cite":true,"duration_ms":46073,"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":"This paper argues that flux-limited diffusion radiation hydrodynamics and Monte Carlo radiative transfer yield gas temperatures agreeing to about 10 percent in the planet-hosting region of a protoplanetary disk, with systematic regional…","keywords":["radiation hydrodynamics","Monte Carlo radiative transfer","protoplanetary disks","planet-disk interaction","accreting planets","flux-limited diffusion","synthetic observations","temperature comparison"],"falsifier":"Recompute the Monte Carlo temperatures for the same 18 snapshots with fully iterative, temperature-dependent opacities and compare zone classifications cell by cell; if substantial populations of cells cross the sublimation boundary, the fixed-opacity premise is the thing to blame, whereas if the temperatures stay within the reported $|\\delta T| \\lesssim 0.1$ bounds, the premise is exonerated.","tokens_in":32613,"feed_emoji":"🪐","tokens_out":7103,"duration_ms":61672,"temperature":0.7,"pith_summary":"The paper is trying to establish how much a fast radiation-hydrodynamics (RHD) model can be trusted when its temperatures are recalculated with a more complete Monte Carlo radiative transfer (MCRT) scheme for accreting planets in protoplanetary disks. Across 18 models spanning two planet masses, three grid resolutions, and three accretion timescales, the two temperature distributions match within about 10 percent in the high-resolution region, measured by the antisymmetric difference $\\delta T = (T_{\\rm RHD}-T_{\\rm MC})/(T_{\\rm RHD}+T_{\\rm MC})$. The agreement is not uniform: the RHD model runs warmer in the dense midplane and colder at the photosphere and in the far outer disk. The paper attributes these systematic differences to angle- and frequency-averaging and missing scattering in the diffusion treatment, and to missing compressional heating and internal-energy advection in the Monte Carlo step. This matters because synthetic images of forming planets are only as reliable as the temperature maps they are built from.","feed_headline":"RHD and Monte Carlo models agree to ~10 percent near planets","feed_subtitle":"Systematic biases traced to flux-angle averaging and missing compression heating; synthetic flux maps shift by up to 100 percent regionally.","key_machinery":"The comparison is carried by a non-iterative opacity-matching step: in each Monte Carlo cell the dust-to-gas ratio is lowered so that the cell's Rosseland mean opacity equals the value $\\kappa_R(T_{\\rm RHD})$ already used in the RHD simulation, allowing both solvers to work with the same optical properties without iterating temperature-dependent opacities. On top of this sits the antisymmetric relative difference $\\delta T = (T_{\\rm RHD}-T_{\\rm MC})/(T_{\\rm RHD}+T_{\\rm MC})$, which maps the agreement region by region. The authors verify a posteriori that 98–100 percent of cells fall in the same dust-sublimation zone under both temperature fields, so the fixed-opacity premise does not drive the result.","core_discovery":"The central claim is that, for local high-resolution radiation-hydrodynamical models of a 10 or 300 Earth-mass planet embedded in a global protoplanetary disk, flux-limited diffusion plus stellar ray tracing and a Monte Carlo radiative-transfer post-processing step give essentially the same gas temperatures: for all cells in the high-resolution region and all tested resolutions and accretion parameters, $|\\delta T| \\lesssim 0.1$–$0.13$, with the planet's Hill region on average $0.01$–$0.1$ colder in the RHD run. The discrepancies that remain are systematic and regional rather than Monte Carlo noise. The RHD model is too warm in the optically thick midplane because the angle-averaged diffusion treatment lets the inward stellar flux and the outward cooling flux interact unphysically; it is too cool just above the photosphere and in the outer disk because frequency-averaged opacities and the neglect of scattering suppress heating there. The MCRT model, in turn, cannot represent negative energy sources such as local compressional work and internal-energy advection, which shows up in spiral wakes. The paper concludes that neither method alone is the \"true\" temperature and that the path to closer agreement runs through multigroup or angularly discretized transport on the hydro side and MCRT schemes that can handle negative source terms.","pith_inferences":["If the artificial flux-flux interaction at the $\\tau=1$ surface is the main midplane bias, then multigroup or angularly discretized transport should push RHD temperatures down toward MCRT values specifically in the densest regions; the paper names such schemes as the next step without testing them here.","The missing compressional work in MCRT implies that post-processed temperatures are least trustworthy exactly where planet-driven spirals compress gas, so time-dependent or negative-source-term MCRT would likely raise inferred spiral-wake temperatures for massive planets.","A practical consequence the paper does not spell out: synthetic observations meant to constrain planet mass from gap or spiral morphology should be computed from MCRT temperatures, while hydrodynamical conclusions can keep using RHD temperatures, because the two uses have different error budgets."],"forward_implications":["In viscosity-dominated optically thick regions and irradiation-dominated optically thin layers the two methods agree closely, so existing RHD temperatures are reliable there without Monte Carlo post-processing.","Switching from RHD to MCRT temperatures changes visual-wavelength flux maps negligibly, but NIR thermal dust flux rises by roughly 160 percent and submillimeter total flux changes by up to about 15 percent, with regional submillimeter differences exceeding 100 percent.","Raising resolution from N1 to N3 does not remove the systematic temperature bias; submillimeter flux similarity improves only in limited regions, such as the inner high-resolution zone and the circumplanetary region for the 300 Earth-mass planet.","Because neither solver contains all the relevant physics, the roughly 10 percent agreement represents a floor on temperature accuracy for this kind of accreting-planet model, not a convergence certificate."],"supporting_citations":[{"why":"Supplies the flux limiter that defines the FLD closure, the angle-averaging approximation whose bias is central to the temperature discrepancies.","marker":"Levermore & Pomraning (1981)"},{"why":"Supplies the irradiation ray-tracing and implicit diffusion scheme used in the RHD runs.","marker":"Kuiper et al. (2010)"},{"why":"Provides the MCRT code that carries the post-processing temperature and flux calculations.","marker":"Ober et al. (2015)"},{"why":"Provides the continuous-absorption estimator that lets MCRT compute dust temperatures from photon energy deposition.","marker":"Lucy (1999)"},{"why":"Provides the immediate re-emission scheme used to speed up the MCRT temperature calculation.","marker":"Bjorkman & Wood (2001)"},{"why":"Supplies precomputed photon paths in optically thick media, enabling the large-scale MCRT simulations of the disks.","marker":"Krieger & Wolf (2020)"},{"why":"Sets the evaporation-temperature prescription shared by the RHD and MCRT opacity treatments.","marker":"Isella & Natta (2005)"},{"why":"Provides the Rosseland opacity floor used when dust evaporates, one of the identified sources of RHD-MCRT discrepancy.","marker":"Bell & Lin (1994)"}],"fun_headline_variants":["RHD and Monte Carlo temps agree within ~10% near planets","Radiation solver choice shifts disk flux maps by up to 100%","Both disk temperature models miss key heating effects","MCRT post-processing checks RHD: 10% typical, 40% worst","Flux-angle averaging and missing compression heat cause gaps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison rests on the assumption that fixing each Monte Carlo cell's opacity to the value from the hydro run does not push cells across the dust-sublimation boundary; if it did, the reported temperature differences would be biased, and the paper's own check finds 98–100 percent zone agreement after the fact.","fun_headline_variants_meta":{"raw":{"variants":["RHD and Monte Carlo temps agree within ~10% near planets","Radiation solver choice shifts disk flux maps by up to 100%","Both disk temperature models miss key heating effects","MCRT post-processing checks RHD: 10% typical, 40% worst","Flux-angle averaging and missing compression heat cause gaps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000822,"raw_usage":{"total_tokens":3697,"prompt_tokens":1144,"completion_tokens":2553,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":760,"completion_tokens_details":{"reasoning_tokens":2464}},"tokens_in":760,"tokens_out":2553,"duration_ms":20278,"temperature":1.0,"reasoning_tokens":2464,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:53:36.744955+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Monte Carlo temperatures for the same 18 snapshots with fully iterative, temperature-dependent opacities and compare zone classifications cell by cell; if substantial populations of cells cross the sublimation boundary, the fixed-opacity premise is the thing to blame, whereas if the temperatures stay within the reported $|\\delta T| \\lesssim 0.1$ bounds, the premise is exonerated.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the flux limiter that defines the FLD closure, the angle-averaging approximation whose bias is central to the temperature discrepancies."},{"cited_title":"2010, A&A, 511, A81","cited_arxiv_id":null,"evidence_quote":"Supplies the irradiation ray-tracing and implicit diffusion scheme used in the RHD runs."},{"cited_title":"L., & Klahr, H","cited_arxiv_id":null,"evidence_quote":"Provides the MCRT code that carries the post-processing temperature and flux calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the immediate re-emission scheme used to speed up the MCRT temperature calculation."},{"cited_title":"& Wolf, S","cited_arxiv_id":null,"evidence_quote":"Supplies precomputed photon paths in optically thick media, enabling the large-scale MCRT simulations of the disks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Rosseland opacity floor used when dust evaporates, one of the identified sources of RHD-MCRT discrepancy."}],"review_version":1}