Issue |
A&A
Volume 694, February 2025
|
|
---|---|---|
Article Number | A52 | |
Number of page(s) | 31 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202451780 | |
Published online | 31 January 2025 |
Monte Carlo post-processing for radiation hydro simulations of accreting planets in protoplanetary disks
1
Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel,
Leibnizstraße 15,
24118
Kiel,
Germany
2
Max Planck Institut für Astronomie,
Königstuhl 17,
69117
Heidelberg,
Germany
★ Corresponding authors; akrieger@astrophysik.uni-kiel.de; klahr@mpia.de
Received:
3
August
2024
Accepted:
2
December
2024
This paper is part of a series investigating the observational appearance of planets accreting from their nascent protoplanetary disk (PPD). We evaluate the differences between gas temperature distributions determined in our radiation hydrodynamical (RHD) simulations and those recalculated via post-processing with a Monte Carlo (MC) radiative transport (RT) scheme. Our MCRT simulations were performed for global PPD models, each composed of a local 3D high-resolution RHD model embedded in an axisymmetric global disk simulation. We report the level of agreement between the two approaches and point out several caveats that prevent a perfect match between the temperature distributions with our respective methods of choice. Overall, the level of agreement is high, with a typical discrepancy between the RHD and MCRT temperatures of the high-resolution region of only about 10 percent. The largest differences were found close to the disk photosphere, at the transition layer between optically dense and thin regions, as well as in the far-out regions of the PPD, occasionally exceeding values of 40 percent. We identify several reasons for these discrepancies, which are mostly related to general features of typical radiative transfer solvers used in hydrodynamical simulations (angle- and frequency-averaging and ignored scattering) and MCRT methods (ignored internal energy advection and compression and expansion work). This provides a clear pathway to reduce systematic temperature inaccuracies in future works. Based on MCRT simulations, we finally determined the expected error in flux estimates, both for the entire PPD and for planets accreting gas from their ambient disk, independently of the amount of gas piling up in the Hill sphere and the used model resolution.
Key words: accretion, accretion disks / hydrodynamics / radiative transfer / methods: numerical / protoplanetary disks / planet-disk interactions
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.