Optical depth in polarised Monte Carlo radiative transfer

Abstract

Context. The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the concept of optical depth. Aims. We investigate in detail the relation between optical depth and the optical properties and density of the attenuating medium in polarised Monte Carlo radiative transfer codes that take into account dichroic extinction. Methods. Based on solutions for the radiative transfer equation, we discuss the optical depth scale in polarised radiative transfer with spheroidal grains. We compare the dichroic optical depth to the extinction and total optical depth scale. Results. In a dichroic medium, the optical depth is not equal to the usual extinction optical depth, nor to the total optical depth. For representative values of the optical properties of dust grains, the dichroic optical depth can differ from the extinction or total optical depth by several ten percent. A closed expression for the dichroic optical depth cannot be given, but it can be derived efficiently through an algorithm that is based on the analytical result corresponding to elongated grains with a uniform grain alignment.Conclusions. Optical depth is more complex in dichroic media than in systems without dichroic attenuation, and this complexity needs to be considered when generating random free path lengths in Monte Carlo radiative transfer simulations. There is no benefit in using approximations instead of the dichroic optical depth.

Linear spectropolarimetric analysis of fairall 9 with VLT/FORS2

“…In order to investigate various effects of different BLR geometry and kinematics on the polarization spectra, we performed three simulations using radiative transfer codes 6 (Goosmann & Gaskell 2007;Marin et al 2012Marin et al , 2015Rojas Lobos et al 2018;Marin 2018) and 7 (Baes et al 2019;Camps & Baes 2020, and references therein). We assumed that the BLR emission is being scattered by the inner part of the dusty torus.…”

SKIRT 9: Redesigning an advanced dust radiative transfer code to allow kinematics, line transfer and polarization by aligned dust grains

“…This means that the polarization state of radiation changes as it moves through the medium and that the extinction by the medium depends on the local polarization state of the radiation. As a result, the calculations to move a photon packet along its path become significantly more complicated (see, e.g., Mishchenko 1991;Baes et al 2019).…”

Thermal Emission and Scattering by Aligned Grains: Plane-Parallel Model and Application to Multiwavelength Polarization of the HL Tau Disk

“…Monte Carlo techniques have been a powerful tool to treat scattering for the three dimensional structure of disks. However, the complete treatment of scattering including aligned grains is notoriously difficult and most of the work has been limited to spherical or randomly aligned grains (see e.g., Dullemond et al 2012, Steinacker et al 2013, Baes et al 2019. Given that the dust in protoplanetary disks can be fairly geometrically thin (e.g., Dutrey et al 2017;Villenave et al 2020) and especially for HL Tau (Pinte et al 2016), we can use a plane-parallel slab to capture the essence of the problem.…”