Accelerating dust temperature calculations with graphics-processing units

When calculating the infrared spectral energy distributions (SEDs) of galaxies in radiation-transfer models, the calculation of dust grain temperatures is generally the most time-consuming part of the calculation. Because of its highly parallel nature, this calculation is perfectly suited for massively parallel general-purpose Graphics Processing Units (GPUs). This paper presents an implementation of the calculation of dust grain equilibrium temperatures on GPUs in the Monte-Carlo radiation transfer code sunrise, using the CUDA API. The GPU can perform this calculation 55 times faster than the 8 CPU cores, showing great potential for accelerating calculations of galaxy SEDs.

[1]  Cambridge,et al.  GPUs for data processing in the MWA , 2009 .

[2]  Ari Laor,et al.  Spectroscopic constraints on the properties of dust in active galactic nuclei , 1993 .

[3]  M. Wolff,et al.  The DIRTY Model. II. Self-consistent Treatment of Dust Heating and Emission in a Three-dimensional Radiative Transfer Code , 2000, astro-ph/0011576.

[4]  P. Jonsson.  sunrise: polychromatic dust radiative transfer in arbitrary geometries , 2006, astro-ph/0604118.

[5]  Adolf N. Witt,et al.  The DIRTY Model. I. Monte Carlo Radiative Transfer through Dust , 2000, astro-ph/0011575.

[6]  B. Draine,et al.  Temperature fluctuations in interstellar grains. I. Computational method and sublimation of small grains , 1989 .

[7]  B. Whitney,et al.  PANCHROMATIC SPECTRAL ENERGY DISTRIBUTIONS OF DUSTY GALAXIES WITH RADISHE. I. PREDICTIONS FOR HERSCHEL: CORRELATING COLORS WITH GALACTIC ENERGY SOURCES , 2007, 0711.4361.

[8]  Todd L. Veldhuizen,et al.  Arrays in Blitz++ , 1998, ISCOPE.

[9]  Robert G. Belleman,et al.  High Performance Direct Gravitational N-body Simulations on Graphics Processing Units , 2007, ArXiv.

[10]  J. Weingartner,et al.  Dust Grain-Size Distributions and Extinction in the Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud , 2001 .

[11]  S. Bianchi.  Dust extinction and emission in a clumpy galactic disk - An application of the radiative transfer code TRADING , 2008, 0807.2946.

[12]  B. Draine,et al.  Infrared Emission from Interstellar Dust. IV. The Silicate-Graphite-PAH Model in the Post-Spitzer Era , 2006, astro-ph/0608003.

[13]  Simon Portegies Zwart,et al.  SAPPORO: A way to turn your graphics cards into a GRAPE-6 , 2009, ArXiv.

[14]  Linhua Jiang,et al.  Modeling the Dust Properties of z ~ 6 Quasars with ART2—All-Wavelength Radiative Transfer with Adaptive Refinement Tree , 2007, 0706.3706.

[15]  Eric B. Ford,et al.  Parallel algorithm for solving Kepler’s equation on Graphics Processing Units: Application to analysis of Doppler exoplanet searches , 2008, 0812.2976.