Accelerating Dust Temperature Calculations with Graphics Processing Units

Abstract

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 69 times faster than the 8 CPU cores, showing great potential for accelerating calculations of galaxy SEDs.

7 pages, 2 figures, accepted to New Astronomy. Minor updates to text and performance based on feedback from referee

References in corpus (7)

• Infrared Emission from Interstellar Dust. IV. The Silicate-Graphite-PAH Model in the Post-Spitzer Era
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg
• Modeling the Dust Properties of z ~ 6 Quasars with ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree
• Parallel Algorithm for Solving Kepler's Equation on Graphics Processing Units: Application to Analysis of Doppler Exoplanet Searches
• Panchromatic Spectral Energy Distributions of Dusty Galaxies with RADISHE. I. Predictions for Herschel: Correlating Colors with Galactic Energy Sources
• High-Resolution Panchromatic Spectral Models of Galaxies including Photoionisation and Dust

Cited by in corpus (34)

• Dusty Star-Forming Galaxies at High Redshift
• A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law
• Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3
• Galaxy Properties from the Ultra-violet to the Far-Infrared: Lambda-CDM models confront observations
• What Does a Submillimeter Galaxy Selection Actually Select? The Dependence of Submillimeter Flux Density on Star Formation Rate and Dust Mass
• Submillimetre galaxies in a hierarchical universe: number counts, redshift distribution, and implications for the IMF
• Galaxy Morphology and Star Formation in the Illustris Simulation at z=0
• Three-Dimensional Dust Radiative Transfer
• The CO-H2 Conversion Factor in Disc Galaxies and Mergers
• A Theory for the Excitation of CO in Star Forming Galaxies
• How to distinguish starbursts and quiescently star-forming galaxies: The `bimodal' submillimetre galaxy population as a case study
• Should we believe the results of UV-mm galaxy SED modelling?
• Deep Learning Identifies High-z Galaxies in a Central Blue Nugget Phase in a Characteristic Mass Range
• Diverse Structural Evolution at z > 1 in Cosmologically Simulated Galaxies
• The evolution of galaxy shapes in CANDELS: from prolate to oblate
• Distinguishing Mergers and Disks in High Redshift Observations of Galaxy Kinematics
• NIHAO XIII: Clumpy discs or clumpy light in high redshift galaxies?
• Automated Distant Galaxy Merger Classifications from Space Telescope Images using the Illustris Simulation
• Powderday: Dust Radiative Transfer for Galaxy Simulations
• Star Formation and Clumps in Cosmological Galaxy Simulations with Radiation Pressure Feedback
• The Origin of the Microlensing Events Observed Towards the LMC and the Stellar Counterpart of the Magellanic Stream
• Stellar Masses of Giant Clumps in CANDELS and Simulated Galaxies Using Machine Learning
• Computational advances in gravitational microlensing: a comparison of CPU, GPU, and parallel, large data codes
• Cosmological Implications of a Stellar Initial Mass Function that Varies with the Jeans Mass in Galaxies
• The Evolution of Primordial Binary Open Star Clusters: Mergers, Shredded Secondaries and Separated Twins
• The Cosmic Evolution of the IMF Under the Jeans Conjecture with Implications for Bottom-Heavy Ellipticals
• GPU-based Monte Carlo dust radiative transfer scheme applied to AGN
• PyCOOL - a Cosmological Object-Oriented Lattice code written in Python
• Analysing Astronomy Algorithms for GPUs and Beyond
• Why is the Milky Way X-factor Constant?
• The nature of giant clumps in high-z discs: a deep-learning comparison of simulations and observations
• Knudsen number dependence of 2D single-mode Rayleigh-Taylor fluid instabilities
• Modelling the spectral energy distribution of galaxies: introducing the artificial neural network
• Using hybrid GPU/CPU kernel splitting to accelerate spherical convolutions