Please use this identifier to cite or link to this item:
|Title:||Benchmarking the calculation of stochastic heating and emissivity of dust grains in the context of radiative transfer simulations||Authors:||Camps, Peter
Fitzgerald, Michael P.
|Issue Date:||2015||Journal:||ASTRONOMY & ASTROPHYSICS||Number:||580||First Page:||A87||Abstract:||Context. Thermal emission by stochastically heated dust grains (SHGs) plays an important role in the radiative transfer (RT) problem for a dusty medium. It is therefore essential to verify that RT codes properly calculate the dust emission before studying the effects of spatial distribution and other model parameters on the simulated observables. <BR /> Aims: We define an appropriate problem for benchmarking dust emissivity calculations in the context of RT simulations, specifically including the emission from SHGs. Our aim is to provide a self-contained guide for implementors of such functionality and to offer insight into the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. <BR /> Methods: The benchmark problem definition includes the optical and calorimetric material properties and the grain size distributions for a typical astronomical dust mixture with silicate, graphite, and PAH components. It also includes a series of analytically defined radiation fields to which the dust population is to be exposed and instructions for the desired output. We processed this problem using six RT codes participating in this benchmark effort and compared the results to a reference solution computed with the publicly available dust emission code DustEM. <BR /> Results: The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range 3 μm ≤ λ ≤ 1000 μm. <BR /> Conclusions: We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of SHGs. This work can serve as a reference for implementors of dust RT codes, and it will pave the way for a more extensive benchmark effort focusing on the RT aspects of the various codes.||Acknowledgments:||We thank the authors of DustEM ( Compiègne et al. 2011 ) for making their code publicly available. This work fits in the CHARM framework (Contemporary physical challenges in Heliospheric and AstRophysical Models), a Phase-VII Interuniversity Attraction Pole (IAP) program organized by BELSPO, the BELgian federal Science Policy Office. T.L. acknowledges the support from the Swedish National Space Board (SNSB). G.N. is supported by the Leverhulme Trust research project grant RPG-2013-418. M.J. acknowledges the support of the Academy of Finland Grant No. 250741. J.S. acknowledges support from the ANR (SEED ANR-11-CHEX-0007-01).||URI:||http://hdl.handle.net/20.500.12386/23133||URL:||https://www.aanda.org/articles/aa/abs/2015/08/aa25998-15/aa25998-15.html||ISSN:||0004-6361||DOI:||10.1051/0004-6361/201525998||Bibcode ADS:||2015A&A...580A..87C||Fulltext:||open|
|Appears in Collections:||1.01 Articoli in rivista|
Show full item record
checked on Sep 30, 2020
checked on Sep 30, 2020
Items in DSpace are published in Open Access, unless otherwise indicated.