SPHRAY: a smoothed particle hydrodynamics ray tracer for radiative transfer
Abstract
We introduce the publicly available code SPHRAY, a smoothed particle hydrodynamics (SPH) ray tracer designed to solve the 3D, timedependent, radiative transfer equation for cosmological density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) raytracing scheme that does not interpolate the SPH particles on to a grid but instead integrates directly through the SPH kernels. Given an arbitrary (series of) SPH density field(s) and a description of the sources of ionizing radiation, the code will calculate the nonequilibrium ionization and temperature state of hydrogen (HI, HII) and helium (HeI, HeII, HeIII). The sources of radiation can include pointlike objects, diffuse recombination radiation, and a background field from outside the computational volume. The MC raytracing implementation allows for the quick introduction of new physics and is parallelizationfriendly. A quick Axis Aligned Bounding Box (AABB) test taken from computer graphics applications allows for the acceleration of the raytracing component. We present the algorithms used in SPHRAY and verify the code by performing the test problems detailed in the recent Radiative Transfer Comparison Project of Iliev et. al. The source code for SPHRAY and example SPH density fields are made available on a companion website (http://www.sphray.org).
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 June 2008
 DOI:
 10.1111/j.13652966.2008.13212.x
 arXiv:
 arXiv:0802.3698
 Bibcode:
 2008MNRAS.386.1931A
 Keywords:

 radiative transfer;
 methods: numerical;
 cosmology: theory;
 Astrophysics
 EPrint:
 17 pages, 16 figures, submitted to MNRAS, comments welcome. source code, high res. figures and examples can be found at http://www.sphray.org