Issue |
A&A
Volume 666, October 2022
|
|
---|---|---|
Article Number | A101 | |
Number of page(s) | 14 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202244521 | |
Published online | 14 October 2022 |
Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission
Sterrenkundig Observatorium, Universiteit Gent,
Krijgslaan 281 S9,
9000
Gent, Belgium
e-mail: maarten.baes@ugent.be
Received:
16
July
2022
Accepted:
11
August
2022
Context. The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however.
Aims. We explore alterations to the standard photon packet life cycle in Monte Carlo radiative transfer that allow the treatment of net stimulated emission without loss of generality or efficiency.
Methods. We present the explicit absorption technique that allows net stimulated emission to be handled efficiently. It uses the scattering rather than the extinction optical depth along a photon packet’s path to randomly select the next interaction location, and offers a separate, deterministic treatment of absorption. We implemented the technique in a special-purpose Monte Carlo code for a two-stream 1D radiative transfer problem and in the fully featured 3D code SKIRT, and we studied its overall performance using quantitative statistical tests.
Results. Our special-purpose code is capable of recovering the analytical solutions to the two-stream problem in all regimes, including the one of strong net stimulated emission. The implementation in SKIRT is straightforward, as the explicit absorption technique easily combines with the variance reduction and acceleration techniques already incorporated. In general, explicit absorption tends to improve the efficiency of the Monte Carlo routine in the regime of net absorption.
Conclusions. Explicit absorption allows the treatment of net stimulated emission in Monte Carlo radiative transfer, it interfaces smoothly with other variance reduction and acceleration techniques, and it tends to improve the efficiency of the simulations in the net absorption regime. We recommend to always include this new technique in Monte Carlo radiative transfer.
Key words: radiative transfer / methods: numerical
© M. Baes et al. 2022
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.